Skip to main content

A novel film–pore–surface diffusion model to explain the enhanced enzyme adsorption of corn stover pretreated by ultrafine grinding

Abstract

Background

Ultrafine grinding is an environmentally friendly pretreatment that can alter the degree of polymerization, the porosity and the specific surface area of lignocellulosic biomass and can, thus, enhance cellulose hydrolysis. Enzyme adsorption onto the substrate is a prerequisite for the enzymatic hydrolysis process. Therefore, it is necessary to investigate the enzyme adsorption properties of corn stover pretreated by ultrafine grinding.

Results

The ultrafine grinding pretreatment was executed on corn stover. The results showed that ultrafine grinding pretreatment can significantly decrease particle size [from 218.50 μm of sieve-based grinding corn stover (SGCS) to 17.45 μm of ultrafine grinding corn stover (UGCS)] and increase the specific surface area (SSA), pore volume (PV) and surface composition (SSA: from 1.71 m2/g of SGCS to 2.63 m2/g of UGCS, PV: from 0.009 cm3/g of SGCS to 0.024 m3/g of UGCS, cellulose surface area: from 168.69 m2/g of SGCS to 290.76 m2/g of UGCS, lignin surface area: from 91.46 m2/g of SGCS to 106.70 m2/g of UGCS). The structure and surface composition changes induced by ultrafine grinding increase the enzyme adsorption capacity from 2.83 mg/g substrate of SGCS to 5.61 mg/g substrate of UGCS. A film–pore–surface diffusion model was developed to simultaneously predict the enzyme adsorption kinetics of both the SGCS and UGCS. Satisfactory predictions could be made with the model based on high R 2 and low RMSE values (R 2 = 0.95 and RMSE = 0.16 mg/g for the UGCS, R 2 = 0.93 and RMSE = 0.09 mg/g for the SGCS). The model was further employed to analyze the rate-limiting steps in the enzyme adsorption process. Although both the external-film and internal-pore mass transfer are important for enzyme adsorption on the SGCS and UGCS, the UGCS has a lower internal-pore resistance compared to the SGCS.

Conclusions

Ultrafine grinding pretreatment can enhance the enzyme adsorption onto corn stover by altering structure and surface composition. The film–pore–surface diffusion model successfully captures features on enzyme adsorption on ultrafine grinding pretreated corn stover. These findings identify wherein the probable rate-limiting factors for the enzyme adsorption reside and could, therefore, provide a basis for enhanced cellulose hydrolysis processes.

Background

Lignocellulosic biomass, such as crop residues, is the only renewable and sustainable resource that can be stored and transported. The annual yield of crop residues is abundant according to the FAO Statistics (2013), which indicated that approximately 12 % of the world’s land area is used for crop production [1]. Among the main crop residues, corn stover is one of the most favorable bioethanol feedstocks because of its wide geographic distribution and high cellulose content. The bioethanol conversion of corn stover has attracted the interest of scientists around the world [2, 3].

For the conversion of lignocellulosic biomass to bioethanol, the key bottleneck is the initial conversion of biomass to sugars. It is well known that lignocellulosic biomass, in its native form, is recalcitrant to hydrolysis with cellulase enzyme systems in the biochemical conversion process. To overcome biomass recalcitrance and improve cellulose accessibility, many chemical pretreatment methods (acid [4], alkali [5], ammonia fiber explosion [2] and so on [6]) were employed. However, these chemical pretreatment methods generate highly toxic effluents and cause negative impacts on the environment. Mechanical comminution is an environmentally friendly pretreatment that can alter the degree of polymerization, crystallinity degree, porosity and specific surface area of lignocellulosic biomass and, thus, enhance cellulose hydrolysis [7]. Most previous studies on the mechanical comminution pretreatment of lignocellulosic biomass were usually carried out by chipping (10–30 mm), grinding and milling (0.2–2 mm) [79].

Recently, ultrafine grinding (approximately 25 μm) technology, which can achieve a small particle size, large specific surface area, and high chemical activity [10], was also sporadically explored in the field of lignocellulose pretreatment. For example, Silva et al. investigated the effects of grinding processes on the enzymatic degradation of wheat straw [11]. The results showed that the ultrafine grinding pretreatment significantly enhanced enzymatic hydrolysis yield up to 10-fold as compared with coarsely grinding. Although some properties, such as the particle size and cellulose crystallinity, had been characterized to explain the hydrolysis mechanism after the ultrafine grinding pretreatment, some intrinsic properties, such as the adsorption kinetics, should be further investigated.

Cellulase adsorption onto the substrate via the binding domain is a prerequisite step for the enzymatic hydrolysis process and directly affects the enzymatic hydrolysis yield of lignocellulosic biomass [12, 13]. Thus, an adequate description of the adsorption step is indispensable for understanding and optimizing hydrolysis reaction, especially for that after the ultrafine grinding pretreatment. It is well known that ultrafine grinding increases the available specific surface area/pore volume [14] and, thus, improves the exposure level of the cellulose-binding domain, which is closely related to the cellulase adsorption kinetics. However, the cellulase adsorption kinetics of lignocellulosic biomass after the ultrafine grinding pretreatment has never been reported until now. Previous experimental and modeling studies on the cellulase adsorption of lignocellulosic biomass mainly focused on those pretreated by chemical methods, such as acid [15], hydrothermal [13], organosolv [13], and SO2-catalyzed steam explosion [16]. These studies commonly characterized cellulase adsorption by the Langmuir isotherm model, which describes the relationship between the amount of enzyme protein binding with substrate and the amount of enzyme protein free in solution after attaining equilibrium adsorption [17]. The Langmuir isotherm model can evaluate the maximum adsorption capacity of the substrate under different enzyme loadings, but it is not capable of expressing the adsorption kinetics of cellulase along with the adsorption time. The adsorption kinetics can be used to better understand the rate-controlling step of the mass transfer involved in the adsorption process. From a mechanistic viewpoint, the adsorption of cellulase onto lignocellulosic biomass can include three consecutive steps: the external diffusion of cellulase from bulk solution across the liquid film surrounding the solid biomass particles, internal diffusion of cellulase through the biomass particles by pore volume diffusion and surface diffusion, and the adsorption of cellulose molecules onto the biomass particles at the active sites (Fig. 1).

Fig. 1
figure 1

