Substrate reactivity as a function of the extent of reaction in the enzymatic hydrolysis of lignocellulose

In an effort to better understand the role of the substrate in the rapid fall off in the rate of enzymatic hydrolysis of cellulose with conversion, substrate reactivity was measured as a function of conversion. These measurements were made by interrupting the hydrolysis of pretreated wood at various degrees of conversion; and, after boiling and washing, restarting the hydrolysis in fresh buffer with fresh enzyme. The comparison of the restart rate per enzyme adsorbed with the initial rate per enzyme adsorbed, both extrapolated back to zero conversion, provides a measurement of the substrate reactivity without the complications of product inhibition or cellulase inactivation. The results indicate that the substrate reactivity falls only modestly as conversion increases. However, the restart rate is still higher than the rate of the uninterrupted hydrolysis, particularly at high conversion. Hence we conclude that the loss of substrate reactivity is not the principal cause for the long residence time required for complete conversion. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 56: 650-655, 1997.     

Kinetics of the enzymatic hydrolysis of cellulose: 1. A mathematical model for a batch reactor process

A mathematical model for enzymatic cellulose hydrolysis, based on experimental kinetics of the process catalysed by a cellulase [see 1,4-(1,3;1,4)-β-d-glucan 4-glucanohydrolase, EC 3.2.1.4] preparation from Trichoderma longibrachiatum has been developed. The model takes into account the composition of the cellulase complex, the structural complexity of cellulose, the inhibition by reaction products, the inactivation of enzymes in the course of the enzymatic hydrolysis and describes the kinetics of d-glucose and cellobiose formation from cellulose. The rate of d-glucose formation decelerated through the hydrolysis due to a change in cellulose reactivity and inhibition by the reaction product, d-glucose. The rate of cellobiose formation decelerated due to inhibition by the product, cellobiose, and inactivation of enzymes adsorbed on the cellulose surface. Inactivation of the cellobiose-producing enzymes as a result of their adsorption was found to be reversible. The model satisfactorily predicts the kinetics of d-glucose and cellobiose accumulation in a batch reactor up to 70–80% substrate conversion on changing substrate concentration from 5 to 100 g l^?1and the concentration of the enzymic preparation from 5 to 60 g l^?1.     

Thermal inactivation and product inhibition of Aspergillus terreus CECT 2663 alpha-L-rhamnosidase and their role on hydrolysis of naringin solutions

The kinetics of thermal inactivation of A. terreus alpha-rhamnosidase was studied using the substrate p-nitrophenyl alpha-L-rhamnoside between 50 degrees C and 70 degrees C. Up to 60 degrees C the inactivation of the purified enzyme was completely reversible, but samples of crude or partially purified enzyme showed partial reversibility. The presence of the product rhamnose, the substrate naringin, and other additives reduced the reversible inactivation, maintaining in some cases full enzyme activity at 60 degrees C. A mechanism for the inactivation process, which permitted the reproduction of experimental results, was proposed. The products rhamnose (inhibition constant, 2.1 mM) and prunin (2.6 mM) competitively inhibited the enzyme reaction. The maximum hydrolysis of supersaturated naringin solution, without enzyme inactivation, was observed at 60 degrees C. Hydrolysis of naringin reached 99% with 1% naringin solution, although the hydrolysis degree of naringin was only 40% due to products inhibition when the initial concentration of flavonoid was 10%. The experimental results fitted an equation based on the integrated Michaelis-Menten's, including competitive inhibition by products satisfactorily.     

Heat-stable protease from Pseudomonas fluorescens T16: purification by affinity column chromatography and characterization.

A heat-stable extracellular protease from Pseudomonas fluorescens T16, a psychrotroph, was purified by affinity column chromatography on a carbobenzoxy-D-phenylalanine-triethylene tetramine-Sepharose-4B column. The purified enzyme is a monomer with a molecular weight of 38,905 +/- 2,000. In an analytical ultracentrifuge, the Schlieren profile revealed a single symmetrical peak. The sedimentation coefficient was estimated to be 3.93S. Alpha-casein was the preferred substrate, with a Km of 0.05 mM. Heating crude enzyme and purified enzyme in buffer at 50, 90, and 120 degrees C resulted in a rapid initial loss of more than 50% of the initial activity followed by a gradual inactivation which exhibited first-order kinetics. The activation energy for the hydrolysis of casein was calculated to be 3.2 kcal/mol (13.4 kJ/mol).     

Heat-stable protease from Pseudomonas fluorescens T16: purification by affinity column chromatography and characterization

