Catalase improves saccharification of lignocellulose by reducing lytic polysaccharide monooxygenase-associated enzyme inactivation

Objectives Efficient enzymatic saccharification of plant cell wall material is key to industrial processing of agricultural and forestry waste such as straw and wood chips into fuels and chemicals. Results Saccharification assays were performed on steam-pretreated wheat straw under ambient and O2-deprived environments and in the absence and presence of a lytic polysaccharide monooxygenase (LPMO) and catalase. A kinetic model was used to calculate catalytic rate and first-order inactivation rate constants of the cellulases from reaction progress curves. The addition of a LPMO significantly (P < 0.01, Student’s T test) enhanced the rate of glucose release from 2.8 to 6.9 h−1 under ambient O2 conditions. However, this also significantly (P < 0.01, Student’s T test) increased the rate of inactivation of the enzyme mixture, thereby reducing the performance half-life from 65 to 35 h. Decreasing O2 levels or, strikingly, the addition of catalase significantly reduced (P < 0.01, Student’s T test) enzyme inactivation and, as a consequence, higher efficiency of the cellulolytic enzyme cocktail was achieved. Conclusion Oxidative inactivation of commercial cellulase mixtures is a significant factor influencing the overall saccharification efficiency and the addition of catalase can be used to protect these mixtures from inactivation. Electronic supplementary material The online version of this article (doi:10.1007/s10529-015-1989-8) contains supplementary material, which is available to authorized users.


Introduction
Cellulose is the planet's most abundant biopolymer and is found in plant cell walls in a strong network with other polysaccharides and the aromatic polymer lignin. Overcoming the recalcitrance of the plant cell wall (and of cellulose in particular) to depolymerisation remains a target for optimisation and effective use of lignocellulosic biomass as an industrial feedstock. As several full-scale factories producing lignocellulosic ethanol have come online recently around the world, the importance of guiding the industry towards optimal process conditions cannot be overestimated.
Various factors that negatively impact the hydrolysis efficiency of cellulase enzymes on pretreated plant cell wall material have been reported (reviewed by Yang et al. 2001). These include substrate-related factors, such as the presence of hemicellulose, cellulose crystallinity, degree of polymerization and increasing recalcitrance as a function of cellulose conversion, and enzyme-related factors, including product inhibition, inhibition from pretreatmentderived xylooligomers and phenols, inactivation and non-productive binding. However, the relative importance of these factors is insufficiently understood due in part to the variety of substrates, enzymes and assay conditions employed in these studies.
Lytic polysaccharide monooxygenases (LPMOs), such as AA9, carry out oxidative cleavage of glycosidic bonds and boost the activity of glycosyl hydrolases. However, in the absence of cellulose and in the presence of O 2 and a reducing agent, AA9 will release superoxide as a product of each redox cycle (Kjaergaard et al. 2014). Superoxide is spontaneously transformed into H 2 O 2 that reacts with transition metals such as iron and copper to form highly toxic hydroxyl radicals and other reactive oxygen species. Catalase (EC 1.11.1.6) fulfils a crucial role in biology as it terminates chains of radical chemistry by the dismutation of H 2 O 2 into water and O 2 .
The importance of LPMO and O 2 -associated inactivation of cellulase mixtures is shown here for the first time, in addition to the ability of catalase to mitigate this process.

Materials and methods
Unless otherwise stated, all reagents were laboratory grade. Steam-pretreated wheat straw was purchased at Lund University (kindly prepared by Mats Galbe). All commercial enzymes were obtained from Novozymes A/S. The aldose oxidase used in this study was cloned from Microdochium nivale (Genbank Accession BD103535) and expressed from Fusarium venenatum (Xu et al. 2001). Cloning and expression of the Thermoascus aurantiacus (Ta) AA9 (Accession ABW56451) was performed as previously described (Dotson et al. 2007). Catalase from Thermoascus aurantiacus (Accession DD046677) was expressed in Aspergillus niger and purified as described previously (Haruhiko and Sadaji 2004). The b-glucosidase from Aspergillus fumigatus (Af) (Accession EAL88289) was cloned and expressed from Aspergillus oryzae as described elsewhere (Teter et al. 2008).