The schematic illustration of mass transfer for enzyme adsorption on a porous substrate. The adsorption of cellulase onto lignocellulosic biomass includes three consecutive steps: a the external diffusion of cellulase from bulk solution across the liquid film surrounding the solid biomass particles, b internal diffusion of cellulase through the biomass particles by pore volume diffusion and surface diffusion, and c the adsorption of cellulose molecules onto the biomass particles at the active sites

This study first investigated the enzyme adsorption kinetics of ultrafine grinding pretreated corn stover. Then, a filmporesurface diffusion model was developed to explain the enzyme adsorption kinetics. The rate-limiting steps in the adsorption process were further investigated. To our knowledge, this is the first work in the literature to reveal the enzyme adsorption behavior of corn stover pretreated by ultrafine grinding, thus identifying wherein the probable rate-limiting factors for the enzyme adsorption reside and could, therefore, provide a basis for enhanced cellulose hydrolysis processes.

Results and discussion

Carbohydrates and lignin content

The cellulose, hemicellulose, and lignin content of the sieve-based grinding corn stover (SGCS) and ultrafine grinding corn stover (UGCS) is listed in Table 1. The amounts are similar to the previously reported values. For example, Li et al. observed 34.9 % glucan, 21.7 % xylan, and 20.5 % lignin in corn stover [18]. Saha et al. observed 37.0 % cellulose, 28.9 % hemicellulose, and 21.2 % lignin in corn stover [19]. It was also shown that there are no significant differences in the carbohydrates and lignin content of both substrates. Previous studies have noted that the chemical components have important adsorption interactions with enzyme molecules, although the enzyme adsorption of lignin is considered a nonproductive one [15]. Even if the two substrates present the similar contents, the grinding may affect the surface composition and, thus, change the adsorption capacity/affinity of the enzyme for the substrate. This is further corroborated by surface composition measurement for the two substrates. The surface areas of cellulose and lignin, which are two dominant components in the cellulase adsorption [20, 21], were measured by determining the maximum adsorption capacity of the dyes Congo Red [22] and Azure B [23] on the substrates, respectively (Table 1). The cellulose surface area of the UGCS (290.76 m2/g) was almost twofold higher than that of SGCS (168.69 m2/g). Compared with the lignin surface area of the SGCS (91.46 m2/g), that of the UGCS (106.7 m2/g) also moderately increased. These results indicated that the substrate pretreated by ultrafine grinding can induce more exposure of the surface composition (especially for cellulose), which will be favorable to enzyme adsorption.

Table 1 Chemical content and structural properties of SGCS and UGCS

Particle size distribution and morphology

Figure 2 shows the particle size distributions of both the SGCS and UGCS. The particle size distribution was characterized by the median diameter (d 50) and the span defined by (d 90d 10)/d 50, where d 10, d 50 and d 90 represent the 10th, 50th and 90th percentiles of the total volume, respectively [11]. The median sizes (d 50) of the UGCS and SGCS were 17.45 μm and 218.50 μm, respectively. The spans of the UGCS and SGCS were 2.72 and 2.93, respectively. The smaller span value indicated a more uniform size distribution. Severe vibration ball milling under the ultrafine grinding condition destroyed the fiber structure and, thus, achieved significant particle size reduction and unified particle size distribution. The ultrafine grinding of crop residues was also reported by several studies. For example, Silva et al. investigated the median particle sizes and particle size distribution spans of wheat straw under the operating conditions of ball milling and jet milling [11]. Ball milling reduced the particle size from 270 to 16 μm over a 0–240 h period. The span first increased to more than 5 during the first 120 h and then decreased to 2.5 at the end of the 120 h. Jet milling reduced the median particle size of wheat straw from 107 to 22 μm and was much more rapid (85 min) than ball milling. A previous study by our team also explored the ultrafine grinding of wheat straw by 8 h of vibration ball milling and reported ultrafine wheat straw powder with a median size of 17.0 μm and a span of 4.0 [24]. Compared with previous studies, our study produced ultrafine powder of corn stover in a shorter time (30 min), which indicated less energy consumption.

Fig. 2
figure 2

Particle size distribution of the UGCS and SGCS. These data were determined by a laser diffraction particle size analyzer

Specific surface area (SSA) and pore volume (PV) distribution

The SSA and PV distribution of the SGCS and UGCS is listed in Table 1 and Fig. 3. The SSA of the UGCS was approximately 1.5-fold higher than that of the SGCS (Table 1). Although the values between the SSA and the surface composition areas were uncomparable due to different measured methods [22], their similar increased trends for the UGCS indicated that the ultrafine grinding pretreatment significantly affects substrate structure and surface composition. The PV of the UGCS was approximately threefold higher than that of the SGCS (Fig. 3a). The UGCS had a wider pore volume distribution (2–300 nm) than the SGCS (2–50 nm) based on differential curves of the pore volume distribution (Fig. 3b), which indicated that mesopores and macropores existed in the UGCS. The SSA and PV properties are important parameters for the conversion of lignocellulosic biomass to biofuels and are often useful to ascertain whether the comminution pretreatment technology is useful or not [25]. Commonly, the comminution pretreatment can enhance the SSA of lignocellulosic biomass. This is because drastic milling to the straw can destroy the structure of the lignocellulose, disorganizing the tightly ordered fibers and exposing more enzyme bonding sites [26, 27]. Piccolo et al. found that the SSA increased by more than 60 % after ball milling compared to untreated wheat straw samples [28]. Furthermore, the SSA is highly sensitive to the particle size of lignocellulosic biomass. Zhang et al. reported a linear correlation of the SSA with particle size for pan-milling cellulose powder [29]. The SSA is not only related to the particle size, but is also strongly related to the PV of the lignocellulosic biomass. The surface area of the substrate can be divided into an interior surface area, reflected by the biomass porosity, and an exterior surface area, largely determined by the particle size. Compared with the sieve-based grinding pretreatment, the ultrafine grinding pretreatment can produce more significant changes in the internal-pore structure, and these changes are mainly responsible for the enzymatic adsorption and hydrolysis of biomass [30, 31]. The size of a cellulase is approximately 5.1 nm [32] and, hence, only those pores larger than 5.1 nm are accessible to enzyme. The pore accessible to enzyme is correlated with the enzyme diffusion resistance and adsorption rate [32, 33]. Compared with the SGCS, the UGCS has a higher volume fraction of pores larger than 5.1 nm in diameter (Fig. 3a).

Fig. 3
figure 3

Pore volume distribution of the UGCS and SGCS as a function of pore diameter. These data were determined by liquid nitrogen adsorption experiments. a Cumulative pore volume versus pore diameter; b Differential pore volume (dV/dD) versus pore diameter

Equilibrium adsorption

The Langmuir isotherm model agreed well with the equilibrium adsorption data of both the SGCS and UGCS (Fig. 4) based on their statistical parameters (R 2 ≥ 0.90, RMSE ≤ 0.20 mg/g). The Langmuir parameters, including the maximum adsorption capacity (q m), affinity constant (K a) and bonding strength (S = q m × K a), are listed in Table 2. A number of previous studies carried out the cellulase equilibrium adsorption of lignocellulosic biomass and also observed robust adaptability of the Langmuir model [34]. For example, Machado et al. investigated the adsorption characteristics of cellulase on Avicel, pretreated sugarcane bagasse, and lignin [13]. Langmuir model isotherms were chosen to compare the kinetic properties of these various enzyme-substrate systems. Qi et al. explored cellulase adsorption of two different pretreated wheat straws and proposed a good fit to the cellulase adsorption data by the Langmuir adsorption isotherm [35]. It is difficult to directly compare the Langmuir parameters of this study to those of previous studies for different combinations for enzyme, substrate, and temperature. Zhang and Lynd collected Langmuir parameters for the cellulase adsorption of lignocellulosic biomass and observed wide variations [36]. However, the Langmuir parameters of both the SGCS and UGCS in this study can be directly compared because of the same experimental conditions. The results showed that the q m (5.61 mg/g) and K a (11.5 mL/mg) values obtained for the UGCS were much higher than those (q m = 2.83 mg/g, K a = 6.22 mL/mg) for the SGCS. These results indicated that the substrate pretreated by ultrafine grinding has a stronger adsorption capability of enzyme molecules. The reason for the high cellulase adsorption amount of UGCS may be because the ultrafine grinding achieved significant changes in the intact cellulose–hemicellulose–lignin network. More generated pores, demonstrated by a high SSA and PV distribution, increased the diffusion of enzyme molecules into the substrate. More importantly, more exposed binding sites of the substrate, demonstrated by a high cellulose and lignin surface area, improved the substrate accessibility to cellulase.

Fig. 4
figure 4

Equilibrium adsorption of cellulose to a SGCS and b UGCS. The equilibrium adsorption experiments were performed with different loadings of the cellulase (1.5–10.5 mg/g substrate for celluclast 1.5 L). The cellulase adsorption data were fitted by Langmuir equilibrium isotherm. Error bars represent the standard deviation of the measurements for the bound enzyme

Table 2 Langmuir adsorption isotherm parameters of SGCS and UGCS

Filmporesurface diffusion model for the enzymatic adsorption kinetics

The cellulase adsorption kinetic profiles of the SGCS and UGCS are shown in Fig. 5. Compared with the kinetic data of the SGCS, the adsorption amount of the UGCS at any time was much higher. This may be explained by the changes induced by the ultrafine grinding pretreatment, which yielded high SSA, PV and surface composition areas. On the one hand, high SSA, PV and surface composition values produced a large exposure area of the substrate and, thus, the binding sites of the substrate to the cellulase are also accordingly increased to achieve high enzymatic adsorption capacity. On the other hand, large pore openings produce less restriction and provide efficient adsorption of the enzyme molecules [37]. Wang et al. [11] investigated the cellulase adsorption and cellulose accessibility to the cellulase of the set of pretreated substrates with different pore volume distributions [31]. The authors found that increasing the pore volume in the substrates increases the cellulose accessibility to cellulase, which correlated well with the amounts of adsorbed cellulase.

Fig. 5
figure 5

Comparison of observed and predicted cellulase adsorption kinetics for SGCS and UGCS. The cellulase adsorption kinetic experiments were performed for 2, 5, 10, 20, 30, 60, 90, and 120 min with an enzyme loading of 6 mg/g substrate. The predicted values were obtained by current filmporesurface diffusion model. Error bars represent the standard deviation of the measurements for the absorbed cellulase amount in substrate

The filmporesurface diffusion models were developed to simultaneously predict the cellulase adsorption kinetics of both the SGCS and UGCS (Fig. 5 and Table 3). The model prediction agreed reasonably well with the observed kinetic data based on high R 2 and low RMSE values (R 2 = 0.95 and RMSE = 0.16 mg/g for the UGCS, R 2 = 0.93 and RMSE = 0.09 mg/g for the SGCS). Internal diffusion is an important mass transport process during the enzyme adsorption of corn stover particles and includes pore and surface diffusion. The fitted pore diffusion coefficients (D p) were found to be 9.45 × 10−7 cm2/min for the SGCS and 6.04 × 10−6 cm2/min for the UGCS. The magnitude of D p is affected by pore structure parameters, such as the pore size, porosity, and tortuosity. The ultrafine grinding pretreatment can reduce the pore diffusion resistance by changing these pore structure parameters and then enhance the pore diffusion coefficient. Compared with the surface diffusion coefficient (D s) of the SGCS, that of the UGCS was smaller by several orders of magnitude. Surface diffusion is often described by a hopping mechanism in which the migrating particles are viewed as hopping between distinct, energetically favorable adsorption sites on the surface [38]. When an adsorbed particle obtains a sufficient activation energy, it can overcome the energy barrier between adsorption sites and jump to a neighboring site. Thus, the speed of surface diffusion depends on the bond strength of the attached sorption site and the affinity of the recipient site. The equilibrium adsorption data showed that the UGCS had a higher affinity constant (K a) and bonding strength (S) than the SGCS. This may explain the low D s value of the UGCS. External film diffusion is another mass transfer process that is characterized by external-film transfer coefficients (K L). The fitted K L value of the UGCS was much less than that of the SGCS. The relationship between K L and the particle size is not straightforward. Badruzzaman et al. quantified the arsenate adsorption on granular ferric hydroxide by the film-surface diffusion model and then evaluated the K L dependence on the particle size [39]. The results showed that the obtained K L values did not correlate with the particle radius. The magnitude of K L is affected not only by the adsorbent particle size but also by the adsorbent-solution system hydraulics.

Table 3 The kinetic and statistical parameters of the filmporesurface diffusion model fitting

