Original Research ARTICLE
Nonlinear temperature sensitivity of enzyme kinetics explains canceling effect—a case study on loamy haplic Luvisol
- 1Department of Agricultural Soil Science, University of Göttingen, Göttingen, Germany
- 2Department of Soil Science of Temperate Ecosystems, University of Göttingen, Göttingen, Germany
- 3Institute of Physicochemical and Biological Problem in Soil Science, Russian Academy of Science, Pushchino, Russia
The temperature sensitivity of enzymes responsible for organic matter decomposition in soil is crucial for predicting the effects of global warming on the carbon cycle and sequestration. We tested the hypothesis that differences in temperature sensitivity of enzyme kinetic parameters Vmax and Km will lead to a canceling effect: strong reduction of temperature response of catalytic reactions. Short-term temperature response of Vmax and Km of three hydrolytic enzymes responsible for decomposition of cellulose (β-glucosidase, cellobiohydrolase) and hemicelluloses (xylanase) were analyzed in situ from 0 to 40°C. The apparent activation energy varied between enzymes from 20.7 to 35.2 kJ mol−1 corresponding to the Q10 values of the enzyme activities of 1.4–1.9 (with Vmax-Q10 1.0–2.5 and Km-Q10 0.94–2.3). Temperature response of all tested enzymes fitted well to the Arrhenius equation. Despite that, the fitting of Arrhenius model revealed the non-linear increase of two cellulolytic enzymes activities with two distinct thresholds at 10–15°C and 25–30°C, which were less pronounced for xylanase. The nonlinearity between 10 and 15°C was explained by 30–80% increase in Vmax. At 25–30°C, however, the abrupt decrease of enzyme-substrate affinity was responsible for non-linear increase of enzyme activities. Our study is the first demonstrating nonlinear response of Vmax and Km to temperature causing canceling effect, which was most strongly pronounced at low substrate concentrations and at temperatures above 15°C. Under cold climate, however, the regulation of hydrolytic activity by canceling in response to warming is negligible because canceling was never observed below 10°C. The canceling, therefore, can be considered as natural mechanism reducing the effects of global warming on decomposition of soil organics at moderate temperatures. The non-linearity of enzyme responses to warming and the respective thresholds should therefore be investigated for other enzymes, and incorporated into Earth system models to improve the predictions at regional and global levels.
The temperature sensitivity of soil organic matter (SOM) decomposition has attracted significant interest because of its importance in the global carbon cycle (C) and potential feedbacks to global warming (Schimel, 1995; Davidson and Janssens, 2006; De Graaff et al., 2014). Temperature sensitivity of reaction rates (e.g., SOM decomposition) is commonly described by the Arrhenius equation. This is based on the energy required to initiate the reaction, termed the activation energy (Ea), (Von Lützow and Kögel-Knabner, 2009). For chemical reactions, the Arrhenius equation predicts exponentially increasing reaction rates with increasing temperature, assuming constant values of activation energy (Kirschbaum, 2000; Craine et al., 2010; Craine and Gelderman, 2011). In nature, however, the decomposition of SOM is mediated by extracellular enzymes, produced by microorganisms (Allison et al., 2010; Glanville et al., 2012; Zimmermann and Bird, 2012; Van Gestel et al., 2013). Therefore, deviations from Arrhenius behavior can occur as a consequence of temperature sensitivity of enzyme systems, through enzyme denaturation and proteolysis, for example (Bennett et al., 2008; Maire et al., 2013; Goyal et al., 2014) or by temperature-accelerated desorption of immobilized enzymes (Nannipieri et al., 1996; Nielsen et al., 2006). Besides that, changes in the temperature dependency of microbial communities may cause expression of various set of isoenzymes (i.e., an enzyme with the same function but a different structure) or changes in enzyme conformation (Bradford, 2013). This, in turn, affects apparent activation energy and temperature sensitivity of enzyme-mediated reactions in situ in soil.
Besides enzyme properties, the rate of enzyme-mediated reactions is dependent on substrate amount, which can strongly vary in time and space considering heterogeneity of soil microhabitats (Cheng et al., 1996; Ekschmitt et al., 2005; Ruamps et al., 2013). Therefore, not only the Arrhenius equation, but the relationship between rate of enzyme-mediated reaction and substrate needs to be considered in the context of temperature sensitivity of SOM decomposition. The rate of catalytic reactions is a saturating function of substrate concentration and is described by the Michaelis-Menten relationship (Michaelis and Menten, 1913). Both parameters of the Michaelis-Menten equation—maximum enzyme activity (Vmax) and the half-saturation constant (Km), which is an intrinsic feature of an enzyme system reflecting the affinity of the enzyme for the substrate (Khalili et al., 2011; Bradford, 2013)—are temperature sensitive (Davidson and Janssens, 2006; Davidson et al., 2006; Stone et al., 2012). As both Vmax and Km values (respectively, in the numerator and denominator of the Michaelis-Menten equation) usually increase with temperature (Stone et al., 2012), a canceling effect (absence or strong reduction of response of the enzyme to temperature) can occur for the resulting enzymatic reaction rate (Berry and Raison, 1981; Atkin and Tjoelker, 2003; Davidson et al., 2006). The canceling effect may be pronounced when the substrate concentration is lower than or close to Km and if both Km and Vmax have similar temperature sensitivities (Larionova et al., 2007; Gershenson et al., 2009). Ecological and evolutionary processes in microbial communities could also reduce the temperature sensitivity of enzymes by reducing Vmax and increasing Km (Hochachka and Somero, 2002). It should be noted that response of pure and isolated enzyme are strongly different compared to the soil enzyme, for instance, the catalytic properties of the soil enzymes (e.g., substrate affinity) are much higher compared to the bacterial and fungal enzymes originated from pure cultures (Skálová et al., 2005; Tischer et al., 2015).
Relative temperature responses are commonly compared by a Q10 index—the change of a reaction rate with a temperature increase of 10°C (Birgander et al., 2013). Q10 values of 2–3 are commonly assumed for respiration rates in the temperature interval 10–20°C (Davidson and Janssens, 2006). However, enzyme activities have been found to be less temperature sensitive than compared to respiration, with Q10 values < 2 (Browman and Tabatabai, 1978; Tabatabai, 1994; Koch et al., 2007), and lower Q10 values are often observed at higher temperatures (Tjoelker et al., 2001; Xu and Qi, 2001). Recently, a cross-latitudinal study for enzymes involved in the C cycle demonstrated Vmax-Q10 values ranging from 1.5 to 2.3 and Km-Q10 values ranging from 0.90 to 1.9 (German et al., 2012). When Vmax and Km cancel each other out (i.e., the Q10 of catalytic reaction ~1), decomposition is restricted by temperature sensitivity of a bottle-neck process that produces available substrate, e.g., by decomposition of recalcitrant or stabilized SOM (Ågren and Wetterstedt, 2007). The canceling effect is usually more pronounced when substrate concentration is lower or close to Km and if both Km and Vmax have similar temperature sensitivities (Larionova et al., 2007; Gershenson et al., 2009). Accordingly, the canceling effect can be an important factor controlling the “actual” temperature sensitivity of organics decomposition in soils (Von Lützow and Kögel-Knabner, 2009). Despite theoretical predictions (Davidson and Janssens, 2006), there is still a lack of experimental data on the occurrence of canceling as dependent on temperature range, substrate amount and enzyme.
This study was designed to test how the canceling effect of three enzymes involved in the C cycle changes with temperature. We hypothesized that the temperature sensitivity of Vmax and Km differ both: (1) in the specific Q10, and (2) in the temperature ranges with maximal changes of enzyme activity and substrate affinity. The consequence of this hypothesis is that the increase in the rate of enzymatic reaction with temperature is not continuous, but may have some thresholds with a much stronger response. Using the modified Arrhenius equation we tested the response of enzyme systems to temperature, e.g., thermal denaturation or changes in enzyme affinity to substrate.
The modern global change models simulate the responses of soil C pools to about 4.8°C increase over the next 100 years (IPCC, 2007; Wieder et al., 2013). Because the main C input in soil is ongoing by above and belowground litter (Hasibeder et al., 2014), it is crucial to understand how enzymes involved in litter decomposition will respond to warming to better predict the links between C input, SOM formation and global warming (German et al., 2011). Thus, soil (haplic Luvisol) was incubated over a temperature range of 0–40°C (with 5° steps) during 5 days to determine the activities of three enzymes: β-glucosidase and cellobiohydrolase, which are responsible for degrading cellulose (Mganga et al., 2015; Wang et al., 2015), and xylanase, which degrades xylooligosaccharides (short xylene chain) into xylose and is thus responsible for breaking down hemicelluloses (Chen et al., 2012).
Material and Methods
Soil samples were taken from the top 10 cm of the Ap horizon of an arable loamy haplic Luvisol located on a terrace plain of the Leine River north-west of the city of Göttingen (Holtensen, Germany) with mean annual temperature 5–15°C. The properties of the soil were: pH 6.5; 5.8% sand, 87.2 % silt, 5.8 % clay; 12.6 g kg−1 C, 1.3 g kg−1 N, 1.4 g cm−3 bulk density, and 30 % water field capacity. The samples were kept cold (~4°C) during transportation to the laboratory. Then the samples were frozen at −20°C until pre-incubation.
Temperature Incubation Experiment
Three enzymes involved in the C cycle were investigated over a temperature range of 0, 5, 10, 15, 20, 25, 30, 35, and 40°C. Nine climate chambers (SBS C120) were used to regulate the temperature with a precision of < 1°C. At each temperature four replicates were incubated. The frozen samples were thawed at 4°C for 1 day and then pre-incubated at 20°C for 14 days before the start of 5 days incubation. Because freezing may affect enzymatic activities (Lee et al., 2007; Stone et al., 2012), all samples were frozen similarly, and they were pre-conditioned after thawing. We therefore assume that this pretreatment corresponded to snow thaw in spring and that the freezing effect was negligible after the pre-incubation.
The kinetics of enzyme activities was assayed using fluorogenically labeled substrates. Three types of fluorogenic substrates based on 4-methylumbelliferone (MUF) were used to assess enzymatic activities; 4-methylumbelliferyl-β-D-glucoside (MUF-G) to detect β-glucosidase activity; 4-methylumbelliferyl-β-D-cellobioside (MUF-C) to detect cellobiohydrolase activity; and 4-methylumbelliferyl-β-D-xylopyranoside (MUF-X) to detect xylanase activity. All substrates and chemicals were purchased from Sigma (Germany).
MUF-substrates were dissolved in 2-methoxyethanol (Hoppe, 1983). Pre-dissolved MUF-substrates were further diluted with sterile universal buffer [MES (C6H13NO4SNa0.5)] to obtain the desired concentrations; we determined enzyme activities in a range of substrate concentrations from low to high (0, 10, 20, 30, 40, 50, 100, 200 μmol g−1soil). Saturation concentrations of fluorogenic substrates were determined in preliminary experiments. Suspensions of 0.5 g soil (dry weight equivalent) with 50 mL water were prepared separately for each of 4 incubated replicates using low-energy sonication (40 J s−1 output energy) for 2 min (Stemmer et al., 1998; Koch et al., 2007). Fifty microliters of soil suspension was added to 50 μL buffer (pH:6.5) and 100 μL of each substrate solution in a 96-well microplate (Puregrade, Germany). The empty microplates as well as working solutions were pre-incubated at the exact temperature of the assay before the experiment started. The calibration solutions were prepared using soil suspension (50 μL) and MUF to obtain a series of concentrations 0–1.2 mM (Ali et al., 2015). Fluorescence was measured in microplates at an excitation wavelength of 355 nm and an emission wavelength of 460 nm, slit width of 25 nm, with a Victor3 1420-050 multi label counter (Perkin Elmer, USA). Activity of each enzyme was determined in each soil sample at exact temperature during 3 h. No enzyme assay took longer than 2 min. After each measurement (i.e., after 30 min, 1 h, 2 h, and 3 h) the microplates were promptly returned to the climate chambers.
Enzyme activities were expressed as MUF release in nmol per g dry soil per hour (nmol g−1 h−1). In addition for all four incubation replicates, the assay of each enzyme at each substrate concentration was performed in three analytical replicates (12 wells in the microplate). Due to high-temperature effects, the chemical reactions of the substrates are considerable. As thermal hydrolysis of MUF-phosphate could occur at temperatures of 40–65°C (Spohn and Kuzyakov, 2013), we checked possible temperature effects on the chemical decomposition of the three MUF-substrates, but no significant effects were detected in the range 0–40°C.
Enzyme Kinetics, Temperature Sensitivity, and Statistical Analyses
The Michaelis-Menten equation used to determine parameters of enzyme activity (V),
Where, Vmax is the maximal rate of enzymatic activity at a given temperature; Km is the half saturation constant, or the substrate concentration at ½Vmax, and [S] is the concentration of the substrate (Michaelis and Menten, 1913; Segel, 1975). All parameters were modeled with the non-linear regression routine of STATISTICA. The t-test (pairwise differences) was applied to distinguish the significant differences for each neighboring pair of independent variables (mean values of Vmax at 0, 5, 10, 15, 20, 25, 30, 35, 40°C), then the temperatures with insignificant differences in Vmax were selected in one temperature-range group (Boone et al., 1998; Melillo et al., 2002). Thus, based on the t-test, the enzyme activities were not significantly different within each categorical group, but they were significantly different from the other groups. Thereafter, One-way analysis of variance (ANOVA) at a probability level of p < 0.05, was applied to distinguish the significant differences in Km values between categorical temperature-range groups (Table S1). Homogeneity of variance and normality of the values was tested by the Leven's test and Shapiro-wilk's W-test.
We used the conventional Q10 function (2) to examine variation in temperature sensitivity and express temperature responses of each enzyme kinetic parameter (i.e., Km or Vmax).
where, R is the rate of a process or a value of kinetic parameters and T is temperature. For unification and comparison of canceling effects at 5and 10°C increment, we used function (3):
The activation energy was calculated according to the classical Arrhenius equation (Equation 4):
where, k is the reaction rate constant; A is the frequency of molecular collisions; Ea is the required activation energy in Joules per mole; R is the gas constant (8.314 J mol−1 K−1) and T is the temperature in Kelvin. In addition to the classical Arrhenius equation (Equation 4), the modified Arrhenius equation (Equation 5) was applied to estimate possible effects of enzyme denaturation on apparent Ea that contributed to enzymes reaction rate (Flickinger and Drew, 1999).
Here, ΔS°and ΔH° are the entropy and enthalpy of denaturation, respectively, and υ is the observed rate of biological processes. Activity of β-glucosidase demonstrated similar temperature sensitivity to other enzymes (see Results) but the strongest pattern of both grouping for Vmax and of non-linearity in temperature response for Km. We assumed, therefore, that any effect of thermal denaturation would be strongest for β-glucosidase as compared with the other tested enzymes. Hence, we applied the modified Arrhenius equation (Equation 7) to β-glucosidase activity, assuming intercept and absolute values of ΔS° and ΔH° determined in previous studies, using differential scanning calorimetry: ΔS° = 86 J mol−1 K−1 and ΔH° = 108 kJ mol−1, (Zoldák et al., 2004; Goyal et al., 2014).
The parameters A and Ea of Equations (4, 5) were fitted by the Marquardt minimization method using ModelMaker version 3.0.3 software (ModelMaker, 1997). The fitting quality and model efficiency were compared by the coefficient of determination (R2)
Enzyme Responses to Temperature
The activities of the three enzymes responded positively to increasing temperature with widest range of activity for β-glucosidase (from 69.9 to 362.8) and the narrowest one for xylanase (from 9.2 to 14.5) (Figure 1). The Q10 values for reaction rates were always >1, with the average range of 1.4–1.9. The Michaelis-Menten kinetics (enzyme activity as a function of substrate concentration) demonstrated stepwise increases in response to temperature across the whole range of substrate concentrations, i.e., 10–200 μmol g−1 (Figure 1). The temperature sensitivity of Vmax and Km decreased with increasing temperature and varied between enzymes, corresponding to Vmax-Q10 values of 1.0–2.5 and Km-Q10 values of 0.94–2.3 (Table 1).
Figure 1. Michaelis-Menten kinetics of β-glucosidase (top), Cellobiohydrolase (middle), and Xylanase (bottom) as a response to temperature increasing from 0 to 40°C with 5° steps. Values are means of four replications (± SE). Curves present fitting of Michaelis-Menten kinetics by non-linear regression. The fitted Vmax and Km values are presented in Table 1. Grouping (I–III) was formed using the t-test based on significant differences between neighboring Vmax values.
Table 1. The Q10 values (± SE of four replicates) of Vmax and Km of three hydrolytic enzymes measured at nine temperatures.
Temperature sensitivity of all tested enzymes fit well to the Arrhenius equation (Figure 2) with the Ea values ranged from 20.7 to 35.3 kJ mol−1 (Table 2). However, the fitting of Vmax for all enzymes, especially for β-glucosidase and cellobiohydrolase, to the classical Arrhenius model (Eq. 4) demonstrated deviations within the temperature range 15–30°C (Figure 2, green solid line). The pairwise comparison t-test considering temperature increase with 5 degree resolution revealed three groups of stepwise increase of enzyme activities (Figure 1). The accuracy of such a grouping, however, was enzyme-dependent and was less pronounced for enzymes with overall low temperature sensitivity, e.g., xylanase activity differed between lowest and highest temperature (40°C temperature increase) only by a factor of three (Figure 1). Accordingly, the “jumping” pattern was more obvious for β-glucosidase compared to the two other enzymes.
Figure 2. Fitting of Arrhenius model to experimental Vmax values of β-glucosidase kinetics (top), cellobiohydrolase, and xylanase (bottom), as a function of temperature from 0 to 40°C (green solid line). The fitting quality was not improved by application of modified Arrhenius model (red dash line and Table S2).
Table 2. The varying Ea-values which were obtained by fitting the Arrhenius equation for three hydrolytic enzymes studied.
To prove whether the deviations from Arrhenius function can be explained by enzyme denaturation, we applied the modified Arrhenius equation (5) considering thermal degradation of the enzymes. Fitting to the modified Arrhenius equation (5) demonstrated negligible decreases in catalytic activity, implying an absence of denaturation. Similar Ea values were obtained for the whole temperature range with Equations (5, 6) (Table S2). Furthermore, fitting to Equation (5) did not reduce the deviations of experimental data from the model predictions, suggesting that the revealed grouping (or jumping effect) cannot be explained by thermal degradation of the enzymes. Statistical index R2 demonstrated similar model efficiency when classical Arrhenius equation were used for fitting as compared with modified Arrhenius equation (Table S2).
We further checked whether the jumping in temperature response of catalytic reactions could be explained by changes in temperature sensitivity of Km. However, the Km values were relatively constant without significant differences between 0 and 25°C. Thus, the low-temperature “threshold” was due to the sudden 30% increase in Vmax over a temperature range of 10–15°C and was unaffected by Km values (Table 1 and Figure 3). The temperature effect on Km, however, revealed distinct threshold with a significant increase in the affinity of all enzymes to the substrate at temperatures above 30°C (Table 1, Figure 3).
Figure 3. The Km of β-glucosidase, cellobiohydrolase and xylanase plotted as a function of temperature. Significant differences of Km values from 25 to 30°C for all enzymes were confirmed by ANOVA (p < 0.05), (See Table S1).
The canceling effect was pronounced at a temperature increase by 5°C and it was enzyme-specific (Figure 4). Generally, the temperature ranges with observable canceling were wider at substrate concentrations below 50 μmol g−1soil (Figures 1, 5). At substrate concentrations above 50 μmol g−1soil, canceling was detected within narrow temperature ranges (25–35°C for cellulose-degrading enzymes and 5–25°C for hemicellulolytic enzymes). Independent of substrate concentration, canceling was never observed at temperatures below 10°C or above 35°C (Figure 5).
Figure 4. Presence of the canceling effect in the area below the dashed red line (1.2). Symbols represent the values of dividing the reaction rates with 5–10°C resolution for β-glucosidase (top), cellobiohydrolase (middle), xylanase (bottom). Values of Q10 ≤ 1 indicate full canceling (below blue line) and all values exceeding 1.2, showing the absence of a canceling (e.g., empty blue symbols). Canceling effect is obvious at moderate (15–25°C) and elevated (25–40°C) temperature ranges and mainly at lower substrate concentrations. The area below the dashed red line was arbitrarily accepted as relevant presence of canceling effect if the relative increase of reaction rates is below 20% (means = 1.20 on the Y axis). This threshold considers the uncertainties of the data: the mean SE of each relative unit at the exact concentration of substrate for each enzyme. The error bars were deliberately omitted for clarity.
Figure 5. Canceling effect as a function of temperature and substrate concentration. Independent of substrate concentration, canceling was never observed at temperatures below 10°C. Canceling effect significantly increased at the intermediate and elevated temperature ranges, especially at low substrate concentration. The canceling effect occurred at warm temperatures at all substrate concentrations for cellobiohydrolase and xylanase.
The canceling effect for 10°C temperature resolution was especially pronounced at low substrate concentration (10–50 μmol g−1soil, Figure 5). With one exception for β-glucosidase: at the increase from 15 to 25°C, a canceling effect was detected at the large concentration range of 10–100 μmol g−1 soil (Figure 4). Thus, pronounced canceling effects were evident over temperature increments of both 5 and 10°C, but at different thresholds of substrate concentration. With a 10°C increment, the canceling effect was significant only at substrate concentrations below 25 μmol g−1 soil. With 5°C resolution, however, the canceling was observed even at saturating substrate concentrations.
Enzyme Responses to Temperature
The Q10 values of the enzyme activities varied from 1.5 to 1.9 within the low temperature range and decreased to 1.4 with increasing temperature. This decrease is in accordance with theoretical predictions (Davidson and Janssens, 2006) and experimental observations of reduced Q10 values at increased temperature (Xu and Qi, 2001). Our results were in line with common range of Ea for the substrates of three tested enzymes, which varied from 30 to 57 kJ mol−1 (Blagodatskaya et al., 2014). Temperature sensitivity of Vmax and Km for the three tested enzymes, with the highest Vmax-Q10 = 2.5 for cellobiohydrolase and the highest Km-Q10 = 2.3 for xylanse (Table 1), was in line with studies of German et al. (2012) and Stone et al. (2012). Notably, higher temperature sensitivity of Km for xylanase suggests that hemicellulose degradation may be more affected by temperature than the other components of plant residues (Stone et al., 2012).
Despite satisfactory approximation of enzyme activities by the Arrhenius equation the fitting demonstrated remarkable deviations within the temperature range 15–30°C. Such deviations were related to the temperature sensitivity of the enzymes and were more pronounced for enzymes with wider range of temperature response (β-glucosidase, cellobiohydolase). Deviations of experimental data from both classic and modified Arrhenius models have been observed previously for parameters of microbial activity, e.g., for microbial respiration and growth (Ratkowsky et al., 1983; Pietikäinen et al., 2005; Sand-Jensen et al., 2007). Furthermore, an inconsistent response of enzyme activity to gradual temperature increase was also observed (but not discussed) for β–glucosidase in soils of Alaska and Costa Rica (Figure 1 in German et al., 2012). Such deviations from Arrhenius equation have been partly explained by protein denaturation at temperatures above 35°C (Daniel et al., 1996; Goyal et al., 2014) or by experimental artifacts arising from changes of substrate pools over time that are more significant at higher temperatures (Sand-Jensen et al., 2007). In our study, however, the fitting quality was not improved when denaturation was considered by the modified Arrhenius model, suggesting a negligible effect of thermal denaturation in the temperature range 0–40°C (Goyal et al., 2014). Besides that, enzymes immobilization by clays or entrapment by organic matter (Allison, 2006; Burns et al., 2013) indirectly alter enzymes kinetics (Nannipieri and Gianfreda, 1998; Nannipieri et al., 2011) and increase stability against thermal denaturation and proteolysis (Burns, 1982; Stotzky, 1986; Sarkar et al., 1989; Nannipieri et al., 1996; Nielsen et al., 2006; Maire et al., 2013). Since, thermal denaturation of enzymes cannot explain the deviations from Arrhenius function in our study, other factors, such as, enzyme immobilization, changes in microbial community structure or shifts in community function (resulting in expression of iso-enzymes) altering temperature sensitivity of the enzyme kinetic parameters (Vmax and Km), must be responsible for this behavior (Davidson and Janssens, 2006).
The parameters of Michaelis-Menten kinetics responded differently to temperature increase at low and high temperature thresholds because the sensitivities of Vmax and Km do not always correspond to each other (Stone et al., 2012; Tischer et al., 2015). The Vmax values raised by 30–100% at the temperature threshold between 10 and 15°C. Such an increase in Vmax was not accompanied by significant changes in the Km values. Thus, the enzyme activity (Vmax) strongly increased, but the enzyme systems (Km) themselves remained unaffected. The absence of significant trends in enzyme affinities (Km) within the range of 0–25°C shows that microorganisms in the studied Luvisol at temperate climate are adapted to diurnal temperature changes and produce a similar enzyme sets below 25°C. The deviations from Arrhenius behavior at 10–15°C therefore cannot be explained by the temperature sensitivity of Km. Thus, the raised enzyme activity between 10 and 15°C, could indicate the increase in enzyme production due to the shift in domination of mesophilic microbial group with optimal growth temperature at 15–25°C. Such microorganisms can occur and produce similar extracellular enzymes at both psychrophilic and mesophilic temperature ranges. However, their activity and enzymes production increases under optimal temperature resulting in local non-linear temperature response. This was supported by the absence of a canceling effect at low temperature (0–10°C) and by its occurrence at temperatures above 15°C indicating weaker temperature sensitivity of mesophilic microorganisms. A two-fold increase of Km between 25 and 30°C indicated strong reduction of enzyme affinities to the substrate due to production of iso-enzymes (Hochachka and Somero, 2002). Thus, “jumping” temperature sensitivity of soil enzymes at 30°C is related to the shift toward temperature-dependent physiological groups producing enzyme isoforms with similar functions, but different temperature sensitivities (Khalili et al., 2011; Bradford, 2013). The uncommon high temperatures affect not only the active microorganisms, but also activate dormant microorganisms (Birgander et al., 2013) expressing the enzymes beneficial in thermophilic range (Hochachka and Somero, 2002; Pepper et al., 2011; Yumoto, 2013). Therefore, the major shift in the temperature response of soil microbial activity above 30°C was expected, because these temperatures are uncommon under a temperate climate with annual soil temperature of 8–10°C (Bárcenas-Moreno et al., 2009).
The observed temperature thresholds with strong temperature sensitivity were clearly obvious for β–glucosidase, the enzyme responsible for decomposition of cellulose. These thresholds were also observable for other enzymes involved in C cycle responsible for the decomposition of cellulose and hemicelluloses. However, we expect that the thresholds, canceling effects and mechanisms of sensitivity may be different for the enzymes of other cycles: N, P, S, as well as the broad group of oxidative enzymes, like peroxidases and phenoloxidases. Furthermore, to generalize the conclusions based on one soil type from temperate climate, more soils from various climate zones need to be tested. Therefore, we need more mechanistic work, in situ studies along with the studies of pure and isolated enzymes from a range of habitats to verify the findings on temperature sensitivity of hydrolytic reactions and to relate them with functional and phylogenic structure of microbial community as well as with temperature responses of specific proteins.
The canceling effect was not significant at lower temperatures (0–10°C), but its significance increased at 15–25°C for all enzymes tested (Figure 6). This finding supports the theoretical predictions that Vmax and Km can cancel each other upon temperature increase (Davidson et al., 2006) especially at the intermediate and elevated temperatures that favors mesophilic microorganisms. Thus, reduced temperature sensitivity of hydrolytic reactions and a smoothed response to warming could be predicted within the mesophilic temperature range for similar soils under temperate climate. Canceling occurrence at moderate temperatures was accompanied by deviations of reaction rates from Arrhenius function (Figure 6). This was mainly due to strong increase of Q10−Vmax. Next deviation from Arrhenius function at 25–30°C was mainly explained by increase of Q10−Km. The occurrence of a canceling effect for cellobiohydrolase and xylanase at warm temperatures (30–35°C) at all substrate concentrations suggested slower decomposition of the main components of plant residues (cellulose and hemicelluloses) than predicted by the Arrhenius equation, even at increased net primary production (Figure 6). As a result of canceling, the decomposition of recalcitrant C (Ågren and Wetterstedt, 2007; Larionova et al., 2007; Gershenson et al., 2009) and specifically of cellulose and hemicellulose (based on our study) will be slowed down. The canceling therefore can be considered a natural mechanism reducing the consequences of global warming for microbial decomposition of soil organics predicted in temperate ecosystems (Tang and Riley, 2014). Under cold climate, however, the regulation of hydrolytic activity by canceling in response to warming could be of minor relevance given that canceling was never observed at temperatures below 10°C (Figure 5). As this conclusion has been done for soil from temperate ecosystem, the relevance of canceling mechanism responsible for weaker soil carbon-climate feed-backs needs to be proven for the soils with contrasting properties (e.g., texture, structure, pH, C content, etc.) in the range of climate zones, e.g., in boreal and tropical environments. Consideration of sensitivity jumping and temperature thresholds can significantly improve the predictions of Earth system models for C cycling at regional and global levels, especially if the response of broad range of functional enzymes is evaluated. To our knowledge, this is the first study explored the canceling effect across a complete range of substrate concentrations in soil for the set of hydrolytic enzymes at a temperature resolution of 5°C increments (nine temperature levels). This study confirms previous studies on “apparent” Km (Tabatabai, 1994; Nannipieri and Gianfreda, 1998; Marx et al., 2005) and provides a basis for implementing both Vmax and Km in models of soil organic matter decomposition under global warming (German et al., 2012). In particular, because Km and Vmax vary independently as temperature changes (Stone et al., 2012), both parameters of the Michaelis-Menten kinetics need to be considered in global earth system models, especially at temperature thresholds for physiological groups of soil microorganisms.
Figure 6. Occurrence and potential outcome of “a canceling effect” in soil of temperate ecosystem as affected by temperature and substrate concentration. Canceling effect is more pronounced at warm temperatures (red arrows) and at low (thick arrow) than at high (thin arrow) substrate concentration. Canceling does not occur at low temperatures (dotted red line). Canceling occurrence at moderate temperatures was accompanied by deviations of reaction rates (green circles) from Arrhenius function (green line). This was mainly due to strong increase of QVmax10. Next deviation from Arrhenius function at warm temperatures was mainly explained by increase of QKm10. As a result of canceling, the decomposition of: (1) recalcitrant C (Ågren and Wetterstedt, 2007; Larionova et al., 2007; Gershenson et al., 2009) and (2) cellulose and hemicellulose (based on our study) will be slowed down. Therefore, at low substrate concentrations and warm temperatures, canceling effect can be considered as natural mechanism mitigating accelerated SOM decomposition as a consequence of global warming.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors wish to thank Ingrit Ostermeyer for support in the laboratory. We gratefully acknowledge the German Academic Exchange Service (DAAD) for supporting BR. The contribution of EB was supported by the Russian Science Foundation (project o 14-14-00625). The authors thank referees and the editor for constructive comments. This work was supported by Georg-August University of Goettingen.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2015.01126
Ali, R. S., Ingwersen, J., Demyan, M. S., Funkuin, Y. N., Wizemann, H. D., Kandeler, E., et al. (2015). Modelling in situ activities of enzymes as a tool to explain seasonal variation of soil respiration from agro-ecosystems. Soil Biol. Biochem. 81, 291–303. doi: 10.1016/j.soilbio.2014.12.001
Bárcenas-Moreno, G., Gómez-Brandón, M., Rousk, J., and Bååth, E. (2009). Adaption of soil microbial communities to temperature: comparison of fungi and bacteria in a laboratory experiment. Glob. Change Biol. 15, 2950–2957 doi: 10.1111/j.1365-2486.2009.01882.x
Bennett, M. P., Drake, C., and Mitchell. (2008). Regulation of membrane proteins by dietary lipids: effects of cholesterol and docosahexaenoic acid acyl chain-containing phospholipids on rhodopsin stability and function. Biophys. J. 95, 1206–1216. doi: 10.1529/biophysj.107.122788
Berry, J. A., and Raison, J. K. (1981). “Responses of macrophytes to temperature,” in Physiological Plant Ecology I. Responses to the Physical Environment, eds O. L. Lange, P. S. Nobel, C. B. Osmond, and H. Ziegler (Berlin: Springer-Verlag), 277–388.
Birgander, J., Reischke, S., Jones, D. L., and Rousk, J. (2013). Temperature adaptation of bacterial growth and 14C-glucose mineralisation in a laboratory study. Soil Biol. Biochem. 65, 294–303. doi: 10.1016/j.soilbio.2013.06.006
Blagodatskaya, Å., Zheng, X., Blagodatsky, S., Wiegl, R., Dannenmann, M., Butterbach-Bahl, K., et al. (2014). Oxygen and substrate availability interactively control the temperature sensitivity of CO2 and N2O emission from soil. Biol. Fertil. Soils 50, 775–783. doi: 10.1007/s00374-014-0899-6
Burns, R. G., DeForest, J. L., Marxsen, J., Sinsabaugh, R. L., Stromberger, M. E., Wallenstein, M. D., et al. (2013). Soil enzymes in a changing environment: current knowledge and future directions. Soil Biol. Biochem. 58, 216–234. doi: 10.1016/j.soilbio.2012.11.009
Chen, R., Blagodatskaya, E., Senbayram, M., Blagodatsky, S., Myachina, O., Dittert, K., et al. (2012). Decomposition of biogas residues in soil and their effects on microbial growth kinetics and enzyme activities. Biomass Bioenergy, 45, 221–229. doi: 10.1016/j.biombioe.2012.06.014
Craine, J. M., and Gelderman, T. M. (2011). Soil moisture controls on temperature sensitivity of soil organic carbon decomposition for a mesic grassland. Soil Biol. Biochem. 43, 455–457. doi: 10.1016/j.soilbio.2010.10.011
Craine, J. M., Spurr, R., McLauchlan, K. K., and Fierer, N. (2010). Landscape-level variation in temperature sensitivity of soil organic carbon decomposition. Soil Biol. Biochem. 42, 373–375. doi: 10.1016/j.soilbio.2009.10.024
Davidson, E. A., Janssens, I. A., and Luo, Y. (2006). On the variability of respiration in terrestrial ecosystems: moving beyond Q10. Glob. Change Biol. 12, 154–164. doi: 10.1111/j.1365-2486.2005.01065.x
De Graaff, M. A., Throop, H. L., Verburg, P. S. J., Arnone, J. A., and Campos, X. (2014). A synthesis of climate change and vegetation cover effects on biogeochemical cycling in shrub dominated drylands. Ecosystems 17, 931–945. doi: 10.1007/s10021-014-9764-6
Ekschmitt, K., Liu, M., Vetter, S., Fox, O., and Wolters, V. (2005). Strategies used by soil biota to overcome soil organic matter stability – why is dead organic matter left over in the soil? Geoderma 128, 167–176. doi: 10.1016/j.geoderma.2004.12.024
German, D. P., Marcelo, K. R. B., Stone, M. M., and Allison, S. D. (2012). The Michaelis–Menten kinetics of soil extracellular enzymes in response to temperature: a cross-latitudinal study. Glob. Change Biol. 18, 1468–1479. doi: 10.1111/j.1365-2486.2011.02615.x
German, D. P., Weintraub, M. N., Grandy, A. S., Lauber, C. L., Rinkes, Z. L., Allison, S. D., et al. (2011). Optimization of hydrolytic and oxidative enzyme methods for ecosystem studies. Soil Biol. Biochem. 43, 1387–1397. doi: 10.1016/j.soilbio.2011.03.017
Gershenson, A., Bader, N. E., and Cheng, W. (2009). Effects of substrate availability on the temperature sensitivity of soil organic matter decomposition. Glob. Change Biol. 15, 176–183. doi: 10.1111/j.1365-2486.2008.01827.x
Glanville, H., Hill, P. W., Maccarone, L., Golyshin, P., Murphy, D., Jones, D. L., et al. (2012). Temperature and water controls on vegetation, microbial dynamics, and soil carbon and nitrogen fluxes in a high Arctic tundra ecosystem. Funct. Ecol. 26, 1366–1380. doi: 10.1111/j.1365-2435.2012.02056.x
Goyal, M., Chaudhuri, T. K., and Kuwajima, K. (2014). Irreversible Denaturation of maltodextrin glucosidase studied by differential scanning calorimetry, circular dichroism, and turbidity measurements. PLoS ONE 9:115877. doi: 10.1371/journal.pone.0115877
Hasibeder, R., Fuchslueger, L., Richter, A., and Bahn, M. (2014). Summer drought alters carbon allocation to roots and root respiration in mountain grassland New Phytol. 205, 1117–1127. doi: 10.1111/nph.13146
Hoppe, H. G. (1983). Significance of exoenzymatic activities in the ecology of brackish water-measurements by means of methylumbelliferyl-substrates. Mar. Ecol. Prog. Ser. 11:299. doi: 10.3354/meps011299
Karhu, K., Auffret, M. D., Dungait, J. A. J., Hopkins, D. W., Prosser, J. I., Singh, B. K., et al. (2014). Temperature sensitivity of soil respiration rates enhanced by microbial community response. Nature 513, 81–84. doi: 10.1038/nature13604
Karhu, K., Fritze, H., Tuomi, M., Vanhala, P., Spetz, P., Kitunen, V., et al. (2010). Temperature sensitivity of organic matter decomposition in two boreal forest soil profiles. Soil Biol. Biochem. 42, 72–82. doi: 10.1016/j.soilbio.2009.10.002
Khalili, B., Nourbakhsh, F., Nili, N., Khademi, H., and Sharifnabi, B. (2011). Diversity of soil cellulase isoenzymes is associated with soil cellulase kinetic and thermodynamic parameters. Soil Biol. Biochem. 43, 1639–1648. doi: 10.1016/j.soilbio.2011.03.019
Kirschbaum, M. U. F. (1995). The temperature-dependence of soil organic-matter decomposition, and the effect of global warming on soil organic-C storage. Soil Biol. Biochem. 27, 753–760. doi: 10.1016/0038-0717(94)00242-S
Koch, O., Tscherko, D., and Kandeler, E. (2007). Temperature sensitivity of microbial respiration, nitrogen mineralization, and potential soil enzyme activities in organic alpine soils: temperature sensitivity in alpine soils. Global Biogeochem. Cycles 21:GB4017. doi: 10.1029/2007gb002983
Larionova, A. A., Yevdokimov, I. V., and Bykhovets, S. S. (2007). Temperature response of soil respiration is dependent on concentration of readily decomposable C. Biogeosciences 4, 1073–1081. doi: 10.5194/bg-4-1073-2007
Lee, Y. B., Lorenz, N., Dick, L. K., and Dick, R. P. (2007). Cold storage and pretreatment incubation effects on soil microbial properties. Soil Sci. Soc. Am. J. 71, 1299–1305. doi: 10.2136/sssaj2006.0245
Maire, V., Alvarez, G., Colombet, J., Comby, A., Despinasse, R., Dubreucq, E., et al. (2013). An unknown oxidative metabolism substantially contributes to soil CO2 emissions. Biogeosciences 10, 1155–1167. doi: 10.5194/bg-10-1155-2013
Marx, M. C., Kandeler, E., Wood, M., Wermbter, N., and Jarvis, S. C. (2005). Exploring the enzymatic landscape: distribution and kinetics of hydrolytic enzymes in soil particlesize fractions. Soil Biol. Biochem. 37, 35–48. doi: 10.1016/j.soilbio.2004.05.024
Melillo, J. M., Steudler, P. A., Aber, J. D., Newkirk, K., Lux, H., Bowles, F. P., et al. (2002). Soil warming and carbon-cycle feedbacks to the climate system. Science 298, 2173–2176. doi: 10.1126/science.1074153
Mganga, K. Z., Razavi, B. S., and Kuzyakov, K. (2015). Microbial and enzymes response to nutrient additions in soils of Mt. Kilimanjaro region depending on land use. Eur. J. Soil Biol. 69, 33–40. doi: 10.1016/j.ejsobi.2015.05.001
Nannipieri, P., and Gianfreda, L. (1998). “Kinetics of enzyme reactions in soil environments,” in Structure and Surface Reactions of Soil Particles, eds P. M. Huang, N. Senesi, and J. Buffle (New York, NY: John Wiley & Sons), 449–479.
Nannipieri, P., Sequi, P., and Fusi, P. (1996). “Humus and enzyme activity,” in Humic Substances in Terrestrial Ecosystems, ed A. Piccolo (Amsterdam: Elsevier), 293–328 doi: 10.1016/B978-044481516-3/50008-6
Nielsen, K. M., Calamai, L., and Pietramellara, G. (2006). “Stabilization of extracellular DNA and proteins by transient binding to various soil components,” in Nucleic Acids and Proteins in Soil. eds P. Nannipieri and K. Smalla (Heidelberg: Springer), 141–188.
Pietikäinen, J., Marie, P., and Bååth, E. (2005). Comparison of temperature effects on soil respiration and bacterial and fungal growth rates. FEMS Microbiol. Ecol. 52, 49–58. doi: 10.1016/j.femsec.2004.10.002
Ratkowsky, D. A., Lowry, R. K., McMeekin, T. A., Stokes, A. N., and Chandler, R. E. (1983). Model for bacterial culture growth rate throughout the entire biokinetic temperature range. J. Bacteriol. 154, 1222–1226.
Ruamps, L. S., Nunan, N., Pouteau, V., Leloup, J., Raynaud, X., Roy, V., et al. (2013). Regulation of soil organic C mineralisation at the pore scale. FEMS Microbiol. Ecol. 86, 26–35. doi: 10.1111/1574-6941.12078
Sand-Jensen, K. A. J., Pedersen, N. L., and Søndergaard, M. (2007). Bacterial metabolism in small temperate streams under contemporary and future climates. Freshw. Biol. 52, 2340–2353. doi: 10.1111/j.1365-2427.2007.01852.x
Skálová, T., Dohnálek, J., Spiwok, V., Lipovová, P., Vondráčková, E., Petroková, H., et al. (2005). Cold-active β-galactosidase from arthrobacter sp. C2-2 forms compact 660 kDa hexamers: crystal structure at 1.9 Aresolution. J. Mol. Biol. 353, 282–294. doi: 10.1016/j.jmb.2005.08.028
Spohn, M., and Kuzyakov, Y. (2013). Distribution of microbial- and root-derived phosphatase activities in the rhizosphere depending on P availability and C allocation – Coupling soil zymography with 14C imaging. Soil Biol. Biochem. 67, 106–113. doi: 10.1016/j.soilbio.2013.08.015
Stone, M. M., Weiss, M. S., Goodale, C. L., Adams, M. B., Fernandez, I. J., German, D. P., et al. (2012). Temperature sensitivity of soil enzyme kinetics under N-fertilization in two temperate forests. Glob. Change Biol. 18, 1173–1184. doi: 10.1111/j.1365-2486.2011.02545.x
Stotzky, G. (1986). “Influence of soil mineral colloids and metabolic processes, growth adhesion, and ecology of microbes and viruses,” in Interactions of Soil Minerals with Natural Organics and Microbes Vol. 17, eds M. Huang and M. Schnitzer (Madison, WI: Soil Science Society of America), 305–428.
Tabatabai, M. (1994). “Soil enzymes,” in Methods of Soil Analysis Part 2: Microbiological and Biochemical Properties, eds R. W. Weaver, S. Angle and P. Bottomley (Madison, WI: Soil Science Society of America), 775–833.
Tischer, A., Blagodatskaya, E., and Hamer, U. (2015). Microbial community structure and resource availability drive the catalytic efficiency of soil enzymes under land-use change conditions. Soil Biol. Biochem. 89, 226–237. doi: 10.1016/j.soilbio.2015.07.011
Tjoelker, M. G., Oleksyn, J., and Reich, P. B. (2001). Modelling respiration of vegetation: evidence for a general temperature-dependent Q10. Glob. Change Biol. 7, 223–230. doi: 10.1046/j.1365-2486.2001.00397.x
Van Gestel, N. C., Reischke, S., and Bååth, E. (2013). Temperature sensitivity of bacterial growth in a hot desert soil with large temperature fluctuations. Soil Biol. Biochem. 65, 180–185. doi: 10.1016/j.soilbio.2013.05.016
Wang, M., Lai, G. L., Nie, Y., Geng, S., Liu, L., Zhu, B., et al. (2015). Synergistic function of four novel thermostable glycoside hydrolases from a long-term enriched thermophilic methanogenic digester. Front. Microbiol. 6:509. doi: 10.3389/fmicb.2015.00509
Xu, M., and Qi, Y. (2001). Budgets of soil erosion and deposition for sediments and sedimentary organic carbon across the conterminous United States. Glob. Biogeochem. Cycles 15, 687–696. doi: 10.1029/2000GB001365
Zimmermann, M., and Bird, M. I. (2012). Temperature sensitivity of tropical forest soil respiration increase along an altitudinal gradient with ongoing decomposition. Geoderma, 187–188, 8–15. doi: 10.1016/j.geoderma.2012.04.015
Zoldák, G., Zubrik, A., Musatov, A., Stupák, M., and Sedlák, E. (2004). Irreversible thermal denaturation of glucose oxidase from Aspergillus niger is the transition to the denatured state with residual structure. J. Biol. Chem. 279, 47601–47609. doi: 10.1074/jbc.M406883200
Keywords: carbon cycle, Michaelis-Menten kinetics, Arrhenius function, soil enzymes, temperature sensitivity, canceling effect, activation energy
Citation: Razavi BS, Blagodatskaya E and Kuzyakov Y (2015) Nonlinear temperature sensitivity of enzyme kinetics explains canceling effect—a case study on loamy haplic Luvisol. Front. Microbiol. 6:1126. doi: 10.3389/fmicb.2015.01126
Received: 05 August 2015; Accepted: 28 September 2015;
Published: 14 October 2015.
Edited by:Paul Bodelier, Netherlands Institute of Ecology (NIOO-KNAW), Netherlands
Reviewed by:Michael S. Strickland, Virginia Tech, USA
Kirsten Hofmockel, Iowa State University, USA
Copyright © 2015 Razavi, Blagodatskaya and Kuzyakov. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Evgenia Blagodatskaya, firstname.lastname@example.org