Skip to main content

Role of surface tryptophan for peroxidase oxidation of nonphenolic lignin



Despite claims as key enzymes in enzymatic delignification, very scarce information on the reaction rates between the ligninolytic versatile peroxidase (VP) and lignin peroxidase (LiP) and the lignin polymer is available, due to methodological difficulties related to lignin heterogeneity and low solubility.


Two water-soluble sulfonated lignins (from Picea abies and Eucalyptus grandis) were chemically characterized and used to estimate single electron-transfer rates to the H2O2-activated Pleurotus eryngii VP (native enzyme and mutated variant) transient states (compounds I and II bearing two- and one-electron deficiencies, respectively). When the rate-limiting reduction of compound II was quantified by stopped-flow rapid spectrophotometry, from fourfold (softwood lignin) to over 100-fold (hardwood lignin) lower electron-transfer efficiencies (k 3app values) were observed for the W164S variant at surface Trp164, compared with the native VP. These lignosulfonates have ~20–30 % phenolic units, which could be responsible for the observed residual activity. Therefore, methylated (and acetylated) samples were used in new stopped-flow experiments, where negligible electron transfer to the W164S compound II was found. This revealed that the residual reduction of W164S compound II by native lignin was due to its phenolic moiety. Since both native lignins have a relatively similar phenolic moiety, the higher W164S activity on the softwood lignin could be due to easier access of its mono-methoxylated units for direct oxidation at the heme channel in the absence of the catalytic tryptophan. Moreover, the lower electron transfer rates from the derivatized lignosulfonates to native VP suggest that peroxidase attack starts at the phenolic lignin moiety. In agreement with the transient-state kinetic data, very low structural modification of lignin, as revealed by size-exclusion chromatography and two-dimensional nuclear magnetic resonance, was obtained during steady-state treatment (up to 24 h) of native lignosulfonates with the W164S variant compared with native VP and, more importantly, this activity disappeared when nonphenolic lignosulfonates were used.


We demonstrate for the first time that the surface tryptophan conserved in most LiPs and VPs (Trp164 of P. eryngii VPL) is strictly required for oxidation of the nonphenolic moiety, which represents the major and more recalcitrant part of the lignin polymer.


Removal of the highly recalcitrant lignin polymer is a key step for the natural recycling of plant biomass in land ecosystems, and a central issue for the industrial use of cellulosic feedstocks in the sustainable production of fuels, chemicals and different materials [13]. White biotechnology must contribute to the development of lignocellulose biorefineries by providing tailor-made microbial and enzymatic biocatalysts enabling “greener” and more efficient biotransformation routes for the complete use of both polysaccharides and lignin as the main biomass constituents [4, 5].

The so-called white-rot basidiomycetes (due to the whitish color of delignified wood) are the main lignin degraders in Nature [6]. The process has been described as an “enzymatic combustion” [7] and would involve peroxidases of the lignin peroxidase (LiP), manganese peroxidase (MnP) and versatile peroxidase (VP) families, together with other oxidoreductases [6, 8]. After some controversy in the past [9], the most recent evidence on the involvement of peroxidases in lignin degradation comes from the availability of massive sequencing tools applied to fungal genomes. The analysis of basidiomycete genomes shows the presence of the above ligninolytic peroxidase genes in the genomes of all typical white-rot (ligninolytic) basidiomycetes sequenced to date, and their absence from all the brown-rot (cellulolytic) basidiomycete genomes [1014].

Among the three peroxidase families LiP, first reported from Phanerochaete chrysosporium [15], and VP, described later from Pleurotus eryngii [16, 17], have attracted the highest interest since they are able to degrade nonphenolic model compounds representing the main substructures in lignin (such as β-O-4′ alkyl-aryl ethers) [1820] by single-electron abstraction forming an aromatic cation radical [21], and subsequent Cα–Cβ bond cleavage [22] (while MnP would act on the minor phenolic units). From the discovery of LiP, the huge number of biochemical and molecular biology studies on these enzymes generally used simple aromatic substrates, such as veratryl (3,4-dimethoxybenzyl) alcohol [2325], and similar studies using the real lignin substrate are extremely rare [26].

A landmark in lignin biodegradation studies was the identification of a solvent-exposed peroxidase residue, Trp171 in P. chrysosporium LiP (isoenzyme H8) [27, 28] and Trp164 in P. eryngii VP (isoenzyme VPL) [29], as the responsible for oxidative degradation of nonphenolic lignin model compounds by long-range electron transfer (LRET) from the protein surface to the heme cofactor of the H2O2-activated enzyme. This single-electron transfer generates a reactive tryptophanyl radical [30, 31], whose exposed nature would enable direct oxidation of the lignin polymer. Recently, the authors have shown that removal of this aromatic residue lowers in different extents the electron transfer from technical lignins (partially phenolic softwood and hardwood water-soluble lignosulfonates) to the peroxide-activated VP transient states (the so-called compounds I and II, CI and CII) [32, 33].

To clarify the role of the surface tryptophan residue in phenolic/nonphenolic lignin degradation, stopped-flow reactions of the above VP and the corresponding tryptophan-less variant are performed in the present study using native (underivatized) and permethylated/acetylated (nonphenolic) softwood and hardwood lignosulfonates as enzyme substrates, together with lignosulfonate steady-state treatments analyzed by size-exclusion chromatography (SEC) and heteronuclear single quantum correlation (HSQC) two-dimensional nuclear magnetic resonance (2D-NMR).


Transient kinetics of VP and its W164S variant: native lignins

Peroxidase catalytic cycle includes two-electron activation of the resting enzyme by H2O2 yielding CI, which is reduced back via CII with one-electron oxidation of two substrate molecules (Additional file 1: Figure S1a). These three enzyme forms present characteristic UV–visible spectra (Additional file 1: Figure S1b, c) that enable to calculate the kinetic constants for CI formation and CI/CII reduction (see “Methods” section).

The transient-state kinetic constants for the reaction of native lignosulfonates with H2O2-activated wild-type recombinant (hereinafter native) VP and its W164S mutated variant were obtained by stopped-flow rapid spectrophotometry, showing CII reduction as the rate-limiting step [34]. In the reactions of native VP CI and CII (Fig. 1a; Additional file 1: Figure S2a, d, continuous lines) relatively similar apparent second-order rate constants (k 2app and k 3app) were obtained for the two lignosulfonates (top of Tables 1, 2) (k 1app for CI formation by H2O2 being 3460 ± 70 s−1 mM−1). The main difference was in the CII reduction dissociation constant (K D3), which was tenfold lower for hardwood than softwood lignosulfonate indicating a higher affinity for the former lignin. Softwood lignosulfonate did not saturate native VP for CI reduction (Additional file 1: Figure S2a, d, red continuous line) and only a k app value can be provided.

Fig. 1
figure 1

Kinetics of CII reduction by native (a), acetylated (b) and permethylated (c) softwood (LSS, red) and hardwood (LSH, blue) lignosulfonates: Native VP (continuous line) vs W164S variant (dashed line). Stopped-flow reactions were carried out at 25 °C in 0.1 M tartrate (pH 3). The lignosulfonate concentrations (here and in Additional file 1: Figure S2) refers to the lignosulfonate basic phenylpropanoid unit. Means and 95 % confidence limits are shown

Table 1 CI reduction constants by acetylated, methylated and native (softwood and hardwood) lignosulfonates: Native VP vs W164S variant
Table 2 CII reduction constants by acetylated, methylated and native (softwood and hardwood) lignosulfonates: Native VP vs W164S variant

In the W164S variant (whose no-saturation kinetic traces are included in Fig. 1a; Additional file 1: Figure S2a, d, dashed lines) substitution of the catalytic tryptophan resulted in impaired oxidation of both lignosulfonates (bottom of Tables 1, 2). The strongest effect was with the hardwood lignosulfonate, where the k 2app and rate-limiting k 3app values experienced ninefold and 125-fold decreases, respectively.