To measure the relative importance of external-film mass transfer to internal-pore mass transfer within the two substrates, the Biot number (B i) was used as an indicator and calculated as B i = K L × R/D e, where K L is the external-film transfer coefficient, R is the particle radius, and D e is the effective diffusion coefficient in the internal pore. Traegner and Suidan suggested that the external-film mass transfer is the rate-controlling step for B i < 1, both the external-film and internal-pore mass transfer are the rate-controlling steps for 1 ≤ B i ≤ 100, and the internal-pore mass transfer is the rate-controlling step for B i > 100 [40]. The calculated B i values for the SGCS and UGCS were 92.81 and 39.57, respectively. Hence, these B i values indicated that both the external-film and internal-pore mass transfer were important for cellulase adsorption on the SGCS and UGCS. However, the smaller B i value of the UGCS indicated a lower internal-pore resistance within the UGCS. The finding showed that the ultrafine grinding pretreatment significantly decreases the particle size and improves the pore diffusion properties, such as the pore size, porosity, pore volume and pore openings, resulting in less internal-pore resistance within the UGCS.

Conclusions

The ultrafine grinding pretreatment was executed on corn stover. The results showed that the ultrafine grinding pretreatment can significantly decrease the particle size (from 218.5 μm of SGCS to 17.45 μm of UGCS) and increase the specific surface area (SSA), pore volume (PV) and surface composition (SSA: from 1.71 m2/g of SGCS to 2.63 m2/g of UGCS, PV: from 0.009 cm3/g of SGCS to 0.024 m3/g of UGCS, cellulose surface area: from 168.69 m2/g of SGCS to 290.76 m2/g of UGCS, lignin surface area: from 91.46 m2/g of SGCS to 106.70 m2/g of UGCS). The structure and surface composition changes induced by ultrafine grinding increase the enzyme adsorption capacity from 2.83 mg/g substrate of SGCS to 5.61 mg/g substrate of UGCS. A filmporesurface diffusion model was developed to simultaneously predict the enzyme adsorption kinetics of both the SGCS and UGCS. The model provided satisfactory predictions based on high R 2 and low RMSE values (R 2 = 0.95 and RMSE = 0.16 mg/g for the UGCS, R 2 = 0.93 and RMSE = 0.09 mg/g for the SGCS). The model was further employed to analyze the rate-limiting steps in the enzyme adsorption process. Although both external-film and internal-pore mass transfer are important for the enzyme adsorption on the SGCS and UGCS, the UGCS has a lower internal-pore resistance compared to the SGCS because the ultrafine grinding pretreatment significantly decreased the particle size and improved the pore diffusion properties such as the pore size, porosity, pore volume and pore openings. These findings identify wherein the probable rate-limiting factors for the enzyme adsorption reside and could, therefore, provide a basis for enhanced cellulose hydrolysis processes.

Methods

Samples and enzyme preparation

Corn stover was collected in 2013 from the Shangzhuang agronomy farm of the China Agricultural University, located in Beijing, China. The corn stover was air dried and milled to coarse particle size (approximately 1–2 cm). Then, it was dried in a forced-air oven at 45 °C for 48 h and milled to a size less than 1 mm in an RT-34 hammer mill (Rong Tsong Precision Technology Co., Taiwan). The milled material was sieved by a JH-300A sieve shaker fitted with a 40-mesh screen to obtain the SGCS samples (Jiahe Machinery Co., Henan province, China). Then, 400 g of powder was further milled using a CJM-SY-B ultrafine vibration grind mill to obtain the UGCS samples (Taiji Ring Nano Products Co., Hebei, China). The corn stover powder was mixed with ZrO2 balls (6–10 mm diameter) in a 1:2 volume ratio for 0.5 h, and the instrument temperature was controlled below 30 °C. All powders obtained were sealed in PVC plastic bags at room temperature before use in all experiments. The celluclast 1.5 L (cellulase) was purchased from Sigma-Aldrich (St. Louis, MO, USA), and the protein content is 36.7 mg/mL.

Analysis of the surface areas of cellulose and lignin

The surface areas of cellulose and lignin were measured according to the literature [21]. The surface areas of cellulose and lignin on the SGCS and UGCS were analyzed by determining the monolayer adsorption maximum of Congo Red (Direct Red 28) [22] and Azure B [23], respectively. Of each adsorption, 100 mg dry material was weighed in 25 mL conical flask; 10 mL of the dye (Congo Red in 30 mM phosphate buffer at pH 6 and Azure B in 50 mM Na-phosphate buffer at pH 7) was added to the conical flasks and incubated for 24 h on a shaker at 200 rpm. Congo Red adsorption was performed at 60 °C and Azure B adsorption at 25 °C. After incubation, the liquid fraction was separated by centrifugation and the supernatant was filtered through a 0.45 μm PTFE filter. The residual dye concentration and reference solutions were determined spectrophotometrically (Congo Red at 498 nm and Azure B at 647 nm) and the amounts of adsorbed dye were calculated. The adsorption experiments were performed in duplicate using Congo Red concentrations of 4, 2, 1, 0.25, 0.05 and 0 g/L and Azure B concentrations of 2, 1, 0.5, 0.25, 0.1 and 0 g/L. The parameters of the adsorption isotherm were fitted to the Langmuir isotherm in MATLAB (Mathworks, Natick, MA, USA). Then, the cellulose surface area was calculated per dry material from the adsorption maximum with 1 g of the adsorbed dye representing 1055 m2 surface [22]. And the surface area of the lignin was obtained from the maximum adsorption capacity and the area (1.297 m2/mg) covered by Azure B [23].

Particle size determination

The particle size distribution was measured using an LS230 laser diffraction particle size analyzer (Beckman Coulter Inc., Miami, FL, USA). The particle measurement range is from 0.375 μm to 2000 μm. Before measurement, the samples were dispersed with distilled water to form a uniform liquid suspension and then were poured into the measurement instrument with ultrasound. An LS v3.29 system based on the Fraunhofer mode was used to measure the particle size.

Specific surface area, pore size and pore volume distribution determination

The specific surface area, pore size and pore volume distribution of the SGCS and UGCS were measured with the Autosorb-iQ porosity analyzer (Quantachrome Instruments, FL, USA). The samples were degassed at 80 °C for 7 h and then cooled in the presence of nitrogen gas under −195 °C, allowing the nitrogen gas to condense on the surfaces and within the pores. The specific surface area was calculated using the Brunauer–Emmett–Teller (BET) model [41], which relates the gas pressures to the volume of gas adsorbed. The pore volume distribution with respect to the pore size was estimated using the Barrett–Joyner–Halenda (BJH) model [42].

Enzyme adsorption kinetic experiments

The enzyme adsorption kinetic experiments were conducted for the SGCS and UGCS with an enzyme loading of 6 mg protein/g substrate, which was among the usual enzyme hydrolysis loading capacity. The lignocellulose substrate-binding studies were performed in centrifuge tubes (10 mL) with a sodium citrate buffer (0.05 M, pH 4.8) using a 1 % (w/v) substrate concentration and incubated for 2, 5, 10, 20, 30, 60, 90, and 120 min in a shaking water bath at 4 °C to avoid hydrolysis. Every experiment was run two times, and substrate blanks without enzyme and enzyme blanks without substrate were also analyzed. After incubation, all samples were centrifuged for 3 min in a refrigerated centrifuge at 6000 rpm. The supernatant was filtered and used to determine the free enzyme by measuring the protein concentration in the supernatant using the Bradford assay by Coomassie brilliant blue dye [43]. The bound enzyme was calculated by subtracting the free enzyme concentration from the initial enzyme concentration loaded.

Equilibrium enzyme adsorption experiments

Different loadings of the enzyme (1.5–10.5 mg/g substrate for celluclast 1.5 L) were performed and incubated for 2 h under the same condition mentioned above. The bound enzyme concentration calculated was correlated with the free enzyme concentration using the following Langmuir equilibrium isotherm:

$$q_{\text{b}} = \frac{{q_{\text{m}} K_{\text{a}} C_{\text{f}} }}{{ 1+ K_{\text{a}} C_{\text{f}} }}$$
(1)

where q b is the equilibrium amount of solid-phase bound enzyme (mg protein/g substrate), q m is the maximum solid-phase bound capacity (mg protein/g substrate), K a is the affinity constant (mL/mg protein), and C f is the equilibrium concentration of free enzyme in solution (mg protein/mL).

The Langmuir adsorption constants (K a and q m) of the SGCS and UGCS were obtained by nonlinear regression using MATLAB (Mathworks, Natick, MA, USA). The binding strength (S in mL/g substrate), another constant from the Langmuir adsorption isotherm, could be used to estimate the stability of the enzyme bound with substrates. The binding strength can be calculated by S = q m × K a.

Filmporesurface diffusion adsorption model

The adsorption of cellulase onto lignocellulosic biomass involves three consecutive steps: external diffusion of the cellulase from the bulk solution across the liquid film surrounding the solid biomass particles, internal diffusion of the cellulase through the biomass particles by pore volume diffusion and surface diffusion, and the adsorption of cellulose molecules onto the biomass particles at the active sites (Fig. 1). The filmporesurface diffusion adsorption model was proposed based on the following assumptions: (a) the adsorbent particles are spherical; (b) the adsorption rate at an active site is instantaneous; and (c) the solute adsorbed amount on the adsorbent can be represented by the Langmuir isotherm equation.

The rate of mass transfer in the external film surrounding the solid particle is assumed to be directly proportional to the concentration difference in the film. Therefore, the external-film mass transfer is given by

$$V_{\text{L}} \frac{{dC_{\text{L}} }}{dt} = - K_{\text{L}} A(C_{\text{L}} - C_{\text{P,r}} \left| {_{\text{r = R}} } \right.)$$
(2)

where t is the adsorption time, V L is the total volume of the liquid phase, C L is the concentration of enzyme in the liquid phase, C P,r|r=R is the enzyme concentration at the particle surface, K L represents the external-film mass transfer coefficient, and A represents the outer surface area of all the particles, estimated as:

$$A = \frac{3m}{{\rho_{\text{a}} R}}$$
(3)

where m is the mass of all the particles, R is the radius of the particle, and ρ a is the apparent density of the particle, estimated as:

$$\rho_{\text{a}} = \frac{{\rho_{\text{s}} }}{{1 + V_{\text{p}} \rho_{\text{s}} }}$$
(4)

where V p is the pore volume per mass of the particle and ρ s is the solid density, estimated as follows:

$$\rho_{\text{s}} = \frac{1}{{M_{\text{c}} /\rho_{\text{c}} + M_{\text{h}} /\rho_{\text{h}} + M_{\text{l}} /\rho_{\text{l}} + M_{\text{o}} /\rho_{\text{o}} }}$$
(5)

where M c, M h, M l, and M o and ρ c, ρ h, ρ l, and ρ o are the mass percentages on a dry basis and the densities of cellulose, hemicellulose, lignin and other compositions in solid particles, respectively.

Based on the mass balance equation for the adsorption of enzyme with internal-pore diffusion in a spherical particle, the following equation can be obtained:

$$\begin{aligned}\varphi \varepsilon \frac{{\partial C_{\text{P,r}} }}{\partial {t}} + \rho_{\text{a}} \frac{{\partial {q}_{\text{r}} }}{\partial {t}} &= \frac{1}{{r^{2} }}\frac{\partial }{\partial {r}}\left[ {r^{2} {D}_{\text{p}} \frac{{\partial {C}_{\text{P,r}} }}{\partial {r}}} \right] \\ & \quad + \frac{1}{{r^{2} }}\frac{\partial }{\partial {r}}\left[ {r^{2} \rho_{\text{a}} D_{\text{s}} \frac{{\partial {q}_{\text{r}} }}{\partial {r}}} \right] \end{aligned}$$
(6)

where C P,r is the enzyme concentration in the particle pores at position r, \(\varphi\) is the ratio of the accessible pore volume to the enzyme (V pa) to the total pore volume (V p), r is the radial position in the particle, q r is the solid-phase enzyme adsorption amount at position r, D p is the pore diffusion coefficient of the enzyme, D s is the surface diffusion coefficient of the enzyme, and ε is the porosity of the solid particle, estimated as:

$$\varepsilon = \frac{{V_{\text{p}} }}{{V_{\text{p}} + 1/\rho_{\text{s}} }}$$
(7)

As the adsorption step occurs much more rapidly than the mass transfer step in physical adsorption, the pore solution concentration and the solid-phase adsorbed amount can be expressed by the Langmuir isotherm equation:

$$q_{\text{r}} = f(C_{\text{P,r}} ) = \frac{{q_{\text{m}} K{}_{\text{a}}C_{\text{P,r}} }}{{1 + K_{\text{a}} C_{\text{P,r}} }}$$
(8)

