Boosting LPMO-driven lignocellulose degradation by polyphenol oxidase-activated lignin building blocks
© The Author(s) 2017
Received: 26 January 2017
Accepted: 3 May 2017
Published: 10 May 2017
Many fungi boost the deconstruction of lignocellulosic plant biomass via oxidation using lytic polysaccharide monooxygenases (LPMOs). The application of LPMOs is expected to contribute to ecologically friendly conversion of biomass into fuels and chemicals. Moreover, applications of LPMO-modified cellulose-based products may be envisaged within the food or material industry.
Here, we show an up to 75-fold improvement in LPMO-driven cellulose degradation using polyphenol oxidase-activated lignin building blocks. This concerted enzymatic process involves the initial conversion of monophenols into diphenols by the polyphenol oxidase MtPPO7 from Myceliophthora thermophila C1 and the subsequent oxidation of cellulose by MtLPMO9B. Interestingly, MtPPO7 shows preference towards lignin-derived methoxylated monophenols. Sequence analysis of genomes of 336 Ascomycota and 208 Basidiomycota reveals a high correlation between MtPPO7 and AA9 LPMO genes.
The activity towards methoxylated phenolic compounds distinguishes MtPPO7 from well-known PPOs, such as tyrosinases, and ensures that MtPPO7 is an excellent redox partner of LPMOs. The correlation between MtPPO7 and AA9 LPMO genes is indicative for the importance of the coupled action of different monooxygenases in the concerted degradation of lignocellulosic biomass. These results will contribute to a better understanding in both lignin deconstruction and enzymatic lignocellulose oxidation and potentially improve the exploration of eco-friendly routes for biomass utilization in a circular economy.
Fungal carbohydrate converting enzymes are considered important for eco-friendly application in plant biomass degradation. The resulting carbohydrates are sources for the production of biochemicals or biofuels and new enzymatically modified cellulose-based products are envisaged for the future.
Next to carbohydrates, phenolic compounds are also major components of lignocellulosic plant biomass. Phenolic compounds are present either as free molecules or in conjugated form as part of lignin or bound to carbohydrates. Lignin is one of the major constituents of (secondary) plant cell walls, together with the polysaccharides cellulose and hemicellulose. Lignin is composed of the three aromatic monolignol units: coniferyl, sinapyl, and para-coumaryl alcohol. Cellulose consists of β-(1→4)-linked glucosyl chains that interact with each other via hydrogen bonds and van der Waals forces, which results in the formation of crystalline cellulose fibrils. Unlike cellulose, hemicellulose is a heteropolymer varying in its monosaccharide composition and linkages between the monosaccharides. Examples are xylan, mannan, or β-(1→3, 1→4)-linked β-glucan. Hemicellulose interacts with lignin through ester and ether linkages, thereby forming a network that embeds the cellulose microfibrils [1, 2].
Recent studies focused on the function of lytic polysaccharide monooxygenases (LPMOs) have confirmed that these enzymes drive the oxidative degradation of cellulose, and they are considered important for the enzymatic degradation of plant biomass . LPMOs are classified based on their sequence in the Carbohydrate Active enzyme (CAZy; ) database as auxiliary activity (AA) families AA9, AA10, AA11, and AA13. In brief, LPMOs have been reported to oxidize β-(1→4)-linked glucan chains at either the C1- or C4-carbon position or both, resulting in the cleavage of glucan chains [3, 5, 6]. Other LPMOs of these AA families have been described to oxidize the (1→4)-linkage of chitin, xylan, hemicellulose, such as xyloglucan and glucomannan, soluble cellodextrins, and starch [7–11].
In order to oxidize polysaccharides, LPMOs demand electrons to activate molecular oxygen in their copper-containing active site [3, 5, 6]. The electrons can be donated by reducing agents, like low molecular weight compounds (ascorbic acid, gallic acid) or the macromolecule lignin [12–14]. Other ways of providing LPMOs with electrons have been reported, such as flavocytochrome-dependent cellobiose dehydrogenases (CDHs), light-induced pigments, light-driven chemical oxidation of water, or diphenol-regenerating GMC (glucose-methanol-choline-oxidase/dehydrogenase)-oxidoreductases [5, 12, 15, 16]. The mechanistic understanding of electron donation systems is highly relevant in order to enable optimization of LPMO activity and, thereby, plant biomass degradation.
Phenolic compounds, including small molecular weight compounds that serve as lignin building blocks and lignin, are intrinsically present in plant biomass and are natural electron donors for LPMO activity. However, monophenols are not optimal electron donors for LPMO activity, because of their relatively high redox potential . Compounds with a 1,2-benzenediol or a 1,2,3-benzenetriol moiety have, compared to monophenols, a lower redox potential . Their low redox potential enables them to reduce the copper ion in the active site of LPMOs and enhance the LPMO activity .
The enzymatic oxidation of phenolic compounds is a well-known reaction in many natural environments. For example, it causes browning of food products, contributes to taste in tea fermentation, and plays a role in plant biomass decomposition [17–19]. These oxidative systems involve the activity of laccases (EC.18.104.22.168), peroxidases (EC.1.11.1), tyrosinases (EC.22.214.171.124), and catechol oxidases (EC 126.96.36.199) [18–20]. In particular tyrosinases, also often referred to as polyphenol oxidases, are of interest due to their ability to hydroxylate phenolic compounds . This so-called monophenolase activity typically involves the ortho-hydroxylation of monophenols into ortho-diphenols, compounds that comprise a 1,2-benzenediol moiety. Tyrosinases also exhibit diphenolase activity, which is characterized by the oxidation of these ortho-diphenols into ortho-quinones . With respect to biomass degradation, the diphenolase activity of polyphenol oxidases is not conductive to polysaccharides oxidation since ortho-quinones cannot be utilized by LPMOs . Based on their monophenolase activity, tyrosinases have been shown to use non-methoxylated monophenols as substrates, rather than using methoxylated monophenols, that are the predominant structural units of lignin .
In this study, we investigated if polyphenol oxidases, in particular tyrosinases, can enhance LPMO activity. Therefore, the tyrosinase MtPPO7 from Myceliophthora thermophila C1 was used, which showed activity towards phenolic compounds and was obtained from DuPont Industrial Biosciences. In addition, we used the commercially available tyrosinase AbPPO from the white button mushroom Agaricus bisporus. Both tyrosinases were each incubated with the previously described MtLPMO9B in the presence of various plant phenolic compounds and cellulose [13, 24]. We found that MtPPO7 is highly active towards methoxylated monophenols, including monomeric lignin building blocks. This activity can strongly boost the release of oxidized gluco-oligosaccharides as catalyzed by MtLPMO9B and thereby the degradation of cellulose. In addition, we found a strong correlation between genes encoding MtPPO like proteins and AA9 LPMOs in fungal genomes of 336 Ascomycota and 208 Basidiomycota.
MtPPO7 addition improves cellulose oxidation by MtLPMO9B
RAC incubated with MtLPMO9B only and with AbPPO or MtPPO7 addition in the presence of 21 reducing agents
Activity [%] MtLPMO9Ba
Activity [%] MtLPMO9B + AbPPOb
Activity [%] MtLPMO9B + MtPPO7b
As expected, from our previous research , monophenols (group Ia) are less efficient electron donors for MtLPMO9B, with the exception of 4-hydroxybenzoic acid (no. 1). The low MtLPMO9B activity towards RAC resulted in a very small amount of released non-oxidized and C1-oxidized gluco-oligosaccharides (Fig. 2; Table 1; Additional file 1: Figure S1). The addition of MtPPO7 to RAC incubated with MtLPMO9B did not increase the electron donor capacities of this group Ia compounds, except for para-coumaric acid (no. 3). In contrast, the addition of MtPPO7 to RAC incubated with MtLPMO9B in the presence of methoxylated monophenols (group Ib, compounds with a 1-hydroxy, 2-methoxy moiety) led in all cases to a high increase (up to ±75 times, no. 5) in the release of non-oxidized and C1-oxidized gluco-oligosaccharides (Fig. 2; Additional file 1: Figure S1). Similarly, the oxidative cleavage of RAC incubated with MtLPMO9B increased, when MtPPO7 was added to group Ic compounds, which comprise a 1-hydroxy-2,6-dimethoxy moiety. For group IIa, IIb, and III phenolic compounds, the addition of MtPPO7 to incubations of RAC with MtLPMO9B led to a moderate increase of non-oxidized and C1-oxidized gluco-oligosaccharides (4–39%), with one exception (no. 19) (Fig. 2; Table 1).
Incubations performed with AbPPO showed deviating results from corresponding experiments with MtPPO7. First, unlike MtPPO7, addition of AbPPO to non-methoxylated monophenols (group Ia) resulted in increased levels of non-oxidized and C1-oxidized oligosaccharides during the incubation of RAC with MtLPMO9B. This different activity was expected on the basis of the known hydroxylating capacity of AbPPO towards non-methoxylated monophenolic compounds. Especially, the addition of AbPPO in the presence of para-coumaric acid (no. 3) increased MtLPMO9B-catalyzed degradation of RAC compared to the incubation without AbPPO (Fig. 2; Table 1). In contrast to MtPPO7, AbPPO addition to group Ib, Ic, II, and III compounds improved the MtLPMO9B-catalyzed RAC degradation relatively moderate (5–361%, Table 1) or even decreased the RAC degradation by 1–44% (Fig. 2; Table 1).
MtLPMO9B-mediated cellulose oxidation in the presence of MtPPO7 using time-dependent incubations
Activity of MtPPO7 towards methoxylated phenolic compounds
The most striking observation from the experiments described above is the conversion of methoxylated phenolic compounds by MtPPO7 into products, which enhance the oxidative degradation of cellulose by MtLPMO9B. The fate of these phenolic compounds upon MtPPO7 incubation was further studied in the absence of MtLPMO9B using UHPLC-UV-MSn.
MtPPO7 was active towards all eight methoxylated monophenols (group Ib) tested. Based on UV, the concentrations of group Ib compounds decreased during MtPPO7 incubation between 28 and 67% compared to the same compounds incubated without MtPPO7 (Fig. 4). Three types of MtPPO7 reactions were observed based on mass differences between substrate and products formed (Additional file 4: Table S1): (A) hydroxylation (m/z +16; e.g., no. 5 and 7), (B) decarboxylation (m/z −44, e.g., no. 8), and (C) demethylation (m/z −14, e.g., no. 11). Especially the hydroxylation of methoxylated phenolic compounds is a key reaction, since the products formed comprise a second hydroxyl group. The formed ortho-diphenols are known to be efficient electron donors for LPMOs. Decarboxylation and demethylation occurred either in the presence of MtPPO7 or as a result of polymerization reactions in both the presence and absence of MtPPO7. Several group Ib compounds spontaneously formed dimers in the absence of MtPPO7 (no. 8 or no. 11). These dimers were almost absent when group Ib compounds were incubated with MtPPO7 (Additional file 4: Table S1). Similar to group Ib compounds, the concentrations of compounds comprising a 1-hydroxyl-2,6-dimethoxy moiety (group Ic) incubated with MtPPO7 decreased, based on UV, between 6% and 64% compared to the incubation without MtPPO7 (Fig. 4; Additional file 4: Table S1). Again, masses indicating decarboxylation (−44 Da) and demethylation (−14 Da) reactions were formed. MtPPO7 was also able to convert group IIb compounds (up to 87% substrate conversion, no. 20) that comprise a 1,2-dihydroxy-3-methoxy moiety (Fig. 4). The reactions observed were decarboxylation (m/z −44, e.g., no. 19) and dimerization (e.g., no. 20) of group IIb compounds based on masses formed (Additional file 4: Table S1).
In general, most phenolic compounds (such as no. 17 and no. 21) that were incubated with MtPPO7 formed insoluble complexes, which likely result from polymerization reactions of ortho-quinones formed. These insoluble complexes were not determined by UHPLC-UV-MSn. The complexation of ortho-quinones resulted, obviously, in a decrease of the UV signal of the substrates (Fig. 4; Additional file 5: Figure S4).
Conversion of guaiacol and 3-methylcatechol by MtPPO7
For a better discrimination between monophenolase and diphenolase activity, MtPPO7 conversion of guaiacol and 3-methylcatechol was monitored over a period of 24 h by UHPLC-UV-MSn. The conversion of guaiacol by MtPPO7 resulted in the initial formation of 3-methoxycatechol (Additional file 4: Table S1; Additional file 5: Figure S4, Additional file 6: Figure S5). Further reactions resulted in the formation of brown pigments indicating that 3-methoxycatechol was oxidized into ortho-quinones, which are likely to polymerize and form insoluble complexes. Other masses determined by UHPLC-UV-MSn indicated the polymerization of guaiacol and compounds originating from the oxidation of guaiacol by MtPPO7 into trimers (m/z 399, 401, 415). Based on the masses detected (Additional file 4: Table S1; Additional file 5: Figure S4), a scheme is presented of possible reaction pathways of guaiacol occurring during MtPPO7 incubation (Additional file 6: Figure S5).
The oxidation of guaiacol by MtPPO7 mainly occurred between 4 and 8 h of incubation, whereas 3-methylcatechol showed to be oxidized by MtPPO7 within the first 2 h. During the incubation of 3-methylcatechol with MtPPO7 pink pigments are formed, which precipitate after sample centrifugation. Based on MSn, masses of 3-methylcatechol (m/z 123) and masses indicating the dimerization of 3-methylcatechol (m/z 245) were detected upon incubation of both 3-methylcatechol with MtPPO7 and 3-methylcatechol alone (Additional file 4: Table S1). No detectable amounts of new products were formed during the incubation of 3-methylcatechol with MtPPO7 compared to the incubation of 3-methylcatechol only (Additional file 4: Table S1).
Structural model of MtPPO7
Genome-wide analysis of AA9 LPMO, AbPPO, and MtPPO7
Principal component analysis (PCA) was performed on all 336 Ascomycota and 208 Basidiomycota. In addition, we used the numbers which describe the presence of AA9 LPMO, MtPPO7s, and AbPPOs encoding genes in each fungal species.
MtPPO7 enhances cellulose oxidation by MtLPMO9B
Many PPOs, tyrosinases in particular, have been described to oxidize monophenols and various phenolic compounds comprising a 1,2-benzenediol or 1,2,3-benzenetriol moiety. However, none of these tyrosinases show a high activity towards methoxylated compounds (group Ib) such as ferulic acid (no. 8) [33–36]. Importantly, MtPPO7 hydroxylates methoxylated phenolic compounds (group Ib) and thereby improves the activity of MtLPMO9B up to 75 times (no. 5, Table 1). On the other hand, AbPPO addition to RAC incubated with MtLPMO9B in the presence non-methoxylated monophenols also resulted in a significant increase of the MtLPMO9B activity (up to 99 times). Especially the activity of MtPPO7 towards methoxylated phenolic compounds is of high relevance, as these compounds are abundant as structural lignin units and, therefore, intrinsically present during plant biomass degradation. Here, we only determined the effect of polyphenol oxidases on the LPMO-mediated cellulose degradation in the presence of phenolic compounds that are known as lignin building blocks but we did not, for example, investigate the impact of these enzymes in the presence of washed or unwashed pretreated biomass. In comparison to the substrate used in this thesis, we expect a lower effect of MtPPO7-like enzymes on the LPMO-mediated substrate oxidation using pretreated biomass, since this biomass contains a mixture of multiple potential electron-donating compounds, which were described in an earlier study .
A recent study described that H2O2, in addition to or instead of O2 (Fig. 1), acts as a co-substrate for LPMOs, while reducing agents are still needed to activate the active site copper . Once copper is activated, LPMOs oxidize substrates under the use of H2O2. The latter can be formed by the reduction of O2 through reducing agents . Considering these recent findings, it can only be hypothesized how MtPPO7 increased the MtLPMO9B-mediated cellulose oxidation, which was investigated in our research. It is possible that the oxidation of the phenolic compounds led to the formation of compounds that have either an enhanced reducing efficiency on MtLPMO9B or led to an enhanced H2O2 generation which increased the MtLPMO9B activity.
Activity of MtPPO7 towards phenolic compounds
Based on our results, we conclude that MtPPO7 improves the activity of MtLPMO9B by two main reactions. First, MtPPO7 hydroxylates methoxylated monophenols at the ortho-position (monophenolase activity) and forms compounds comprising a 1,2-benzenediol moiety. These compounds have, compared to methoxylated monophenols, a lower redox potential and are known to be good electron donors for LPMOs [12–14]. Secondly, MtPPO7 inhibits the formation of non-enzymatic coupling reactions, which occurred during the incubation of several cinnamic acid derivatives in the absence of MtPPO7, but were not formed when MtPPO7 was present during the incubations (Additional file 4: Table S1). Non-enzymatic coupling reactions of cinnamic acid derivatives of group I, II, or III result in the formation of bulky phenolic polymers. The aliphatic acrylic acid (prop-2-enoic acid) group present in cinnamic acid derivatives constitutes an elongation of the conjugated aromatic ring. This group takes part in polymerization reactions, which can be caused by radical formation due to the presence of copper (2.5 µM) [18, 38]. However, it remains unclear, how MtPPO7 prevents (reduces) the polymerization of cinnamic acid derivatives.
Structural characteristics of MtPPO7
The structural model of MtPPO7 shows a four-helix bundle architecture with the presence of six histidine residues coordinating the two copper ions in the active site, which is typical for PPO-like tyrosinases and catechol oxidases (Fig. 5). We clearly showed that the preference of MtPPO7 towards methoxylated phenolic compounds and the low activity towards diphenols is different from reported fungal or plant catechol oxidases and tyrosinases (Fig. 4). Here, we further discuss the structure–function relationship of these PPOs in comparison with the generated model of MtPPO7. However, it should be noted that the structural model of MtPPO7 was generated based on a catechol oxidase from Aspergillus oryzae (AoCO4) and both proteins share an overall amino acid sequence identity of only 38% . It has already been shown that marginal structural differences determine whether PPOs exhibit a stronger monophenolase or diphenolase activity . Based on our model, it is possible that structural features discussed below deviate compared to the crystal structure which is, however, not available yet.
Firstly, the presence of a thioether bridge between a cysteine and histidine is assumed to be important for the monophenolase and diphenolase activity. This bond has been described to be present in many fungal PPOs, such as AbPPO3/4 from A. bisporus, NcPPO from Neurospora crassa, and also in plant catechol oxidases, such as IbCO from Ipomoea batatas [26, 44, 45]. However, similar to AoCO4 from A. oryzae , MtPPO7 does not form this thioether bond (Fig. 5). Secondly, MtPPO7 shows the presence of the PPO typical ‘gate residue’ (Phe300 in MtPPO7, Phe261 in IbCO, and Val299 in AoCO4). It has been hypothesized that a bulky ‘gate residue,’ such as a phenylalanine, limits the accessibility of monophenols to the copper A ion (CuA) and, therefore, might be a limiting factor for the monophenolase activity [44, 46]. However, in AoCO4 the bulky phenylalanine is replaced by Val299 and still AoCO4 shows a rather low activity towards monophenols [31, 42]. Our data suggest that Phe300 in MtPPO7 does not limit the accessibility for (methoxylated) monophenols. Thirdly, it also has recently been proposed that the presence of certain amino acid residues next to the CuB-coordinating histidines indicates whether a tyrosinase or catechol oxidase exhibits a stronger monophenolase and a weak diphenolase activity or a weak monophenolase and a strong diphenolase activity . However, none of the described amino acid residues at that specific position (Pro, Asn, Glu, and Gln) is present in MtPPO7, which has a tyrosine (Tyr268) at the corresponding position. In MtPPO7 a tyrosine (Tyr268) occupies. Whether Tyr268 is responsible for the observed specificity towards methoxylated monophenols can only be hypothesized at this moment. Fourthly, MtPPO7 shares with tyrosinases and catechol oxidases the characteristic ‘tyrosine motif’ (Tyr-X-Tyr/Phe or Tyr/Phe-X-Tyr; in MtPPO7 residues 399-401). This motif is highly conserved among plants and fungi [31, 42, 47]. Fifthly, MtPPO7 does not contain the ‘YG motif’ (Gly-Tyr motif), which is a typical feature of fungal tyrosinases and catechol oxidases . In conclusion, based on substrate specificity (Fig. 1) and modeled structure (Fig. 5), we do not classify MtPPO7 as a typical tyrosinase or catechol oxidase.
Based on genome analysis, we found a positive correlation between genes encoding AA9 LPMOs and MtPPO7-like enzymes in Ascomycota and Basidiomycota, which has not been shown previously. Interestingly, this correlation strengthens the evidence that fungi benefit from the concerted activity of AA9 LPMOs and MtPPO7-like enzymes in nature, which is in agreement with major findings obtained from the experimental data of this work. The proportion of annotated AA9 LPMOs within Ascomycota and Basidiomycota (82 and 86%, respectively) is about 5–10% lower with the findings of a recent study (92% in Ascomycota and Basidiomycota), which may result from the different species chosen (Fig. 6) . The two PPOs (MtPPO7 and AbPPO) used in this work share only a low sequence identity and, therefore, are considered to belong to different subgroups of polyphenol oxidases. It is important to realize that active tyrosinases, such as the mushroom tyrosinase AbPPO3, have been described to consist of large (H) and small (L) subunits [45, 48]. We focused for the genome analysis on genes encoding the protein of the catalytic H subunit, which contains the highly conserved binuclear copper site. The function of the L subunit is so far not known . All purified enzymes used in this work are secreted. However, we did not differentiate between secreted and non-secreted fungal AA9 LPMOs and PPOs in the genome analysis. It may be that most of the enzymes detected here take part in the external fungal metabolism, as can be deducted from the high abundance of the annotated MtPPO7s and AbPPOs throughout Ascomycota and Basidiomycota. We, therefore, specifically selected genomes of plant biomass degrading fungi that contain at least ten genes encoding known cellulose-degrading enzymes per fungus and found again a similar correlation between genes encoding AA9 LPMOs, AbPPOs, and MtPPO7s (Fig. 7) [32, 49]. In contrast to MtPPO7s, genes encoding AbPPOs showed a weak correlation with AA9 LPMOs in Ascomycota and Basidiomycota (Fig. 7), especially for the selected cellulose-degrading Basidiomycota (Fig. 7c). One possible explanation is that a cellulose-rich environment is associated with the abundance of lignin as seen in soft- and hardwoods. Still, there may be considerable variation in the substrate specificity and specific activity of PPOs within the MtPPO7- or AbPPO-like gene families. The latter seems to be consistent with the low sequence identity among PPOs and, in addition, by the diverse substrate specificities of PPOs that have only marginal differences in their amino acid sequence [31, 43, 50].
For the first time, we demonstrate the importance of the coupled action of different monooxygenases in the concerted degradation of plant biomass. We demonstrated that MtPPO7 is particularly active towards methoxylated phenolic compounds that are the predominant structural units of lignin. This feature distinguishes MtPPO7 from the well-known mushroom tyrosinase AbPPO, which stimulates the LPMO activity via its ability to hydroxylate non-methoxylated monophenols. However, its strong diphenolase activity limits the applicability of this tyrosinase for producing electron-donating capacity for LPMOs. In addition, we established that genes encoding MtPPO7-like enzymes and AA9 LPMOs are highly correlated throughout Ascomycota and Basidiomycota, suggesting that AA9 LPMOs benefit from the activity of MtPPO7-like enzymes in nature. Further understanding in both lignin deconstruction and enzymatic lignocellulose oxidation will lead to more eco-friendly routes for biomass utilization, which is seen as a prerequisite for a circular economy.
Enzyme expression, production, and purification
MtPPO7 (UniProt: KX772412) was over-expressed in the homologous Myceliophthora thermophila C1 strain. A low protease/low (hemi-)cellulase-producing M. thermophila C1 strain was used to produce MtPPO7 [51, 52]. The MtPPO7-containing culture broth was fractionated to obtain a pure MtPPO7 preparation. This preparation was provided by DuPont Industrial Biosciences. MtLPMO9B was expressed and purified as described in .
The pure MtPPO7 fraction was analyzed by LC–mass spectrometry confirming the presence of the MtPPO7 by ‘The Scripps Research Institute’ (San Diego, CA, USA).
Purification and identification of mushroom tyrosinase
Tyrosinase from the edible button mushroom Agaricus bisporus was purified from a commercial enzyme preparation (Sigma-Aldrich, Steinheim, Germany) as described previously . The purified enzyme preparation (referred to as AbPPO) was shown to contain the isoforms PPO3 and PPO4 .
Cellulose substrate and reducing agents
Regenerated amorphous cellulose (RAC) was prepared from Avicel PH-101 as described previously [7, 53]. All reducing agents used throughout this study (Table 1) were purchased from Sigma-Aldrich (Steinheim, Germany).
Incubation conditions for MtLPMO9B, MtPPO7, and AbPPO
Regenerated amorphous cellulose (1–2 mg mL−1, see Figure captions) was dissolved in 50 mM ammonium acetate buffer (pH = 5.0), with or without addition of reducing agents (final concentration of 2 mM). The standard enzyme concentrations of MtLPMO9B, MtPPO7, and AbPPO used in this work were 5.0, 5.0, and 0.7 µg protein mg−1 substrate, respectively. All samples were incubated for 20 h at 50 °C in a head-over-tail rotator in portions of 1 mL total volume (Stuart rotator, Bibby Scientific, Stone, UK) at 20 rpm. Supernatants of all incubations were analyzed by HPAEC and RP-UHPLC-UV-ESI–MSn.
Oligosaccharides were analyzed by high-performance anion-exchange chromatography (HPAEC) with pulsed amperometric detection (PAD) using a HPAEC system (ICS-5000, Dionex, Sunnyvale, CA, USA) as described previously .
Supernatants of all incubations were subjected to an Accela reversed-phase high-performance liquid chromatography (RP-UHPLC) system coupled to electron spray ionization mass spectrometry (Thermo Scientific, San Jose, CA, USA). Injected samples (5 μL) were separated using an Acquity C18 column (2.1 × 150 mm, 1.7 μm particle size) with an Acquity UHPLC Shield RP18 Vanguard guard column (2.1 × 5 mm, 1.7 μm particle size). Both columns were purchased from Waters (Milford, MA, USA). Gradient elution with eluent A (H2O + 1% (v/v) acetonitrile + 0.1% (v/v) HOAc) and eluent B (acetonitrile + 0.1% (v/v) HOAc) was performed according to the following steps: From 0 to 17.7 min a linear gradient from 5 to 60% B; from 17.7 to 21.7 min, isocratic 100% B, and from 21.7 to 26 min, isocratic 5% B. The flow rate and the injection volume were 0.4 mL min−1 and 5 μL, respectively. The column temperature was set to 30 °C and the photodiode array detector was operated in the range of 200–400 nm.
Samples were further analyzed using an LTQ-Velos mass spectrometer (Thermo Scientific) equipped with a ESI–MS. Data were collected over a m/z range of 90–1500 in both negative (NI) and positive (PI) modes. The collision energy was set to 35%.
An alignment was made of the amino acid sequence of MtPPO7 and the amino acid sequence of catechol oxidase from Aspergillus oryzae (AoCO4), which scored highest in a Blast search using the MtPPO7 sequence against the Protein Data Bank (38% amino acid identity). Using this alignment and the available structure of AoCO4 (PDB id: 4J3P) as template, structural models were obtained for MtPPO7 using the Modeller program version 9.16. One hundred comparative models were generated, after which the model with lowest corresponding DOPE score was selected for inspection and image generation using Pymol (Pymol, The PyMOL Molecular Graphics System, version 188.8.131.52 Schrödinger, LLC, New York, NY, USA).
Fungal genomes were obtained from the JGI MycoCosm portal . In total, protein sequence annotations of 336 Ascomycota and 208 Basidiomycota genomes were used. BLAST databases for those protein sequences were created using ‘makeblastdb’ program in BLAST + v2.2.30 . The protein BLAST was performed separately using ten MtPPO7 sequences [PDB id: 4j3p and closely related MtPPO7s: Q2UNF9, Aspergillus oryzae (strain ATCC 42149); G2QC95, Myceliophthora thermophila (ATCC 42464); Q2H7I7, Chaetomium globosum (ATCC 6205); G0SFX8, Chaetomium thermophilum (DSM 1495); L7IAQ4, Magnaporthe oryzae (strain Y34); L7JMT9, Magnaporthe oryzae (strain P131); G4N2I5, Magnaporthe oryzae (strain 70-15); A0A084GCK1, Scedosporium apiospermum; A0A0C4DYF2, Magnaporthiopsis poae (ATCC 64411); J3P591, Gaeumannomyces graminis var. tritici (strain R3-111a-1)], ten AA9 LPMO sequences (PDB id: 4d7u and closely related LPMOs: Q7SHI8, Neurospora crassa (strain ATCC 24698); G2QCJ3, Myceliophthora thermophila (strain ATCC 42464); F7W1P4, Sordaria macrospora (strain ATCC MYA-333); G2RB73, Thielavia terrestris (strain ATCC 38088); Q2H8N9, Chaetomium globosum (strain ATCC 6205); G0S408, Chaetomium thermophilum (strain DSM 1495); F8MLY8, Neurospora tetrasperma (strain FGSC 2508); T0L448, Colletotrichum gloeosporioides (strain Cg-14); A0A0H4K9X4 and A0A1C9CXI0, Myceliophthora thermophila C1); and four AbPPOs (AbPPO1, AbPPO2, AbPPO3, and AbPPO4 [55, 56]) as query sequences. Resulting sequences below E-value cut-off of 0.001 with query coverage above 60% for AA9 LPMOs, 65% for MtPPO7s, and 85% for AbPPOs were considered for further analysis. Selection of cellulase-rich Ascomycota and Basidiomycota was based on the presence of at least 10 genes encoding cellulose-degrading enzymes, which are classified in the CAZy database as glycosyl hydrolase families GH1, GH3, GH5, GH6, GH7, GH12, GH45. The GH gene families were selected based on Kubicek et al. . Previous data  were used to determine the number of annotated genes encoding cellulose-degrading enzymes. Based on this selection, 27 Ascomycota and 23 Basidiomycota species were selected for the Pearson correlation analysis (Fig. 7; Additional file 7: Table S2, Additional file 8: Table S3). All the statistical analyses were performed in R .
MF, MJK, JV, SWAH, DW, J-PV, and MAK designed research; MF performed research and analyzed data; SKM and MF performed genome analysis; AHW and MF performed structural modeling; and MF, SKM, AHW, WJHB, HG, and MAK wrote the paper. All authors read and approved the final manuscript.
This research was funded by The Graduate School VLAG (Advanced Studies in Food Technology, Agrobiotechnology, Nutrition and Health Sciences) of Wageningen University & Research (Wageningen, The Netherlands) and by a VICI Grant (865.14.001) of the Netherlands Organization for Scientific Research (NWO) to D.W.
The authors declare that they have no competing interests.
Availability of supporting data
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lam TBT, Kadoya K, Iiyama K. Bonding of hydroxycinnamic acids to lignin: ferulic and p-coumaric acids are predominantly linked at the benzyl position of lignin, not the β-position, in grass cell walls. Phytochemistry. 2001;57(6):987–92. doi:10.1016/S0031-9422(01)00052-8.View ArticleGoogle Scholar
- Takahashi N, Koshijima T. Ester linkages between lignin and glucuronoxylan in a lignin-carbohydrate complex from beech (Fagus crenata) wood. Wood Sci Technol. 1988;22(231–241):231–41. doi:10.1007/BF00386018.View ArticleGoogle Scholar
- Vaaje-Kolstad G, Westereng B, Horn SJ, Liu Z, Zhai H, Sorlie M, et al. An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Science. 2010;330(6001):219–22. doi:10.1126/science.1192231.View ArticleGoogle Scholar
- Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014:42;D490–5. doi:10.1093/nar/gkt1178.View ArticleGoogle Scholar
- Phillips CM, Beeson WT, Cate JH, Marletta MA. Cellobiose dehydrogenase and a copper-dependent polysaccharide monooxygenase potentiate cellulose degradation by Neurospora crassa. ACS Chem Biol. 2011;6(12):1399–406. doi:10.1021/cb200351y.View ArticleGoogle Scholar
- Quinlan RJ, Sweeney MD, Lo Leggio L, Otten H, Poulsen JC, Johansen KS, et al. Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. Proc Natl Acad Sci USA. 2011;108(37):15079–84. doi:10.1073/pnas.1105776108.View ArticleGoogle Scholar
- Frommhagen M, Sforza S, Westphal AH, Visser J, Hinz SW, Koetsier MJ, et al. Discovery of the combined oxidative cleavage of plant xylan and cellulose by a new fungal polysaccharide monooxygenase. Biotechnol Biofuels. 2015;8:101. doi:10.1186/s13068-015-0284-1.View ArticleGoogle Scholar
- Vu VV, Beeson WT, Span EA, Farquhar ER, Marletta MA. A family of starch-active polysaccharide monooxygenases. Proc Natl Acad Sci USA. 2014;111(38):13822–7. doi:10.1073/pnas.1408090111.View ArticleGoogle Scholar
- Isaksen T, Westereng B, Aachmann FL, Agger JW, Kracher D, Kittl R, et al. A C4-oxidizing lytic polysaccharide monooxygenase cleaving both cellulose and cello-oligosaccharides. J Biol Chem. 2014;289(5):2632–42. doi:10.1074/jbc.M113.530196.View ArticleGoogle Scholar
- Forsberg Z, Rohr AK, Mekasha S, Andersson KK, Eijsink VG, Vaaje-Kolstad G, et al. Comparative study of two chitin-active and two cellulose-active AA10-type lytic polysaccharide monooxygenases. Biochemistry. 2014;53(10):1647–56. doi:10.1021/bi5000433.View ArticleGoogle Scholar
- Agger JW, Isaksen T, Varnai A, Vidal-Melgosa S, Willats WG, Ludwig R, et al. Discovery of LPMO activity on hemicelluloses shows the importance of oxidative processes in plant cell wall degradation. Proc Natl Acad Sci USA. 2014;111(17):6287–92. doi:10.1073/pnas.1323629111.View ArticleGoogle Scholar
- Kracher D, Scheiblbrandner S, Felice AKG, Breslmayr E, Preims M, Ludwicka K, et al. Extracellular electron transfer systems fuel cellulose oxidative degradation. Science. 2016. doi:10.1126/science.aaf3165.Google Scholar
- Frommhagen M, Koetsier MJ, Westphal AH, Visser J, Hinz SWA, Vincken J-P, et al. Lytic polysaccharide monooxygenases from Myceliophthora thermophila C1 differ in substrate preference and reducing agent specificity. Biotechnol Biofuels. 2016;9(1):1–17. doi:10.1186/s13068-016-0594-y.View ArticleGoogle Scholar
- Westereng B, Cannella D, Wittrup Agger J, Jørgensen H, Larsen Andersen M, Eijsink VGH, et al. Enzymatic cellulose oxidation is linked to lignin by long-range electron transfer. Scientific Rep. 2015;5:18561. doi:10.1038/srep18561.View ArticleGoogle Scholar
- Bissaro B, Forsberg Z, Ni Y, Hollmann F, Vaaje-Kolstad G, Eijsink VGH. Fueling biomass-degrading oxidative enzymes by light-driven water oxidation. Green Chem. 2016;18:5357–66. doi:10.1039/c6gc01666a.View ArticleGoogle Scholar
- Cannella D, Mollers KB, Frigaard NU, Jensen PE, Bjerrum MJ, Johansen KS, et al. Light-driven oxidation of polysaccharides by photosynthetic pigments and a metalloenzyme. Nat Commun. 2016;7:11134. doi:10.1038/ncomms11134.View ArticleGoogle Scholar
- Gao Z-J, Liu J-B, Xiao X-G. Purification and characterisation of polyphenol oxidase from leaves of Cleome gynandra L. Food Chem. 2011;129(3):1012–8. doi:10.1016/j.foodchem.2011.05.062.View ArticleGoogle Scholar
- Dec J, Haider K, Bollag J-M. Decarboxylation and demethylation of naturally occurring phenols during coupling reactions and polymerization. Soil Sci. 2001;166(10):660–71.View ArticleGoogle Scholar
- Kuijpers TFM, van Herk T, Vincken J-P, Janssen RH, Narh DL, van Berkel WJH, et al. Potato and mushroom polyphenol oxidase activities are differently modulated by natural plant extracts. J Agric Food Chem. 2014;62(1):214–21. doi:10.1021/jf4043375.View ArticleGoogle Scholar
- Beloqui A, Pita M, Polaina J, Martínez-Arias A, Golyshina OV, Zumárraga M, et al. Novel polyphenol oxidase mined from a metagenome expression library of bovine rumen: biochemical properties, structural analysis, and phylogenetic relationships. J Biol Chem. 2006;281(32):22933–42. doi:10.1074/jbc.M600577200.View ArticleGoogle Scholar
- Lerch K. Neurospora tyrosinase: structural, spectroscopic and catalytic properties. Mol Cell Biochem. 1983;52(2):125–38. doi:10.1007/bf00224921.View ArticleGoogle Scholar
- Konrad L. Tyrosinase: molecular and active-site structure. Enzymatic browning and its prevention. ACS Symposium Series: American Chemical Society; 1995. p. 64–80.Google Scholar
- Vanholme R, Demedts B, Morreel K, Ralph J, Boerjan W. Lignin biosynthesis and structure. Plant Physiol. 2010;153(3):895–905. doi:10.1104/pp.110.155119.View ArticleGoogle Scholar
- Kuijpers TFM, Narváez-Cuenca C-E, Vincken J-P, Verloop AJW, van Berkel WJH, Gruppen H. Inhibition of enzymatic browning of chlorogenic acid by sulfur-containing compounds. J Agric Food Chem. 2012;60(13):3507–14. doi:10.1021/jf205290w.View ArticleGoogle Scholar
- Kuijpers TFM, Gruppen H, Sforza S, van Berkel WJH, Vincken J-P. The antibrowning agent sulfite inactivates Agaricus bisporus tyrosinase through covalent modification of the copper-B site. FEBS J. 2013;280(23):6184–95. doi:10.1111/febs.12539.View ArticleGoogle Scholar
- Mauracher SG, Molitor C, Al-Oweini R, Kortz U, Rompel A. Latent and active abPPO4 mushroom tyrosinase cocrystallized with hexatungstotellurate(VI) in a single crystal. Acta Crystallographica Sect D. 2014;70(9):2301–15. doi:10.1107/S1399004714013777.View ArticleGoogle Scholar
- Mauracher SG, Molitor C, Al-Oweini R, Kortz U, Rompel A. Crystallization and preliminary X-ray crystallographic analysis of latent isoform PPO4 mushroom (Agaricus bisporus) tyrosinase. Acta Crystallographica Secti F Struct Biol Commun. 2014;70(Pt 2):263–6. doi:10.1107/s2053230x14000582.View ArticleGoogle Scholar
- Wu J, Chen H, Gao J, Liu X, Cheng W, Ma X. Cloning, characterization and expression of two new polyphenol oxidase cDNAs from Agaricus bisporus. Biotechnol Lett. 2010;32(10):1439–47. doi:10.1007/s10529-010-0329-2.View ArticleGoogle Scholar
- Grigoriev IV, Nikitin R, Haridas S, Kuo A, Ohm R, Otillar R, et al. MycoCosm portal: gearing up for 1000 fungal genomes. Nucleic Acids Res. 2014;42(Database issue):D699–704. doi:10.1093/nar/gkt1183.View ArticleGoogle Scholar
- Berka RM, Grigoriev IV, Otillar R, Salamov A, Grimwood J, Reid I, et al. Comparative genomic analysis of the thermophilic biomass-degrading fungi Myceliophthora thermophila and Thielavia terrestris. Nat Biotechnol. 2011;29(10):922–7. doi:10.1038/nbt.1976.View ArticleGoogle Scholar
- Gasparetti C, Faccio G, Arvas M, Buchert J, Saloheimo M, Kruus K. Discovery of a new tyrosinase-like enzyme family lacking a C-terminally processed domain: production and characterization of an Aspergillus oryzae catechol oxidase. Appl Microbiol Biotechnol. 2010;86(1):213–26. doi:10.1007/s00253-009-2258-3.View ArticleGoogle Scholar
- Zhao Z, Liu H, Wang C, Xu J-R. Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genom. 2013;14(1):1–15. doi:10.1186/1471-2164-14-274.View ArticleGoogle Scholar
- Selinheimo E, NiEidhin D, Steffensen C, Nielsen J, Lomascolo A, Halaouli S, et al. Comparison of the characteristics of fungal and plant tyrosinases. J Biotechnol. 2007;130(4):471–80. doi:10.1016/j.jbiotec.2007.05.018.View ArticleGoogle Scholar
- Zhou P, Smith NL, Lee CY. Potential purification and some properties of Monroe apple peel polyphenol oxidase. J Agric Food Chem. 1993;41(4):532–6. doi:10.1021/jf00028a004.View ArticleGoogle Scholar
- Paul B, Gowda LR. Purification and characterization of a polyphenol oxidase from the seeds of field bean (Dolichos lablab). J Agric Food Chem. 2000;48(9):3839–46. doi:10.1021/jf000296s.View ArticleGoogle Scholar
- Zou Y, Hu W, Jiang A, Ma K. Partial purification and characterization of a novel extracellular tyrosinase from Auricularia auricula. Appl Biochem Biotechnol. 2013;172(3):1460–9. doi:10.1007/s12010-013-0638-8.View ArticleGoogle Scholar
- Bissaro B, Rohr AK, Skaugen M, Forsberg Z, Horn SJ, Vaaje-Kolstad G, et al. Fenton-type chemistry by a copper enzyme: molecular mechanism of polysaccharide oxidative cleavage. bioRxiv. 2016. doi:10.1101/097022.Google Scholar
- Rosazza JPN, Huang Z, Dostal L, Volm T, Rousseau B. Review: biocatalytic transformations of ferulic acid: an abundant aromatic natural product. J Ind Microbiol. 1995;15(6):457–71. doi:10.1007/bf01570016.View ArticleGoogle Scholar
- Schweigert N, Zehnder AJB, Eggen RIL. Chemical properties of catechols and their molecular modes of toxic action in cells, from microorganisms to mammals. Environ Microbiol. 2001;3(2):81–91. doi:10.1046/j.1462-2920.2001.00176.x.View ArticleGoogle Scholar
- Raymond J, Rakariyatham N, Azanza JL. Purification and some properties of polyphenoloxidase from sunflower seeds. Phytochemistry. 1993;34(4):927–31. doi:10.1016/S0031-9422(00)90689-7.View ArticleGoogle Scholar
- Lee CY, Smith NL, Pennesi AP. Polyphenoloxidase from DeChaunac grapes. J Sci Food Agric. 1983;34(9):987–91. doi:10.1002/jsfa.2740340915.View ArticleGoogle Scholar
- Hakulinen N, Gasparetti C, Kaljunen H, Kruus K, Rouvinen J. The crystal structure of an extracellular catechol oxidase from the ascomycete fungus Aspergillus oryzae. JBIC, J Biol Inorg Chem. 2013;18(8):917–29. doi:10.1007/s00775-013-1038-9.View ArticleGoogle Scholar
- Solem E, Tuczek F, Decker H. Tyrosinase versus catechol oxidase: one asparagine makes the difference. Angewandte Chemie Int Edit. 2016;55(8):2884–8. doi:10.1002/anie.201508534.View ArticleGoogle Scholar
- Klabunde T, Eicken C, Sacchettini JC, Krebs B. Crystal structure of a plant catechol oxidase containing a dicopper center. Nat Struct Biol. 1998;5(12):1084–90. doi:10.1038/4193.View ArticleGoogle Scholar
- Ismaya WT, Rozeboom HJ, Weijn A, Mes JJ, Fusetti F, Wichers HJ, et al. Crystal structure of Agaricus bisporus mushroom tyrosinase: identity of the tetramer subunits and interaction with tropolone. Biochemistry. 2011;50(24):5477–86. doi:10.1021/bi200395t.View ArticleGoogle Scholar
- Decker H, Schweikardt T, Tuczek F. The first crystal structure of tyrosinase: all questions answered? Angewandte Chemie Int Edit. 2006;45(28):4546–50. doi:10.1002/anie.200601255.View ArticleGoogle Scholar
- Marusek CM, Trobaugh NM, Flurkey WH, Inlow JK. Comparative analysis of polyphenol oxidase from plant and fungal species. J Inorg Biochem. 2006;100(1):108–23. doi:10.1016/j.jinorgbio.2005.10.008.View ArticleGoogle Scholar
- Strothkamp KG, Jolley RL, Mason HS. Quaternary structure of mushroom tyrosinase. Biochem Biophys Res Commun. 1976;70(2):519–24. doi:10.1016/0006-291X(76)91077-9.View ArticleGoogle Scholar
- Kubicek CP, Kubicek EM. Enzymatic deconstruction of plant biomass by fungal enzymes. Curr Opin Chem Biol. 2016;35:51–7. doi:10.1016/j.cbpa.2016.08.028.View ArticleGoogle Scholar
- Goldfeder M, Kanteev M, Adir N, Fishman A. Influencing the monophenolase/diphenolase activity ratio in tyrosinase. Biochimica et Biophysica Acta Proteins Proteom. 2013;1834(3):629–33. doi:10.1016/j.bbapap.2012.12.021.View ArticleGoogle Scholar
- Visser H, Joosten V, Punt PJ, Gusakov AV, Olson PT, Joosten R, et al. Development of a mature fungal technology and production platform for industrial enzymes based on a Myceliophthora thermophila isolate, previously known as Chrysosporium lucknowense C1. Ind Biotechnol. 2011;7:214–23. doi:10.1089/ind.2011.0003.View ArticleGoogle Scholar
- Punt PJ, Burlingame RP, Pynnonen CM, Olson PT, Wery J, Visser J, Heinrich, et al. Chrysosporium lucknowense protein production system. Patent WO/2010/107303; 2010.Google Scholar
- Zhang YHP, Cui J, Lynd LR, Kuang LR. A transition from cellulose swelling to cellulose dissolution by o-phosphoric acid: evidence from enzymatic hydrolysis and supramolecular structure. Biomacromol. 2006;7(2):644–8. doi:10.1021/bm050799c.View ArticleGoogle Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST + : architecture and applications. BMC Bioinform. 2009;10(1):1–9. doi:10.1186/1471-2105-10-421.View ArticleGoogle Scholar
- Wichers HJ, Recourt K, Hendriks M, Ebbelaar CEM, Biancone G, Hoeberichts FA, et al. Cloning, expression and characterisation of two tyrosinase cDNAs from Agaricus bisporus. Appl Microbiol Biotechnol. 2003;61(4):336–41. doi:10.1007/s00253-002-1194-2.View ArticleGoogle Scholar
- N-y Li, W-m Cai, Q-l Jin, Q-p Qin, F-l Ran. Molecular cloning and expression of polyphenoloxidase genes from the mushroom, Agaricus bisporus. Agric Sci China. 2011;10(2):185–94. doi:10.1016/S1671-2927(09)60305-9.View ArticleGoogle Scholar
- Kubicek CP, Starr TL, Glass NL. Plant cell wall–degrading enzymes and their secretion in plant-pathogenic fungi. Ann Rev Phytopathol. 2014;52(1):427–51. doi:10.1146/annurev-phyto-102313-045831.View ArticleGoogle Scholar
- Team RC. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016. https://www.R-project.org.