Transient kinetics of VP and its W164S variant: nonphenolic lignins

The residual reduction of W164S CI and CII in the above experiments could be due to the presence of more easily oxidizable phenolic units. Using NMR after sample acetylation, the lignosulfonate phenolic content was estimated as ~20–30 % of lignin units. Methylation was optimized using pyrolysis–gas chromatography/mass spectrometry (Py-GC/MS) to follow the reaction progress (Additional file 1: Figure S3) till complete derivatization (of both phenolic and alcoholic hydroxyls), as shown by NMR after secondary acetylation (Fig. 2).

Fig. 2
figure 2

Lignosulfonate permethylation: 1H-NMR analysis after secondary acetylation confirming the previous complete methylation of softwood lignosulfonate (b) compared with the untreated sample (a). Regions of phenolic and alcoholic acetates are indicated

Then, new transient-state kinetic constants were calculated for the derivatized (nonphenolic) lignosulfonates. Figure 1b, c (and Additional file 1: Figure S2b/e, c/f) show the kinetic traces for the acetylated and methylated lignosulfonates, respectively, whose CI and CII reduction constants are included in Tables 1 and 2, respectively. With these nonphenolic lignins no strong difference between CI and CII reduction rates was observed, in contrast with native lignosulfonate where CII reduction is clearly the rate-limiting step. In most native VP reactions (continuous lines), saturation kinetics was observed (except for CI reduction by methylated softwood lignosulfonate) and only a k 2app value can be provided. The opposite tendency was found for the W164S variant (dashed line) where saturation was more rarely observed.

For native VP, lignin methylation (and in lower extent acetylation) significantly decreased CI reduction (Additional file 1: Figure S2, left) resulting in 20–40-fold lower k 2app values, while CII reduction was much less affected (Fig. 1). However, for the W164S variant, similar decreases in both CI and CII reduction were observed, resulting in 25–45-fold lower k app for the methylated samples.

When the effect of W164S mutation on the nonphenolic lignin constants was considered (bottom of Tables 1, 2), small decreases in CI reduction were observed (similar to those obtained with native lignins). However, for reduction of W164S CII strong k 3app decreases with respect to native VP were observed (up to 350-fold for the methylated hardwood lignosulfonate). More importantly, the previously observed reduction of W164S CII by native softwood lignosulfonate (Fig. 1a, red dashed line) disappeared when the acetylated or methylated samples were evaluated as W164S substrates (Fig. 1b, c red dashed lines, respectively).

Steady-state treatment of native lignin with VP and its W164S variant

In addition to the above stopped-flow reactions, the effect of the enzymatic treatments was also analyzed by SEC and 2D-NMR spectroscopy during steady-state reactions.

Native VP significantly modified the molecular-mass distribution and main peak (Mp) of softwood and especially of hardwood lignosulfonates (green continuous lines in Fig. 3a, b, respectively), with respect to the controls (red and blue lines), revealing a clear polymerization tendency in the latter case (Mp of ~20,000 Da compared with ~5500 Da in the control) and the disappearance of a broad shoulder around 11 mL elution volume (~6800 Da) in the former case. More importantly, the W164S variant only caused a very limited modification in the molecular-mass distribution of the two lignins, in agreement with its low kinetic constants for rate-limiting CII reduction. Such modification included a modest displacement of Mp (to 6500 Da) in hardwood lignosulfonate and a slight decrease of the softwood lignosulfonate shoulder (dashed lines).

Fig. 3
figure 3

SEC profiles of softwood (a) and hardwood (b) lignosulfonates treated for 24 h with native VP and its W164S variant and control without enzyme, and sulfonated polystyrene standards (c). Lignosulfonate samples (12 g L−1) after a 24-h treatment with 1.2 µM native VP (green line) and its W164S variant (dashes) in presence of 9.5 mM H2O2, and the corresponding softwood (red) and hardwood (blue) lignosulfonate controls without enzyme, were analyzed in a Superdex-75 column using 0.15 M NaOH as eluent (0.5 mL·min−1) and detection at 280 nm. Sulfonated polystyrenes (Mp 78,400, 29,500, 10,200 and 4210 Da, from left to right) were used as molecular mass standards in c (arrow shows the excluded blue dextran elution volume)

Then, the structural modifications of guaiacyl (G) and syringyl-guaiacyl (S–G) lignins from softwood and hardwood, respectively, were analyzed by 2D-NMR (Fig. 4). The main lignosulfonate units and side-chain interunit linkages are shown in Fig. 4g (no p-hydroxyphenyl units were detected). Both sulfonated (A) and non-sulfonated (A) β-O-4′ substructures were found in the control lignins, together with less abundant (non-sulfonated) phenylcoumaran (B) and resinol (C) substructures (Fig. 4a, d).

Fig. 4
figure 4

HSQC NMR spectra of softwood (a–c) and hardwood (d–f) lignosulfonates treated for 24 h with native VP (b, e) and its W164S variant (c, f) and control without enzyme (a, d), and formulae of the main structures identified (g). Signals correspond to 13C-1H correlations at the different positions of lignin native/α-oxidized/α-sulfonated syringyl (red signals) and guaiacyl (green signals) units, α-sulfonated/non-sulfonated side chains in β-O-4′ (blue signals), phenylcoumaran (cyan signals), and resinol (purple signals) substructures, and methoxyls (orange signal) (gray, unassigned signals). Same amount of sample (40 mg before enzymatic treatment) and DMSO-d 6 (0.75 mL) were used for all the spectra, which were normalized to the same intensity of the DMSO signal (not shown) for comparison. List of signalsCH ppm): 53.2/3.46, Cβ/Hβ in phenylcoumarans (B β ); 53.4/3.00, Cβ/Hβ in resinols (C β ); 55.5/3.66, C/H in methoxyls (MeO); 59.4/3.4 and 3.72, Cγ/Hγ in β–O–4′ (A γ ); 61.1/4.00, Cγ/Hγ in sulfonated β–O–4´ (A γ ); 65.6/3.93, Cα/Hα in sulfonated β–O–4′ linked to a G-unit (A α(G) ); 67.2/4.02, Cα/Hα in sulfonated β–O–4′ linked to a S-unit (A α(S) ); 70.8/4.16 and 3.77, Cγ/Hγ in β-β′ resinols (C γ ); 71.1/4.72, Cα/Hα in β–O–4′ linked to a G-unit (A α(G) ); 71.5/4.85, Cα/Hα in β–O–4′ linked to a S-unit (A α(S) ); 79.3/4.91, Cβ/Hβ in sulfonated β–O–4′ linked to a G unit (A β(G) ); 80.9/4.67, Cβ/Hβ in sulfonated β–O–4′ linked to a S unit (A β(S) ); 83.3/4.24, Cβ/Hβ in β–O–4′ linked to a G unit (A β(G) ); 84.9/4.59, Cα/Hα in β-β′ resinols (C α ); 85.7/4.08, Cβ/Hβ in β–O–4′ linked to a S unit (A β(S) ); 86.7/5.41, Cα/Hα in phenylcoumarans (B α ); 103.8/6.68, C2/H2 and C6/H6 in syringyl units (S 2,6 ); 106.2/7.29, C2/H2 and C6/H6 in α-oxidized syringyl units (S’ 2,6 ); 108.0/6.68, C2/H2 and C6/H6 in sulfonated syringyl units (S 2,6 ); 114.0/6.60 and 114.3/6.87, C2/H2 and C5/H5 in guaiacyl units (G 2 /G 5 ); and 122.8/6.75, C6/H6 in guaiacyl units (G 6 ) (minor, and largely overlapping, signals of C2/H2, C5/H5 and C6/H6 correlations in non-sulfonated guaiacyl units would appear at 110.7/6.93, 114.2/6.65 and 118.6/6.79 ppm, respectively; not shown). Three additional aromatic signals in the treated samples, at 126.1/7.14, 127.7/7.21 and 128.9/7.22 ppm, were assigned to protein (phenylalanine residues in the added enzyme)