A heat-stable extracellular protease from Pseudomonas fluorescens T16, a psychrotroph, was purified by affinity column chromatography on a carbobenzoxy-D-phenylalanine-triethylene tetramine-Sepharose-4B column. The purified enzyme is a monomer with a molecular weight of 38,905 +/- 2,000. In an analytical ultracentrifuge, the Schlieren profile revealed a single symmetrical peak. The sedimentation coefficient was estimated to be 3.93S. Alpha-casein was the preferred substrate, with a Km of 0.05 mM. Heating crude enzyme and purified enzyme in buffer at 50, 90, and 120 degrees C resulted in a rapid initial loss of more than 50% of the initial activity followed by a gradual inactivation which exhibited first-order kinetics. The activation energy for the hydrolysis of casein was calculated to be 3.2 kcal/mol (13.4 kJ/mol).     

Acetic acid formation in Escherichia coli fermentation

Theoretical analysis of cellulase product inhibition (by cellobiose and glucose) has been performed in terms of the mathematical model for enzymatic cellulose hydrolysis. The analysis showed that even in those cases when consideration of multienzyme cellulase system as one enzyme (cellulase) or two enzymes (cellulase and beta-glucosidase) is valid, double-reciprocal plots, usually used in a product inhibition study, may be nonlinear, and different inhibition patterns (noncompetitive, competitive, or mixed type) may be observed. Inhibition pattern depends on the cellulase binding constant, enzyme concentration, maximum adsorption of the enzyme (cellulose surface area accessible to the enzyme), the range in which substrate concentration is varied, and beta-glucosidase activity. A limitation of cellulase adsorption by cellulose surface area that may occur at high enzyme/substrate ratio is the main reason for nonlinearity of double-reciprocal plots. Also, the results of calculations showed that material balance by substrate, which is usually neglected by researchers studying cellulase product inhibition, must be taken into account in kinetic analysis even in those cases when the enzyme concentration is rather low. (c) 1992 John Wiley & Sons, Inc.     

Kinetics of the enzymatic hydrolysis of cellulose: 2. A mathematical model for the process in a plug-flow column reactor

The kinetics of enzymatic cellulose hydrolysis in a plug-flow column reactor catalysed by cellulases [see 1,4-(1,3;1,4)-β-d-glucan 4-glucanohydrolase, EC 3.2.1.4] from Trichoderma longibrachiatum adsorbed on cellulose surface have been studied. The maximum substrate conversion achieved was 90–94%. The possibility of enzyme recovery for a reactor of this type is discussed. A mathematical model for enzymatic cellulose hydrolysis in a plug-flow column reactor has been developed. The model allows for the component composition of the cellulase complex, adsorption of cellulases on the substrate surface, inhibition by reaction products, changes in cellulose reactivity and the inactivation of enzymes in the course of hydrolysis. The model affords a reliable prediction of the kinetics of d-glucose and cellobiose formation from cellulose in a column reactor as well as the degree of substrate conversion and reactor productivity with various amounts of adsorbed enzymes and at various flow rates.     

Origin of initial burst in activity for Trichoderma reesei endo-glucanases hydrolyzing insoluble cellulose

The kinetics of cellulose hydrolysis have long been described by an initial fast hydrolysis rate, tapering rapidly off, leading to a process that takes days rather than hours to complete. This behavior has been mainly attributed to the action of cellobiohydrolases and often linked to the processive mechanism of this exo-acting group of enzymes. The initial kinetics of endo-glucanases (EGs) is far less investigated, partly due to a limited availability of quantitative assay technologies. We have used isothermal calorimetry to monitor the early time course of the hydrolysis of insoluble cellulose by the three main EGs from Trichoderma reesei (Tr): TrCel7B (formerly EG I), TrCel5A (EG II), and TrCel12A (EG III). These endo-glucanases show a distinctive initial burst with a maximal rate that is about 5-fold higher than the rate after 5 min of hydrolysis. The burst is particularly conspicuous for TrCel7B, which reaches a maximal turnover of about 20 s(-1) at 30 °C and conducts about 1200 catalytic cycles per enzyme molecule in the initial fast phase. For TrCel5A and TrCel12A the extent of the burst is 2-300 cycles per enzyme molecule. The availability of continuous data on EG activity allows an analysis of the mechanisms underlying the initial kinetics, and it is suggested that the slowdown is linked to transient inactivation of enzyme on the cellulose surface. We propose, therefore, that the frequency of structures on the substrate surface that cause transient inactivation determine the extent of the burst phase.     

