Comparative kinetic analysis of two fungal β-glucosidases
© Chauve et al; licensee BioMed Central Ltd. 2010
Received: 19 November 2009
Accepted: 11 February 2010
Published: 11 February 2010
The enzymatic hydrolysis of cellulose is still considered as one of the main limiting steps of the biological production of biofuels from lignocellulosic biomass. It is a complex multistep process, and various kinetic models have been proposed. The cellulase enzymatic cocktail secreted by Trichoderma reesei has been intensively investigated. β-glucosidases are one of a number of cellulolytic enzymes, and catalyze the last step releasing glucose from the inhibitory cellobiose. β-glucosidase (BGL1) is very poorly secreted by Trichoderma reesei strains, and complete hydrolysis of cellulose often requires supplementation with a commercial β-glucosidase preparation such as that from Aspergillus niger (Novozymes SP188). Surprisingly, kinetic modeling of β-glucosidases lacks reliable data, and the possible differences between native T. reesei and supplemented β-glucosidases are not taken into consideration, possibly because of the difficulty of purifying BGL1.
A comparative kinetic analysis of β-glucosidase from Aspergillus niger and BGL1 from Trichoderma reesei, purified using a new and efficient fast protein liquid chromatography protocol, was performed. This purification is characterized by two major steps, including the adsorption of the major cellulases onto crystalline cellulose, and a final purification factor of 53. Quantitative analysis of the resulting β-glucosidase fraction from T. reesei showed it to be 95% pure. Kinetic parameters were determined using cellobiose and a chromogenic artificial substrate. A new method allowing easy and rapid determination of the kinetic parameters was also developed. β-Glucosidase SP188 (Km = 0.57 mM; Kp = 2.70 mM) has a lower specific activity than BGL1 (Km = 0.38 mM; Kp = 3.25 mM) and is also more sensitive to glucose inhibition. A Michaelis-Menten model integrating competitive inhibition by the product (glucose) has been validated and is able to predict the β-glucosidase activity of both enzymes.
This article provides a useful comparison between the activity of β-glucosidases from two different fungi, and shows the importance of fully characterizing both enzymes. A Michaelis-Menten model was developed, including glucose inhibition and kinetic parameters, which were accurately determined and compared. This model can be further integrated into a cellulose hydrolysis model dissociating β-glucosidase activity from that of other cellulases. It can also help to define the optimal enzymatic cocktails for new β-glucosidase activities.
Discovered four decades ago by Reese and Mandels [1, 2], the cellulolytic enzyme system secreted by the filamentous fungi Trichoderma reesei (synonym Trichoderma viride) is the initial parent of most fungal strains for industrial cellulase production. The hydrolysis step converting cellulose to glucose is recognized as the major limiting step in the development of biological processes for production of biofuels from lignocellulosic raw materials because of the low efficiency of cellulases and their cost.
The enzymatic hydrolysis of cellulose involves three types of cellulases (cellobiohydrolases, endoglucanases and β-glucosidases) working in synergy . Endoglucanases (EC 126.96.36.199) randomly cleave the β-1,4 glycosidic linkages of cellulose; cellobiohydrolases (EC 188.8.131.52) attack cellulose chain ends to produce the constitutive unit of cellulose, cellobiose (a dimer of glucose linked by a β-1,4 glycosidic bond); and β-glucosidases (EC 3.2.1 21) hydrolyse cellobiose into two molecules of glucose. In fact, enzymatic hydrolysis of cellulose is a six-step complex process, the last step being a homogenous catalysis reaction involving the action of β-glucosidase on cellobiose. Cellobiose is a strong inhibitor of both cellobiohydrolases and endocellulases, and the β-glucosidase action can reduce its effect. In addition, the produced glucose also inhibits cellulolysis, although to a lesser extent . Unexpectedly, the amount of β-glucosidase-1 (BGL1) generated by T. reesei hyperproducing strains represents a very low percentage of the total secreted proteins [3, 5]. This very low level of β-glucosidase activity  often limits the amount of this enzyme in commercial cellulase preparations. This limitation can be alleviated either by overexpressing β-glucosidase in T. reesei or by adding extra β-glucosidase from other sources [7, 8]. Supplementing the native T. reesei enzymatic cocktail with β-glucosidase from other fungi is often performed to avoid inhibition of cellobiose in standardized hydrolysis test . The most common commercial β-glucosidase preparation (Novozymes SP188; Novo Nodisk A/S, Bagsvaerd, Danmark) is produced by Aspergillus niger.
Kinetic parameters of β-glucosidase from Aspergillus niger and Trichoderma reesei
Novozymes S188 (A. niger)
Sigma (A. niger)
20 to 80
A. niger (Culture)
Cellobiase 250L (A. niger)
C + S
T. reesei Rut C30
C - UC - NC + S
T. reesei (Culture)
T. reesei QM9414
At this stage, a comparison of β-glucosidases from A. niger and T. reesei is required to validate a model of cellulose hydrolysis using cellulase mixtures supplemented or not with β-glucosidase. To our knowledge, both enzymes have been previously compared only once , and without any previous purification, which can greatly affect the values of the kinetic constants determined and the representativeness of the deduced models. We focused our study on the parameters important for future industrial utilization of the enzymes in a biofuel production process, especially in quantifying the effects of cellobiose, glucose and temperature. Because β-glucosidase from T. reesei represents a very low proportion of the cellulolytic mixture secreted, a procedure for its purification was developed. In addition, a specific method suitable for a rapid reaction time course was also used. This led to the first consistent comparative kinetic analysis between A. niger and T. reesei β-glucosidase on their natural and artificial substrates, and can be further used for modeling addition of extra β-glucosidase.
Enzyme samples, substrates, buffer and reagents
A commercial β-glucosidase preparation (SP188; Novozymes), microcrystalline cellulose (Avicel PH101), cellobiose, p-nitrophenol and ρ-nitrophenyl-β-D-glucopyranoside (pNPG) used in this study were purchased from Sigma-Aldrich (Sigma-Aldrich, Lyon, France). Glucose was obtained from Prolabo. All buffer components and salts used were reagent grade and obtained from Sigma Aldrich or GE Healthcare (GE Healthcare, Saclay, France) for enzyme purification.
T. reesei β-glucosidase was obtained from a crude enzyme preparation using the hypercellulolytic mutant strain CL847 . Enzymatic productions were performed in 2.5 L working volume fermentors using two steps: (i) growth on glucose,(ii) fed batch with a sugar mixture containing glucose, lactose and xylose as the carbon source, as previously described [5, 14]. At the end of the fermentation, the mycelia were removed by centrifugation from the culture medium at 20,000 g and the clear medium containing the enzymes secreted by T. reesei (40 mg proteins per mL) was further used for purification of β-glucosidase.
Preparation of β-glucosidase from SP188 and T. reesei CL847
Fast protein liquid chromatography (Akta FPLC System; GE Healthcare) was used for the following purification experiments.
For biochemical characterization of the β-glucosidase, the SP188 commercial preparation was desalted to remove the ingredients present in the industrial formulation. For this purification step, a desalting column (Hi-trap; GE Healthcare) previously equilibrated with 25 mM imidazole pH 7.6 was used.
The purified T. reesei β-glucosidase used in this study was prepared as follows. A sample (0.625 mL) of the cell-free medium containing 25 mg of extracellular proteins was mixed for 30 minutes at 4°C with 20 mL of 50 mM citrate buffer pH 4.8 containing 5% Avicel PH101. After centrifugation (6000 g for 10 minutes at 4°C), the supernatant was recovered; this step was repeated once using the same amount of fresh Avicel. This allows the adsorption of proteins containing a cellulose binding domain (> 80% of the total protein content of the initial crude sample). The resulting supernatant was concentrated using a centrifuge (Vivaspin 20; Sartorius stedim, Aubagne, France) with a molecular weight cut-off of 50 kDa, allowing separation between β-glucosidase and any remaining low molecular weight proteins. This β-glucosidase-enriched sample was then desalted as described for SP188. The desalted sample was then injected into a column (10/100 GL; MonoQ) previously equilibrated with 25 mM imidazole pH 7.6 (linear flow rate of 0.7 mL/min) using FPLC. Because β-glucosidase has a high apparent isoelectric point (pI) value (close to 8), highly purified β-glucosidase was recovered by collecting the non-retained fraction. The pH of the β-glucosidase fraction was adjusted to pH 4.8 using 1 M HCl for storage. When necessary, qualitative analysis of the sample was performed using one-dimensional (1D) gel electrophoresis according to Laemmli , and protein concentration was measured according to Lowry , using bovine serum albumin as the protein standard.
β-glucosidase activity was determined by measuring the amount of glucose and p-nitrophenol released, respectively, from cellobiose and p-nitrophenyl-β-D-glucopyranoside (pNPG) used as substrates. Native substrate assays (duplicate) involved incubating the substrate at concentrations from 0.23 to 16.4 mM in 50 mM citrate buffer pH 4.8 and enzyme solution (in a range from 2.5 to 10 mg/L) at 50°C with a time interval of 5 minutes and a total period of 20 minutes. The final volume of the assays was 100 μL and they were carried out in 1.5 mL Eppendorf tubes. The reaction was stopped by heating at 95°C for 5 min and cooled to 4°C. pNPG triplicate assays were performed using the same protocol as above but using a substrate range from 0.15 to 5 mM and and enzyme concentration of 1.25 to 5 mg/L. Reactions were prepared in 96-well PCR plates, and incubations were performed in a thermal cycler (PCR i-Cycler; Bio-Rad Laboratories, Marnes la coquette, France). Enzymatic activity was stopped by adding 100 μL of a 2% Na2CO3 solution.
For both pNPG and cellobiose substrates, inhibition parameters were determined using the same concentrations as above and glucose concentration up to 1 g/L. Inhibition effect of higher glucose concentrations (up to 60 g/L), was tested with pNPG as substrate.
The thermal stability was investigated by heating the enzyme to 50°C, removing aliquots every 24 h for 72 h in total, and assaying for kinetic measurement with 5 mM pNPG. Results showed no loss of activity within 72 h for both enzymes at 50°C (data not shown). Six additional experiments, in the validation domain of the model, were performed to evaluate the quality of the model prediction. These assays were not used to fit the parameters of the model.
Glucose concentration was determined using a high performance liquid chromatography system (ISC300; Dionex, Jouy En Josas, France) equipped with an amperometric detector. Data acquisition from the detector and determination of retention times and peak areas were performed using Chromeleon software (Dionex). A chromatography column (CarboPac PA1; measuring 250 mm × 4 mm internal diameter) with polystyrene-divinylbenzenzene pellicular resin (particle size 10 μm) and a pressure limit of 2000 psi was used. Purified water (18 MΩ; Milli-Q Plus Water Purification System; Millipore, Bedford, MA, USA) was used to prepare all solutions. Eluent 1 was 200 mM NaOH (Fischer Scientific, Illkirch, France) and eluent 2 was water (flow-rate: 1.0 mL/min). The detector operated in the amperometric mode using a gold working electrode and an Ag/AgCl pH reference electrode. Glucose, xylose and cellobiose solutions at 10 mg/L were used as external standards.
p-Nitrophenol was quantified by spectrophotometry at 405 nm (Multiskan Ascent Scanner; Thermo Labsystem, Issy les moulineaux, France) using 96-well plates. A range of p-nitrophenol solutions from 0 to 0.3 mM was used for calibration.
Km values were determined from the mean value of rates (in duplicates or triplicates) at six concentrations of substrate (PNPG or cellobiose), and Ki values were determined by the same procedure with four concentrations of glucose.
Results and discussion
Purification of β-glucosidases
Purification of T. reesei β-glucosidase may appear challenging because it represents < 2% of the secretome of hyperproducing strains [5, 17], which mainly contains cellobiohydrolases and endoglucanases. The whole purification procedure is detailed in Figure 1. The method we used does not involve any fixation of the target enzyme onto a column, so that it preserves its biochemical integrity. It takes advantage of the capacity of processive and non-processive cellulases to adsorb onto cellulosic polymers via their cellulose binding domains , which are absent in T. reesei BGL1. With this first purification step, 80% of the proteins contained in the original crude were trapped. Taking into account the theoretical and apparent molecular weight of β-glucosidase (close to 75 and 117 kDa, respectively), the enriched sample was concentrated using a 50 kDa cut-off membrane. Qualitative analysis was performed using 1D gel electrophoresis of the recovered fraction. We were able to detect a clear major band corresponding to the β-glucosidase and two other remaining proteins (estimated purity yield close to 60%) with apparent molecular weights of 100 kDa and 20 kDa. Using the pI of β-glucosidase, the concentrated sample was finally purified using a monoQ anionic column equilibrated at pH 7.4, to trap the last contaminating proteins in the sample. As described by Herpoël-Gimbert et al. , quantitative analysis of the resulting β-glucosidase fraction revealed a degree of purity of 95%.
Full kinetic model using β-glucosidase from A. niger
where [S]0 is the initial substrate concentration (mM), [P]0 is the initial product concentration (mM), [E]0 is the initial enzyme concentration (mM, with Mβ-glucosidase = 117 kDa ), Km is the M-M constant (mM), Kp is the inhibition constant (mM), kcat is the first order kinetic constant (per minute), and rini is the initial rate of hydrolysis (mM/min)
Kinetic parameters of the A. niger β-glucosidase preparation using a Michaelis-Menten model
K m , mM
K p competitive inhibition, mM
K m /K p
where k cat : first order kinetic constant (per min), A is the pre-exponential factor, E a is the activation energy (J/mol), R is the universal gas constant (8.314 J/mol/K), and T is the temperature (°K)
Additional experimental assays were performed to evaluate the quality of our model at temperatures between 40 and 60°C, at cellobiose concentrations between 0.5 and 16 mM, at glucose concentrations between 0 and 5 mM, and at enzyme concentrations between 2.5 and 10 mg/L. For each condition tested, the model predicts glucose conversion rate with high accuracy (4% error). A more important relative error was observed when very low enzyme concentrations, correlated with very low conversion level, were used.
Comparison of β-glucosidase from A. niger and T. reesei
Comparison of β-glucosidase from A. niger and from Trichoderma reesei for the Michaelis-Menten model
k cat, per min
1582 ± 79
5276 ± 264
1897 ± 95
2445 ± 107
K m , mM
K p competitive inhibition, mM
K m /K p
This work contributes to a more reliable kinetic model of cellulolysis by the determination of kinetic parameters of purified β-glucosidase of T. reesei using an M-M model. In addition, this model is also suitable for the prediction of a commercial β-glucosidase preparation obtained from A. niger (SP188) which is frequently added into reactions to achieve complete enzymatic hydrolysis of pretreated lignocellulosic raw materials. As it integrates the complete set parameters of β-glucosidase, this model can also be used for the prediction of the action of β-glucosidase from other microorganisms. It also includes the inhibition effect of the final product (glucose) and the influence of temperature. The comparison of β-glucosidase activities on its natural substrate (cellobiose) and on a chromophoric substrate (pNPG) shows that (i) β-glucosidase affinity is higher for pNPG, (ii) the global hydrolysis activity is higher for T. reesei β-glucosidase and (iii) the Km/KP ratios are similar for both substrates, confirming the competitive inhibition mechanism. Thus, it validates the use of pNPG as a substrate for the determination of the kinetic parameters, which can be easily determined using the 96-well plate method described here.
Comparing the β-glucosidases of T. reesei and A. niger was previously performed by Grous et al.  but they used non-purified enzymes, which can introduce some inaccuracy in the values of parameters. Using the innovative purification protocol based on the adsorption capacity of enzymes onto cellulose, we successfully purified β-glucosidase from T. reesei despite its low proportion in the crude preparation. Our results show that the presence of glucose does not equally affect β-glucosidase activities from T. reesei or A. niger. Glucose inhibition can have a deleterious effect in a bioethanol production process, in which high ethanol and accordingly high sugar concentrations are expected.
As already mentioned, this β-glucosidase model is the first part of a global kinetic model of the enzymatic hydrolysis of cellulose, dissociating the activity of β-glucosidase as a final homogeneous reaction and the hydrolysis of the cellulose by endo- and exo-glucanases in a heterogeneous phase. It can be used both with supplemented and non-supplemented cocktails, and allows prediction of the cellulose saccharification
MC is a recipient of an ADEME doctoral fellowship. We acknowledge the financial support given by FUTUROL project. We are also grateful to Sandra Menir for technical assistance, and to Antoine Margeot and Serge Perez for helpful discussions.
- Reese ET, Mandels M: Enzymatic degradation. Cellulose and Cellulose Derivatives. Edited by: Bikales NM, Segal L. 1971, New York: Wiley Interscience, 1079-1094.Google Scholar
- Mandels M, Reese ET: Induction of cellulose in Trichoderma viride as influenced by carbon source and metals. J Bact. 1957, 73: 269-278. 10.1002/path.1700730133.View ArticleGoogle Scholar
- Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS: Microbial cellulose utilization: Fundamentals and biotechnology. Microbiol Mol Biol Rev. 2002, 66: 506-577. 10.1128/MMBR.66.3.506-577.2002.View ArticleGoogle Scholar
- Dekker RFH: Kinetic, inhibition, and stability properties of a commercial beta-D-glucosidase (cellobiase) preparation from aspergillus-niger and its suitability in the hydrolysis of lignocellulose. Biotechnold Bioeng. 1986, 28: 1438-1442. 10.1002/bit.260280918.View ArticleGoogle Scholar
- Herpoël-Gimbert I, Margeot A, Dolla A, Jan G, Molle D, Lignon S, Mathis H, Sigoillot JC, Monot F, Asther M: Comparison secretome analysis of two Trichoderma reesei RUT-C30 and CL847 hypersecretory strains. Biotechnol Biofuels. 2008, 1: 1-18. 10.1186/1754-6834-1-18.View ArticleGoogle Scholar
- Kubicek CP: Release of carboxymethylcellulose and b-glucosidase from cell walls of Trichoderma reesei. Eur J Appl Microbial Biotechnol. 1981, 13: 226-231. 10.1007/BF00500103.View ArticleGoogle Scholar
- Xin Z, Yinbo Q, Peiji G: Acceleration of ethanol production from paper mill waste fiber by supplementation with [beta]-glucosidase. Enzyme Microb Technol. 1993, 15: 62-65. 10.1016/0141-0229(93)90117-K.View ArticleGoogle Scholar
- Kumar R, Singh S, Singh OV: Bioconversion of lignocellulosic biomass: biochemical and molecular perspectives. J Ind Microbiol Biotechnol. 2008, 35: 377-391. 10.1007/s10295-008-0327-8.View ArticleGoogle Scholar
- Xiao ZZ, Zhang X, Gregg DJ, Saddler JN: Effects of sugar inhibition on cellulases and beta-glucosidase during enzymatic hydrolysis of softwood substrates. Appl Biochem Biotechnol. 2004, 113-16: 1115-1126. 10.1385/ABAB:115:1-3:1115.View ArticleGoogle Scholar
- Zhang YHP, Lynd LR: Toward an aggregated understanding of enzymatic hydrolysis of cellulose: Noncomplexed cellulase systems. Biotechnol Bioeng. 2004, 88: 797-824. 10.1002/bit.20282.View ArticleGoogle Scholar
- Bansal P, Hall M, Realff MJ, Lee JH, Bommarius AS: Modeling cellulase kinetics on lignocellulosic substrates. Biotechnol Adv. 2009, 27: 833-48. 10.1016/j.biotechadv.2009.06.005.View ArticleGoogle Scholar
- Grous W, Converse A, Grethlein H, Lynd L: Kinetics of cellobiose hydrolysis using cellobiase composites from Trichoderma reesei and Aspergillus niger. Biotechnol Bioeng. 1985, 27: 463-470. 10.1002/bit.260270411.View ArticleGoogle Scholar
- Durand H, Tiraby G, Pourquié J: Génétique des micro-organismes industriels. 9. Ème Colloque De La Société Française De Microbiologie. 1984, 39-50.Google Scholar
- Pourquié J: Biochemistry and Genetics of cellulose degradation. FEMS Symposium No. 43. 1988, 43: 71-86.Google Scholar
- Laemmli UK: Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970, 227: 680-685. 10.1038/227680a0.View ArticleGoogle Scholar
- Lowry OH, Rosebrough NJ, Farr AL, Randall RJ: Protein measurement with folin phenol reagent. J Biol Chem. 1951, 193: 265-75.Google Scholar
- Himmel ME, Adney WS, Fox JW, Mitchell DJ, Baker JO: Isolation and characterization of 2 forms of beta-d-glucosidase from Aspergillus niger. Appl Biochem Biotechnol. 1993, 39: 213-225. 10.1007/BF02918991.View ArticleGoogle Scholar
- Medve J, Stahlberg J, Tjerneld F: Isotherms for adsorption of cellobiohydrolase I and II from Trichoderma reesei on microcrystalline cellulose. Appl Biochem Biotechnol. 1997, 66: 39-56. 10.1007/BF02788806.View ArticleGoogle Scholar
- Seidle HF, Marten I, Shoseyov O, Huber RE: Physical and kinetic properties of the family 3 beta-glucosidase from Aspergillus niger which is important for cellulose breakdown. Protein J. 2004, 23: 11-23. 10.1023/B:JOPC.0000016254.58189.2a.View ArticleGoogle Scholar
- Decker CH, Visser J, Schreier P: beta-glucosidases from five black Aspergillus species: Study of their physico-chemical and biocatalytic properties. J Agric Food Chem. 2000, 48: 4929-4936. 10.1021/jf000434d.View ArticleGoogle Scholar
- Calsavara LPV, De Moraes FF, Zanin GM: Modeling cellobiose hydrolysis with integrated kinetic models. Appl Biochem Biotechnol. 1999, 77-9: 789-806. 10.1385/ABAB:79:1-3:789.View ArticleGoogle Scholar
- Chen HZ, Hayn M, Esterbauer H: Purification and Characterization of 2 Extracellular Beta-Glucosidases from Trichoderma-Reesei. Biochim Biophys Acta. 1992, 1121: 54-60.View ArticleGoogle Scholar
- Chirico WJ, Brown RD: Purification and Characterization of a beta-glucosidase from Trichoderma reesei. Eur J Biochem. 1987, 165: 333-341. 10.1111/j.1432-1033.1987.tb11446.x.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.