After 24-h treatment with native VP (Fig. 4b, e) the intensity of the signals of the different aromatic lignin units (G, S, S and S ) and side-chain interunit linkages (A, A, B and C) decreased simultaneously, maintaining similar linkage percentages. However, the methoxyl numbers per unit increased up to ~ twofold. In the hardwood lignosulfonate, this was accompanied by higher abundance of Cα-oxidized syringyl units (S′) with respect to total syringyl units, while the S/G ratio also increased (from 2.0 in the control to 3.5 in the 24-h treated sample). Concerning side-chain signals, only those of the main sulfonated β-O-4′ substructures (A α , A β and A γ ) remained in the softwood lignosulfonate, while those of phenylcoumaran (B), resinol (C) and β-O-4′ (A) non-sulfonated side chains disappeared. In contrast, signals of sulfonated (A) and non-sulfonated β-O-4′ (A) and resinol (C) side chains could be observed in the hardwood lignosulfonate, albeit with low intensities. More interestingly, in the lignosulfonates treated for 24 h with the W164S variant (Fig. 4c, f) only minor changes in the aliphatic/aromatic HSQC signals were observed (spectra with similar intensities of most signals, and only slight increases of methoxyl content and S/G ratio compared with the control).

Steady-state treatment of nonphenolic vs native lignosulfonates

With the purpose of further investigating lignosulfonate modification by VP, including the observed small changes by the W164S variant, derivatized (nonphenolic) lignosulfonates were treated in new steady-state experiments.

The native VP was able to modify the nonphenolic lignosulfonates but the changes in the molecular-mass distribution (Additional file 1: Figure S4, green continuous line) and molecular structure of lignins (Additional file 1: Figure S5b, e) were modest, compared with those observed for the native (partially phenolic) lignosulfonates (Fig. 3a, b, green continuous line, and Fig. 4b, e, respectively). These modifications include lower-intensity signals in the NMR spectra of nonphenolic hardwood lignosulfonate (the S′ signal being the exception) and displacement of the Mp in the SEC profile, while lower changes were observed for the nonphenolic softwood lignosulfonate.

In contrast, the SEC profiles of the W164S-treated (green dashed lines) and control (red and blue lines) lignosulfonates were very similar (Additional file 1: Figure S4), and the same was observed in the 2D-NMR spectra (Additional file 1: Figure S5 c/f, a/d, respectively). The above results indicate that the phenolic lignin moiety: (i) promoted lignosulfonate modification under steady-state conditions; and (ii) was also responsible for the (small) modification of native lignin by the W164S variant.

Comparison with P. chrysosporium LiP

The transient state kinetic constants for reduction of P. chrysosporium LiP CII (the rate-limiting step in catalytic cycle) by native and methylated (nonphenolic) softwood and hardwood lignosulfonates were obtained (Additional file 1: Figure S6). Only the hardwood lignosulfonate (blue continuous line) showed saturation kinetics, with k 3 and K D3 values in the same order of those obtained for VP. In contrast, the k 3app value for softwood lignosulfonate was over one order of magnitude lower than found for VP. Methylation reduced the electron transfer rate (k 3app) and for the hardwood lignosulfonate this decrease was much higher than found for VP.

SEC and 2D-NMR spectroscopy of lignosulfonates treated with LiP revealed the same modification trends of the VP treatments. The shoulder (~11 mL) in softwood lignosulfonate SEC (Additional file 1: Figure S7a, red) was reduced without Mp displacement (green line), while in hardwood lignosulfonate (Additional file 1: Figure S7b) the Mp was displaced (~7500 Da) but the polymerization degree was lower than found with VP. Most signals in the 2D-NMR spectra of softwood (Additional file 1: Figure S8a–c) and hardwood (Additional file 1: Figure S8d–f) lignosulfonates showed slightly reduced intensities after LiP treatment, although the decreases were much lower than found with VP. The only exception was the (Cα-oxidized) S′ units increasing after the LiP treatment, as also found for VP. Moreover, the main changes were produced during the first 3 h, while with VP they progressed during the whole treatment.


Etherified lignin from phenolic monolignols

Although lignin precursors (monolignols) are phenolic (4-hydroxycinnamylic) alcohols, the resulting polymer is basically nonphenolic, since β-O-4′ and other ether interunit linkages involving the C4 hydroxyl (as found also in phenylcoumaran and 4-O-5′ structures) are predominant [35, 36]. In this way, only average 15 % phenylpropanoid units in native (natural) lignins would be phenolic [37], with the highest contents in grasses and conifers [38, 39]. This phenolic content increases in technical (industrial) lignins since cleavage of some interunit ether linkages is always produced in chemical pulping [40]. However, while kraft pulping (the most widespread process nowadays) is largely based on lignin (ether linkage) breakdown releasing a highly phenolic (up to 70 % or more) lignin, the sulfite process is based on lignin solubilization (after sulfonation) resulting in only moderate increases of the phenolic content [38, 41]. These two lignosulfonate characteristics (water solubility and phenolic content more similar to native lignin) were exploited in the present study to estimate electron transfer rates between (native and derivatized) lignin and ligninolytic peroxidases (including a mutated VP variant).

Electron transfer as seen from the peroxidase side (stopped-flow data)

Accurate kinetic constants for lignin degradation (by basidiomycete peroxidases) are difficult to be obtained under steady-state conditions due to the impossibility to follow lignin oxidation during short incubation periods (for maximal enzyme activity). However, the corresponding electron transfer rates can be precisely estimated (from the “peroxidase side”) by following the reduction of the H2O2-activated enzyme transient states (CI and CII) by lignin, using rapid spectrophotometry in single-turnover reactions under stopped-flow conditions.

Transient-state kinetic constants for P. chrysosporium LiP reduction by in vitro synthesized lignin (dehydrogenation polymer, DHP) had been reported (with k 3app ~60 s−1 mM−1) [26], the differences with the LiP constants obtained here being most probably related to differences in lignin preparations and solubility limitations using DHP. Also, DHP has a significant phenolic content [42] that will affect electron-transfer estimation, as shown here for lignosulfonates. Moreover, no mutated variants were included in these LiP studies [26] and, therefore, the catalytic residue/s remained unidentified. The first evaluation of several (three) possible LRET pathways for peroxidase oxidation of lignin was reported for P. eryngii VP [29] showing that only the pathway initiated at Trp164, homologous to LiP Trp171 [27], was operative. The VP and LiP site-directed mutagenesis studies used VA as a simple model for nonphenolic lignin. Other nonphenolic compounds (from dimers to tetramers) including the lignin most frequent linkages were used in subsequent studies [18, 20, 28, 43] but site-directed mutagenesis studies using the lignin polymer as substrate have been only recently reported, as discussed below.

Using water-soluble lignosulfonates, we estimated the reduction constants of P. eryngii VP transient states and, unexpectedly, some reduction of both CI and CII was observed for the W164S variant lacking the putative catalytic residue [32]. In the present study, we compared the transient-state kinetic constants of P. eryngii VP (and its W164S variant) and P. chrysosporium LiP on native (20–30 % phenolic) and nonphenolic (derivatized) softwood and hardwood lignosulfonates. With this purpose, samples were methylated with methyl iodide [44], which has advantages with respect to other methylating agents applied to lignosulfonates [45, 46].