Differentiating Eq. (8) yields:

$$dq_{\text{r}} = f^{\prime } (C_{\text{P,r}} )dC_{\text{P,r}}$$
(9)

Substituting Eq. (9) into Eq. (6) gives:

$$\left[ {\varphi \varepsilon + \rho_{\text{a}} f^{\prime } (C_{\text{P,r}} )} \right]\frac{{\partial C_{\text{P,r}} }}{\partial t} = \frac{1}{{r^{2} }}\frac{\partial }{\partial r}\left[ {r^{2} D_{\text{e}} \frac{{\partial C_{\text{P,r}} }}{\partial r}} \right]$$
(10)

where D e is the effective diffusion coefficient in the internal pore, given as:

$$D_{\text{e}} = D_{p} + f^{\prime}(C_{\text{P,r}} )\rho_{\text{a}} D_{\text{s}}$$
(11)

Substituting \(f^{'} (C_{\text{P,r}} ) = \frac{{q_{\text{m}} K_{\text{a}} }}{{(1 + K_{\text{a}} C_{\text{P,r}} )^{2} }}\) into Eq. (10) gives:

$$\begin{aligned} \left( {\varphi \varepsilon + \frac{{\rho_{\text{a}} q_{\text{m}} K_{\text{a}} }}{{(1 + K_{\text{a}} C_{\text{P,r}} )^{2} }}} \right)\frac{{\partial C_{\text{P,r}} }}{\partial t} &= \hfill \left( {D_{\text{p}} + \frac{{D_{\text{s}} \rho_{\text{a}} q_{\text{m}} K_{\text{a}} }}{{(1 + K_{\text{a}} C_{\text{P,r}} )^{2} }}} \right)\frac{{\partial^{2} C_{\text{P,r}} }}{{\partial r^{2} }} \\ & \quad- \frac{{2q_{\text{m}} K_{\text{a}}^{2} D_{\text{s}} \rho_{\text{a}} }}{{(1 + K_{\text{a}} C_{\text{P,r}} )^{3} }}\left( {\frac{{\partial C_{\text{P,r}} }}{\partial r}} \right)^{2} + \frac{2}{r}\left( {D_{\text{p}} + \frac{{D_{\text{s}} \rho_{\text{a}} q_{\text{m}} K_{\text{a}} }}{{(1 + K_{\text{a}} C_{\text{P,r}} )^{2} }}} \right)\frac{{\partial C_{\text{P,r}} }}{\partial r} \hfill \\ \end{aligned}$$
(12)

The average enzyme adsorption amount in the solid particles (q a) is given by

$$q_{\text{a}} = \frac{{\int_{0}^{R} {q_{\text{r}} \times } 4\pi r^{2} dr}}{{\frac{4}{3}\pi R^{3} }} = \frac{3}{{R^{3} }}\int_{0}^{R} {q_{\text{r}} r^{2} } dr$$
(13)

The initial and boundary conditions are listed as:

$$t = 0\begin{array}{*{20}c} {} & {C_{\text{L}} = C_{0} } \\ \end{array} \begin{array}{*{20}c} {} & {\begin{array}{*{20}c} {C_{\text{p,r}} \left| {_{0 \le r \le R} } \right. = 0} & {} \\ \end{array} } \\ \end{array}$$
(14)
$$\begin{aligned} t > 0 & {K_{\text{L}} (C_{\text{L}} - C_{\text{P,r}} \left| {_{\text{r = R}} } \right. )= D_{\text{e}} \frac{{\partial C_{\text{P,r}} }}{\partial r}\text{|}_{\text{r = R}} } \\&\quad {\frac{{\partial C_{\text{P,r}} }}{\partial r}\text{|}_{\text{r = 0}} = 0} \end{aligned}$$
(15)

The filmporesurface diffusion model can be numerically solved by combining Eqs. (215). The C P,r values can be used to calculate the average enzyme adsorption amount in the solid particles (q a) according to Eqs. (8) and (13). The predicted q a values were compared with the observed values and were used to estimate the model parameters. The model parameters (K L, D p, and D s) were simultaneously fitted to all experimental data using a custom-written program in MATLAB (Mathworks, Natick, MA, USA). Table 4 provides a description of these symbols.

Table 4 The model symbols

Abbreviations

RMSE:

root mean squared error

SGCS:

sieve-based grinding corn stover

UGCS:

ultrafine grinding corn stover

PV:

pore volume

SSA:

specific surface area

References

  1. Council WE, Wales RiEa. World Energy Resources 2013 Survey. 2013. p. 7.1–7.24.

  2. Gao XD, Kumar R, Singh S, Simmons BA, Balan V, Dale BE, et al. Comparison of enzymatic reactivity of corn stover solids prepared by dilute acid, AFEX (TM), and ionic liquid pretreatments. Biotechnol Biofuels. 2014;7:71. doi:10.1186/1754-6834-7-71.

    Article  Google Scholar 

  3. Chen L, Zhang H, Li J, Lu M, Guo X, Han L. A novel diffusion-biphasic hydrolysis coupled kinetic model for dilute sulfuric acid pretreatment of corn stover. Bioresour Technol. 2015;177:8–16. doi:10.1016/j.biortech.2014.11.060.

    Article  CAS  Google Scholar 

  4. Shen J, Wyman CE. A novel mechanism and kinetic model to explain enhanced xylose yields from dilute sulfuric acid compared to hydrothermal pretreatment of corn stover. Bioresour Technol. 2011;102(19):9111–20. doi:10.1016/j.biortech.2011.04.001.

    Article  CAS  Google Scholar 

  5. Cabrera E, Munoz MJ, Martin R, Caro I, Curbelo C, Diaz AB. Alkaline and alkaline peroxide pretreatments at mild temperature to enhance enzymatic hydrolysis of rice hulls and straw. Bioresour Technol. 2014;167:1–7. doi:10.1016/j.biortech.2014.05.103.

    Article  CAS  Google Scholar 

  6. Pierre G, Sannier F, Goude R, Nouviaire A, Maache-Rezzoug Z, Rezzoug SA, et al. Evaluation of thermomechanical pretreatment for enzymatic hydrolysis of pure microcrystalline cellulose and cellulose from Brewers’ spent grain. J Cereal Sci. 2011;54(3):305–10. doi:10.1016/j.jcs.2011.06.004.

    Article  CAS  Google Scholar 

  7. Sun Y, Cheng JY. Hydrolysis of lignocellulosic materials for ethanol production: a review. Bioresour Technol. 2002;83(1):1–11. doi:10.1016/S0960-8524(01)00212-7.

    Article  CAS  Google Scholar 

  8. Zhu JY, Wang GS, Pan XJ, Gleisner R. Specific surface to evaluate the efficiencies of milling and pretreatment of wood for enzymatic saccharification. Chem Eng Sci. 2009;64(3):474–85. doi:10.1016/j.ces.2008.09.026.

    Article  CAS  Google Scholar 

  9. Adapa P, Tabil L, Schoenau G. Grinding performance and physical properties of non-treated and steam exploded barley, canola, oat and wheat straw. Biomass Bioenergy. 2011;35(1):549–61. doi:10.1016/j.biombioe.2010.10.004.

    Article  CAS  Google Scholar 

  10. Chen HZ. Biotechnology of lignocellulose: theory and practice. Dordrecht: Springer; 2014.

    Book  Google Scholar 

  11. Silva GG, Couturier M, Berrin JG, Buleon A, Rouau X. Effects of grinding processes on enzymatic degradation of wheat straw. Bioresour Technol. 2012;103(1):192–200. doi:10.1016/j.biortech.2011.09.073.

    Article  CAS  Google Scholar 

  12. Gan Q, Allen SJ, Taylor G. Kinetic dynamics in heterogeneous enzymatic hydrolysis of cellulose: an overview, an experimental study and mathematical modelling. Process Biochem. 2003;38(7):1003–18. doi:10.1016/S0032-9592(02)00220-0.

    Article  CAS  Google Scholar 

  13. Machado DL, MoreiraNeto J, da Cruz Pradella JG, Bonomi A, Rabelo SC, daCosta AC. Adsorption characteristics of cellulase and β-glucosidase on Avicel, pretreated sugarcane bagasse, and lignin. Biotechnol Appl Biochem. 2015. doi:10.1002/bab.1307.

    Google Scholar 

  14. Ji GY, Gao CF, Xiao WH, Han LJ. Mechanical fragmentation of corncob at different plant scales: impact and mechanism on microstructure features and enzymatic hydrolysis. Bioresour Technol. 2016;205:159–65. doi:10.1016/j.biortech.2016.01.029.

    Article  CAS  Google Scholar 

  15. Zheng Y, Pan ZL, Zhang RH, Jenkins BM. Kinetic modeling for enzymatic hydrolysis of pretreated creeping wild ryegrass. Biotechnol Bioeng. 2009;102(6):1558–69. doi:10.1002/bit.22197.

    Article  CAS  Google Scholar 

  16. Börjesson J, Engqvist M, Sipos B, Tjerneld F. Effect of poly(ethylene glycol) on enzymatic hydrolysis and adsorption of cellulase enzymes to pretreated lignocellulose. Enzyme Microb Tech. 2007;41(1–2):186–95. doi:10.1016/j.enzmictec.2007.01.003.

    Article  Google Scholar 

  17. Bansal P, Hall M, Realff MJ, Lee JH, Bommarius AS. Modeling cellulase kinetics on lignocellulosic substrates. Biotechnol Adv. 2009;27(6):833–48. doi:10.1016/j.biotechadv.2009.06.005.

    Article  CAS  Google Scholar 

  18. Li C, Cheng G, Balan V, Kent MS, Ong M, Chundawat SPS, et al. Influence of physico-chemical changes on enzymatic digestibility of ionic liquid and AFEX pretreated corn stover. Bioresour Technol. 2011;102(13):6928–36. doi:10.1016/j.biortech.2011.04.005.

    Article  CAS  Google Scholar 

  19. Saha BC, Qureshi N, Kennedy GJ, Cotta MA. Enhancement of xylose utilization from corn stover by a recombinant Escherichia coli strain for ethanol production. Bioresour Technol. 2015;190:182–8. doi:10.1016/j.biortech.2015.04.079.

    Article  CAS  Google Scholar 

  20. Kumar R, Wyman CE. An improved method to directly estimate cellulase adsorption on biomass solids. Enzyme Microb Tech. 2008;42(5):426–33. doi:10.1016/j.enzmictec.2007.12.005.

    Article  CAS  Google Scholar 

  21. Pihlajaniemi V, Sipponen MH, Kallioinen A, Nyyssola A, Laakso S. Rate-constraining changes in surface properties, porosity and hydrolysis kinetics of lignocellulose in the course of enzymatic saccharification. Biotechnol Biofuels. 2016;9:18. doi:10.1186/s13068-016-0431-3.

    Article  Google Scholar 

  22. Wiman M, Dienes D, Hansen MAT, van der Meulen T, Zacchi G, Liden G. Cellulose accessibility determines the rate of enzymatic hydrolysis of steam-pretreated spruce. Bioresour Technol. 2012;126:208–15. doi:10.1016/j.biortech.2012.08.082.

    Article  CAS  Google Scholar 

  23. Sipponen MH, Pihlajaniemi V, Littunen K, Pastinen O, Laakso S. Determination of surface-accessible acidic hydroxyls and surface area of lignin by cationic dye adsorption. Bioresour Technol. 2014;169:80–7. doi:10.1016/j.biortech.2014.06.073.

    Article  CAS  Google Scholar 

  24. Yang Y, Ji G, Xiao W, Han L. Changes to the physicochemical characteristics of wheat straw by mechanical ultrafine grinding. Cellulose. 2014;21(5):3257–68. doi:10.1007/s10570-014-0381-5.

    Article  Google Scholar 

  25. Barakat A, Mayer-Laigle C, Solhy A, Arancon RAD, de Vries H, Luque R. Mechanical pretreatments of lignocellulosic biomass: towards facile and environmentally sound technologies for biofuels production. RSC Adv. 2014;4(89):48109–27. doi:10.1039/c4ra07568d.

    Article  CAS  Google Scholar 

  26. Fougere JD, Lynch M, Zhao J, Zheng Y, Li K. Impact of mechanical downsizing on the physical structure and enzymatic digestibility of pretreated hardwood. Energy Fuels. 2014;28(4):2645–53. doi:10.1021/ef5001387.

    Article  CAS  Google Scholar 

  27. Lin ZX, Huang H, Zhang HM, Zhang L, Yan LS, Chen JW. Ball milling pretreatment of corn stover for enhancing the efficiency of enzymatic hydrolysis. Appl Biochem Biotechnol. 2010;162(7):1872–80. doi:10.1007/s12010-010-8965-5.

    Article  CAS  Google Scholar 

  28. Piccolo C, Wiman M, Bezzo F, Lidén G. Enzyme adsorption on SO2 catalyzed steam-pretreated wheat and spruce material. Enzyme Microb Tech. 2010;46(3–4):159–69. doi:10.1016/j.enzmictec.2009.11.007.

    Article  CAS  Google Scholar 

  29. Zhang W, Liang M, Lu C. Morphological and structural development of hardwood cellulose during mechanochemical pretreatment in solid state through pan-milling. Cellulose. 2007;14(5):447–56. doi:10.1007/s10570-007-9135-y.

    Article  CAS  Google Scholar 

  30. Huang RL, Su RX, Qi W, He ZM. Understanding the key factors for enzymatic conversion of pretreated lignocellulose by partial least square analysis. Biotechnol Prog. 2010;26(2):384–92. doi:10.1002/btpr.324.

    CAS  Google Scholar 

  31. Wang QQ, He Z, Zhu Z, Zhang YHP, Ni Y, Luo XL, et al. Evaluations of cellulose accessibilities of lignocelluloses by solute exclusion and protein adsorption techniques. Biotechnol Bioeng. 2012;109(2):381–9. doi:10.1002/bit.23330.

    Article  CAS  Google Scholar 

  32. Grethlein HE. The effect of pore-size distribution on the rate of enzymatic-hydrolysis of cellulosic substrates. Nat Biotechnol. 1985;3(2):155–60. doi:10.1038/nbt0285-155.

    Article  CAS  Google Scholar 

  33. Kumar R, Wyman CE. Cellulase adsorption and relationship to features of corn stover solids produced by leading pretreatments. Biotechnol Bioeng. 2009;103(2):252–67. doi:10.1002/Bit.22258.

    Article  CAS  Google Scholar 

  34. Pierre G, Maache-Rezzoug Z, Sannier F, Rezzoug SA, Maugard T. High-performance hydrolysis of wheat straw using cellulase and thermomechanical pretreatment. Process Biochem. 2011;46(11):2194–200. doi:10.1016/j.procbio.2011.09.002.

    Article  CAS  Google Scholar 

  35. Qi B, Chen X, Su Y, Wan Y. Enzyme adsorption and recycling during hydrolysis of wheat straw lignocellulose. Bioresour Technol. 2011;102(3):2881–9. doi:10.1016/j.biortech.2010.10.092.

    Article  CAS  Google Scholar 

  36. Zhang YHP, Lynd LR. Toward an aggregated understanding of enzymatic hydrolysis of cellulose: noncomplexed cellulase systems. Biotechnol Bioeng. 2004;88(7):797–824. doi:10.1002/bit.20282.

    Article  CAS  Google Scholar 

  37. Shih-Yuan Chen Y-TC, Jey-Jau Leec, Soofin Cheng. Tuning pore diameter of platelet SBA-15 materials with short mesochannels for enzyme adsorption. J Mater Chem. 2011;21:5693–703.

    Article  Google Scholar 

  38. Medved I, Cerny R. Surface diffusion in porous media: a critical review. Microporous Mesoporous Mat. 2011;142(2–3):405–22. doi:10.1016/j.micromeso.2011.01.015.

    Article  CAS  Google Scholar 

  39. Badruzzaman M, Westerhoff P, Knappe DRU. Intraparticle diffusion and adsorption of arsenate onto granular ferric hydroxide (GFH). Water Res. 2004;38(18):4002–12. doi:10.1016/j.watres.2004.07.007.

    Article  CAS  Google Scholar 

  40. Traegner UK, Suidan MT. Parameter evaluation for carbon adsorption. J Environ Eng. 1989;115(1):109–28. doi:10.1061/(ASCE)0733-9372(1989)115:1(109).

    Article  CAS  Google Scholar 

  41. Brunauer S, Emmett PH, Teller E. Adsorption of gases in multimolecular layers. J Am Chem Soc. 1938;60:309–19. doi:10.1021/Ja01269a023.

    Article  CAS  Google Scholar 

  42. Barrett EP, Joyner LG, Halenda PP. The determination of pore volume and area distributions in porous substances 1 Computations from nitrogen isotherms. J Am Chem Soc. 1951;1:373–80. doi:10.1021/Ja01145a126.

    Article  Google Scholar 

  43. Bradford MM. Rapid and sensitive method for quantitation of microgram quantities of protein utilizing principle of protein-dye binding. Anal Biochem. 1976;72(1–2):248–54. doi:10.1006/abio.1976.9999.

    Article  CAS  Google Scholar 

  44. Ehrnrooth EML. Change in pulp fiber density with acid-chlorite delignification. J Wood Chem Technol. 1984;4(1):91–109. doi:10.1080/02773818408062285.

    Article  CAS  Google Scholar 

  45. Richard TL, Veeken AHM, de Wilde V, Hamelers HVM. Air-filled porosity and permeability relationships during solid-state fermentation. Biotechnol Prog. 2004;20(5):1372–81. doi:10.1021/bp0499505.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

HYZ performed the experiments, analyzed the experimental results, and wrote the manuscript. LJC conceived this study, analyzed the experimental results, proposed the model, and wrote the manuscript. MSL and JBL assisted in the experimental analysis. LJC and LJH coordinated the overall study. All authors read and approved the manuscript.

Acknowledgements

The authors would like to thank the reviewers for positive criticism to improve the quality of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Availability of supporting data

Not applicable.

Consent for publication

Not applicable.

Ethical approval and consent to participate

Not applicable.

Funding

This research was supported by the National Natural Science Foundation of China (Project No. 31571569), the Beijing Nova Program (Project No. Z131105000413056), the Beijing Youth Talent Plan Program in University (Project No. YETP0317), and the National High-level Personnel of Special Support Program for Outstanding Young Talents (Project No. Zutingzi 2015-48).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Longjian Chen.

Rights and permissions

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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, H., Chen, L., Lu, M. et al. A novel film–pore–surface diffusion model to explain the enhanced enzyme adsorption of corn stover pretreated by ultrafine grinding. Biotechnol Biofuels 9, 181 (2016). https://doi.org/10.1186/s13068-016-0602-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-016-0602-2

Keywords