The effect of histidine residue modification on tyrosinase activity and conformation: inhibition kinetics and computational prediction

We found that the histidine chemical modification of tyrosinase conspicuously inactivated enzyme activity. The substrate reactions with diethylpyridinecarbamate showed slow-binding inhibition kinetics (K(I) = 0.24 +/- 0.03 mM). Bromoacetate, as another histidine modifier, was also applied in order to study inhibition kinetics. The bromoacetate directly induced the exposures of hydrophobic surfaces following by complete inactivation via ligand binding. For further insights, we predicted the 3D structure of tyrosinase and simulated the docking between tyrosinase and diethylpyridinecarbamate. The docking simulation was shown to the significant binding energy scores (-3.77 kcal/mol by AutoDock4 and -25.26 kcal/mol by Dock6). The computational prediction was informative to elucidate the role of free histidine residues at the active site, which are related to substrate accessibility during tyrosinase catalysis.     

Kinetic modeling of cellulose hydrolysis with first order inactivation of adsorbed cellulase

Enzymatic hydrolysis involves complex interaction between enzyme, substrate, and the reaction environment, and the complete mechanism is still unknown. Further, glucose release slows significantly as the reaction proceeds. A model based on Langmuir binding kinetics that incorporates inactivation of adsorbed cellulase was developed that predicts product formation within 10% of experimental results for two substrates. A key premise of the model, with experimental validation, suggests that V(max) decreases as a function of time due to loss of total available enzyme as adsorbed cellulases become inactivated. Rate constants for product formation and enzyme inactivation were comparable to values reported elsewhere. A value of k(2)/K(m) that is several orders of magnitude lower than the rate constant for the diffusion-controlled encounter of enzyme and substrate, along with similar parameter values between substrates, implies a common but undefined rate-limiting step associated with loss of enzyme activity likely exists in the pathway of cellulose hydrolysis.     

Kinetic studies of enzymatic hydrolysis of insoluble cellulose: (II). Analysis of extended hydrolysis times

The kinetics of enzymatic hydrolysis of pure insoluble cellulose by means of unpurified culture filtrate of Trichoderma reesei was studied, emphasizing the kinetic characteristics associated with the extended hydrolysis times. The changes in the hydrolysis rate and extent of soluble protein adsorption during the progress of reaction, either apparent or intrinsic, were investigated. The hydrolysis rate declined drastically during the initial hours of hydrolysis. The factors causing the reduction in the hydrolysis rate were examined; these include the transformation of cellulose into a less digestible form and product inhibition. The structural transformation can be partially explained by changes in the crystallinity index and surface area. The product inhibition was caused by the deactivation of the adsorbed soluble protein by the products, which essentially represents the so-called "un-competitive" inhibition. The kinetics of beta-glucosidase were also studied. The result has shown that the action of beta-glucosidase is competitively inhibited by glucose. It has been found that the integrated form of the initial rate expression cannot be used in predicting the progress of reaction because the digestibility of cellulose changes drastically as the hydrolysis proceeds, and that the rate expression for enzymatic hydrolysis of cellulose cannot be simplified or approximated by resorting to the pseudo-steady-state assumption. A mechanistic kinetic model of cellulose hydrolysis should include the following major influencing factors: (1)mode of action of enzyme, (2) structure of cellulose, and (3) mode of interaction between the enzyme and cellulose molecules.     

A mechanistic model for enzymatic saccharification of cellulose using continuous distribution kinetics II: cooperative enzyme action, solution kinetics, and product inhibition

The projected cost for the enzymatic hydrolysis of cellulosic biomass continues to be a barrier for the commercial production of liquid transportation fuels from renewable feedstocks. Predictive models for the kinetics of the enzymatic reactions will enable an improved understanding of current limitations, such as the slow-down of the overall conversion rate, and may point the way for more efficient utilization of the enzymes in order to achieve higher conversion yields. A mechanistically based kinetic model for the enzymatic hydrolysis of cellulose was recently reported in Griggs et al. (2011) (Part I). In this article (Part II), the enzyme system is expanded to include solution-phase kinetics, particularly cellobiose-to-glucose conversion by β-glucosidase (βG), and novel adsorption and product inhibition schemes have been incorporated, based on current structural knowledge of the component enzymes. Model results show cases of cooperative and non-cooperative hydrolysis for an enzyme system consisting of EG(I) and CBH(I). The model is used to explore various potential rate-limiting phenomena, such as substrate accessibility, product inhibition, sterically hindered enzyme adsorption, and the molecular weight of the cellulose substrate.     