First, we found that lignin methylation and acetylation—introducing ether (as found in nonphenolic lignin) and ester linkages at the phenolic hydroxyls, respectively—significantly lower the electron transfer rates, indicating that the phenolic units are easier to be oxidized by the enzyme. The above correlated with the lower lignin modification after steady-state treatment discussed below. Preferential degradation of the phenolic lignin moiety had been described after fungal decay by P. eryngii [47]. In spite of the above decrease of electron transfer rates, the constants for VP CI and CII reduction by the nonphenolic lignosulfonates (k 2app 100–320 and k 3app 80–600 s−1 mM−1) are much higher than reported for veratryl alcohol (k 2app 2.8 and k 3app 1.3 s−1 mM−1) [48]. This is mainly due to lower K D revealing that VP is more efficient binding polymeric lignin than simple aromatics. Moreover, although LiP is better reduced by veratryl alcohol [49, 50] than VP, its reduction constants by nonphenolic lignosulfonates are worst that found for VP, indicating that VP is more efficient than LiP abstracting electrons from nonphenolic lignin (under the present experimental conditions). This correlates with the significantly higher lignosulfonate modification found after VP treatment.

Second, and more importantly, we demonstrated that the solvent-exposed catalytic tryptophan (Trp164 of P. eryngii VP) is required for oxidizing the main nonphenolic lignin moiety, since CII reduction is practically absent in the W164S mutated variant. This is shown by both transient-state kinetic constants (50–60 fold lower k 3app values for nonphenolic than native lignin) and SEC and 2D-NMR results. Since they have a similar phenolic moiety, differences between CII reduction by the two native lignosulfonates could be related to the smaller size of the monomethoxylated units in softwood lignin, enabling contact and direct electron transfer to the heme cofactor at the main access channel. In contrast, we found that nonphenolic lignin can reduce the CI of the W164S variant, although with only 20–25 % efficiency compared with native VP. The above suggests that in native VP catalytic cycle (Additional file 1: Figure S1a) the Trp164 radical is required for nonphenolic lignin oxidation at the CII level (VP-IIB) while at the CI level both the porphyrin radical (VP-IA) and the Trp164 radical (VP-IB) would be able to oxidize nonphenolic lignin.

Additional aspects of lignin modification as shown by SEC and 2D-NMR

2D-NMR spectroscopy represents the state-of-the-art technology for structural characterization of lignins [5153], with broad application to lignin-engineered transgenic plants for biorefineries [54, 55]. This technique has been also used to study delignification of lignocellulosic feedstocks by fungal laccases in the presence of redox mediators [56, 57]. In a recent study, the authors used for the first time 2D-NMR to demonstrate lignosulfonate degradation by VP [32, 33].

After assigning the main signals of sulfonated and non-sulfonated lignin structures, their 2D-NMR spectra (normalized to the same amount of sample at the beginning of treatment and the same solution volume in the NMR tubes) showed (i) from small to large decreases in the intensity of the above signals and (ii) variable structural modifications of lignins, during their steady-state treatment (the extent of the above changes is clearly illustrated in the difference spectra of softwood and hardwood lignosulfonates—treated samples minus their controls—included as Additional file 1: Figure S9, S10, respectively). In laccase-mediator treatment of lignosulfonates, the decrease of HSQC signals was mainly due to the condensation reactions giving rise to quaternary (unprotonated) carbons [58]. However, degradation of lignin aromatic (and aliphatic) structures is produced during VP treatment, as shown by 13C NMR spectroscopy [32]. Unexpectedly, VP caused a stronger modification than LiP, resulting in the disappearance (or strong decline) of lignin signals. The observed increase of methoxyls (per aromatic unit) suggests the formation of non-aromatic methoxyl-containing (e.g. muconate type) structures [59]. The relative abundance of (Cα-oxidized) S′ units also increased in the treated lignins, as previously reported for the lignin-degrading laccase-mediator system [57, 60]. Such oxidation is among the first reactions in lignin biodegradation.

In contrast with the above results using native (unmodified) peroxidase, the VP variant lacking surface Trp164 only caused a modest modification of the NMR spectra, confirming that its lignin-degrading ability is largely associated to the presence of this surface residue. Moreover, when derivatized lignosulfonates were treated with the Trp164-less variant, the spectra were superimposable to those of the enzyme-less controls, demonstrating that this catalytic residue is strictly required for degradation of the nonphenolic lignin.

In addition to the structural modification revealed by 2D-NMR, the SEC profiles revealed repolymerization of a part of the products from lignin degradation by VP, resulting in residual lignins with increased molecular masses. This behavior, which is due to the coupling tendency of phenoxy and other aromatic radicals already reported in early “ligninase” studies [61], has been described for other oxidoreductases [6264], being especially remarkable in laccase-mediator treatments [58].


Data from stopped-flow (single turnover) analyses and steady-state treatments (the latter analyzed by SEC and 2D-NMR) of native and derivatized (nonphenolic) lignosulfonates unambiguously demonstrate that: (i) the minor phenolic moiety of lignin is preferentially degraded by ligninolytic VP; and (ii) a solvent exposed tryptophan residue (conserved in both VPs and LiPs) is required for electron transfer between the nonphenolic lignin and the H2O2 activated enzyme.


Enzyme production

Native VP from P. eryngii (mature protein-coding sequence of isoenzyme VPL2, GenBank AF007222) and its W164S mutated variant [29] were produced in Escherichia coli and in vitro activated as reported elsewhere [65]. The mature protein-coding sequence of P. chrysosporium LiP-H8 (GenBank Y00262) was also produced in E. coli and in vitro activated [66, 67].

The recombinant enzymes were purified by anion-exchange chromatography (Resource Q column, GE Healthcare, Uppsala, Sweden) using a 0–0.3 M NaCl gradient (2 mL min−1, 20 min) in 1 mM CaCl2-containing 10 mM tartrate, pH 5.5 (for VP and its W164S variant), or succinate, pH 6 (for LiP). The R z (A410/A280 ~4) values were indicative of the purity of the enzymes, and the electron absorption spectra confirmed the correct folding and cofactor incorporation.

Native and derivatized softwood and hardwood lignins

Two water-soluble sulfonated lignins were used in this study: softwood (Picea abies) and hardwood (Eucalyptus grandis) lignosulfonates kindly provided by G. E. Fredheim (Borregaard AS, Sapsborg, Norway). The lignosulfonate samples were dialyzed in 10 mM EDTA, 50 mM Tris (pH 8) with the aim of removing Mn2+ traces (which reduce H2O2-activated VP), and then in Milli-Q water.

Lignosulfonates (50 mg) were acetylated in a 50-mL pear-shaped flask with 3 mL of a pyridine-acetic anhydride (1:1, v/v) solution, stirring for 24 h at room temperature. Then, 10 mL of aqueous methanol (50 %) were added and the mixture was evaporated to dryness under vacuum. The solvent treatment was repeated three times with toluene (3 × 10 mL), and once with methanol (10 mL). Finally, the acetylated lignosulfonates (60–65 mg) were dried at 50 °C overnight. Acetylated lignosulfonates were used as enzyme substrate, and for estimation of phenolic and alcoholic hydroxyl content by NMR, as described below.

For lignosulfonates O-methylation with methyl iodide [44, 68], ~65 mg of sample were dissolved in 10 mL of dimethylsulfoxide (DMSO), methyl iodide (1 mL) and finely powdered NaOH (1 g) were added, and the mixture was vigorously vortexed for 10 min. Then, additional NaOH (300 mg) and methyl iodide (1 mL) were added, the mixture was stirred for 1 h, and the reaction quenched by adding 10 mL of water and adjusting the pH below 7 with 1 M HCl. The methylated lignosulfonates (45–55 mg) were dialyzed, concentrated under vacuum and freeze-dried.

Enzyme (transient-state) kinetics

Reduction of peroxidase CI and CII in 0.1 M tartrate (pH 3) by softwood and hardwood lignosulfonates (native and derivatized samples) was followed in a stopped-flow rapid spectrophotometry equipment (Bio-Logic, Claix, France) with a three-syringe module (SFM300) synchronized to a diode array detector (J&M, Essingen, Germany), and BioKine software.

