DOI: https://doi.org/10.21203/rs.3.rs-36005/v1
Background: Highly efficient enzymatic saccharification of pretreated lignocellulose is a primary key step in achieving lignocellulosic biorefinery. Cellobiohydrolase I (Cel7A) secreted by Trichoderma reesei is an industrially used cellulase possessing carbohydrate binding module 1 (TrCBM1) as the C-terminal domain. Non-productive binding of TrCBM1 to lignin significantly decreases enzymatic saccharification efficiency and enhance cost of biomass conversion due to required additional enzymes. Understanding of the interaction mechanism between lignin and TrCBM1 is essentially required to realize cost-effective biofuels production, but the binding sites in lignin have not been clearly elucidated.
Results: Three types of 13C-labeled b-O-4 lignin oligomer models were synthesized and characterized. The 2D 1H-13C HSQC spectra of the 13C-labeled lignin models exhibited that 13C-labels were correctly incorporated in the (1) aromatic rings and b positions, (2) a positions, and (3) methoxy groups, respectively. The TrCBM1 binding sites in lignin were analyzed by observing NMR chemical shift perturbations (CSPs) using the synthetic 13C-labeled b-O-4 lignin oligomer models. Obvious CSPs were observed in signals from the aromatic regions in oligomers bound to TrCBM1, whereas perturbations in the signals from aliphatic regions and methoxy groups were insignificant. This indicated that hydrophobic interactions and p–p stacking were dominating factors in non-productive binding. The synthetic lignin models have two configurations whose terminal units were differently aligned and donated C(I) and C(II). The C(I) ring showed remarkable perturbation compared with C(II), which indicated that binding of TrCBM1 is evidently affected by configuration of lignin models. Long-chain lignins (DP 4.16–4.70) clearly bound to TrCBM1. Interactions with short-chain lignins (DP 2.64–3.12) were insignificant, indicating that a DP greater than 4 was necessary for TrCBM1 binding.
Conclusion: The CSP analysis using 13C-labeled b-O-4 lignin oligomer models enabled us to identify TrCBM1 binding sites in lignin at the atomic level. This specific interaction analysis will lead to new molecular design of cellulase having controlled affinity to cellulose and lignin for cost-effective biorefinery process.
The enzymatic saccharification of lignocellulose is a key process for the sustainable manufacture of green chemicals and biofuels [1]. Trichoderma reesei is a filamentous fungus that is widely used for the production of commercially available cellulolytic enzyme cocktails. Cellobiohydrolase I (Cel7A) accounts for up to 60% of the cellulases secreted by T. reesei, and carbohydrate-binding module 1 (TrCBM1) (Fig. 1) is connected to the C-terminus of the Cel7A catalytic domain by a highly glycosylated linker [2]. TrCBM1 enhances enzymatic activity by binding to cellulose [3]. TrCBM1 also has a strong affinity for lignin, so lignin significantly inhibits the enzymatic saccharification of pretreated lignocellulose [4]. Non-productive binding of TrCBM1 to lignin should thus be suppressed to achieve efficient enzymatic saccharification. However, the mechanism of interaction is not entirely understood at the molecular level.
Binding of TrCBM1 to lignin is affected by various saccharification conditions, such as temperature [6], pH [7], and the substrate concentration [8]. The chemical properties of lignin have a critical impact on its cellulase binding affinity. Many pretreatments increase the number of phenolic OH groups in lignin and its degree of condensation, which enhance binding between cellulase and lignin [9, 10]. Electrostatic repulsion due to large numbers of aliphatic OH groups and negatively charged functionalities, such as carboxyl and sulfone groups, interferes with cellulase binding [9, 11]. Softwood lignin is mainly synthesized via the radical coupling of guaiacyl monomers, whereas hardwood lignin is synthesized from both guaiacyl (G) and syringyl (S) monomers. The ratio of syringyl to guaiacyl units, denoted by S/G, is one of the most important factors governing the physicochemical properties of lignified plant cell walls. The S/G ratio affects non-productive binding, but it is not the decisive factor for adsorption. Guo et al. found that a low S/G ratio corresponded to a high adsorption capacity [12]. Yang et al. reported that organosolv lignin isolated from softwood lodgepole pine had a higher adsorption affinity for commercial cellulase than lignin from hardwood poplar [13]. However, genetic engineering studies have suggested the opposite, and comparisons of enzymatic hydrolysis in alfalfa cultivars and Eucalyptus mutants with high and low S/G ratios have yielded inconsistent results [14, 15]. Although the relationship between the chemical structure of lignin and its binding affinity for cellulolytic enzymes and CBMs has been researched extensively, evidence for the identification of binding atoms in whole lignin molecules is entirely lacking. This is because no effective method for identifying the interactive sites in lignin, a heterogenous polymer, is available.
NMR chemical shift perturbation (CSP) analysis is a powerful method for the identification of substrate binding sites in proteins [16]. We previously performed CSP analysis to identify the binding amino acid residues in TrCBM1 with softwood and hardwood milled wood lignin (MWL) using 15N-labeled TrCBM1 [17]. The results suggested that the aromatic rings in lignin participated in interactions with amino acid residues in TrCBM1, because the flat plane surface including Y5, Y31, and Y32 in TrCBM1 was a main interaction site. However, the binding sites in lignin have not been characterized at the molecular level. In this study, we synthesized 13C-labeled β-O-4 lignin oligomer model compounds with different labelling positions and performed CSP analysis to identify the TrCBM1 binding sites in the lignin models. Pure TrCBM1 without a catalytic domain and linker was expressed and purified and then added to NMR sample tubes containing the 13C-labeled lignin models to monitor perturbations in signals from the binding sites in the model compounds. NMR analysis of the 13C-labeled lignin model compounds provided the first direct evidence for the identification of binding atoms in the linear lignin chains. The results were in good agreement with our previous binding site analysis of the protein counterpart, TrCBM1 with MWL [17]. Herein, we provide direct evidence to reveal the interaction sites in the β-O-4 lignin substructure that bind to TrCBM1.
Expression and purification of TrCBM1
Escherichia coli BL21 (DE3) was used to express a histidine (His) tag-TrCBM1-green fluorescent protein (GFP) fusion protein. TrCBM1 was obtained after the addition of enterokinase and thrombin to cleave the His tag and GFP, respectively. We previously performed sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) and 2D 1H-15N HSQC NMR to analyze the molecular mass and conformation of 15N-labeled TrCBM1 [17]. In this work, we performed SDS-PAGE and MALDI-TOF-MS to characterize unlabeled TrCBM1. Pure TrCBM1 was a single protein with a molecular mass of 5195.8 Da. The MALDI-TOF-MS spectrum of TrCBM1 and a full-length SDS-PAGE gel are shown in Figure S1 of the Additional file.
In recent work, we synthesized the unlabeled β-O-4 lignin oligomer model compound 4 (Fig. 2) and fully characterized it by performing 1D 1H NMR, 1D 13C NMR, 2D 1H-13C HSQC, 2D 1H-13C HMBC, 2D 1H-13C long-range heteronuclear single quantum multiple bond correlation (LR-HSQMBC), 2D 1H-1H TOCSY, 2D 1H-1H EXSY and 2D 1H-1H ROESY [18]. We applied the protocol to synthesize three 13C-labeled β-O-4 lignin oligomer model compounds. 13C-labels were placed in the aromatic rings and β positions of compound 4(Arβ), the α positions of compound 4(α), and methoxy groups in compound 4(m) using 13C-labeled vanillins and t-butyl-2-bromoacetate. The 13C-vanillins 1(Ar), 1(α), and 1(m) were individually refluxed in acetone with unlabeled or 13C-labeled t-butyl-2-bromoacetate to afford 13C-t-butoxycarbonylmethyl vanillins 2(Arβ), 2(α), and 2(m), respectively. The 13C-t-butoxycarbonylmethyl vanillins were dissolved in anhydrous tetrahydrofuran (THF) and polymerized using lithium diisopropylamide (LDA). The ester groups in polymerized 13C-oligomers 3(Arβ), 3(α), and 3(m) were reduced with NaBH4 to obtain 13C-labeled β-O-4 lignin oligomer model compounds 4(Arβ), 4(α), and 4(m), respectively.
The 2D 1H-13C HSQC spectra of the 13C-labeled lignin models contained signals that corresponded to the 1H-13C correlations (Fig. 3), which was evidence that 13C was incorporated at the designated positions. No HSQC signals attributable to side products were observed in the NMR spectra, which indicated that high-purity lignin models were prepared. HSQC signals at δC/δH 85.9/4.52–4.64, 72.6/4.06–4.08, and 113.2–123.6/6.81–7.06 ppm in the HSQC spectra of lignin model 4(Arβ) with 13C-labeling in the β positions and aromatic rings (Fig. 3a) were assigned to the β positions, Terminal C1, and aromatic regions, respectively. The HSQC spectra of lignin model 4(α) (Fig. 3b) contained three signals at δC/δH 74.7–75.0/4.92–4.98, 74.8–75.0/4.72–4.76, and 66.2/4.43–4.49 ppm, which were assigned to α positions in the A, B, and C rings, respectively. A single signal was observed at 58.1–58.3/3.57–3.85 ppm in the HSQC spectra of lignin model 4(m) with 13C-labeled methoxy groups (Fig. 3c).
The long- and short-chain lignin models were separated via silica gel chromatography and subjected to NMR and binding analysis. Size exclusion chromatography (SEC) revealed that the 13C-labeled lignin models had narrow molecular weight distributions (Figure S2, Additional file). The degree of polymerization (DP) was calculated for each model compound from its weight-average molecular weight (Mw). The DPs of the long-chain lignin models ranged from 4.16 to 4.70, whereas those of the short-chain models ranged from 2.64 to 3.12 (Table 1). Although differences between the Mws of the long- and short-chain models were small, their NMR spectra (Fig. 3d) and SEC profiles (Figure S2, Additional file) were clearly distinct. Therefore, the long- and short-chain lignin models could be used to evaluate the effects of Mw on molecular interactions with TrCBM1.
| Lignin models | Chain length | Mn | Mw | Mw/Mn | DPa |
|---|---|---|---|---|---|
| Compound 4(Arβ) (13C-labeled aromatic rings and β positions) | Long | 843 | 923 | 1.095 | 4.70 |
| Short | 480 | 520 | 1.083 | 2.64 | |
| Compound 4(α) (13C-labeled α positions) | Long | 750 | 918 | 1.224 | 4.67 |
| Short | 531 | 614 | 1.156 | 3.12 | |
| Compound 4(m) (13C-labeled methoxy groups) | Long | 629 | 818 | 1.301 | 4.16 |
| Short | 482 | 587 | 1.218 | 2.98 | |
| Compound 4 (unlabeled) | Long | 882 | 964 | 1.092 | 4.91 |
| Short | 552 | 560 | 1.014 | 2.85 | |
| a DP: Degree of polymerization calculated from the Mw and theoretical molecular mass of the lignin model. | |||||
Analysis of interactions between lignin models and TrCBM1
Interactions between the lignin models and TrCBM1 were evaluated via CSP analysis using 2D 1H-13C HSQC NMR. To analyze the binding positions in lignin model compounds 4(Arβ), 4(α), and 4(m) with a high degree of sensitivity, the compounds were separated into high and low molecular mass fractions. Assignment of the 13C-lignin model HSQC signals was accomplished for all of the carbon atoms and non-exchangeable protons. Almost all of the signals from the A, B, and C rings were assigned separately, because the tops of the peaks were distinguishable [18]. We previously found that the lignin models have two diastereomers. Although these A and B rings have similar alignments, the C ring are differently aligned and designated as C(I) and C(II) [18]. Full assignment of the peaks in the lignin model spectra enabled us to analyze the interactions of carbon atoms and protons at the atomic level.
The 2D 1H-13C HSQC spectra of the lignin models (50 µM) in the presence and absence of TrCBM1 were superimposed (Fig. 4). Several HSQC signals acquired in the presence of 200 µM TrCBM1 showed obvious perturbation, particularly those of the aromatic regions in the long-chain lignin spectra. Most of these signals exhibited larger perturbations when the concentration of TrCBM1 was increased to 350 µM (Fig. 4b). In contrast, there was no significant perturbation of NMR signals from the aliphatic regions and methoxy groups in the presence of TrCBM1 (Fig. 4c). The results clearly showed that the aromatic regions in the long-chain lignin models were the primary sites of interaction with TrCBM1. However, no distinct perturbations were observed in signals from the methoxy groups or any of the aromatic and aliphatic regions when the short-chain lignin models were used for CSP analysis (Figures S3, Additional file). This indicated that the length of the lignin chains was an important factor in TrCBM1 binding.
The change in the chemical shift (Δδ) of each signal was calculated by subtracting the chemical shift value in the spectrum of the lignin model from that of the model recorded in the presence of TrCBM1. The Δδ values obtained in the 1H and 13C-axis are plotted in Fig. 5. Δδ values were large only for the signals from the aromatic regions in the long-chain lignin models, whereas the Δδ values of signals from the aliphatic regions and methoxy groups in both the long- and short-chain lignin models were small. The Δδ values of the A 2, C 2(I), B 5, C 5(I), C 5(II), A 6, C 6(I), and C 6(II) signals from the aromatic regions of the long-chain lignin models were all larger than 0.006 ppm on the 1H-axis, whereas the Δδ values of the A 2 and C 5(I) signals on the 13C-axis were greater than 0.05 ppm (Fig. 5a). For the short-chain lignin models, only the Δδ values of the B 5 signals on the 1H-axis exceeded 0.006 ppm (Fig. 5b). Line broadening was observed in the B 5 signals of the long-chain lignin models in the presence of TrCBM1 (350 µM). The observed line broadening could be attributed to both the on and off rates of complex formation and the multiple binding states of lignin due to non-specific binding [17].
The unlabeled long- and short-chain lignin models of compound 4 were used to evaluate TrCBM1 binding affinity according to the Langmuir adsorption model. The Mws and molecular weight distributions of the unlabeled lignin models (Table 1) were nearly identical to these of the 13C-labeled lignin models (Figure S2, Additional file). The TrCBM1–His tag fusion protein was used instead of TrCBM1 without a His tag to conduct the adsorption experiments, because the soluble lignin models could not be separated from the TrCBM1 protein via centrifugation. After incubating the sample solutions containing the lignin models and TrCBM1–His tag at 50 °C for 1 h, cOmplete His Tag Purification Resin (Roche, Basel, Switzerland) was added to the solutions. The unbound lignin models in the supernatant could then be separated from the bound lignin models, which were adsorbed to the precipitated TrCBM1–His tag bound to the His tag resin. We confirmed that all of the TrCBM1–His tag would adsorb to the His tag resin by performing control experiments. We did not observe unwanted binding between the His tag resin and lignin models (Figures S6 and S7, Additional file). The concentrations of the unbound lignin models in the supernatants were determined to calculate the adsorption parameters summarized in Table 2. Although the long- and short-chain lignin models had similar Γmax values, the KL of the long-chain lignin model was eight times higher than the KL of the short-chain lignin model. This result was consistent with the results of the NMR interaction analysis. The percentages of the lignin models that bound to TrCBM1 were calculated using the KL values. In the CSP analysis, 84.6% and 91.4% of the long-chain lignin models were bound to TrCBM1 at TrCBM1 concentrations of 200 and 350 µM, respectively. At TrCBM1 concentrations of 200 and 350 µM, 33.0% and 46.6% of the short-chain lignin models, respectively, were bound to TrCBM1. The chain length was thus an essential factor in binding between TrCBM1 and the lignin chains, which contained β-O-4 linkages exclusively. We also found that strong adsorption required a DP above 4.
| Lignin model | Langmuir affinity constant, KL (mL/mg) | Amount of adsorption at saturation, Γmax (µg/mg) |
|---|---|---|
| Long-chain | 36.26 | 8.77 |
| Short-chain | 4.47 | 9.15 |
Understanding the mechanism of interaction between cellulase and lignin is essential for the efficient enzymatic saccharification of lignocelluloses. T. reesei is the most important industrial cellulase-producing filamentous fungus. Cel7A is the most abundant secreted cellulolytic enzyme and contains a catalytic domain and TrCBM1. Non-productive binding between TrCBM1 and lignin decreases the efficiency of enzymatic saccharification of pretreated plant biomass, but the interactive sites in lignin have not been previously elucidated.
The flat plane surface of TrCBM1 binds to both lignin and cellulose, and contains the three tyrosine residues (Y5, Y31, and Y32) and neighboring amino acid residues in the underside of TrCBM1 (Fig. 1) [7, 19]. T17, V18, and T24 residues are present in the cleft, which is located on the opposite side of the flat plane surface. This site also interacts with lignin and cellulosic substrates [17]. Although the TrCBM1 binding sites have been extensively investigated via point mutation analysis and NMR, evidence to identify the lignin atoms that participate in TrCBM1 binding is lacking. Therefore, we conducted CSP analysis to investigate interactions between 13C-labeled β-O-4 lignin oligomer compounds and TrCBM1. Interactive sites in the long-chain lignin model were mapped in Fig. 6, which are based on Δδ values observed in CSP analysis. Interactions between the aromatic rings and TrCBM1 were obvious, whereas the aliphatic regions and methoxy groups were not major binding sites. This suggests that TrCBM1 adsorption on lignin occurs via hydrophobic interactions and π–π stacking of the aromatic rings in lignin and the three tyrosine residues on TrCBM1. Rahikainen et al. discussed the importance of TrCBM1 hydrophobicity. They reported that a Y32A TrCBM1 mutant had a lower association constant than wild-type TrCBM1, whereas a Y32W mutant increased the lignin and cellulose binding affinities of TrCBM1 [7]. The hydrophobicity of lignin results in a significant amount of non-productive binding with cellulase [13]. Therefore, hydrophobic interaction can reasonably be interpreted as a dominant driving force behind non-productive TrCBM1 binding by lignin.
Interestingly, there is a variety of Δδ in the aromatic rings of lignin model shown in Fig. 6. The interaction patterns of terminal units on the A and C(I) rings exhibited similar interaction patterns although the A 5 position was not evaluated due to its extremely weak HSQC signal intensity. Internal units, B ring, showed remarkable interactions with TrCBM1 only at the B 5 position, suggesting that the interaction patterns of the terminal and internal units of the lignin model differed. We previously reported two configuration of lignin models whose terminal units were differently aligned and designated them C(I) and C(II) [18]. We demonstrated that the TrCBM1 binding behavior depended significantly on the configurations of the lignin model. The molecular alignment of C(I) apparently resulted in preferential binding to TrCBM1. We are thus the first to reveal that interactive sites in the lignin chains are significantly influenced by their molecular configurations.
Hydrogen bonding and electrostatic interactions are not negligible in non-productive TrCBM1 binding. Phenolic OH groups in lignin promote hydrogen bonding to cellulolytic enzymes [10, 20]. Hydrophilic amino acid residues in TrCBM1, including Q7 and T17 in addition to main chain of H4 and I11 were previously shown to participate in binding with lignin via hydrogen bonds and electrostatic interactions [17]. In our CSP analysis, aromatic rings in the lignin models were the primary sites of interaction with TrCBM1. This was attributed mainly to hydrophobic interactions. However, the results could also be interpreted in terms of hydrogen–π interactions between the aromatic rings in lignin and OH groups on amino acid residues [21, 22]. OH groups in the three tyrosine residues on the flat plane surface of TrCBM1 and the OH group in the T17 residue in the cleft were thought to donate hydrogens for hydrogen–π interactions. Although CSP analysis indicated that interactions at the aliphatic positions were insignificant compared with those at the aromatic positions, signals from the α positions had slightly higher Δδ values than signals from the other aliphatic positions in the lignin models (Fig. 5). This result could be attributed to hydrogen bonding between OH groups near the α positions of the lignin model and TrCBM1.
It is uncertain whether there is a correlation between the molecular weight of lignin and non-productive binding [13, 23]. This is because a decrease in the molecular weight of lignin changes its hydrophobicity and OH content, which are the predominant factors affecting non-productive binding. Changes in these structural features of lignin have disturbed us to estimate the effect of the molecular weight of lignin on cellulase binding. Our CSP results clearly indicated that the long-chain lignin model had a stronger affinity for TrCBM1 than the short-chain lignin model (Table 2). Mattinen et al. found that cellohexaose adsorbed to TrCBM1, whereas short cellooligosaccharides, such as cellobiose and cellotriose, did not bind to TrCBM1 [24]. Based on these results and our observations, stacking on the flat plane surface of TrCBM1 and interactions with the cleft required lignin and cellooligosaccharides with long chains. The low affinity of the short-chain lignin model could also be attributed to its conformation. The short-chain lignin model had a folded conformation in 90% water [18], and the folded short lignin chain could not cover the full length of the three tyrosine residues on the flat plane surface of TrCBM1. Therefore, we concluded that β-O-4 lignin chains with DPs above 4 were essential for strong adsorption of TrCBM1.
Highly efficient enzymatic saccharification is hindered by the non-productive binding of lignin to cellulolytic enzymes. Until now, the binding mechanism was not fully understood. We synthesized 13C-labeled β-O-4 lignin oligomer model compounds to identify the TrCBM1 binding sites in lignin via NMR CSP analysis. Signals from the aromatic regions in the lignin models exhibited obvious perturbations, whereas signals from the aliphatic regions and methoxy groups were not significantly perturbed. These results suggested that hydrophobic interactions and π–π stacking were the principal forces driving interactions between aromatic rings in the lignin models and tyrosine residues in TrCBM1. TrCBM1 bound differently with terminal and internal units in the lignin models. In addition, the binding patterns associated with the C(I) and C(II) terminal alignments differed. This indicated that binding of the lignin models to TrCBM1 was strongly affected by the molecular configuration. Perturbation of signals from the long-chain lignin models (DP 4.16–4.70) was obvious due to their strong binding affinity relative to that of the short-chain lignin models (DP 2.64–3.12). This indicated that a chain length greater than DP4 was necessary for strong interactions between lignin and TrCBM1. This is the first study to characterize the interactive sites in a lignin model compound at the atomic level using purified TrCBM1. A detailed understanding of non-productive binding will lead to the establishment of a fundamental theory for the structural alteration of lignin and enzymes that are not susceptible to unfavorable binding.
E. coli BL21 (DE3) was obtained from Merck (Darmstadt, Germany). Vanillins with 13C-labeled carbonyl group and 13C-labeled aromatic ring were obtained from Cambridge Isotope Laboratories (Tewksbury, MA, USA). Vanillin with 13C-labeled methoxy group and [2-13C]-tert-butyl bromoacetate were purchased from Taiyo Nippon Sanso (Tokyo, Japan). Other laboratory chemicals were obtained from Wako Pure Chemical Industries (Osaka, Japan), nacalai tesque (Kyoto, Japan), and Tokyo Chemical Industry (Tokyo, Japan).
Preparation of TrCBM1
The CBM1 of T. reesei Cel7A (accession number: CAH10320) was expressed and purified as described previously [17]. E. coli BL21 (DE3) was first subjected to heat shock transformation with plasmids for the expression of His-tagged TrCBM1 fused with GFP, hereafter referred to as His tag–TrCBM1–GFP fusion protein. LB medium was inoculated with the transformant, followed by incubation at 37 °C and 200 rpm until the OD600 reached 1.2. Protein expression was induced using 1 mM isopropyl β-1-thiogalactopyranoside, and the bacteria were incubated at 37 °C and 200 rpm for 5 h. Following centrifugation and sonication, His tag–TrCBM1–GFP was isolated from the supernatant and purified via Ni affinity chromatography followed by anion exchange chromatography. The His tag and GFP region were removed using enterokinase (New England Bio Labs, MA, USA) and thrombin (GE Healthcare, IL, USA), respectively. The obtained TrCBM1 was concentrated in 50 mM acetic acid-d4 buffer prepared with D2O using a Vivaspin turbo ultrafiltration device (Sartorius, Göttingen, Germany). We skipped the His tag cleavage process to prepare His tag–TrCBM1 for adsorption analysis. The His tag–TrCBM1 was purified via Ni affinity and anion exchange chromatography.
Synthesis of β- O -4 lignin oligomer model compounds
β-O-4 lignin oligomer models were synthesized as described in our previous work using a modified protocol originally developed by Katahira et al [18, 25]. Vanillin (1) was dissolved in acetone, and the solution was refluxed for 1.5 h in the presence of KI, K2CO3, and t-butyl-2-bromoacetate to obtain t-butoxycarbonylmethyl vanillin (2). The monomer (2) was dissolved in anhydrous THF, and 1.5 M LDA was added dropwise to the solution under a nitrogen atmosphere over 30 min while maintaining the temperature at − 30 °C. The solution was stirred for 30 min at − 30 °C and for an additional 1.5 h at 0 °C–5 °C to obtain the polymerized oligomer (3). Oligomer 3 was dissolved in t-butanol containing NaBH4. Methanol was added to the mixture stepwise over 1 h with constant stirring at 65 °C–70 °C, and the mixture was stirred for 6 h to obtain β-O-4 lignin oligomer model compound 4. The synthesized lignin models were separated into long- and short-chain fractions via silica gel chromatography. Sequential elution was performed using ethyl acetate, 5:1 (v/v) ethyl acetate/methanol, and 2:1 (v/v) ethyl acetate/methanol. Three different 13C-labeled lignin model compounds were synthesized for CSP analysis. Vanillin with 13C-labeled aromatic rings (1(Ar)) and [2-13C]-tert-butyl bromoacetate were used as starting materials to synthesize a lignin model with 13C-labeled aromatic rings and β positions (4(Arβ)). Similarly, vanillins with 13C-labeled carbonyl and methoxy group were used to synthesize a lignin model with 13C-labeled α positions (4(α)) and a model with 13C-labeled methoxy groups (4(m)). The lignin models were characterized using 2D 1H-13C HSQC NMR and SEC. SEC was performed using three TSKgel SuperMultipore HZ-M columns (Tosho, Tokyo, Japan) on a Shimadzu instrument equipped with an LC-20AD pump and an SPD-M20A diode array detector (Shimadzu, Kyoto, Japan). THF was used for elution at a flow rate of 0.35 mL/min at 40 °C.
All NMR spectra were recorded at 298 K on a Bruker Avance III HD 600 spectrometer equipped with a cryogenic probe and a Z-gradient (Bruker BioSpin, MA, USA). The instrument was controlled using Bruker Topspin NMR software version 3.5. For CSP analysis, 50 µM of the 13C-labeled lignin models with 13C-labeled aromatic rings and β positions (4(Arβ)), α positions (4(α)), and methoxy groups (4(m)) were individually dissolved in 50 mM acetic acid-d4 buffer, which was prepared with D2O (pD 5.0) and 10% (v/v) DMSO-d6. Each NMR sample had a volume of 250 µL and contained 20 µM 2,2-dimethyl-2-silapentane-5-sulfonic acid as an internal standard. To identify the TrCBM1 binding sites in the 13C-labeled lignin models, changes in chemical shift (ΔδC,H, ppm) were calculated by comparing chemical shift values in the edited 1H-13C HSQC spectra of the 13C-labeled lignin models in the presence and absence of 200 and 350 µM TrCBM1. CSP analysis of both the long- and short-chain lignin models was performed to evaluate the effect of molecular weight on TrCBM1 binding to the lignin models. Assignment of the lignin model signals was based on our previous report [18].
Adsorption experiment of lignin models with TrCBM1
The Langmuir adsorption isotherm model was applied to evaluate the TrCBM1 binding affinities of the lignin models. Each unlabeled long- and short-chain lignin model was dissolved in 50 mM acetic acid buffer (pH 5.0) containing 3,200 µg/mL of TrCBM1 with a His tag (His tag–TrCBM1). The lignin models were prepared in 1.5 mL micro-tubes at concentrations of 10, 20, 30, 50, 80, 100, 150, and 200 µM. The total volume of each solution was 50 µL. The sample solutions were incubated at 50 °C for 60 min with constant shaking at 1,000 rpm using a Comfort thermomixer (Eppendorf, Hamburg, Germany). To each solution, 100 µL of a 50% (v/v) suspension of cOmplete His Tag Purification Resin (Roche, Basel, Switzerland) in the same buffer was added. The mixtures were resuspended by vortexing for 10 s and centrifuged at 500 g for 10 s. The free lignin model concentrations in the supernatants were determined by measuring absorbance at 280 nm using a BioSpec-nano spectrophotometer (Shimadzu, Kyoto, Japan). The quantity of each lignin model bound to TrCBM1 was calculated by subtracting the amount of the free lignin model compound from the initially loaded amount. Experiments were conducted in duplicate, and the average values were used to calculate the Langmuir affinity constants using Eq. (1).
ΓC = Γmax[KLC/(1 + KLC)], (1)
where ΓC is the amount of the bound lignin model, Γmax is the amount of the lignin model bound to His tag–TrCBM1 at saturation, KL is the Langmuir affinity constant, and C is the concentration of the free lignin model in the supernatant. A solution containing only His tag–TrCBM1 without a lignin model was used as a blank.
Chemical shift perturbation; DMSO:Dimethyl sulfoxide; DP:Degree of polymerization; EXSY:Exchange spectroscopy; G:Guaiacyl; GFP:Green fluorescent protein; His tag:Histidine tag; HMBC:Heteronuclear multiple bond correlation; HSQC:Heteronuclear single-quantum correlation; LDA:Lithium diisopropylamide; LR-HSQMBC:Long-range heteronuclear single quantum multiple bond correlation; MALDI-TOF-MS:Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry; Mw:Weight-average molecular weight; MWL:Milled wood lignin; NMR:Nuclear Magnetic Resonance; ROESY:Rotating-frame overhauser effect spectroscopy; S:Syringyl; SDS-PAGE:Sodium dodecyl sulfate-polyacrylamide gel electrophoresis; SEC:Size exclusion chromatography; THF:Tetrahydrofuran; TOCSY:Total correlation spectroscopy.
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Availability of data and materials
All data generated or analyzed during this study are included in this manuscript and its Additional files.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
Y.T., T.N., M.K., T.W. designed the experiments. Y.T. performed the experiments. Y.T., K.K., T.N. analyzed and interpreted the NMR data. Y.T., T.W., T.N., M.K. wrote the paper.
All authors read and approved the final manuscript.
Acknowledgements
We thank Dr. Isamu Katsuyama for supporting chemical synthesis of 13C-labeled lignin oligomer model compounds.
Funding
This work was supported by JSPS KAKENHI Grant Number JP18J20331, a collaboration program of RISH (M2-2) and joint usage/research programs of IAE (ZE30A-36, ZE31A-37).