Enzymatic saccharification of lignocellulose under low O 2 conditions
A disposable 280 L Atmosbag (Sigma Aldrich) was used for conducting saccharification experiments at 20 g scale under controlled O 2 levels. First, the level of dissolved O 2 in the reactions could be decreased to 1 % relative to the concentration in ambient air (as measured by Mettler Toledo O 2 electrode) by flushing the tent with N 2 and degassing all liquids prior to the experiment. These conditions caused control reactions with cellobiose and aldose oxidase to be significantly reduced in a 1 h experiment ( Supplementary Fig. 1). Duplicate samples of steam pretreated wheat straw (Lund University) were hydrolysed at 10 % (w/v) dry matter content (3.6 % initial cellulose) at around pH 5 at 50°C using free-fall mixing (Boekel Scientific, Big Shot III Hybridization oven model 230402) for up to 144 h.
Four doses of Cellic CTec3, 2.8, 5.6, 8.4 and 11.2 mg protein/g cellulose, were tested with and without addition of a constant dose of 0.22 mg Ta catalase/g cellulose. Half of the samples set were tested under standard conditions (without control of O 2 levels) and the other half was tested under a N 2 saturated atmosphere as described above. The pH was measured and adjusted daily and the glucose concentrations in samples withdrawn from the test tubes were determined. pH was adjusted using 1 M KOH and 75 mg Lactrol per kg slurry was added to prevent contamination by lactic acid bacteria.
The effect of LPMO activity was determined in a similar experimental set-up with enzyme blends consisting of 80 % Celluclast 1.5L, 10 % A. fumigatus b-glucosidase, and 10 % Thermoascus aurantiacus AA9 or BSA (96 % Sigma) by protein mass. The total protein dose in each experiment was 8.4 mg protein/g cellulose.
For the experiment with added H 2 O 2 , lids for the reaction tubes with a septum were used. A defined amount of H 2 O 2 was added to the reaction mixture using Hamilton syringes after about 18 h of incubation at which time the reaction mixtures were liquefied. H 2 O 2 was added four times at intervals of about 3.5 h during which the tubes were held at 50°C.

Determination of free sugar
The concentration of glucose was determined by HPLC analysis using an Aminex HPX 87H column, 5 mM H 2 SO 4 isocratic eluent and refractive index detection. Fractional cellulose conversion to glucose concentrations were calculated by dividing the glucose concentrations measured by HPLC at each time point by the initial concentration of glucose equivalents in the substrate.

Kinetic modelling
A two-stage kinetic model ( Fig. 1) was used to analyze cellulose hydrolysis progress curves generated under ambient and reduced O 2 conditions. The first stage in the model involves the conversion of cellulose to cellobiose (panel a). As described, all of the components in the commercial cellulases tested here, except for b-glucosidase(s), contributed to the formation of cellobiose. Similarly, only b-glucosidase, which accounted for either 5 or 10 % of the protein mass of the commercial enzyme preparations tested here, catalyzed the conversion of cellobiose to glucose (panel b). Both stages were assumed to follow Michaelis-Menten kinetics. It was also assumed that the cellulases in Stage 1 were subject to competitive inhibition by glucose (K G a ) and cellobiose (K G2 ) and subject to first-order inactivation (k i ) to the inactive form, E*. b-glucosidase, that was subject to competitive inhibition by glucose (K G b ). The mechanism of cellulose hydrolysis has been described using Michealis-Menten kinetics with product inhibition by others (reviewed by Sousa et al. 2011). This model was employed because it provided a suitable framework for interpreting the progress curves generated here under a limited range of substrate concentrations and enzyme doses. Additional complexities, such as enzyme adsorption onto cellulose (Scheiding et al. 1984, Philippidis et al. 1992, Nidetzky and Steiner 1993 and increasing recalcitrance of the substrate (Drissen et al. 2007) reported by others were not included in this model but are discussed below. The rate equations for the first (Eq. 1) and second (Eq. 2) stages of the model are shown below.
where, S is the substrate cellulose (g/l); t is time (h); k s is the cellulase catalytic rate constant (h -1 ); E a is the concentration of active cellulase (g/l); K s is the cellulase Michaelis-Menten constant (g/l); G is glucose (g/l); K G a is the competitive glucose inhibition constant for cellulase (g/l); G2 is cellobiose (g/l); K G2 is the competitive cellobiose inhibition constant for cellulase (g/l) where, k cat is the b-glucosidase catalytic rate constant (h -1 ); Bg is the concentration of b-glucosidase (g/l); K M is the b-glucosidase Michaelis-Menten constant (g/l); K G b is the competitive glucose inhibition constant for b-glucosidase (g/l) Inactivation of cellulases under the assay conditions described above was assumed to be first-order and is described in Eq. 3. A cellulase performance half-life (t 1/2 ) was calculated according to Eq. 4. bglucosidase activity was assumed to be stable.
where, k i is the cellulase inactivation rate constant (h -1 ) The differential equations were applied to each data set using a 4th order Runge-Kutta numerical integration using Microsoft Excel. The model was used to fit all doses of enzyme tested under a given set of experimental conditions by varying, k s and k i . Optimal values of k s and k i were determined simultaneously using the Excel Solver by minimization of least squres. Other parameters in the kinetic model were fixed to the values shown in Table 1.

Results and discussion
The effect of O 2 deprivation Steam-pretreated wheat straw was saccharified using the commercial enzyme product Cellic CTec3 at four different dosages and the fractional conversion of cellulose to glucose plotted as a function of time (Fig. 2). Progress curves of the data obtained for the ambient assay condition are shown in Fig. 2a. The initial rate of glucose release was, as expected, positively correlated with the enzyme dosage used and approximately 90 % cellulose conversion to glucose was observed by 144 h using the highest enzyme dose. It also appeared that the glucose release data over time resulting from each enzyme dose in Fig. 2a approached different asymptotic values, indicating that substrate depletion alone did not account for the decrease in glucose release over time. When the conversion data were plotted as a function of enzyme dose multiplied by time (Et) as shown in Fig. 2b, it was clear that the conversion data did not superimpose but rather deviated from a common Et relationship. This observation is consistent with time-dependent inactivation of the enzyme (Selwyn 1965) and kinetic modelling of these data indicated that the cellulase half-life was 39 h under these conditions (Table 2). When CTec3 was tested under limited O 2 conditions ( Fig. 2, panels e, f), the cellulase half-life increased significantly (P \ 0.01, Student's T test) to the extent that no significant level of enzyme inactivation was observed. Interestingly, the cellulase catalytic rate constant measured under O 2 -limited conditions   These results indicate that oxidative reactions harmful to the cellulolytic enzyme cocktail take place in the reaction mixture. Cellic CTec3 contains significant amounts of AA9 which requires O 2 for activity and previous reports have pointed to the need for O 2 to be available for optimal saccharification efficiency (Cannella andJorgensen 2014, Müller et al. 2015). In line with these reports, O 2 deprivation (Fig. 2e, f) significantly reduced the catalytic rate constant. However, a marked increase in enzyme half-life was also observed. The combined effect is that reactions conducted in ambient conditions saccharified faster; however, given sufficient incubation time, the highest glucose release was observed in the low O 2 regime ( Supplementary Fig. 2).
The rate of glucose production decreases markedly as cellulose hydrolysis proceeds and that this cannot be accounted for simply by substrate depletion. Product inhibition, substrate transformation, non-productive binding and enzyme inactivation have been proposed to account for this phenomenon. Importantly, the time-dependent enzyme inactivation observed above is distinct from conversion-dependent phenomena, which follow a common Et relationship. While such effects were not addressed in our kinetic model, any conversion-dependent phenomena, such as increasing substrate recalcitrance (Desai and Converse 1997), Drissen et al. 2007), would be captured in the cellulase product inhibition constant (K G a ), which was fixed for modelling purposes, along with true glucose inhibition effects. A Langmuir-like dependence of glucose production rate as a function of cellulase dose has been reported by others (Philippidis et al. 1993), Drissen et al. 2007). However, the observation here that progress curves measured under O 2 -limited conditions follow a common Et relationship indicates that the doses used did not reach the maximum binding capacity of the substrate. If this were the case, then a deviation from the Et relationship in which higher doses fall below the trend established by lower doses would be observed. In all likelihood, if even higher cellulase doses were tested, such a deviation in the Et relationship would be observed at some point. However, the absence of such an observation here supports the omission of a Langmuir-type dependency on cellulase catalytic rate from this model.

The effect of catalase addition
Catalase function is essential to control oxidative stress in living systems, where H 2 O 2 is a by-product of mitochondrial electron transport, b-oxidation of fatty acids and photorespiration (Kocabas et al. 2008). In addition to the biotic (enzymatic) generation of H 2 O 2 , this compound is also an unavoidable reactive oxygen species produced by abiotic (non-enzymatic) routes during aerobic degradation of organic matter (Remucal and Sedlak 2011). H 2 O 2 is either directly or indirectly harmful to enzymes and inhibits the activity of cellulases (Reese and Mandels 1980). The effect of adding the thermostable catalase from Thermoascus aurantiacus was therefore tested. Catalase 10.5 ± 0.6 a 0 ± 1.3 a n/a a The model was fitted to each progress curves shown in Fig. 2 by varying k s and k i . Values shown are best fit values and their associated standard deviations. Statistically significant differences relative to the ambient air and no catalase control, as determined using Student's T test, are indicated a P value \ 0.01 relative to value of equivalent parameter measured under ambient air conditions without catalase When catalase was added (Fig. 2 c and d) under ambient air conditions, the deviation in the Et relationship decreased substantially. Modelling of these data indicated that the addition of catalase significantly (P \ 0.05, Student's T test) increased the enzyme half-life to 75 h. No significant effect on the cellulase catalytic rate constant was observed, implying that catalase had no net effect on the rate of cellulase catalysis under these conditions. Addition of catalase under O 2 -deprived conditions (panel g and h) had no detectable effects.

LPMO activity is important for inactivation
In order to probe the impact of LPMO activity on the O 2 -and catalase-related inactivation observed above, enzymes mixtures with and without supplemental AA9 were tested. Celluclast, which is a wild-type Trichoderma reesei whole enzyme mixture, was used as an LPMO-poor cellulase cocktail. Mixtures were prepared containing 80 % Celluclast by protein mass and 10 % b-glucosidase to ensure rapid and complete conversion of cellobiose to glucose, thereby avoiding product inhibition of the cellobiohydrolases and simplifying product quantification. The remaining 10 % of the mix comprised either the Ta AA9 or an inert protein in the form of bovine serum albumin.
The progress curves for these experiments are shown in Fig. 3 and values for the kinetic parameters shown in Table 3. The k s and t 1/2 of the Celluclast ? b-glucosidase were 2.8 h -1 and 65 h, respectively, under ambient air (panel a). Addition of the AA9 increased the catalytic rate to 6.9 h -1 and resulted in an enzyme mixture with a lower t 1/2 (35 h), showing that the AA9 is closely associated with both the rate of catalysis and the rate of inactivation of the cellulase mixture. Inactivation of the cellulase mixture containing AA9 was reduced significantly (P value \0.01, Student's T test) when catalase was added while catalase addition, as before, had no significant effect on k s . Similarly, catalase addition had no effect when added simply to a blend of Celluclast ?b-glucosidase, indicating that catalase increases the stability of a cellulase mixture only when an AA9 is present.
Similar experiments done under O 2 -deprived conditions are shown in panel b. The k s and t 1/2 of the mixture without Ta AA9 and catalase under O 2deprived conditions (2.8 h -1 and 78 h, respectively) were similar to those measured under ambient conditions. Addition of AA9 under O 2 -limiting conditions increased the catalytic rate constant of the enzyme mixture from 2.8 to 4.6 h -1 . However, this k s was lower than that measured under ambient conditions (6.9 h -1 ) indicating that the reduced O 2 levels limited the catalytic activity of the AA9. No significant changes in t 1/2 were observed upon addition of AA9 or catalase under O 2 -limited conditions. This further demonstrates that the effects of O 2 here, both on the catalytic rate and inactivation of cellulase mixtures, are AA9-dependent.

Addition of H 2 O 2 hampers cellulolytic activity
The addition of T. aurantiacus catalase in the experiments above clearly reduced inactivation of a cellulase mixture containing AA9 under ambient O 2 conditions. Yet catalase addition had no apparent effect on catalytic rate, implying catalase interrupts O 2 -dependent radical reactions mediated by the AA9 that are damaging to one or more components of the enzyme mixture but that are distinct from the AA9 catalytic activity on cellulose. To better understand the effects H 2 O 2 may have in this system, CTec3 was incubated with pretreated wheat straw under both ambient and O 2 -deprived conditions with and without different doses of H 2 O 2 . These progress curves are shown in Fig. 4. The lowest dose, 0.01 mg H 2 O 2 /g slurry, did not have any impact on glucose release in the ambient condition [ Fig. 4a, blue (control) and red markers). In contrast, this dose of H 2 O 2 had a strong negative effect under the N 2 condition as it led to a reduction in conversion of approximately 10 % (Fig. 4c). At 0.01 mg H 2 O 2 /g slurry, the amount of H 2 O 2 added at each injection is only 1/10 of the concentration previously found to reduce the activity of Trichoderma cellulases by 50 % in assays using Avicel as the substrate (Reese and Mandels 1980). Still, the final yield of glucose is severely reduced to a level below the condition with ambient air ( Supplementary  Fig. 2). The addition of catalase fully prevented the loss of cellulolytic activity ( Supplementary Fig. 3).
The robustness of the system was tested by adding higher dosages of H 2 O 2 with or without catalase present (Fig. 4). In the absence of catalase, the addition of two doses of 1 mg H 2 O 2 /g had a catastrophic effect on glucose release regardless of the level of O 2 tested here (Fig. 4 a, b, black lines). The data show that the fixed concentration of catalase used in this study can detoxify both the 0.01 and the 0.1 mg added H 2 O 2 /g whereas the highest concentration 1 mg/g is too high for efficient detoxification in the O 2-deprived reactions ( Supplementary Fig. 3).
The formation reaction paths of H 2 O 2 are dependent on physical parameters such as temperature, dry matter content and on chemical parameters such as pH, concentration of transition metal and importantly the availability of O 2 . H 2 O 2 will spontaneously react with a large number of chemical species in the reaction mixture and the absolute concentration may therefore be relatively low compared to the flux. Indeed, it is highly likely that it is the hydroxyl radical derived from H 2 O 2 rather than H 2 O 2 itself that leads to oxidative inactivation of enzymes. A transient concentration in the range of 0.01 mg/g of slurry under O 2 deprivation seems to have a larger impact than what is seen in ambient air but more sophisticated analyses are required to determine the O 2 -dependent flux in H 2 O 2 . In conclusion, oxidative inactivation of commercial cellulase mixtures is a significant factor influencing the overall saccharification efficiency and the addition of catalase can successfully be used to protect cellulases mixtures from inactivation.
Acknowledgments This work was supported by a grant from the ERA-IB Initiative, which for the Danish partners was funded by the Danish Council for Strategic Research (Grant No. 12-134922). Lars Østergård, Novozymes A/S is thanked for good discussions and advice and Dr. Hai Yu Ren is thanked for her contributions to identifying the beneficial effect of catalase on biomass saccharification.
Funding KSJ is from the European Union's Seventh Framework Programme for research, technological development and demonstration under grant agreement no. 608473. This paper reflects only the author's views and the Union is not liable for any use that may be made of the information contained herein. Supplementary Fig. 1-The reduction in cellobiose concentration through oxidation by aldose oxidase either in ambient air or under *1 % O 2 tension relative to ambient air. Supplementary Fig. 2 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. CTec3 was incubated with pretreated wheat straw in the presence of 0 (blue circles), 0.01 (red triangles), 0.1 (green squares) and 1 (black diamonds) mg H 2 O 2 /g slurry added four times (as indicated by the arrows) after the cellulosic material was liquefied. Effects of H 2 O 2 were tested under ambient air (panel a) and after purging with N 2 (panel b). For this experiment lids for the reaction tubes with a septum were used and the indicated amount of H 2 O 2 was added to the reaction mixture using Hamilton syringes. The first addition was made after about 18 h of incubation at which time the reaction mixtures were liquefied. H 2 O 2 was added four times at intervals of about 3.5 h during which the tubes continued incubation at 50°C