CI reduction was studied by mixing the enzyme (1 µM final concentration) with H2O2 (1 µM final concentration) for 0.6 s, resulting in CI formation. Next, different amounts of lignosulfonate (5–350 µM final concentration) in 0.1 M (final concentration) tartrate (pH 3) were added, and the reactions were followed at 416 nm (isosbestic point of VP CII and resting state). CII reduction was studied by mixing a solution of enzyme and ferrocyanide (both at 1 µM final concentration) with H2O2 at equimolar ratio. The solution was aged for 6 s, and CII formation was achieved. Then, different amounts of lignosulfonate (5–350 µM final concentration) in 0.1 M (final concentration) tartrate (pH 3) were added, and the reaction was followed at 406 nm (Soret maximum of resting VP and LiP). The lignin concentrations in these and other experiments were referred to the basic phenylpropanoid unit in softwood and hardwood lignosulfonates.

All kinetic traces exhibited single-exponential character from which pseudo first-order rate constants (k 2obs and k 3obs for CI and CII reduction, respectively) were calculated. Plots of k 2obs and k 3obs vs substrate concentration fitted to linear or hyperbolic models. From those kinetics that fitted to a linear model apparent second-order rate constants (k 2app and k 3app for CI and CII reduction, respectively) were obtained. Plots of k obs vs substrate concentration that fitted to a Michaelis–Menten model yielded dissociation constants of the CI-lignin and CII-lignin complexes (K D2 and K D3, respectively) and first-order rate constants (k 2 and k 3, respectively). The corresponding apparent second-order rate constants, k 2app (k 2/K D2) and k 3app (k 3/K D3), were calculated with the equation: k obs = (k/K D)[S]/(1 + [S]/K D), where [S] indicates substrate concentration.

Lignin treatment under steady-state conditions

Lignosulfonates (12 g L−1) were treated with VP, its W164S variant, and LiP (all 1.2 µM concentration, added in two doses at the beginning and after 6 h of reaction) and H2O2 (9.5 mM, final concentration, added continuously over 24 h with a syringe pump) in 50 mM phosphate (pH 5), at 25 °C, and samples were taken after different times (3, 12 and 24 h). Control treatments were performed under the same conditions but in the absence of enzyme. Although VP and LiP show the highest activity at pH 3 (as used in stopped-flow experiments) the above long-term lignosulfonate treatments were performed at pH 5 (to maintain the enzyme active during the whole incubation period) after preliminary experiments where treatments at pH 3.5 and 5 were compared.

SEC analyses

Changes in the molecular-mass distribution of lignosulfonates after 24-h peroxidase treatment and controls were analyzed by SEC using a Superdex-75 column (HR-10/30, 3000–70,000/100,000 Da range; GE Healthcare) with 0.15 M NaOH as the mobile phase, at a flow rate of 0.5 mL·min−1, and UV (280 nm) detection. Blue dextran (Serva, Heindelberg, Germany) was used to determine the exclusion volume of the column, and a kit of sulfonated polystyrenes sodium salt standards with Mp in the 4210–976,000 Da range (PSS, Mainz, Germany) was used for calibration and mass determination (Ve/Vo vs Log[Mp], where Ve and Vo are the elution and void volumes respectively).

NMR analyses

Samples after different times (3, 12 and 24 h) of native and derivatized lignosulfonate treatment and the corresponding controls were freeze-dried for NMR analyses. Solution NMR spectra, including 1H-NMR and HSQC 2D-NMR, were recorded at 25 °C on an AVANCE III 500 MHz instrument (Bruker) equipped with a cryogenically cooled 5 mm TCI gradient probe with inverse geometry. The lignosulfonate samples (40 mg initial weight, before treatments) were dissolved in 0.75 mL of deuterated DMSO-d 6. The central solvent peak was used as the internal reference (at δCH 39.5/2.49 ppm), and the other signals were normalized to the same intensity of the DMSO signals (since the same DMSO volume and initial amount of sample was used in all the cases).

The HSQC experiment used Bruker’s “hsqcetgpsisp.2” adiabatic pulse program with spectral widths from 0 to 10 ppm (5000 Hz) and from 0 to 165 ppm (20,625 Hz) for the 1H and 13C dimensions. The number of transients was 64, and 256 time increments were always recorded in the 13C dimension. The 1 J CH used was 145 Hz. Processing used typical matched Gaussian apodization in the 1H dimension and squared cosine-bell apodization in the 13C dimension. Prior to Fourier transformation, the data matrices were zero-filled to 1024 points in the 13C dimension. Signals were assigned by literature comparison [32, 51, 58, 6972].

In the aromatic region of the spectrum, the C2–H2, C5–H5 and C6–H6 correlation signals were integrated to estimate the amount of lignins and the S/G ratio. In the aliphatic oxygenated region, the signals of methoxyls, and Cβ–Hβ (or Cα–Hα) correlations in the side chains of sulfonated and non-sulfonated β-O-4′, phenylcoumaran and resinol substructures were integrated. The intensity corrections introduced by the adiabatic pulse program permits to refer the latter integrals to the previously obtained number of lignin units.

The percentage of phenolic structures was calculated by referring the phenolic acetate signal in the HSQC 2D-NMR spectra (at 20.5/2.23 ppm) to the total number of lignin aromatic units (G + S + S′). To overcome differences in coupling constants of aliphatic and aromatic 13C-1H couples, the latter was estimated from the intensity of the methoxyl signal, taking into account the S/G ratio of the sample, and the number of methoxyls of G and S units [73].



compound I (of peroxidase catalytic cycle)


compound II (of peroxidase catalytic cycle)


dehydrogenation polymer (in vitro synthesized lignin)




ethylenediaminetetraacetic acid


guaiacyl (lignin unit)


heteronuclear single-quantum correlation

k 2 and k 3 :

first-order rate constants for CI and CII reduction, respectively

k 2app and k 3app :

apparent second-order rate constants for CI and CII reduction, respectively

K D2 and K D3 :

equilibrium dissociation constants for CI and CII reduction, respectively

k obs :

pseudo-first-order rate constant


lignin peroxidase


long-range electron transfer


main peak (in SEC)


nuclear magnetic resonance


pyrolysis-gas chromatography/mass spectrometry


syringyl (lignin unit)


size-exclusion chromatography


