Subscriber access provided by UNIVERSITY OF SASKATCHEWAN LIBRARY
Elucidating the interactive impacts of substrate-related properties on lignocellulosic biomass digestibility: A sequential analysis Liyuan Chai, Mingren Liu, Xu Yan, Xunqiang Cheng, Tingzheng Zhang, Mengying Si, Xiao-Bo Min, and Yan Shi ACS Sustainable Chem. Eng., Just Accepted Manuscript • DOI: 10.1021/ acssuschemeng.8b00592 • Publication Date (Web): 18 Mar 2018 Downloaded from http://pubs.acs.org on March 19, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
1
Elucidating the interactive impacts of substrate-related properties on
2
lignocellulosic biomass digestibility: A sequential analysis
3 4
Liyuan Chai1, 2, Mingren Liu 1, Xu Yan1, 2, Xunqiang Cheng1, Tingzheng Zhang1,
5
Mengying Si1, Xiaobo Min1, 2, Yan Shi*1, 2
6
1
7 8 9
School of Metallurgy and Environment, Central South University, Changsha 410083, China
2
Chinese National Engineering Research Center for Control & Treatment of Heavy Metal Pollution, Changsha 410083, China
10
*To whom correspondence should be addressed. Y. Shi, E-mail:
[email protected] 11
(Y. Shi); Fax: +86-0731-88710171; Tel: +86-0731-88830875
12
L. Chai, M. Liu, X. Yan, X. Cheng, T. Zhang, M. Si, X. Min, Y. Shi, Mailing address:
13
No.932 South Lushan Road, Changsha Hunan 410083, P.R. China
14 15
Key words: Enzymatic hydrolysis; Lignin; Cellulose; NMR spectroscopy;
16
Lignocellulose biomass; Physicochemical properties; Statistical analysis
1 ACS Paragon Plus Environment
ACS Sustainable Chemistry & Engineering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
17
ABSTRACT
18
A lack of insight into interactive effects among substrate-related factors holds
19
back the determination of dominant factors in efficient sugar conversion. Herein,
20
thirteen factors defining compositional and physicochemical properties of
21
lignocellulose pretreated by dilute acid/base and enzymes were analyzed through an
22
innovative sequence of correlation analysis, principle component analysis, multiple
23
linear regression and multiscale statistical validation. Results showed that the lignin
24
content, cellulose content and O/C ratio principally affected enzymatic hydrolysis.
25
The dominant role was played by the lignin content due to its major recalcitrance
26
providing to biomass and concomitant impacts on surface lignin and porosity
27
properties. The structural features of lignin played a less pronounced role with high
28
lignin content remained. Besides, the sequential analysis revealed different inhibition
29
mechanisms for glucan and carbohydrate conversion, i.e., non-productive binding of
30
enzymes and steric hindrance of lignin, respectively. The established weighing order
31
of interactive factors enlightens more efficient pretreatment strategies.
2 ACS Paragon Plus Environment
Page 2 of 36
Page 3 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
32
INTRODUCTION
33
The efficient bioconversion of lignocellulosic carbohydrates into sugars is
34
essential for the valorization of lignocellulose. 1 To overcome the innate recalcitrance
35
of lignocellulose, different pretreatment strategies have been developed to open the
36
substrate structure for cellulose degradation.2-3 The efforts to elucidate their
37
underlying mechanisms have been made by investigating the concomitant changes in
38
substrate-related factors, such as biomass porosity, lignin/hemicellulose distribution,
39
contents and structural features as well as cellulose crystallinity and accessibility.2-6
40
Unfortunately, the interactive changes in the biomass during pretreatments often lead
41
to conflicting trends in weighting the dominant factors.1-5,
42
influencing mechanisms of these factors thus presents a significant challenge and the
43
development of a strategy that optimizes biomass digestibility remains hindered.
7-12
Interpreting the
44
The interactions among these substrate factors are mainly affected by the
45
distribution, content and composition of lignin and hemicellulose in the highly
46
dynamic lignocellulosic structure.4, 6, 12 According to a popular morphological model,
47
highly ordered cellulose elemental fibrils are cross-linked with hemicellulose and
48
embedded in a non-cellulosic polysaccharide matrix. The lignin, as the plasticizer, is
49
partially bonded with hemicellulose via lignin-carbohydrate complexes.1, 7-8 Therefore,
50
the negative role of lignin in lignocellulose digestion has been long recognized.
51
Biomass pretreatments that reduce the lignin content are thought to improve the
52
saccharification efficiency.3, 5, 9 Contradicting this view, the sugar release from poplar 3 ACS Paragon Plus Environment
ACS Sustainable Chemistry & Engineering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 36
53
samples with high syringyl /guaiacyl (S/G) ratio shows less dependency on lignin
54
content.9 This study reported a tangled interaction between the lignin structural
55
features and the lignin content. The increasing S/G ratio implies less cross-linking in
56
the triaxial structure and abates the negative role of the lignin content9, 13 whereas it
57
was also reported the S/G ratio contributes relatively little to sugar release.
58
Moreover, phenolic compounds of lignin should aggravate the non-productive binding
59
of cellulase but their counteraction effect with alcoholic and carboxylic hydroxyls
60
remains unclear.4, 15 Consequently, no clear picture of lignin inhibition has emerged.
61
Similar obscureness was also displayed on the impacts of hemicellulose. Cellulose
62
accessibility proves linearly proportional to hemicellulose removal.3,
63
cellulase reaction rate is also linearly related to the pore volume of the biomass3, 8, 10
64
but followed by a rapid decline in the conversion rate.16 After the initial hydrolysis of
65
hemicellulose and amorphous cellulose, biomass porosity increased17 but the
66
concomitant increase in crystalline region and lignin content also hinders the
67
enzymatic hydrolysis.18-19 Therefore, all substrate-related factors are mechanically
68
important, but the factor playing the dominant role remains unidentified.
10
14
The initial
69
The statistical analyses of these arguments have been developed over many years.
70
Pairs of linear relationships between factors and biomass digestibility are commonly
71
tested by correlation analysis (CA). Unfortunately, CA alone often leads to conflicting
72
or confusing trends.3, 5, 9-10, 16-19 In particular, CA is susceptible to database diversity; a
73
small sample size or individual data from a specific pretreatment may undermine the 4 ACS Paragon Plus Environment
Page 5 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
74
accuracy of the correlation results and the reliability of further interpretation. Worse
75
still, CA alone cannot handle the tangled interactions among different factors. In some
76
studies, statistical analysis has been complemented with principle component analysis
77
(PCA) which reduces the dimensionality of the problem.5, 20 However, current studies
78
have focused on a limited scale of factors without removing the irrelevant ones.5, 20
79
This reduces the robustness of the total variance interpretation and compromises the
80
PCA performance. For a more integrated perspective, a database providing diverse
81
factors after pretreatments should be combined with a systematic and mechanical
82
statistical method.
83
In this study, thirteen factors were comprehensively collected after dilute
84
acid/base and combined enzyme-chemical pretreatments. These factors included the
85
compositional variances, substrate surface properties, cellulose crystallinity, biomass
86
porosity and structural features of lignin. The latter pretreatment introduced laccase
87
and endo-xylanase for facile bioprocessing and diversified the physicochemical
88
changes.21 To identify the dominant factors of biomass digestibility, the interactions
89
between the thirteen factors were analyzed through a sequential analysis which
90
integrates CA, PCA and stepwise multiple linear regression (MLR).22-23 Analogous to
91
impurity removal in cascade filtration, the method excludes the factors that are
92
irrelevant or subjected to the collinearity problem, thus guaranteeing that no chance
93
correlations or collinearity will undermine the accuracy of the MLR results. To our
94
knowledge, our study is the first sequential study on the complex influence of 5 ACS Paragon Plus Environment
ACS Sustainable Chemistry & Engineering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
95
interacting factors on enzymatic hydrolysis. We thus provide a systematic and reliable
96
approach for optimizing pretreatment strategies.
97
EXPERIMENTAL SECTION
98
Pretreatment and enzymatic hydrolysis
99
Corn stover ( 0.6, p < 0.05, variance inflation factor < 10 (for evaluating
185
the collinearity level) and F-statistic > 15.22 The error in the model during its
186
application was then estimated by both internal validation (leave-one-out method,
187
LOO; bootstrap method, BOOT) and experimental test (Table S11). The
188
cross-validation coefficient Q2 was defined as follows:
189
Q2 =1 −
∑ (yi ŷi )2
(2)
2 ∑ (yi y) i
190
where yi and ŷi were the detected and predicted glucan/carbohydrate conversion,
191
respectively. yi was the average value in dataset. The BOOT was performed on 5000
192
randomly sampled subsets. The Q2 and the root mean square error (RMSE) of LOO
193
2 and BOOT method were also calculated. A model was considered robust if its QLOO
194
2 and QBOOT exceeded 0.5 and RMSE less than 0.3. All statistical analyses were
195
performed in SPSS software (Version 21.0, SPSS Inc., USA), the R programming
196
language (3.4.1, USA) and MATLAB (2014a, MathWorks, USA).
197
RESULTS AND DISCUSSION 10 ACS Paragon Plus Environment
Page 11 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
198
Correlation analysis of thirteen substrate-related properties with potential
199
influence on biomass digestibility
200
Thermochemical pretreatment with dilute acid/base as well as facile bioprocessing
201
with laccase and xylanase digestion were employed in present study, providing
202
diverse changes in the substrate-related factors. 21 For simplicity, we hereafter refer to
203
the dilute acid and acid-enzyme pretreatments as acid pretreatment (a similar
204
terminology is adopted for base pretreatment; Table S1). Quantile−quantile plots of
205
glucan and carbohydrate conversions were performed. All plots were distributed in
206
the diagonal line, confirming the validity of a random normal distribution assumption
207
(Figure S1).22 Figures 1 and 2, respectively, present the matrix of correlation
208
coefficients among factors and the relations between glucan/carbohydrate conversions
209
and the thirteen factors. The data falls into three classes: acid pretreated data, base
210
pretreated data and untreated data.
211
Compositional analysis
212
The acid/base pretreatment induced delignification and hemicellulose removal,
213
with concomitant changes in the physicochemical factors of the biomass substrate.
214
Herein, the lignin content varied from 1.61% to 26.2% and was strongly negatively
215
correlated with both glucan (r = −0.830, P < 0.05) and carbohydrate conversions (r =
216
−0.866, P < 0.05, Figures 1 and 2a). It’s observed that the glucan and carbohydrate
217
conversion of untreated and acid-pretreated samples were less than 50%, whereas 11 ACS Paragon Plus Environment
ACS Sustainable Chemistry & Engineering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
218
those of base-pretreated samples exceeded 50% (Table S2). Lignin forms not only
219
hydrophobic networks that restrict the entry of cellulases and hemicellulases into the
220
polysaccharide1, 3 but also LCCs that non-productively bind with the enzymes.10, 32
221
Residual lignin after pretreatment could further block the progress of cellulase down
222
the cellulose chain.33 Therefore, the lignin content exerted a pronounced adverse
223
effect on the enzymatic digestibility.
224
By contrast, the hemicellulose content was not significantly correlated with glucan
225
(r = 0.163, P > 0.05) or carbohydrate conversions (r = 0.151, P > 0.05, Figures 1 and
226
2b); this finding is consistent with a previous CA of diversely pretreated samples.20
227
Although hemicellulose removal improved the access of cellulases to cellulose, its
228
impact on biomass digestibility might be less important. Hemicellulose obstacles are
229
frequently overcome by adding accessory enzymes such as xylanase in Cellic®
230
CTec2.34 With deacetylation during lignin removal, the inhibition effect of acetyl
231
groups on endo-xylanases is also reduced.
232
hydrolysis of hemicellulose is thus achieved. Due to the low correlation coefficients,
233
the hemicellulose content was excluded from the PCA step in the sequential study.
11, 19, 35
A more efficient enzymatic
234
The cellulose content was positively correlated with glucan (r = 0.636, P < 0.05)
235
and carbohydrate conversion (r = 0.638, P < 0.05, Figures 1 and 2c). Given that the
236
acid/base pretreatment only moderately affects the cellulose itself, cellulose content
237
mainly varies by concomitant mass reduction during lignin and/or hemicellulose
238
removal.
2, 5, 36
Specifically, closed regression coefficients of lignin and hemicellulose 12 ACS Paragon Plus Environment
Page 12 of 36
Page 13 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
239
contents to the increase in cellulose contents were obtained (-0.819 vs. -0.890, Table
240
S9). Therefore, the positive relation between cellulose content and glucan and
241
carbohydrate conversions resulted from a combined effect of lignin and hemicellulose
242
removal.
243
Surface O/C ratio
244
Besides the lignin in bulk substrate, the surface lignin increases the adhesion
245
forces to cellulase, thus aggravating non-productive binding of enzymes.27, 37 Herein,
246
the impact of the surface-lignin coverage was identified from the ratio of surface
247
oxygen (O1s) to carbon (C1s) obtained by XPS.38 Generally, a higher O/C ratio
248
indicates that the reduced-oxygen components (such as lignin, 0.33) are less
249
distributed than the oxygen-rich components (such as cellulose, 0.8339).15 For
250
confirmation, XPS C1s peaks were deconvoluted. The dominant peak at 284.8 ± 0.1
251
eV (C1) was corresponded to the C-C bonds from lignin structure.27 Notably, the O/C
252
ratio was negatively correlated with the content of C1 (r = −0.581, P < 0.05, Table S4),
253
confirming that lowering surface lignin coverage could increase the surface O/C ratio.
254
As presented in Table S4 and Figure 2d, the O/C ratio of lignocellulosic biomass
255
varied from 0.371 to 0.658 and was positively correlated with both glucan (r = 0.783,
256
P < 0.05) and carbohydrate conversion (r = 0.742, P < 0.05), revealing that lower
257
surface lignin coverage contributed to biomass enzymatic digestibility. Indeed, a low
258
surface lignin coverage was highly correlated with lignin content (r= −0.666, P
0.05, respectively; Figures 1 and 2e). In fact, cellulose
271
crystallinity affected by pretreatment is often less straightforward to enzymatic
272
hydrolysis in a complicated biomass structure.2, 40-41 The good correlation between CrI
273
and enzymatic digestibility is often obtained from nearly pure cellulose.42 Therefore,
274
some studies have associated the CrI with actual cellulose content and constructed a
275
new parameter (CrI/cellulose) to assess the actual changes in crystallinity.2,
276
Interestingly, the untreated biomass possessed the lowest CrI value but the highest
277
CrI/cellulose (1.411, Table S3) and the Cr/cellulose values were lower in the
278
base-pretreated samples. The base pretreatment led to cellulose swelling which
14 ACS Paragon Plus Environment
5
Page 15 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
279
decreased the crystallinity and partially transformed the recalcitrant cellulose Iα to
280
amorphous-like cellulose IIII, thus boosting the enzymatic hydrolysis.5,
281
suggests that the total crystallinity is mainly increased with the increase in the
282
cellulose content, although the crystallinity of the cellulose itself deceases.44-45
283
Correspondingly, the correlation between the CrI/cellulose and glucan/carbohydrate
284
conversions was improved (r = −0.548/−0.483, P < 0.05, Figures 1 and 2f). In the
285
next stage of PCA, CrI/cellulose was preserved as the glucan-conversion-relevant
286
factor.
287
Biomass Porosity
7, 43
This
288
Considering that hydrolysis requires intimate contact between the enzymes and a
289
valid polysaccharide surface, biomass porosity should be considered.7 Herein, the
290
biomass porosity was estimated from the specific surface area (SSA), pore volume
291
(PV) and average pore sizes (APS) determined by nitrogen adsorption. Unexpectedly,
292
the glucan and carbohydrate conversions were poorly correlated with SSA (r =−0.218
293
and −0.268, P > 0.05, respectively; Figures 1 and 2g) and with PV (r = −0.270 and
294
−0.302, P > 0.05, respectively; Figures 1 and 2h). We then separately analyzed the
295
acid and base pretreatment data and yielded very different results. In acid pretreatment,
296
the porosity of biomass substrate is significantly promoted by hemicellulose removal
297
and lignin rearrangement.3 In present study, the acid pretreatment increased the SSA
298
from 0.623 m2 g-1 (untreated) to 3.41–7.72 m2 g-1, and the PV from 0.00242 m3 g-1 to
15 ACS Paragon Plus Environment
ACS Sustainable Chemistry & Engineering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 16 of 36
299
0.0147–0.0280 m3 g-1 (Figure 2g, h). The larger SSA and PV indicates a more porous
300
structure, facilitating enzyme access to the embedded cellulosic microfibrils.2,
301
Therefore, the glucan and carbohydrate conversions were strongly positively
302
correlated with SSA (r = 0.861 and 0.754, P < 0.05, respectively) and with PV (r =
303
0.848 and 0.814, P < 0.05, respectively). However, the base-pretreated biomass was
304
less porous (with SSA and PV ranging from 1.22 to 1.97 m2 g-1 and from 0.00335 to
305
0.00715 m3 g-1, respectively; Figure 2g, h). The effect of delignification on cellulose
306
accessibility is often limited.3, 46 By partially reducing the mechanical support of the
307
whole biomass, delignification might induce the partial aggregation of adjacent
308
cellulose microfibrils.2,
309
between glucan conversion/carbohydrate conversion and SSA (r = −0.101/−0.107, P >
310
0.05) and PV (r = −0.129/−0.135, P > 0.05) after the base pretreatment. Notably,
311
although polysaccharide digestion increases the structural porosity, the enzymatic
312
hydrolysis rate still drops.45, 46 As enzymatic hydrolysis slightly increases the CrI,
313
crystalline regions of cellulose cannot tell the whole story.42, 48 Therefore, in addition
314
to enzyme factors (including traffic jams),49 the subsequently increased lignin content
315
might play more important roles. This suggests that biomass porosity does not play
316
the dominant role.17, 19, 46, 50
47
15
Corroborating this view, we observed poor correlations
317
The APS is also considered as a positive factor for enzyme invasion.32 However,
318
the glucan and carbohydrate conversions were significantly negatively correlated with
319
APS (r = −0.634 and −0.561, P < 0.05, respectively; Figures 1 and 2i). The APS 16 ACS Paragon Plus Environment
Page 17 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
320
decreased with the increase of the severity in pretreatments, while adverse trends were
321
observed in both SSA and PV. This indicated more nano-pores (>10nm in present
322
study, consistent with the reported pore size to enable essential accessibility of
323
cellulase protein molecule40, 51) was formed during the ‘peeling away’ of the cell wall
324
components in pretreatments.3, 52 At the most fundamental level, these nano-pores
325
provides more reaction volumes for the synergistical enzymes, thus facilitating
326
enzymatic hydrolysis.53 Consequently, in the next stage of PCA, the APS was
327
preserved in the sequential study, whereas the SSA and PV were excluded.
328
Lignin structural features
329
S/G ratio. Lignin subunits were identified by 2D-HSQC spectra. The S/G ratio varied
330
from 0.59 to 6.39 (Figure 2j; for relative contents of syringyl (S), guaiacy (G) and
331
p-hydroxyphenyl (H) units, see Table S5). The S/G ratio was positively correlated
332
with the glucan and carbohydrate conversions (r = 0.651 and 0.712, P < 0.05,
333
respectively; Figures 1 and 2j). This result revealed that a higher S/G ratio facilitated
334
sugar release, consistent with the trend observed in a large-scale poplar study.9 Due to
335
the open C-5 position, G units not only feature higher affinity to enzymes than S units
336
and thus aggravate the non-productive binding of cellulases.11 They could further lead
337
to more stable 5-5 and β-5 linkages and form a more cross-linked and recalcitrant
338
structure of lignin.13, 54 By contrast, the methoxylated and blocked C-5 position in S
339
units preferentially forms β-β linkages, which contributes to the linearity of lignin.9, 55
17 ACS Paragon Plus Environment
ACS Sustainable Chemistry & Engineering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 36
340
Considering G-rich lignin with lower reactivity is also related to less labile β-O-49,
341
the relative content of lignin linkages was investigated. The S/G ratio was positively
342
correlated with the relative content of β-O-4 (r =0.458, P < 0.05, Table S5), but
343
negatively correlated with the relative content of β-5 (r = −0.413, P < 0.05, Table S5).
344
These results confirmed the negative role of G-rich lignin.
345
The negative influence of lignin content reportedly decreases at higher S/G ratios
346
(≥2.0).9 However, the interaction effect between the lignin content and S/G ratio
347
remains nonquantitative. Therefore, the relative contribution of these two factors in
348
biomass digestibility was additionally weighed by standardised MLR, wherein the
349
coefficients of the lignin content were much higher than those of the S/G ratio (glucan
350
conversion regression: −0.679 vs. 0.263, carbohydrate conversion regression: 0.681 vs.
351
0.323, Table S9). These results implied that enzymatic digestibility might be affected
352
by both the lignin content and composition in untreated biomass, but only dominated
353
by lignin content in pretreated biomass.56 Regardless of the S/G ratio effect, lowering
354
lignin content is preponderant in improving the digestibility.11
355
Hydroxyl Features. The non-productive binding of cellulase is often relevant to the
356
hydroxyl groups of lignin8, which could be quantified by
357
(Figure S3, Table S6). Among the phenolic hydroxyls (PhOH), we detected C5
358
substituted condensed phenolic OH, guaiacyl phenolic OH, catechol type OH and
359
p-Hydroxyl-phenyl-OH.69 Aliphatic (AOH) and carboxylic acid hydroxyls (COH)
360
were summed as a single factor ACOH (aliphatic + carboxylic acid hydroxyls). 4, 57-58 18 ACS Paragon Plus Environment
31
P NMR spectroscopy
Page 19 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
361
The PhOH and ACOH were then summed to give the TOH. Interestingly, the PhOH
362
was significantly negatively correlated with carbohydrate conversion (r = −0.509, P
0.05, Figures 1
364
and 2 m). PhOH in residual lignin could interact with the amide groups in enzymes
365
through hydrogen bond and strengthen the hydrophobic interactions between aromatic
366
amino acids in enzymes and the aromatic rings of PhOH. Thus, phenolic compounds
367
became disruptive in enzymatic functions.28 This result confirmed the negative role of
368
PhOH in cellulase inactivation. However, given the abovementioned structural model
369
of lignocellulose (in which lignin is directly linked to hemicellulose), PhOH should
370
less strikingly affect the glucan conversion59 and was thus excluded from the
371
glucan-conversion-relevant factors in the following PCA stage.
372
Unexpectedly, glucan/carbohydrate conversion was negatively correlated with
373
ACOH (r = −0.500/−0.580, P < 0.05, Figure 2k) and TOH (r =−0.509/−0.632, P
15). Therefore, it was highly unlikely that the results were observed under the
437
null hypothesis.22 All variance inflation factors (VIF) were below 10.0, indicating
438
insignificant correlations among the factors. These results guaranteed reliability of the
439
regression. In the next internal validation by LOO and BOOT, the cross-validation
440
2 2 coefficients were relatively high (QLOO-Glucan = 0.672 and QLOO-Carbohydrate = 0.654). The
441
2 BOOT step also yielded sufficiently high coefficients ( QBOOT-Glucan = 0.692 and
22 ACS Paragon Plus Environment
Page 22 of 36
Page 23 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
442
2 QBOOT-Carbohydrate = 0.698). Their RMSELOO-Glucan, RMSELOO-Carbohydrate, RMSEBOOT-Glucan ,
443
RMSEBOOT-Carbohydrate were 0.109, 0.195, 0.096, 0.171, respectively. The validation
444
coefficients (> 0.5) and RMSEs (< 0.3) confirmed the robustness of the regressions.
445
The experimental validation involved in both mild and severe pretreatments was
446
further conducted on four additionally-tested samples (Table S11). The errors between
447
the predicted and measured sugar conversions were less than 10% and their
448
2 cross-validation coefficient QExt also exceeded 0.5 (Table S11), confirming the
449
validity of the statistical analysis.
450
Eqs. I1 and I2 with sufficient predictability (R2 = 0.784 and 0.832 for glucan and
451
carbohydrate conversion, respectively) can be highlighted for promoting enzymatic
452
digestibility by delignification. However, the major inhibition mechanism of lignin
453
might differ in glucan and carbohydrate conversions. The glucan conversion was
454
mainly affected by O/C ratio and lignin content, indicating that its underlying
455
mechanism involves surface-lignin coverage. Specially, during acid pretreatment,
456
lignin migrates to the biomass surface and aggregates into droplets with a low O/C
457
ratio,61 while the largely dissolved surface lignin in base pretreatment increases the
458
O/C ratio and improves the biomass digestibility.27 Given that a significantly
459
heterogeneous and hydrophobic surface with high lignin content impedes the glucan
460
conversion, we inferred that lignin negatively impacts the glucan conversion mainly
461
through the non-productive binding of cellulase.
462
conversion was mainly influenced by cellulose content besides lignin content. Given
62-63
By contrast, the carbohydrate
23 ACS Paragon Plus Environment
ACS Sustainable Chemistry & Engineering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
463
that cellulose is minimally depolymerised in acid/base pretreatment, the cellulose
464
content is indirectly varied by both lignin and hemicellulose removal. The removal of
465
the recalcitrant components unblocked the initial availability of the biomass substrate.
466
This also implied that, for carbohydrate conversion, the steric hindrance might be the
467
major negative factor. To avoid these two inhibitors, removing lignin to a low content
468
should be the most direct and efficient approach to optimise the amenability of the
469
biomass to enzymatic hydrolysis.
470
CONCLUSION
471
In this work, the substrate-related factors of corn stover pretreated by acid/base
472
and enzymes were comprehensively characterised. A statistically-designed sequence
473
was employed to account for the interaction effects of multivariate and multiscale
474
factors involved in biomass recalcitrance. Among thirteen factors, lignin content,
475
CrI/cellulose, APS, TOH, ACOH, PhOH were significantly negative factors, whereas
476
cellulose content, O/C ratio and S/G ratio played the positive roles. The hemicellulose
477
content, CrI, SSA, and PV proved irrelevant to biomass digestibility. Three
478
representative factors (i.e., lignin content, O/C ratio and cellulose content) were then
479
selected by PCA. Finally, a standardised model was constructed by stepwise MLR,
480
which mechanically interpreted the impacts of these three factors. The equations were
481
determined as Glucan Conversion = 0.225 − 1.38 × Lignin Content + 0.911 × O/C
482
ratio and Carbohydrate Conversion = 0.261 − 3.142 × Lignin Content + 1.488 ×
483
Cellulose Content. Lignin content was systematically and quantitatively proved most 24 ACS Paragon Plus Environment
Page 24 of 36
Page 25 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Sustainable Chemistry & Engineering
484
dominant. The influence of the lignin structural features (S/G ratios and hydroxyl
485
groups) were important but less pronounced when high lignin content remained. More
486
importantly, the sequential analysis revealed different inhibition mechanisms for
487
glucan and carbohydrate conversions, i.e., non-productive binding of cellulase to
488
lignin and steric hindrance of lignin, respectively. Further sequential analysis on
489
larger-scale samples would provide deeper insights into biomass digestibility, thus
490
enabling the optimisation of the bioconversion process.
491
ASSOCIATED CONTENT
492
Supporting Information
493
The Supporting Information is available free of charge on the ACS Publications
494
website. 31
495
Compositional analysis and lignin extraction methods; Q-Q plots, XRD and
496
NMR spectra (Figures S1-3); Pretreatment details, compositional analysis and
497
mass recovery of samples, XPS, CrI and CrI/cellulose, 2D NMR, 31P NMR, PCA
498
results, standardized MLR results and experimental test (Tables S1-11, PDF).
499
AUTHOR INFORMATIN
500
Corresponding Author
501
Yan Shi
502
*E-mail:
[email protected]; Fax: +86-0731-88710171; Tel: +86-0731-88830875
503
Notes
504
The authors declare no competing financial interest.
505
ACKNOWLEDGMENTS
25 ACS Paragon Plus Environment
P
ACS Sustainable Chemistry & Engineering 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 26 of 36
506
This work was supported by key project of National Natural Science Foundation
507
of China (51634010), National Natural Science Foundation of China (31400115,
508
51474102), China Postdoctoral Science Foundation (2017M612594).
509
TABLES AND FIGURES
510
Table 1. Development of regression equations with the factors selected by PCA using
511
stepwise MLR method.
512 #
Equation
R2
F
p
VIF
I1
Glucan Conversion =0.225−1.38×LigninContent + 0.911×O/C ratio
0.784
30.8