Thermal Inactivation and Product Inhibition of Aspergillus terreus CECT 2663 α-L-Rhamnosidase and Their Role on Hydrolysis of Naringin Solutions

The kinetics of thermal inactivation of A. terreus α-rhamnosidase was studied using the substrate p-nitrophenyl α-L-rhamnoside between 50°C and 70°C. Up to 60°C the inactivation of the purified enzyme was completely reversible, but samples of crude or partially purified enzyme showed partial reversibility. The presence of the product rhamnose, the substrate naringin, and other additives reduced the reversible inactivation, maintaining in some cases full enzyme activity at 60°C. A mechanism for the inactivation process, which permitted the reproduction of experimental results, was proposed. The products rhamnose (inhibition constant, 2.1 mM) and prunin (2.6 mM) competitively inhibited the enzyme reaction. The maximum hydrolysis of supersaturated naringin solution, without enzyme inactivation, was observed at 60°C. Hydrolysis of naringin reached 99% with 1% naringin solution, although the hydrolysis degree of naringin was only 40% due to products inhibition when the initial concentration of flavonoid was 10%. The experimental results fitted an equation based on the integrated Michaelis-Menten's, including competitive inhibition by products satisfactorily.     

A kinetic model for the burst phase of processive cellulases

Cellobiohydrolases (exocellulases) hydrolyze cellulose processively, i.e. by sequential cleaving of soluble sugars from one end of a cellulose strand. Their activity generally shows an initial burst, followed by a pronounced slowdown, even when substrate is abundant and product accumulation is negligible. Here, we propose an explicit kinetic model for this behavior, which uses classical burst phase theory as the starting point. The model is tested against calorimetric measurements of the activity of the cellobiohydrolase Cel7A from Trichoderma reesei on amorphous cellulose. A simple version of the model, which can be solved analytically, shows that the burst and slowdown can be explained by the relative rates of the sequential reactions in the hydrolysis process and the occurrence of obstacles for the processive movement along the cellulose strand. More specifically, the maximum enzyme activity reflects a balance between a rapid processive movement, on the one hand, and a slow release of enzyme which is stalled by obstacles, on the other. This model only partially accounts for the experimental data, and we therefore also test a modified version that takes into account random enzyme inactivation. This approach generally accounts well for the initial time course (approximately 1 h) of the hydrolysis. We suggest that the models will be useful in attempts to rationalize the initial kinetics of processive cellulases, and demonstrate their application to some open questions, including the effect of repeated enzyme dosages and the 'double exponential decay' in the rate of cellulolysis.     

A kinetic study of Trichoderma reesei Cel7B catalyzed cellulose hydrolysis

One prominent feature of Trichoderma reesei (Tr) endoglucanases catalyzed cellulose hydrolysis is that the reaction slows down quickly after it starts (within minutes). But the mechanism of the slowdown is not well understood. A structural model of Tr- Cel7B catalytic domain bound to cellulose was built computationally and the potentially important binding residues were identified and tested experimentally. The 13 tested mutants show different binding properties in the adsorption to phosphoric acid swollen cellulose and filter paper. Though the partitioning parameter to filter paper is about 10 times smaller than that to phosphoric acid swollen cellulose, a positive correlation is shown for two substrates. The kinetic studies show that the reactions slow down quickly for both substrates. This slowdown is not correlated to the binding constant but anticorrelated to the enzyme initial activity. The amount of reducing sugars released after 24 h by Cel7B in phosphoric acid swollen cellulose, Avicel and filter paper cellulose hydrolysis is correlated with the enzyme activity against a soluble substrate p-nitrophenyl lactoside. Six of the 13 tested mutants, including N47A, N52D, S99A, N323D, S324A, and S346A, yield ∼15–35% more reducing sugars than the wild type (WT) Cel7B in phosphoric acid swollen cellulose and filter paper hydrolysis. This study reveals that the slowdown of the reaction is not due to the binding of the enzyme to cellulose. The activity of Tr- Cel7B against the insoluble substrate cellulose is determined by the enzyme’s capability in hydrolyzing the soluble substrate.     

Relation between filter paper activity and saccharifying ability of a Trichoderma cellulase system

Summary Tests made to study the relation between filter paper activity and actual saccharifying ability of Trichoderma cellulases show that 30 IU/g of cellulose were sufficient to achieve over 80% hydrolysis of a 25 g/L cellulose suspension in 24 h. With the same enzyme/substrate ratio, but double the concentration of substrate, about 60% hydrolysis was achieved. End- product inhibition is one factor which seriously limits the degree of hydrolysis and therefore the concentration of sugars achievable by enzymatic hydrolysis at high levels of substrate concentration or enzyme/substrate ratio.     

Modeling intrinsic kinetics of enzymatic cellulose hydrolysis

A multistep approach was taken to investigate the intrinsic kinetics of the cellulase enzyme complex as observed with hydrolysis of noncrystalline cellulose (NCC). In the first stage, published initial rate mechanistic models were built and critically evaluated for their performance in predicting time-course kinetics, using the data obtained from enzymatic hydrolysis experiments performed on two substrates: NCC and alpha-cellulose. In the second stage, assessment of the effect of reaction intermediates and products on intrinsic kinetics of enzymatic hydrolysis was performed using NCC hydrolysis experiments, isolating external factors such as mass transfer effects, physical properties of substrate, etc. In the final stage, a comprehensive intrinsic kinetics mechanism was proposed. From batch experiments using NCC, the time-course data on cellulose, cello-oligosaccharides (COS), cellobiose, and glucose were taken and used to estimate the parameters in the kinetic model. The model predictions of NCC, COS, cellobiose, and glucose profiles show a good agreement with experimental data generated from hydrolysis of different initial compositions of substrate (NCC supplemented with COS, cellobiose, and glucose). Finally, sensitivity analysis was performed on each model parameter; this analysis provides some insights into the yield of glucose in the enzymatic hydrolysis. The proposed intrinsic kinetic model parametrized for dilute cellulose systems forms a basis for modeling the complex enzymatic kinetics of cellulose hydrolysis in the presence of limiting factors offered by substrate and enzyme characteristics.     

Comparing the thermodynamic stabilities of a related thermophilic and mesophilic enzyme

Several models have been proposed to explain the high temperatures required to denature enzymes from thermophilic organisms; some involve greater maximum thermodynamic stability for the thermophile, and others do not. To test these models, we reversibly melted two analogous protein domains in a two-state manner. E2cd is the isolated catalytic domain of cellulase E2 from the thermophile Thermomonospora fusca. CenAP30 is the analogous domain of the cellulase CenA from the mesophile Cellulomonas fimi. When reversibly denatured in a common buffer, the thermophilic enzyme E2cd had a temperature of melting (Tm) of 72.2 degrees C, a van't Hoff enthalpy of unfolding (DeltaHVH) of 190 kcal/mol, and an entropy of unfolding (DeltaSu) of 0.55 kcal/(mol*K); the mesophilic enzyme CenAP30 had a Tm of 56.4 degrees C, a DeltaHVH of 107 kcal/mol, and a DeltaSu of 0. 32 kcal/(mol*K). The higher DeltaHVH and DeltaSu values for E2cd suggest that its free energy of unfolding (DeltaGu) has a steeper dependence on temperature at the Tm than CenAP30. This result supports models that predict a greater maximum thermodynamic stability for thermophilic enzymes than for their mesophilic counterparts. This was further explored by urea denaturation. Under reducing conditions at 30 degrees C, E2cd had a concentration of melting (Cm) of 5.2 M and a DeltaGu of 11.2 kcal/mol; CenAP30 had a Cm of 2.6 M and a DeltaGu of 4.3 kcal/mol. Under nonreducing conditions, the Cm and DeltaGu of CenAP30 were increased to 4.5 M and 10.8 kcal/mol at 30 degrees C; the Cm for E2cd was increased to at least 7.4 M at 32 degrees C. We were unable to determine a DeltaGu value for E2cd under nonreducing conditions due to problems with reversibility. These data suggest that E2cd attains its greater thermal stability (DeltaTm = 15.8 degrees C) through a greater thermodynamic stability (DeltaDeltaGu = 6.9 kcal/mol) compared to its mesophilic analogue CenAP30.     

Inactivation of o-diphenoloxidase in the pyrocatechol oxidation reaction

The inactivation kinetics of o-diphenoloxidase isolated from potato tubers was studied in the process of pyrocatechol oxidation. The enzyme when saturated with the substrate is inactivated with the inactivation rate constant kin = 0.5-1.0 min-1; kin depends on the initial concentration of pyrocatechol. The ultimate yield of the enzymic reaction product increases linearly with the initial concentration of the enzyme. Introduction of ethylene-diaminosulphate, a substance which condenses with o-quinones, does not increase the operation stability of o-diphenoloxidase. The data obtained evidence for inactivation of o-diphenoloxidase either at the level of the enzyme-substrate complex or due to bimolecular reaction with the substrate.