versatile peroxidase


  1. Martínez AT, Ruiz-Dueñas FJ, Martínez MJ, del Río JC, Gutiérrez A. Enzymatic delignification of plant cell wall: from nature to mill. Curr Opin Biotechnol. 2009;20:348–57.

    Article  Google Scholar 

  2. Bozell JJ, Petersen GR. Technology development for the production of biobased products from biorefinery carbohydrates-the US Department of Energy’s “Top 10” revisited. Green Chem. 2010;12:539–54.

    Article  CAS  Google Scholar 

  3. Ragauskas AJ, Beckham GT, Biddy MJ, Chandra R, Chen F, Davis MF, Davison BH, Dixon RA, Gilna P, Keller M, Langan P, Naskar AK, Saddler JN, Tschaplinski T, Tuskan GA, Wyman CE. Lignin valorization: improving lignin processing in the biorefinery. Science. 2014;344:1246843.

    Article  Google Scholar 

  4. Pandey A, Hofer R, Larroche C, Taherzadeh M, Nampoothiri M. Industrial biorefineries and white biotechnology. Amsterdam: Elsevier; 2015.

    Google Scholar 

  5. Shahid M, Mohammad F, Chen G, Tang RC, Xing T. Enzymatic processing of natural fibres: white biotechnology for sustainable development. Green Chem. 2016;18:2256–81.

    Article  CAS  Google Scholar 

  6. Martínez AT, Speranza M, Ruiz-Dueñas FJ, Ferreira P, Camarero S, Guillén F, Martínez MJ, Gutiérrez A, del Río JC. Biodegradation of lignocellulosics: microbiological, chemical and enzymatic aspects of fungal attack to lignin. Int Microbiol. 2005;8:195–204.

    Google Scholar 

  7. Kirk TK, Farrell RL. Enzymatic “combustion”: the microbial degradation of lignin. Annu Rev Microbiol. 1987;41:465–505.

    Article  CAS  Google Scholar 

  8. Kersten P, Cullen D. Extracellular oxidative systems of the lignin-degrading basidiomycete Phanerochaete chrysosporium. Fungal Genet Biol. 2007;44:77–87.

    Article  CAS  Google Scholar 

  9. Lewis NG, Yamamoto E. Lignin: occurrence, biogenesis and biodegradation. Annu Rev Plant Phys. 1990;41:455–96.

    Article  CAS  Google Scholar 

  10. Floudas D, Binder M, Riley R, Barry K, Blanchette RA, Henrissat B, Martínez AT, Otillar R, Spatafora JW, Yadav JS, Aerts A, Benoit I, Boyd A, Carlson A, Copeland A, Coutinho PM, de Vries RP, Ferreira P, Findley K, Foster B, Gaskell J, Glotzer D, Górecki P, Heitman J, Hesse C, Hori C, Igarashi K, Jurgens JA, Kallen N, Kersten P, et al. The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes. Science. 2012;336:1715–9.

    Article  CAS  Google Scholar 

  11. Ruiz-Dueñas FJ, Lundell T, Floudas D, Nagy LG, Barrasa JM, Hibbett DS, Martínez AT. Lignin-degrading peroxidases in Polyporales: an evolutionary survey based on ten sequenced genomes. Mycologia. 2013;105:1428–44.

    Article  Google Scholar 

  12. Riley R, Salamov AA, Brown DW, Nagy LG, Floudas D, Held BW, Levasseur A, Lombard V, Morin E, Otillar R, Lindquist EA, Sun H, LaButti KM, Schmutz J, Jabbour D, Luo H, Baker SE, Pisabarro AG, Walton JD, Blanchette RA, Henrissat B, Martin F, Cullen D, Hibbett DS, Grigoriev IV. Extensive sampling of basidiomycete genomes demonstrates inadequacy of the white-rot/brown-rot paradigm for wood decay fungi. Proc Natl Acad Sci USA. 2014;111:9923–8.

    Article  CAS  Google Scholar 

  13. Floudas D, Held BW, Riley R, Nagy LG, Koehler G, Ransdell AS, Younus H, Chow J, Chiniqui J, Lipzen A, Tritt A, Sun H, Haridas S, LaButti K, Ohm RA, Kües U, Blanchette RA, Grigoriev IV, Minto RE, Hibbett DS. Evolution of novel wood decay mechanisms in Agaricales revealed by the genome sequences of Fistulina hepatica and Cylindrobasidium torrendii. Fungal Genet Biol. 2015;76:78–92.

    Article  CAS  Google Scholar 

  14. Nagy LG, Riley R, Tritt A, Adam C, Daum C, Floudas D, Sun H, Yadav JS, Pangilinan J, Larsson KH, Matsuura K, Barry K, LaButti K, Kuo R, Ohm RA, Bhattacharya SS, Shirouzu T, Yoshinaga Y, Martin FM, Grigoriev IV, Hibbett DS. Comparative genomics of early-diverging mushroom-forming fungi provides insights into the origins of lignocellulose decay capabilities. Mol Biol Evol. 2016;33:959–70.

    Article  CAS  Google Scholar 

  15. Tien M, Kirk TK. Lignin-degrading enzyme from the hymenomycete Phanerochaete chrysosporium Burds. Science. 1983;221:661–3.

    Article  CAS  Google Scholar 

  16. Ruiz-Dueñas FJ, Martínez MJ, Martínez AT. Molecular characterization of a novel peroxidase isolated from the ligninolytic fungus Pleurotus eryngii. Mol Microbiol. 1999;31:223–36.

    Article  Google Scholar 

  17. Camarero S, Sarkar S, Ruiz-Dueñas FJ, Martínez MJ, Martínez AT. Description of a versatile peroxidase involved in natural degradation of lignin that has both Mn-peroxidase and lignin-peroxidase substrate binding sites. J Biol Chem. 1999;274:10324–30.

    Article  CAS  Google Scholar 

  18. Kirk TK, Tien M, Kersten PJ, Mozuch MD, Kalyanaraman B. Ligninase of Phanerochaete chrysosporium. Mechanism of its degradation of the non-phenolic arylglycerol β-aryl ether substructure of lignin. Biochem J. 1986;236:279–87.

    Article  CAS  Google Scholar 

  19. Busse N, Wagner D, Kraume M, Czemak P. Reaction kinetics of versatile peroxidase for the degradation of lignin compounds. Am J Biochem Biotechnol. 2013;9:365–97.

    Article  CAS  Google Scholar 

  20. Fernández-Fueyo E, Ruiz-Dueñas FJ, Martínez MJ, Romero A, Hammel KE, Medrano FJ, Martínez AT. Ligninolytic peroxidase genes in the oyster mushroom genome: heterologous expression, molecular structure, catalytic and stability properties and lignin-degrading ability. Biotechnol Biofuels. 2014;7:2.

    Article  Google Scholar 

  21. Kersten PJ, Tien M, Kalyanaraman B, Kirk TK. The ligninase of Phanerochaete chrysosporium generates cation radicals from methoxybenzenes. J Biol Chem. 1985;260:2609–12.

    CAS  Google Scholar 

  22. Cho DW, Parthasarathi R, Pimentel AS, Maestas GD, Park HJ, Yoon UC, Dunaway-Mariano D, Gnanakaran S, Langan P, Mariano PS. Nature and kinetic analysis of carbon-carbon bond fragmentation reactions of cation radicals derived from SET-oxidation of lignin model compounds. J Org Chem. 2010;75:6549–62.

    Article  CAS  Google Scholar 

  23. Martínez AT. Molecular biology and structure-function of lignin-degrading heme peroxidases. Enzyme Microb Technol. 2002;30:425–44.

    Article  Google Scholar 

  24. Hammel KE, Cullen D. Role of fungal peroxidases in biological ligninolysis. Curr Opin Plant Biol. 2008;11:349–55.

    Article  CAS  Google Scholar 

  25. Ruiz-Dueñas FJ, Martínez AT. Microbial degradation of lignin: how a bulky recalcitrant polymer is efficiently recycled in nature and how we can take advantage of this. Microbial Biotechnol. 2009;2:164–77.

    Article  Google Scholar 

  26. Johjima T, Itoh H, Kabuto M, Tokimura F, Nakagawa T, Wariishi H, Tanaka H. Direct interaction of lignin and lignin peroxidase from Phanerochaete chrysosporium. Proc Natl Acad Sci USA. 1999;96:1989–94.

    Article  CAS  Google Scholar 

  27. Doyle WA, Blodig W, Veitch NC, Piontek K, Smith AT. Two substrate interaction sites in lignin peroxidase revealed by site-directed mutagenesis. Biochemistry. 1998;37:15097–105.

    Article  CAS  Google Scholar 

  28. Mester T, Ambert-Balay K, Ciofi-Baffoni S, Banci L, Jones AD, Tien M. Oxidation of a tetrameric nonphenolic lignin model compound by lignin peroxidase. J Biol Chem. 2001;276:22985–90.

    Article  CAS  Google Scholar 

  29. Pérez-Boada M, Ruiz-Dueñas FJ, Pogni R, Basosi R, Choinowski T, Martínez MJ, Piontek K, Martínez AT. Versatile peroxidase oxidation of high redox potential aromatic compounds: site-directed mutagenesis, spectroscopic and crystallographic investigations of three long-range electron transfer pathways. J Mol Biol. 2005;354:385–402.

    Article  Google Scholar 

  30. Pogni R, Baratto MC, Teutloff C, Giansanti S, Ruiz-Dueñas FJ, Choinowski T, Piontek K, Martínez AT, Lendzian F, Basosi R. A tryptophan neutral radical in the oxidized state of versatile peroxidase from Pleurotus eryngii: a combined multi-frequency EPR and DFT study. J Biol Chem. 2006;281:9517–26.

    Article  CAS  Google Scholar 

  31. Smith AT, Doyle WA, Dorlet P, Ivancich A. Spectroscopic evidence for an engineered, catalytically active Trp radical that creates the unique reactivity of lignin peroxidase. Proc Natl Acad Sci USA. 2009;106:16084–9.

    Article  CAS  Google Scholar 

  32. Sáez-Jiménez V, Baratto MC, Pogni R, Rencoret J, Gutiérrez A, Santos JI, Martínez AT, Ruiz-Dueñas FJ. Demonstration of lignin-to-peroxidase direct electron transfer: a transient-state kinetics, directed mutagenesis, EPR and NMR study. J Biol Chem. 2015;290:23201–13.

    Article  Google Scholar 

  33. Sáez-Jiménez V, Baratto MC, Pogni R, Rencoret J, Gutiérrez A, Santos JI, Martínez AT, Ruiz-Dueñas FJ. Demonstration of lignin-to-peroxidase direct electron transfer. A transient-state kinetics, directed mutagenesis, EPR and NMR study. J Biol Chem. 2015;290:30268.

    Article  Google Scholar 

  34. Ruiz-Dueñas FJ, Morales M, Mate MJ, Romero A, Martínez MJ, Smith AT, Martínez AT. Site-directed mutagenesis of the catalytic tryptophan environment in Pleurotus eryngii versatile peroxidase. Biochemistry. 2008;47:1685–95.

    Article  Google Scholar 

  35. Fengel D, Wegener G. Wood: chemistry, ultrastructure, reactions. Berlin: Walter De Gruyter; 1984.

    Google Scholar 

  36. Ralph J, Lundquist K, Brunow G, Lu F, Kim H, Schatz PF, Marita JM, Hatfield RD, Ralph SA, Christensen JH, Boerjan W. Lignins: natural polymers from oxidative coupling of 4-hydroxyphenylpropanoids. Phytochem Rev. 2004;3:29–60.

    Article  CAS  Google Scholar 

  37. Kurek B. Potential applications of fungal peroxidases in the biological processing of wood, lignocelluloses and related compounds. In: Penel C, Gaspar T, Greppin H, editors. Plant peroxidases 1980–1990, topics and detailed literature on molecular, biochemical, and physiological aspects. Geneva: University of Geneva; 1992. p. 139–86.

    Google Scholar 

  38. Lai YZ. Determination of phenolic hydroxyl groups. In: Lin SY, Dence CW, editors. Methods in lignin chemistry. Berlin: Springer-Verlag; 1992. p. 423–34.

    Chapter  Google Scholar 

  39. Camarero S, Bocchini P, Galletti GC, Martínez AT. Pyrolysis-gas chromatography/mass spectrometry analysis of phenolic and etherified units in natural and industrial lignins. Rapid Commun Mass Spectrom. 1999;13:630–6.

    Article  CAS  Google Scholar 

  40. Sixta H, Potthast A, Krotschek AW. Chemical pulping processes. In: Sixta H, editor. Handbook of pulping. Weinheim: Wiley; 2006. p. 109–509.

    Chapter  Google Scholar 

  41. Broussard O, Petit M, Elie N, Baumberger S, Arnaud A, Ducrot PH, Allais F. Monitoring of free phenol content in lignosulfonates by ClO2 titration and UV difference spectroscopy. Holzforschung. 2016. doi:10.1515/hf-2015-0111.

    Google Scholar 

  42. Jacquet G, Pollet B, Lapierre C, Francesch C, Rolando C, Faix O. Thioacidolysis of enzymatic dehydrogenation polymers from p- hydroxyphenyl, guaiacyl, and syringyl precursors. Holzforschung. 1997;51:349–54.

    Article  Google Scholar 

  43. Hammel KE, Tien M, Kalyanaraman B, Kirk TK. Mechanism of oxidative Cα–Cβ cleavage of a lignin model dimer by Phanerochaete chrysosporium ligninase. J Biol Chem. 1985;260:8348–53.

    CAS  Google Scholar 

  44. Ciucanu I, Costello CE. Elimination of oxidative degradation during the per-O-methylation of carbohydrates. J Am Chem Soc. 2003;125:16213–9.

    Article  CAS  Google Scholar 

  45. Ritter DM, Olleman ED, Pennington DE, Wright KA. Functional derivatives of lignin sulfonic acid. J Am Chem Soc. 1950;72:1347–51.

    Article  Google Scholar 

  46. Hassi HY. Chemical characterization of lignosulfonates. Raleigh: NCSU; 1985 (PhD thesis).

    Google Scholar 

  47. Camarero S, Galletti GC, Martínez AT. Preferential degradation of phenolic lignin-units by two white-rot fungi. Appl Environ Microbiol. 1994;60:4509–16.

    CAS  Google Scholar 

  48. Ruiz-Dueñas FJ, Pogni R, Morales M, Giansanti S, Mate MJ, Romero A, Martínez MJ, Basosi R, Martínez AT. Protein radicals in fungal versatile peroxidase: catalytic tryptophan radical in both Compound I and Compound II and studies on W164Y, W164H and W164S variants. J Biol Chem. 2009;284:7986–94.

    Article  Google Scholar 

  49. Khindaria A, Yamazaki I, Aust SD. Veratryl alcohol oxidation by lignin peroxidase. Biochemistry. 1995;34:16860–9.

    Article  CAS  Google Scholar 

  50. Gelpke MDS, Lee J, Gold MH. Lignin peroxidase oxidation of veratryl alcohol: effects of the mutants H82A, Q222A, W171A, and F267L. Biochemistry. 2002;41:3498–506.

    Article  Google Scholar 

  51. Ralph J, Marita JM, Ralph SA, Hatfield RD, Lu F, Ede RM, Peng J, Quideau S, Helm RF, Grabber JH, et al. Solution-state NMR of lignin. In: Argyropoulos DS, editor. Advances in lignocellulosics characterization. Atlanta: Tappi Press; 1999. p. 55–108.

    Google Scholar 

  52. Ralph J, Landucci LL. NMR of lignin. In: Heitner C, Dimmel D, Schmidt J, editors. Lignin and lignans: advances in chemistry. Boca Raton: CRC Press; 2010. p. 137–243.

    Chapter  Google Scholar 

  53. Mansfield SD, Kim H, Lu FC, Ralph J. Whole plant cell wall characterization using solution-state 2D NMR. Nat Protoc. 2012;7:1579–89.

    Article  CAS  Google Scholar 

  54. Ralph J. Perturbing lignification. In: Entwistle K, Harris PJ, Walker J, editors. The compromised wood workshop. Christchurch: University of Canterbury; 2007. p. 85–112.

    Google Scholar 

  55. Vanholme R, Morreel K, Ralph J, Boerjan W. Lignin engineering. Curr Opin Plant Biol. 2008;11:278–85.

    Article  CAS  Google Scholar 

  56. Rico A, Rencoret J, del Río JC, Martínez AT, Gutiérrez A. Pretreatment with laccase and a phenolic mediator degrades lignin and enhances saccharification of Eucalyptus feedstock. Biotechnol Biofuels. 2014;7:6.

    Article  Google Scholar 

  57. Rico A, Rencoret J, del Río JC, Martínez AT, Gutiérrez A. In-depth 2D NMR study of lignin modification during pretreatment of Eucalyptus wood with laccase and mediators. Bioenerg Res. 2015;8:211–30.

    Article  CAS  Google Scholar 

  58. Prasetyo EN, Kudanga T, Ostergaard L, Rencoret J, Gutiérrez A, del Río JC, Santos JI, Nieto L, Jimenez-Barbero J, Martínez AT, Li JB, Gellerstedt G, Lepifre S, Silva C, Kim SY, Cavaco-Paulo A, Klausen BS, Lutnaes BF, Nyanhongo GS, Guebitz GM. Polymerization of lignosulfonates by the laccase-HBT (1-hydroxybenzotriazole) system improves dispersibility. Bioresour Technol. 2010;101:5054–62.

    Article  Google Scholar 

  59. Areskogh D, Li JB, Gellerstedt G, Henriksson G. Structural modification of commercial lignosulphonates through laccase catalysis and ozonolysis. Ind Crops Prod. 2010;32:458–66.

    Article  CAS  Google Scholar 

  60. Gutiérrez A, Rencoret J, Cadena EM, Rico A, Barth D, del Río JC, Martínez AT. Demonstration of laccase-mediator removal of lignin from wood and non-wood plant feedstocks. Bioresour Technol. 2012;119:114–22.

    Article  Google Scholar 

  61. Haemmerli SD, Leisola MSA, Fiechter A. Polymerisation of lignins by ligninases from Phanerochaete chrysosporium. FEMS Microbiol Lett. 1986;35:33–6.

    Article  CAS  Google Scholar 

  62. Hatakka A, Mettälä A, Toikka M, Hortling B, Brunow G. Modification of lignin by laccase and manganese peroxidase. In: Messner K, Srebotnik E, editors. Biotechnology in the pulp and paper industry: advances in applied and fundamental research. Vienna: Facultas-Universitätsverlag; 1996. p. 333–8.

    Google Scholar 

  63. Areskogh D, Li JB, Gellerstedt G, Henriksson G. Investigation of the molecular weight increase of commercial lignosulfonates by laccase catalysis. Biomacromolecules. 2010;11:904–10.

    Article  CAS  Google Scholar 

  64. Zhou HF, Yang DJ, Qiu XQ, Wu XL, Li Y. A novel and efficient polymerization of lignosulfonates by horseradish peroxidase/H2O2 incubation. Appl Microbiol Biotechnol. 2013;97:10309–20.

    Article  CAS  Google Scholar 

  65. Pérez-Boada M, Doyle WA, Ruiz-Dueñas FJ, Martínez MJ, Martínez AT, Smith AT. Expression of Pleurotus eryngii versatile peroxidase in Escherichia coli and optimisation of in vitro folding. Enzyme Microb Technol. 2002;30:518–24.

    Article  Google Scholar 

  66. Doyle WA, Smith AT. Expression of lignin peroxidase H8 in Escherichia coli: folding and activation of the recombinant enzyme with Ca2+ and haem. Biochem J. 1996;315:15–9.

    Article  CAS  Google Scholar 

  67. Fernández-Fueyo E, Ruiz-Dueñas FJ, Miki Y, Martínez MJ, Hammel KE, Martínez AT. Lignin-degrading peroxidases from genome of selective ligninolytic fungus Ceriporiopsis subvermispora. J Biol Chem. 2012;287:16903–16.

    Article  Google Scholar 

  68. Ciucanu I, Kerek F. A simple and rapid method for the permethylation of carbohydrates. Carbohyd Res. 1984;131:209–17.

    Article  CAS  Google Scholar 

  69. Lebo SE, Braaten SM, Fredheim GE, Lutnaes BF, Lauten RA, Myrvold BO, McNally TJ. Recent advances in the characterization of lignosulfonates. In: Hu T, editor. Characterization of lignocellulosic materials. New York: Blackwell Pub.; 2008. p. 189–205.

    Chapter  Google Scholar 

  70. Lutnaes BF, Myrvold BO, Lauten RA, Endeshaw MM. 1H and 13C NMR data of benzylsulfonic acids—model compounds for lignosulfonate. Magn Reson Chem. 2008;46:299–305.

    Article  CAS  Google Scholar 

  71. Rencoret J, Marques G, Gutiérrez A, Nieto L, Santos I, Jiménez-Barbero J, Martínez AT, del Río JC. HSQC-NMR analysis of lignin in woody (Eucalyptus globulus and Picea abies) and non-woody (Agave sisalana) ball-milled plant materials at the gel state. Holzforschung. 2009;63:691–8.

    Article  CAS  Google Scholar 

  72. Magina S, Marques AP, Evtuguin DV. Study on the residual lignin in Eucalyptus globulus sulphite pulp. Holzforschung. 2015;69:513–22.

    Article  CAS  Google Scholar 

  73. Martínez AT, Rencoret J, Marques G, Gutiérrez A, Ibarra D, Jiménez-Barbero J, del Río JC. Monolignol acylation and lignin structure in some nonwoody plants: a 2D NMR study. Phytochemistry. 2008;69:2831–43.

    Article  Google Scholar 

Download references

Authors’ contributions

VSJ and FJRD performed most of the biochemical experimental work and data analysis. JR performed the NMR study and sample derivatization. MARC contributed optimization of lignosulfonate methylation. AG contributed chemical analyses. All authors contributed to the discussion of results, and critically revised the manuscript. VSJ and ATM wrote the paper. All authors read and approved the final manuscript.


We thank Dr. Guro E. Fredheim (Borregaard AS, Sarpsborg, Norway) for providing the lignosulfonate preparations, and Dr. Manuel Angulo (CITIUS, University of Seville) for performing the NMR analyses. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

Competing interests

The authors declare that they have no competing interests.


This work was supported by the INDOX (KBBE-2013-613549 to ATM) and EnzOx2 (H2020-BBI-PPP-2015-RIA-720297 to ATM) EU projects, and the NOESIS (BIO2014-56388-R to FJR-D), BIORENZYMERY (AGL2014-53730-R to AG) and LIGNIN (CTQ2014-60764-JIN to JR) projects of the Spanish Ministry of Economy and Competitiveness (MINECO) co-financed by FEDER funds.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Angel T. Martínez.

Additional information

Verónica Sáez-Jiménez and Jorge Rencoret contributed equally to this work

Additional file


Additional file 1. Additional figures including VP cycle, and additional kinetic, Py-GC/MS, SEC and NMR results. Fig. S1. VP catalytic cycle and CI, CII and resting state electronic absorption spectra. Fig. S2. Kinetics of CI reduction by native, acetylated and permethylated softwood and hardwood lignosulfonates: Native VP vs W164S variant. Fig. S3. Lignosulfonate permethylation: Py-GC/MS of softwood lignosulfonate before and after 1 h methylation with methyl iodide. Fig. S4. SEC profiles of softwood and hardwood nonphenolic lignosulfonates treated for 24 h with native VP and its W164S variant and controls without enzyme. Fig. S5. HSQC NMR spectra of acetylated softwood and hardwood lignosulfonates treated for 24 h with native VP and its W164S variant, and control without enzyme. Fig. S6. Kinetics of reduction of LiP CII by native and permethylated softwood and hardwood lignosulfonates. Fig. S7. SEC profiles of softwood and hardwood lignosulfonates treated for 24 h with native LiP and controls without enzyme. Fig. S8. HSQC NMR spectra of native softwood and hardwood lignosulfonates treated for 3 and 24 h with LiP-H8, and the corresponding controls without enzyme. Fig. S9. Difference spectra of peroxidase-treated softwood lignosulfonates minus their controls. Fig. S10. Difference spectra of peroxidase-treated hardwood lignosulfonates minus their controls.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( 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

Sáez-Jiménez, V., Rencoret, J., Rodríguez-Carvajal, M.A. et al. Role of surface tryptophan for peroxidase oxidation of nonphenolic lignin. Biotechnol Biofuels 9, 198 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: