# Discovery of fungal oligosaccharide-oxidising flavo-enzymes with previously unknown substrates, redox-activity profiles and interplay with LPMOs

Majid Haddad Momeni, Folmer Fredslund, Bastien Bissaro, Olanrewaju Raji, Thu Vuong, Sebastian Meier, Tine Sofie Nielsen, Vincent Lombard, Bruno Guigliarelli, Frédéric Biaso, et al.

► To cite this version:

Majid Haddad Momeni, Folmer Fredslund, Bastien Bissaro, Olanrewaju Raji, Thu Vuong, et al. | Discovery of fungal oligosaccharide-oxidising flavo-enzymes with previously unknown substrates, redox-activity profiles and interplay with LPMOs. Nature Communications, 2021, 12 (1), 10.1038/s41467-021-22372-0. | hal-03234953

HAL Id: hal-03234953

https://hal.inrae.fr/hal-03234953

Submitted on 26 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-021-22372-0

OPEN

<figure><img src="image_3.png" /></figure>

Check for updates

# Discovery of fungal oligosaccharide-oxidising flavo-enzymes with previously unknown substrates, redox-activity profiles and interplay with LPMOs

Majid Haddad Momeni<sup>1,2</sup>, Folmer Fredslund<sup>1,3</sup>, Bastien Bissaro<sup>2</sup>, Olanrewaju Raji<sup>4</sup>, Thu V. Vuong<sup>4</sup>, Sebastian Meier<sup>5</sup>, Tine Sofie Nielsen<sup>1</sup>, Vincent Lombard<sup>6</sup>, Bruno Guigliarelli<sup>7</sup>, Frédéric Biaso<sup>7</sup>, Mireille Haon<sup>2</sup>, Sacha Grisel<sup>2</sup>, Bernard Henrissat<sup>1,6,8,9</sup>, Ditte Hededam Welner<sup>1,3</sup>, Emma R. Master<sup>4,10</sup>, Jean-Guy Berrin<sup>2,10</sup> & Maher Abou Hachem<sup>1,10</sup>

Oxidative plant cell-wall processing enzymes are of great importance in biology and biotechnology. Yet, our insight into the functional interplay amongst such oxidative enzymes remains limited. Here, a phylogenetic analysis of the auxiliary activity 7 family (AA7), currently harbouring oligosaccharide flavo-oxidases, reveals a striking abundance of AA7-genes in phytopathogenic fungi and Oomycetes. Expression of five fungal enzymes, including three from unexplored clades, expands the AA7-substrate range and unveils a cellooligosaccharide dehydrogenase activity, previously unknown within AA7. Sequence and structural analyses identify unique signatures distinguishing the strict dehydrogenase clade from canonical AA7 oxidases. The discovered dehydrogenase directly is able to transfer electrons to an AA9 lytic polysaccharide monooxygenase (LPMO) and fuel cellulose degradation by LPMOs without exogenous reductants. The expansion of redox-profiles and substrate range highlights the functional diversity within AA7 and sets the stage for harnessing AA7 dehydrogenases to fine-tune LPMO activity in biotechnological conversion of plant feedstocks.

<footnote n="1">¹Department of Biotechnology and Biomedicine, Technical University of Denmark, Lyngby, Denmark. ²INRAE, Aix Marseille Univ, Biodiversité et Biotechnologie Fongiques (BBF), Marseille, France. ³The Novo Nordisk Center for Biosustainability, Lyngby, Denmark. ⁴Department of Chemical Engineering and Applied Chemistry, University of Toronto, Toronto, ON, Canada. ⁵Department of Chemistry, Technical University of Denmark, Lyngby, Denmark.</footnote>

<footnote n="6">⁶Architecture et Fonction des Macromolécules Biologiques, UMR 7257 CNRS, USC 1408, Aix Marseille Univ, Marseille, France. ⁷Aix-Marseille Univ, CNRS, UMR7281 Unité de Bioénergétique et Ingénierie des Protéines (BIP), Marseille, France. ⁸INRAE, USC1408 Architecture et Fonction des Macromolécules Biologiques (AFMB), Marseille, France. ⁹Department of Biological Sciences, King Abdulaziz University, Jeddah, Saudi Arabia. ¹⁰Department of Bioproducts and Biosystems, Aalto University, Espoo, Finland. ☉email: jean-guy.berrin@inrae.fr; maha@bio.dtu.dk</footnote>
The involvement of oxidative processes in polysaccharide degradation by fungi has been proposed by the pioneering work of Eriksson et al. in 1974¹. This notion has gained strong support by the recent discovery of lytic polysaccharide monooxygenases (LPMOs) that uniquely catalyse the oxidative cleavage of glycosidic bonds in (semi)crystalline polysaccharides such as starch²⁻⁴, chitin⁵, cellulose⁶⁻⁸ and cellulose-bound hemicelluloses, e.g., xyloglucan⁹ and xylan¹⁰. Besides LPMOs, filamentous fungi co-secrete an impressing diversity of carbohydrate-specific oxidoreductases¹¹. To date, only four fungal oligosaccharide-oxidising enzymes from the auxiliary activity family 7 (AA7) in the Carbohydrate Active enZyme (CAZy) database¹², have been characterized. In addition, oligosaccharide oxidases from plants have been reported¹³,¹⁴, but are currently not assigned into AA7. Our insight is, thus, clearly limited regarding the biological roles as well as the diversity of substrates and redox features within this family. Fungal AA7 enzymes catalyse the oxidation of the reducing end C1-OH in e.g. cellooligosaccharides¹⁵ and lactose¹⁶, xylooligosaccharides¹⁷ as well as chitooligosaccharides¹⁸ to the corresponding lactones. Electrons derived from oligosaccharide oxidation reduce the FAD cofactor that is subsequently re-oxidised via electron transfer to O₂ to generate H₂O₂ (oxidase activity). Notably, comparable dehydrogenase and oxidase activities have been observed for an AA7 enzyme via electron transfer to artificial electron acceptors instead of O₂¹⁶.

The tertiary structures of hitherto described AA7 enzymes¹⁵⁻¹⁷,¹⁹ share a common fold comprising an N-terminal FAD-binding domain (F domain) and a C-terminal substrate-binding domain (S domain). This fold is common within the vanillyl alcohol oxidase (VAO, EC. 1.1.3.38) super family²⁰ that harbours AA7. All previously characterised AA7 oxidases are distinguished by a cysteinyl and histidyl bi-covalently tethered FAD cofactor. By contrast, other VAO family members harbour mostly a monovalently (or less commonly a non-covalently) bound FAD cofactor²⁰.

Reactive oxygen species (ROS), and especially H₂O₂, play important roles in lignocellulose degradation by fungi, but the underpinning molecular details of these roles remain poorly understood¹¹. Although O₂ has long been considered as the co-substrate of LPMOs⁵,⁷, recent findings suggest that H₂O₂ is the more favourable co-substrate during polysaccharide oxidative cleavage²¹⁻²³. LPMO catalysis is mediated by a Cu cofactor⁶ that must be reduced from Cu(II) to Cu(I) to prime the reaction. Enzymatic priming of LPMOs by the modular pyrroloquinoline-quinone-dependent pyranose dehydrogenase²⁴ (CAZy family AA12 dehydrogenase domain appended to an AA8 cytochrome *b* haem domain) and FAD-dependant glucose-methanol-choline (GMC) superfamily of oxidoreductases (AA3), most notably the fungal cellobiose dehydrogenase (CDH) that also possesses a cytochrome *b* haem domain²⁵⁻²⁷, has been reported. The activity on cellooligosaccharides and the transcriptional co-regulation as well as the co-secretion with cellulose active AA9 LPMOs²⁸, justified extensive studies on CDH as a model for direct electron transfer and priming of LPMOs²⁵,²⁷,²⁹. Not all fungi possess the CDH/LPMO pair, suggesting the presence of additional redox partners and mechanisms for LPMO activation. By analogy, the co-secretion of AA7s with LPMOs upon fungal growth on plant biomass³⁰,³¹ prompted us to hypothesize a redox interplay between these two enzyme classes.

Here, we report phylogenetic analyses, combined with the selection and characterisation of five fungal AA7s, three of which belong to previously unexplored clades. We demonstrate activity on four saccharides, previously not reported as AA7 substrates. In-depth analysis of a cellooligosaccharide dehydrogenase with a mono histidyl-tethered FAD highlights clade-dependant redox-profiles

within AA7. Importantly, we unveil direct activation and potentiation of LPMO activity on cellulose by this newly-discovered dehydrogenase. This study provides biochemical, structural and mechanistic insights into AA7 enzymes as components of the fungal redox network secreted during growth on biomass. Our findings suggest a way to tune LPMO activity for enzymatic degradation of major recalcitrant polysaccharides using AA7 dehydrogenases.

## Results

**Clade-dependant variations in the FAD covalent tethering residues in AA7.** To date, the molecular specificity signatures in AA7s have not been unveiled. To explore the sequence diversity in this family, we used the sequence of the previously characterised *Fusarium graminearum* chitooligosaccharide oxidase *FgChitO*¹⁸ as a query to retrieve sequences comprising 470–570 amino acids (aa) from a BLAST search. The sequences (*n* = 1927), originating from fungi, eukaryotic microorganisms and plants, were aligned and curated to generate a phylogenetic tree formed by six clades (Fig. 1a). Clade I, which is the largest (34% of all sequences), is dominated by plant and fungal sequences (Fig. 1b). Indeed, this clade contains plant non-carbohydrate active enzymes, e.g., the berberine bridge enzyme from *Eschscholzia californica* *EcBBE*³² and the monolignol oxidase from *Arabidopsis thaliana* *AtBBE*-like 15³³, as well as the oligo-galacturonide oxidase from *A. thaliana* *AtOGOX1*¹³, although none of these sequences are currently assigned into AA7.

The sequences in clades IIa, III and V are mainly from Ascomycota, whereas mostly Basidiomycota sequences populate clades IV and VI. Remarkably, the majority of the retrieved Ascomycota sequences were from genera known to harbour phytopathogens, e.g., *Fusarium*, *Magnaporthe*, *Colletotrichum*, *Bipolaris*, *Alternaria* or *Botrytis*. The enrichment of phytopathogen sequences is striking in clade II. Thus, fungus-like eukaryotic plant pathogens from Oomycota clustered in clade IIb (90.5% of sequences), whereas the remaining sequences in clade II originate mainly from Ascomycota phytopathogens. All the four previously described AA7 fungal oligosaccharide oxidases clustered in a branch of clade V (Va) (Fig. 1a). The scarce insight into AA7 is evident from the lack of characterized members from (sub)clades II, III, IV, Vb and VI.

Since bi-covalent FAD tethering has been the hallmark of all hitherto characterised AA7s, we analysed the conservation of the cysteine and histidine FAD-tethering residues across the phylogenetic tree. Strikingly, clade II harboured exclusively non-canonical sequences with substitution of the FAD-binding cysteine, histidine or both, while these residues were highly conserved in the other clades (Fig. 1c). We selected five AA7 candidates belonging to different clades (Fig. 1a, Supplementary Table 1) for recombinant expression and functional characterization.

**Identification of a strict oligosaccharide dehydrogenase in AA7.** The selected fungal sequences share 25–39% amino acid sequence identity and originate from *Aspergillus nidulans* (*AnAA7A*, clade I), *Fusarium graminearum* (*FgCelDH7C*, clade IIa and *FgChi7B*, clade Va), *Magnaporthe oryzae* (*MoChi7A*, clade Vb) and *Polyporus brumalis* (*PbChi7A*, clade VI) (Supplementary Table 1). The selected AA7 enzymes were successfully expressed in *Pichia pastoris* and highly pure enzymes were obtained from a single affinity purification. The substrate specificity of each recombinant enzyme was assessed against a panel of 40 compounds including saccharides with a degree of polymerization (DP) 1–4 as well as sugar alcohols and aromatic alcohols (Supplementary Table 2) by monitoring both their oxidase activity (i.e. H₂O₂ production using a peroxidase coupled assay) and dehydrogenase
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Phylogenetic analysis of AA7-like sequences. **a** The phylogenetic tree is based on 1927 sequences. Biochemically characterized enzymes are coloured according to clade, with green circles indicating enzymes from the present study. The PDB entries (in square brackets) are given for available enzyme structures. Clade Va harbours the canonical previously described oligosaccharide oxidases: SsGOOX¹⁵ from *Sarocladium strictum* active on cellooligosaccharides, MnLaO¹⁶ from *Microdochium nivale* active on lactose, TtXyIO¹⁷ from *Thermothelomyces thermophilus* active on xylooligosaccharides and FgChiT0¹⁸ from *Fusarium graminearum* active on chitooligosaccharides. Clade I contains characterised plant enzymes³³,³⁹ that are not assigned in AA7 including the oligogalacturonide oxidase (AtOGOX1)¹³ from *Arabidopsis thaliana*¹³,³³,³⁹. **b** Clade-wise taxonomic distribution of putative AA7 sequences. **c** Clade-wise conservation percentage of the histidine and cysteine FAD-cofactor tethering residues. The number of sequences within each clade is indicated above each bar. The accessions of the sequences in the tree are provided in the Source Data file.</figcaption></figure>

activity using 2,6-dichlorophenolindophenol (DCIP) as an artificial electron acceptor.

FgChi7B (clade Va), MoChi7A (clade Vb) and PbChi7A (clade VI) displayed the highest normalised rates (V₀/E) on chitooligosaccharides, followed by *N*-acetyl glucosamine (GlcNAc) (Supplementary Fig. 1). Importantly, we also discovered oxidase activity against substrates not reported before in AA7. Thus, MoChi7A was active on *N*-acetyl galactosamine (GalNAc), lacto-*N*-biose (LNB) and galacto-*N*-biose (GNB), while PbChi7A oxidised mannooligosaccharides (Supplementary Table 3, Supplementary Fig. 1). The Kₘ values of MoChi7A and PbChi7A

towards GlcNAc were about 20-fold lower compared to FgChi7B (Table 1), highlighting a marked difference in affinity for monosaccharides between these enzymes. No activity was detected for AnAA7A on any of the tested substrates, suggesting that AA7 targets a wider range of substrates than currently reported. The pH-dependence of activity profiles showed the highest relative activities in the pH range 7.0–7.5, except for PbChi7A, which appeared to have a broader profile (Supplementary Fig. 2).

Apart from FgCelDH7C (*vide infra*), the dehydrogenase activity for a given substrate was in the same range or higher
<caption>**Table 1** Apparent kinetic parameters of AA7-catalysed oxidation of saccharides.</caption>

**Oxidase kinetics**

<table><thead><tr><th rowspan="2">Substrate</th><th rowspan="2">Parameter</th><th colspan="3">Enzyme</th></tr><tr><th>MoChi7A</th><th>PbChi7A</th><th>FgChi7B</th></tr></thead><tbody><tr><td rowspan="3">GalNAc</td><td>k<sub>cat</sub> (s<sup>-1</sup>)</td><td>3.7 ± 0.1</td><td>ND</td><td>ND</td></tr><tr><td>K<sub>M</sub> (mM)</td><td>1.69 ± 0.12</td><td>ND</td><td>ND</td></tr><tr><td>k<sub>cat</sub>/K<sub>M</sub> (s<sup>-1</sup> M<sup>-1</sup>)</td><td>2.16 × 10<sup>3</sup></td><td></td><td></td></tr><tr><td rowspan="5">GlcNAc</td><td>k<sub>cat</sub> (s<sup>-1</sup>)</td><td>3.8 ± 0.1</td><td>3.0 ± 0.2<br>3.3 ± 0.2<sup>a</sup></td><td>6.8 ± 0.3</td></tr><tr><td>K<sub>M</sub> (mM)</td><td>0.80 ± 0.05</td><td>0.69 ± 0.20<br>0.50 ± 0.10<sup>*</sup></td><td>15.7 ± 1.13</td></tr><tr><td>k<sub>cat</sub>/K<sub>M</sub> (s<sup>-1</sup> M<sup>-1</sup>)</td><td>4.75 × 10<sup>3</sup></td><td>4.27 × 10<sup>3</sup><br>5.94 × 10<sup>3a</sup></td><td>4.30 × 10<sup>2</sup></td></tr><tr><td>k<sub>cat</sub> (s<sup>-1</sup>)</td><td>4.4 ± 0.1<br>18.3 ± 0.4<sup>a</sup></td><td>3.8 ± 0.04</td><td>2.3 ± 0.03<br>28.7 ± 0.4</td></tr><tr><td>K<sub>M</sub> (mM)</td><td>0.064 ± 0.01<br>0.33 ± 0.03<sup>a</sup></td><td>0.3 ± 0.02</td><td>0.27 ± 0.01<br>1.32 ± 0.07</td></tr><tr><td rowspan="3">Chitobiose</td><td>k<sub>cat</sub>/K<sub>M</sub> (s<sup>-1</sup> M<sup>-1</sup>)</td><td>6.83 ± 10<sup>4</sup><br>5.56 ± 10<sup>4a</sup></td><td>1.27 × 10<sup>4</sup></td><td>8.58 × 10<sup>3</sup><br>2.17 × 10<sup>4a</sup></td></tr><tr><td>k<sub>cat</sub> (s<sup>-1</sup>)</td><td>3.71 ± 0.09</td><td>3.55 ± 0.04</td><td>6.14 ± 0.29</td></tr><tr><td>K<sub>M</sub> (mM)</td><td>0.09 ± 0.01</td><td>2.0 ± 0.1<sup>a</sup><br>0.28 ± 0.01<br>0.75 ± 0.1<sup>a</sup></td><td>0.25 ± 0.03</td></tr><tr><td>Chitotriose</td><td>k<sub>cat</sub>/K<sub>M</sub> (s<sup>-1</sup> M<sup>-1</sup>)</td><td>4.22 × 10<sup>4</sup></td><td>1.25 × 10<sup>4</sup><br>2.67 × 10<sup>3a</sup></td><td>2.44 × 10<sup>4</sup></td></tr></tbody></table>

**Dehydrogenase kinetics of FgCelDH7C.**

<table><thead><tr><th></th><th>Cellobiose</th><th>Cellotriose</th><th>Cellotetraose</th><th>Cellopentaose</th></tr></thead><tbody><tr><td>k<sub>cat</sub> (s<sup>-1</sup>)</td><td>30.8 ± 0.81</td><td>36.1 ± 0.74</td><td>39.27 ± 0.76</td><td>45.28 ± 2.03</td></tr><tr><td>K<sub>M</sub> (mM)</td><td>5.35 ± 0.39</td><td>3.94 ± 0.22</td><td>6.06 ± 0.27</td><td>11.25 ± 0.94</td></tr><tr><td>k<sub>cat</sub>/K<sub>M</sub> (s<sup>-1</sup> M<sup>-1</sup>)</td><td>5757</td><td>9162</td><td>6480</td><td>4025</td></tr></tbody></table>

(up to 12-fold) than the corresponding oxidase activity, and the substrate specificity profiles were mostly similar using both assays (Table 1, Supplementary Fig. 1). By contrast, FgCelDH7C acted as a dehydrogenase with a preference for cellooligosaccharides with kcat values of 32–42 s⁻¹ and the highest catalytic efficiency on cellotriose (Table 1, Supplementary Fig. 1h). The enzyme was also active on glucose, α-(1,4)-glucooligomers (maltooligosaccharide) and lactose, all sharing a reducing end glucosyl unit. Interestingly, the oxidase activity was estimated to be several orders of magnitude lower (V₀/E ≈ 1 – 10 × 10⁻⁵ s⁻¹) (Supplementary Fig. 1g, h). This uniquely higher dehydrogenase/oxidase activity ratio compared to canonical AA7s distinguishes FgCelDH7C as an AA7 dehydrogenase with only trace oxidase activity. Interestingly, FgCelDH7C had a markedly different flavin absorbance spectrum than typical oxidase homologues (Fig. 2a), suggestive of large changes in the FAD chemical environment. The stability of the enzyme was also evaluated to verify the structural integrity (unfolding temperature Tm > 55 °C, Supplementary Fig. 2b). Both MoChi7A, which catalyses efficient oxidation of substrates not reported within AA7, and the strict dehydrogenase FgCelDH7C were selected for further in situ NMR spectroscopy to bring insight into their redox chemistry. Initial experiments identified the chemical shift assignments of substrates, intermediates and products based on two-dimensional ¹H-¹³C NMR analyses (Supplementary Fig. 3). Time-resolved two-dimensional ¹H-¹³C HSQC NMR spectra were acquired to follow the time-course conversion of GlcNAc and GalNAc by MoChi7A, showing the initial formation of a 1,5-pyrano (δ) lactone prior to rearrangement to a 1,4-furano (γ) lactone. The carboxylic acid open form was the main product of GlcNAc

oxidation (due to lactone hydrolysis), whereas the 1,4-furano (γ) lactone accumulated as the main product of GalNAc oxidation (Supplementary Fig. 4). Similarly, assignments were established for the oxidation of cellobiose by FgCelDH7C (Supplementary Fig. 5). A time series of one-dimensional ¹H NMR spectra was used for kinetic analysis of the FgCelDH7C catalysed oxidation of cellobiose in the presence of DCIP that displayed exchange line broadening during the reaction (Supplementary Fig. 6a). Line shape analysis showed that oxidised DCIP was reduced in a slow to intermediate exchange regime (exchange rate ≈ 87 s⁻¹) (Supplementary Fig. 6b). Interestingly, the sharp signals of reduced DCIP were only observed upon stoichiometric conversion of the oxidized DCIP form and hydrolysis of the lactone to the acid (Supplementary Fig. 6a), most likely due to the irreversible hydrolysis breaking the equilibrium. Notably, the oxidase reaction with O₂ as an electron acceptor occurred orders of magnitude slower in the absence of DCIP as electron acceptor (Fig. 2b, c), corroborating the activity data from the coupled peroxidase assay above. Further NMR analyses and activity assays validated the dependence of cellooligosaccharide oxidation on DCIP concentration (Supplementary Fig. 7).

**Structural elements underpinning the dehydrogenase/oxidase activity profiles.** We determined the crystal structures of FgChi7B and FgCelDH7C at resolutions of 2.4 and 1.6 Å, respectively (Supplementary Table 4, Supplementary Fig. 8). Both enzymes share the canonical AA7 fold comprising an FAD-binding F domain and a substrate-binding S domain that is formed by a central β-sheet flanked by α-helices (Supplementary Fig. 9a, b). Moreover, the two enzymes share an aromatic cluster comprising: (i) the catalytic base tyrosine (FgCelDH7C Y454), (ii) a tyrosine/phenylalanine (FgCelDH7C F99) that stacks onto and stabilises the catalytic base, and (iii) an aromatic residue that stacks mainly onto the saccharide unit penultimate to the reducing end (FgCelDH7C F383) (Fig. 3a). Interestingly, the residues of this aromatic cluster are also highly conserved in most sequences in clades IIa and Va (Fig. 3b), highlighting their importance in catalysis and substrate binding.

By contrast, the structure of FgCelDH7C revealed a unique active site architecture, as compared to all hitherto available AA7 structures. The first notable overall structural difference is that a loop (residues Q335-F346) that flanks the active site in FgCelDH7C appears shorter due to a preceding helical segment that folds away from the active site, which provides a less occluded and more solvent-exposed FAD-cofactor compared to FgChi7B and all hitherto reported AA7 structures (Supplementary Fig. 9).

Strikingly, the mode of FAD-anchoring represents a second unique characteristic distinguishing FgCelDH7C from the other AA7 structures. Thus, the isoalloxazine ring of the FAD cofactor in FgChi7B is bi-covalently anchored at the C6 atom to C162 Sγ (6-S-cysteinyl) and the 8α-methyl group to H102 Nδ¹ (8α-N1-histidyl), akin to all previously determined AA7 structures (Supplementary Fig. 9a). By contrast, FgCelDH7C is monocovalently histidyl-tethered due to the substitution of the FAD-tethering cysteine to a glycine (Fig. 3c). Notably, this cysteine is conserved in the majority of AA7 except in clade II sequences that possess a glycine, an alanine or a serine at this position (Fig. 3b, d, Supplementary Fig. 10). A third structural signature of FgCelDH7C, compared to all available AA7 structures, is the substitution of a histidine residue facing the isoalloxazine ring in the FAD-binding domain to a serine (S165) (Fig. 3c).

A fourth different motif in FgCelDH7C entails four arginine residues that protrude into the substrate-binding pocket (Fig. 3c, d), providing a distinctively positive charged milieu compared to other structurally characterised AA7s (Supplementary Fig. 11).

<footnote><sup>a</sup>Apparent kinetic parameters using the dehydrogenase assay. The kinetic parameters are determined from a global fit of the Michaelis-Menten expression to triplicate (n = 3 independent experiments), except for cellopentaose (n = 1 experiment). The data are shown as means ± standard deviations. For cellopentaose the error estimates of the fit to the single data set are shown.</footnote>
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

**Fig. 2** Spectral properties and time-course NMR analysis of the AA7 cellooligosaccharide dehydrogenase *FgCelDH7C*. **a** Spectral comparison of the flavin absorbance of *FgCelDH7C* and the oxidases *FgChi7B* as well as *MoChi7A*, all at 20 µM. The data are means ± standard deviations (n = 3 independent experiments). **b** and **c** Show the time resolved in situ ¹H NMR analyses of the oxidation of cellobiose by *FgCelDH7C* in the presence (b) or absence (c) of 1.3 mM 2,6-dichlorophenolindophenol (DCIP) as an electron acceptor, respectively. The presented spectra in panels b and c are from a single experiment (n = 1). Source data are provided as a Source Data file.

Notably, members of clade IIb share the above signatures with clade IIa, except for the inverted substitution of the FAD-tethering residues. Hence, the cysteine is conserved, but the histidine counterpart is substituted for an arginine (Fig. 1, Supplementary Fig. 10). In conclusion, the present analysis allows the identification of active-site signatures associated with the functional variation in AA7, e.g., activity on oligosaccharides and the dehydrogenase versus the typical oxidase activity.

**The cellooligosaccharide dehydrogenase FgCelDH7C primes and fuels LPMO activity.** To investigate if AA7 enzymes can trigger the activity of LPMOs, we set up AA7-LPMO mediated cellulose degradation assays to monitor the release of native (unmodified) and oxidised cellooligosaccharides from cellulose. The two cellulose-active LPMOs from *Podospora anserina* PaLPMO9E and PaLPMO9H³⁴ were used in these assays. The interplay between *FgCelDH7C* and the C1-oxidising PaLPMO9E was evident from the marked increase in native and C1-oxidised cellooligosaccharides released from phosphoric acid swollen cellulose (PASC) relative to the control reaction, which was fuelled by electrons from ascorbate (Supplementary Fig. 12a). Since *FgCelDH7C* is inactive on cellulose, the release of oligosaccharides in this assay is solely attributed to the PaLPMO9E activity. However, the generation of C1-oxidised cellooligosaccharides, by both *FgCelDH7C* and PaLPMO9E, precluded determining the contribution of each enzyme to the total amount of oxidized species. Therefore, a similar assay was also performed with the C4-oxidising PaLPMO9H (instead of PaLPMO9E) as it allowed attributing both C4 and C1/C4 oxidised cellooligosaccharides (identification based on the previous work³⁴) exclusively to the LPMO activity. This assay with PaLPMO9H resulted in a significant increase in single- and double-oxidized oligosaccharides relative to the ascorbate control (Supplementary Fig. 12b). We performed similar assays on Avicel (higher crystallinity cellulose than PASC), which revealed that the addition of the preferred *FgCelDH7C* substrate cellotetraose (DP4) results in a considerable increase in the level of released cellooligosaccharides in a dose-dependent manner (Supplementary Fig. 13). A similar fuelling of LPMO activity was also observed when the same assay was performed using the previously characterised C4-oxidising AA9 LPMO from *Lentinus similis* (LsAA9A)³⁵,³⁶ (Supplementary Fig. 14). These findings indicate that the ability of *FgCelDH7C* to fuel cellulose active LPMOs is not LPMO-specific.

We further used the CDH from *P. anserina* to benchmark the *FgCelDH7C-PaLPMO9H* system, as CDHs are recognised as key redox partners to AA9 LPMOs²⁶,²⁷,³⁴. The observed release of cello-oligomers was markedly higher when the LPMO reaction was fuelled by a 3-fold lower *FgCelDH7C* concentration compared to *PaCDHB* (Fig. 4, Supplementary Fig. 15). Indeed, the *FgCelDH7C-PaLPMO9H* pair released 1.5- and 2-fold higher C1-C4 double-oxidised and native cellooligosaccharides, respectively, as compared to the *PaCDHB*-fuelled reaction. To evaluate if this interplay was specific for *FgCelDH7C*, we used the promiscuous oxidase *MoChi7A* in a similar cellulose degradation assay. Interestingly, the high level of AA7-fuelled LPMO activity appeared specific to *FgCelDH7C* based on the low amounts of C1- and C1-C4 oxidised species observed in assays where *FgCelDH7C* was replaced with *MoChi7A* that displays low cellooligosaccharides oxidase side-activity (Supplementary Fig. 16).

**Mechanistic insights into the LPMO-AA7 interplay.** To investigate the mechanism of the AA7-LPMO interplay, we analysed the effect of H₂O₂ or the superoxide species (O₂⁻) by performing the reaction in the presence of either horseradish peroxidase (HRP) that converts H₂O₂ to H₂O²¹ or superoxide dismutase
<figure><img src="image_8.png" /><figcaption>**Fig. 3** Active site signatures of AA7 oligosaccharide oxidases versus dehydrogenases. The active sites of the clade IIa discovered dehydrogenase *FgCelDH7C* (PDB: 6YJI) and the canonical clade Va cellooligosaccharide oxidase SsGOOX (in complex with 5-amino-5-deoxy-cellobiono-1,5-lactam ABL, PDB: 2AXR) are shown. **a** The FAD-tethering histidine in addition to an aromatic cluster comprising the tyrosine base catalyst, a phenylalanine/tyrosine and the substrate-stacking aromatic residue are conserved features in fungal oligosaccharide dehydrogenases and oxidases. **b** Sequence logos of patches spanning the structurally similar active site residues shown in (a). **c** Active site differences between clades IIa and Va. **d** Sequence logos of patches spanning the structurally divergent active site residues shown in c from clades IIa and Va. The amino acid numbering of the deposited protein sequences is used in the figure.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 4** Analysis of AA7-LPMO interplay in cellulose degradation. Reactions were performed on Avicel using *FgCelDH7C* and *PaLPMO9H* with subsequent ionic chromatography (HPAEC-PAD) analysis. **a** A representative chromatogram part showing C1-C4 double-oxidised species in the cellulose degradation assay including combinations of Avicel (5 mg mL⁻¹), *FgCelDH7C* (0.4 µM), *PaLPMO9H* (4 µM), *PaCDHB* (1.2 µM), cellotetraose (DP4, 0.8 mM) and ascorbate (Asc., 1 mM) as indicated in the figure. **b** Comparison of the cellulose degradation assay based on the cumulative area under the peaks of native (DP3, DP5 and DP6), C1 oxidised (except DP2 and DP4 which were added as substrates for the CDH and AA7, respectively) and C1-C4 double-oxidised cellooligosaccharides from the reactions in (a). The data in (a) and (b) (n = 3 independent reactions) were generated in NaOAc/NaOH buffer (50 mM, pH 5.2) at 35 °C. The bar plot in (b) shows the means of total peak area (n = 3 independent reactions, each shown as a white circle) with standard deviations. Source data are provided as a Source Data file.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 5** Mechanistic insights into the activation of LPMO by FgCelDH7C. **a, b** Effect of horseradish peroxidase (HRP) and superoxide dismutase (SOD) on the interplay between PaLPMO9H (4 µM) and FgCelDH7C (0.41 µM) in Avicel degradation based on the total area of native, C1 oxidised (C1 ox) and double C1-C4 oxidised (C1-C4ox) cello-oligosaccharides as analysed by HPAEC-PAD. Controls prepared in the absence of FgCelDH7C (1) or in the absence of both SOD and FgCelDH7C (2). The assays were performed for 18 h at 35 °C and terminated using NaOH (0.1 M) prior to the HPAEC-PAD analyses. The total peak area (white circles, *n* = 3 independent experiments) are shown and the bar plots display the means with standard deviations. **c** X-band Electron Paramagnetic Resonance (EPR) spectra of PaLPMO9H–Cu(II) (100 µM) in the presence of cellotriose (DP3, 1 mM) before (blue line) and after (red line) addition of FgCelDH7C (AA7). **d** X-band EPR spectra of PaLPMO9H–Cu(II) (20 µM) in buffer (blue line), in the presence of FgCelDH7C (AA7) pre-reduced with dithionite (red line), or directly fully reduced with dithionite (10 eq., green line). All EPR solutions and experiments were performed under anaerobic conditions. Samples were in 50 mM NaOAc, pH 5.2 and spectra were recorded at 50 K with a 4 mW microwave power and a 30 Gauss modulation amplitude. The data in (c) and (d) are based on a single experiment (*n* = 1). Source data for the a and b panels are provided as a Source Data file.</figcaption></figure>

(SOD) that converts superoxide ions³⁷ to H₂O₂ and O₂, respectively. The addition of increasing concentrations (5–38 nM) of HRP, reduced the AA7-potentiated LPMO activity (Fig. 5a). These results are consistent with H₂O₂, which is supplied by FgCelDH7C and by the side-activity of primed LPMOs, being a preferred LPMO co-substrate. In the presence of increasing concentrations (0.07–0.5 µM) of SOD³⁷, we observed a similar inhibitory trend (Fig. 5b). These findings suggest that O₂⁻ plays a specific role in the FgCelDH7C-fuelled LPMO reaction (Fig. 5b). Indeed, the priming reduction of LPMOs with the superoxide ion has been previously proposed³⁸. In addition to O₂⁻-mediated priming, we also investigated the possibility of direct priming electron transfer between the AA7 and the LPMO using Electron Paramagnetic Resonance (EPR) spectroscopy. The extent of LPMO-Cu(II) reduction to Cu(I) was monitored as a decrease of the Cu(II) EPR signal due to the silence of LPMO-Cu(I) species in EPR²⁵. The addition of FgCelDH7C to PaLPMO9H pre-incubated with

cellotriose (DP3) under anaerobic conditions led to a decrease (35–40%) in PaLPMO9H-Cu(II) signal (Fig. 5c), whereas DP3 alone failed to induce a change in the Cu(II) EPR signal intensity, consistent with the involvement of FgCelDH7C in the observed priming of PaLPMO9H (Fig. 5d). The fact that only partial reduction of the Cu(II) center was observed could be due to kinetic or thermodynamic barriers. We have previously determined the redox potentials of the LPMOs PaAA9E and PaAA9H at pH 5 to +155 mV and +326 mV, respectively²⁵. The redox potentials of canonical AA7 oxidases with a bi-covalently anchored-FAD is about +130 mV, whereas mutants that abolish the cysteinyl-FAD bond exhibited redox potentials in the 60 mV range³⁹,⁴⁰. Accordingly, the reduction of LPMOs by FgCelDH7C is likely to be thermodynamically feasible assuming a redox potential in the same range as the abolished cysteinyl-FAD bond mutants. To test the possibility of a direct electron transfer between FgCelDH7C and PaLPMO9H, we used a strategy
previously developed to analyse electron transfer between redox partner proteins⁴¹. We first established that PaLPMO9H-Cu(II) could be efficiently reduced by the chemical reductant dithionite (Fig. 5d). Then, equimolar solutions of FgCelDH7C and PaLPMO9H were prepared anaerobically. The FgCelDH7C solution was pre-reduced with sub-stoichiometric (80%) amounts of dithionite to avoid excess dithionite in the medium. Identical volumes of the two protein solutions were then mixed which resulted in about 50% reduction of the LPMO active site from Cu (II) to the Cu(I) state (Fig. 5d). Altogether, our data are indicative of a molecular oxygen-independent, direct electron transfer from the AA7-bound FAD to the active site of the LPMO.

Discussion

This study sheds light on an enigmatic family of flavo-enzymes mainly occurring in fungi and fungi-like eukaryotic microorganisms as well as in plants. Our phylogenetic analysis mapped the four hitherto described AA7 fungal oligosaccharide oxidases together with a large set of mostly uncharacterised orthologues into six clades. The majority of clade I sequences share a highly conserved motif involving the catalytic tyrosine base and a stacking tyrosine (or phenylalanine) that is likely to rigidify the position of the catalytic base (Fig. 3a, Supplementary Fig. 10). Curiously, a few atypical enzymes lack the catalytic tyrosine, e.g., *EcBBE* from clade I possess a histidine and an adjacent glutamic acid at the equivalent position in the structure. This is in agreement with an alternative mechanism, involving acid/base catalysis previously proposed for *EcBBE*⁴². Nonetheless, the conservation of the catalytic motif in most sequences supports the conservation of the canonical mechanism across all AA7 clades. These findings expand significantly the sequence inventory that can be considered as AA7.

The substrate aromatic-stacking platform (Phe, Tyr or Trp) (Fig. 3a, Supplementary Fig. 10), is conserved in all previously and presently described oligosaccharides-active AA7s, in addition to a clade I enzyme that is specific for aromatic alcohols³³, indicative of shared substrate-stacking function of the above aromatic platform. By contrast, sequences that lack the aromatic-stacking residue are unlikely to be active on oligosaccharides, e.g., *EcBBE*, which features in alkaloid biosynthesis³². In conclusion, the above described aromatic-stacking platform offers a signature for the identification of putative carbohydrate-active AA7 members, which are likely to have more polar substrate-binding pockets than counterparts active on aromatic substrates. An active site arginine (*FgCelDH7C* R273) (Fig. 3c, d) has been previously proposed as a specificity signature of carbohydrate-active AA7 members¹⁸. This arginine correlates with activity on substrates that lack an *N*-acetyl C2-substituent, e.g., cello-oligosaccharides and xylo-oligosaccharides. A glutamine or smaller amino acids (Fig. 3d) at this position allow the accommodation of *N*-acetylated substrates, e.g., chito-oligosaccharides consistent with our structural analysis of *FgChi7B* and the previously characterised *FgChitO*¹⁸,¹⁹ (Supplementary Fig. 17).

Interestingly, FgCelDH7C offers structural and functional insight into the previously unknown clade II (Fig. 3b, d), which is exclusively populated with atypical sequences of enzymes with potentially mono-histidyl (clade IIa), mono-cysteinyl or non-covalently bound FAD (clade IIb). The lower (<10,000-fold) oxidase activity of FgCelDH7C compared to canonical AA7s correlated with conspicuous changes in its active site architecture. The loss of cysteinyl-FAD is likely to elicit a large decrease in the midpoint redox potential, as observed in *EcBBE*³⁹ (clade I) and *SsGOOX*⁴⁰ (clade Va) when the cysteinyl bond is abolished³⁹. In both cases, the redox potentials decreased from about 130 mV to 53–61 mV accompanied by about 20-fold decrease in oxidase

activity. The changes in the flavin absorbance spectrum due to the loss of the cysteinyl-bond in both enzymes were similar to the FgCelDH7C spectrum (Fig. 2a), which provides a possible spectroscopic feature for lower redox-potential AA7 enzymes.

Oxygen gating and activation at the *re*-side of the isoalloxazine ring in the FAD-domain are also important for the oxidase activity. The entry of O₂ to the active site in flavo-enzymes is favoured through hydrophobic tunnels rather than the solvent-filled substrate-binding pocket⁴³. These tunnels converge to a cavity sterically gated by apolar residues³³,⁴⁴, where the O₂ is positioned and activated for hydride transfer. Strikingly, we unveiled that the milieu of the FAD, particularly on the *re*-side, is considerably different in clade II as compared to canonical AA7s. An eye-catching substitution of a histidine, which is highly conserved in all other clades, to a serine, leucine or valine in clade II is observed (Fig. 6a, b, Supplementary Fig. 10, sequence patch 3). Mutation of a histidine at the equivalent position in the FAD-dependant aryl alcohol oxidase from *Pleurotus eryngii* was shown to reduce the rate of the oxidative half-reaction by an order of magnitude⁴⁵. The histidine was proposed to contribute to the positioning and polarisation of O₂⁴⁵, thereby reducing its activation-free energy barrier. The substitution of this histidine in clade II enzymes is accompanied with other marked changes in the O₂ binding cavity (Fig. 6a, b). These steric and chemical changes are likely to increase the energy barrier for hydride transfer to O₂⁴³,⁴⁴, thereby severely reducing the oxidase activity. On the other hand, the uniquely open and positively charged active site of FgCelDH7C, compared to structural homologues, may favour the accommodation of organic electron acceptors to promote the dehydrogenase oxidative half reaction. Taken together, these findings give insights into the structural elements associated with the distinctively high dehydrogenase/oxidase activity profile of FgCelDH7C and possibly AA7 homologues in clade II.

Clade I plant enzymes catalyse diverse reactions, e.g., in alkaloid synthesis³² and lignin building block synthesis demonstrated for the monolignol oxidase from *Arabidopsis*³³. By contrast, the biological roles of oligosaccharide-specific oxidases have remained an enigma. Both canonical AA7 oligosaccharide-specific oxidases and clade II dehydrogenases are ubiquitous in phytopathogenic fungi and Oomycota plant parasites, responsible for some of the most devastating plant diseases in agriculture. Recently, the activity of oligogalacturonide oxidases (clade I) from *Arabidopsis* has been reported to reduce the elicitor activity of oligogalacturonides that trigger plant immune-response¹³. Interestingly, plants overexpressing the oligogalacturonide oxidase *AtOGOX1* were more resistant to infection by the phytopathogen *Botrytis cinerea*, supporting a role in plant immunity¹³. Fungal LPMOs are another category of oxidoreductases that have been hypothesized to play a role in plant pathogenesis⁴⁶, although this is yet to be demonstrated despite the abundance of LPMO genes in most fungal phytopathogens.

Interestingly, the co-secretion of fungal LPMOs and AA7 enzymes has been observed³⁰,³¹, which inspired us to investigate the interplay between these two enzyme families. A balanced supply of ROS, e.g., H₂O₂ has been shown to be crucial for sustained LPMO activity, due to the affinity (low µM apparent K_M value)²² of LPMOs to H₂O₂ and their high sensitivity to oxidative damage²¹,⁴⁷. We have shown that FgCelDH7C potentiates LPMO-catalysed cellulose degradation in the absence of any exogenous reductant. The quenching of LPMO activity with either an alternative artificial acceptor for the AA7 (Supplementary Fig. 12) or an enzyme that scavenges H₂O₂ (Fig. 4a) is consistent with the proposed peroxygenase LPMO mechanism²¹. In addition, O₂⁻ species that are generated by flavo-enzymes in the presence of O₂⁴⁸, may contribute to LPMO priming
<figure><img src="image_11.png" /><figcaption>**Fig. 6** Oxygen-binding cavity in AA7 oxidoreductases and a schematic model for AA7-LPMO interplay. **a, b** Putative oxygen-binding cavity in typical AA7 oxidases and dehydrogenases, respectively. The substitution of conserved histidine and glycine in oxidases to serine (or other small residues) and valine in dehydrogenases, respectively, is observed. **c** Schematic model of the interplay of cellobiose dehydrogenase (CDH) versus the AA7 celooligosaccharide dehydrogenase with LPMOs during cellulose degradation. Both dehydrogenase classes oxidise celooligosaccharides to the corresponding lactones. The electrons harvested from this oxidation are stored in the FAD-cofactor and subsequently delivered directly to the LPMO in the case of AA7. The transfer of priming electrons in CDH proceeds first from the dehydrogenase domain to the cytochrome *b* haem domain in (closed form). A subsequent large conformational change to the open form is required to expose the haem domain to the LPMO active site for electron transfer. The dotted lines signify the low oxidase side-activity from both classes of dehydrogenases, which generates the H₂O₂ preferred co-substrate to fuel cellulose oxidative degradation by LPMOs. Low levels of H₂O₂ are also generated at the active site of free primed LPMOs, but this is left out from the figure for clarity.</figcaption></figure>

consistent with previous findings²¹,³⁸. To evaluate the specificity of the observed AA7-mediated activation of LPMOs, MoChi7A that possesses a low oxidase activity on celooligosaccharides was used instead of *FgCelDH7C* at equimolar concentrations in a similar assay. Under these conditions, the release of oligomers from cellulose was not observed to the same extent, indicating that the priming and/or the balanced fuelling of LPMO activity is more readily achieved with the AA7 dehydrogenase (Supplementary Fig. 15). Remarkably, the performance of *FgCelDH7C* was comparable with CDH that has been shown to mediate direct transfer of priming electrons from celooligosaccharide oxidation to LPMOs via the cytochrome *b* domain⁴⁹. Our EPR experiments were consistent with direct electron transfer to the LPMO when the FAD-cofactor in *FgCelDH7C* was reduced either with dithionite or cellotriose in the absence of O₂ (Fig. 5c, d). This electron transfer and the interaction with the LPMO is likely to occur at the substrate-binding face (*si*-side) of the FAD, which is consistent with the rather open active site topology of *FgCelDH7C* (Supplementary Fig. 9). By comparison, the electron

transfer from the larger (about 800 aa) and bi-modular CDH to the LPMO requires inter-domain electron transfer from the dehydrogenase to the cytochrome domain and a subsequent large domain movement to the open CDH form⁴⁹,⁵⁰. The deletion of the cytochrome domain impairs the priming of the LPMO, but the LPMO fuelling remains efficient²⁹, likely via the oxidase side-activity of CDH. Similarly, the deletion of the cytochrome domain in the pyrroloquinoline-quinone-dependent pyranose dehydrogenase impaired the fuelling of LPMO activity²⁴. In summary, the one-module AA7 dehydrogenase confers direct priming of LPMOs in contrast to other known LPMO redox partners that require either small organic molecules or an additional cytochrome haem domain to mediate LPMO priming.

The expansion of the substrate range of fungal oligosaccharide-oxidising flavo-enzymes belonging to the AA7 family and the potentiation of LPMO activity, highlight the functional diversity within AA7. These findings set the stage for harnessing other AA7 dehydrogenases to fine-tune the activity of LPMOs that lack known redox partners. The demonstration of a functional and
efficient AA7-LPMO system offers attractive avenues for the
development of commercial enzyme blends targeting recalcitrant
biomass.

**Methods**

**Bioinformatics and phylogenetic analysis.** The sequences of 2200 putative AA7 sequences were retrieved by a BlastP search against the non-redundant protein database interfaced the NCBI server (https://blast.ncbi.nlm.nih.gov/Blast.cgi) using the *FgChitO* sequence (Genbank accession: XP_011325372) as a query and default search settings. The retrieved sequences were filtered by excluding non-eukaryotic proteins and those displaying <25% amino acid (aa) identity and <90% coverage to the query. Two additional AA7-like sequences, which were differentially upregulated in the secretome of *Aspergillus nidulans* grown on starch³¹, and a longer sequence possessing an N-terminal extension of about 250 amino acid residues (aa) (*MoChi7A*) were also included in the sequence inventory. Only sequences comprising 470–570 aa (n = 1927) were aligned using *MAFFT*⁵¹ with default parameters, and curated using Gblock.⁹¹b⁵². The alignment was used to construct a phylogenetic tree using *NGphylogeny*,⁵³,⁵⁴ and rendered by iTOL⁵⁴. The visualization of amino acid conservation was made using WebLogo⁵⁵.

**Production of selected AA7 enzymes.** Based on the phylogenetic analysis above, five selected AA7 gene fragments encoding mature peptides without signal peptides (Supplementary Table 1) were codon optimised for *P. pastoris* and purchased from GENEWIZ (NJ, USA). The gene fragments were cloned within the XbaI and XhoI restriction sites of the pPICZαA vector (Invitrogen, Carlsbad, CA, USA) in frame with the *Saccharomyces cerevisiae* α-mating factor secretion signal and fused to a C-terminal (His)₆ tag. The synthetic gene-bourselling plasmids were propagated in *Escherichia coli* DH5α, linearised with PmE1 and thereafter transformed into competent *P. pastoris* X33 by electroporation following the protocols from the Easy Select Expression System (Invitrogen). Six transformants per construct were screened for production in deep well plates for three days with MeOH addition 3% (v/v) every 24 h. The best-secreting transformants based on SDS-PAGE gel electrophoresis, from the constructs of the selected enzymes (Supplementary Table 1) were chosen for larger-scale production. These clones were grown in shake flasks containing 2 L of BMGY (Pichia expression manual, Invitrogen) containing 1 mL L⁻¹ of Pichia PTM4 trace element solution (2 g L⁻¹ CuSO₄·5H₂O, 3 g L⁻¹, MnSO₄·H₂O, 0.2 g L⁻¹ Na₂MoO₄·2H₂O, 0.02 g L⁻¹ H₃BO₃, 0.5 g L⁻¹ CaSO₄·2H₂O, 0.5 g L⁻¹ CoCl₂, 12.5 g L⁻¹ ZnSO₄·7H₂O, 22 g L⁻¹, FeSO₄·7H₂O, H₂SO₄ 1 mL L⁻¹) and biotin 0.2 g L⁻¹ to an OD₆₀₀ of 2–6, for 16 h at 30°C using an orbital shaker (200 rpm). The cells from each construct were harvested by centrifugation (4000 × g, 10 min, 4°C) and expression was induced by re-suspending the cells into 400 mL of BMMY medium at 20 °C and the culture was continued for 3 days with methanol supplementation to 3% (v/v) every 24 h. The cells were harvested by centrifugation (5,000 × g, 10 min, 4°C) and the pH of the culture supernatants was adjusted to 7.8 using NaOH (2 M), followed by sterile filtration using 0.22 µm filters (Millipore, Burlington, MA, USA). The filtered supernatants were loaded onto 5 mL His Trap HP columns (GE Healthcare, Uppsala, Sweden) connected to an Äkta purifier 100 (GE Healthcare) at 5 mL min⁻¹, equilibrated with buffer A (Tris-HCl 50 mM pH 7.8, NaCl 150 mM, imidazole 10 mM). Non-bound proteins were washed by 10 column volumes (CVs). The bound proteins were eluted with a 34% buffer B (Tris-HCl 50 mM pH 7.8, NaCl 150 mM, imidazole 500 mM) gradient in 10 CVs. Fractions containing recombinant enzymes were pooled, concentrated and buffer exchanged against 50 mM NaOAc, pH 5.2 using a 10 kDa Vivaspin ultra-filtration unit (Sartorius, Göttingen, Germany). Protein purity was assessed using SDS-PAGE analysis and enzyme concentration was determined by measuring A₂₅₀ using a Nanodrop ND-2000 (Thermo Fisher Scientific, Waltham, MA, USA) and the theoretically calculated molar extinction coefficients (Table S1) using the Protparam tool on the EXPASY server (web.expasy.org/proparam/).

**Activity assays and kinetic analyses**

**Initial screening assay.** The oxidase activity assay of the expressed MoChi7A, FgChi7B, FgCelDH7C and AnAA7A was conducted by coupling H₂O₂ production to HRP mediated oxidation of 4-aminoantipyrine (AAP) and 3,5-dichloro-2-hydroxybenzensulfonic acid (DCHBS) (both from Sigma Aldrich, St. Louis, MI, USA) in 96-well microtiter plates. Reactions (200 µL, 0.1 mM AAP, 1 mM DCHBS, 8 U mL⁻¹ HRP (Sigma Aldrich) in 50 mM NaOAc buffer, pH 5.2) were monitored by measurements of A₅₀₅ nm at 35 °C for 20 min with 30 s intervals. The assay was performed in triplicates using 0.1 µM AA7 against 6 carbohydrates, two sugar alcohols (all 2 mM except 1 mM for LNB/GNB) and 6 aromatic alcohols (1 mM) (Supplementary Table 2) by the addition of 100 µL substrate into pre-tempered reaction mixtures in microtiter plates. The substrate screening of PbChi7A was conducted using a modified reaction set up (250 µL, 0.1 mM AAP, 1 mM phenol, 8 U mL⁻¹ HRP (Sigma Aldrich) in 50 mM Tris-HCl pH 8.0) incubated in triplicates at 37 °C for 20 min and measured sequentially at A₅₀₅ nm with 30 s intervals.

**Oxidation kinetics of AA7.** The steady-state apparent kinetic parameters of MoChi7A and FgChi7B (0.65 µM) were determined in an analogous assay, where the reaction mixtures (200 µL) were monitored for 3 min with 20 s measurements intervals. The initial rates for the oxidation of N-acetylglucosamine (GlcNAc) and N-acetylgalactosamine (both in 0.09–12.5 mM) (GalNAc) (only for MoChi7A), chitobiose (0.0–0.5 mM) and chitotriose (0.03–0.6 mM) were determined in triplicates. The kinetic parameters of PbChi7A were determined using the same HRP assay on 12 substrates, including mono- (1–1500 mM), di- (1–250 mM), tri-saccharides (0.05–25 mM) (Supplementary Table 2). The kinetic parameters were determined by fitting the Michaelis-Menten equation in Origin, versions 9.55 and 18 (Northampton, Massachusetts, USA) to the initial rates, determined from the slopes of the linear reaction phases (Table 1). Due to the low oxidase activity of FgCelDH7C, the activity assay was performed at higher enzyme concentrations (150 and 500 µM) in duplicates each and the absorbance was measured for one hour with 1 min intervals to determine the normalised reaction rate (V₀/E).

**Dehydrogenase assay of AA7 activity.** The dehydrogenase activity of the enzymes was also assayed using the redox mediator 2,6-dichlorophenolindophenol (DCIP) (ε₅₂₀ = 6.8 mM⁻¹ cm⁻¹) as a terminal electron acceptor (as opposed to O₂ for the oxidase activity). Reactions (200 µL, 0.1 mM DCIP, 5 mM carbohydrate or 1 mM alcohol substrates, 0.65 µM enzyme in 50 mM HEPES, pH 7.0) were initiated by adding 100 µL carbohydrates or alcohol substrates to 100 µL pre-tempered mixtures containing all other components. The decrease in A₅₃₀ was measured continuously for 20 min with 30 s intervals at 35 °C.

**Temperature and pH stability/optimum measurements.** The pH-dependence of the activity of MoChi7A, FgChi7B and FgCelDH7C (100 nM for all) towards 2 mM of GlcNAc, chitobiose and cellobiose, respectively, was analysed in triplicates using the dehydrogenase assay as above (DCIP as a terminal electron acceptor). The assay was carried out at eight pH values in the range 5.1–9.0 by exchanging the buffer in the screening assay to the universal Britton-Robinson buffer (acetic acid, phosphoric acid, and boric acid, 50 mM of each). The dehydrogenase assay was used to allow for comparison of all three enzymes in the same assay and the initial rates of the oxidation at each pH value we determined from the linear parts of the progress curves as above. Differential scanning calorimetry was employed to examine the conformational stability of FgCelDH7C (1 mg mL⁻¹) at pH values of 5.2 and 7.2 in 50 mM (in acetate and phosphate buffers, respectively) in the temperature range of 15–90 °C at a scanning rate of 1 °C min⁻¹ using a NanoDSC instrument (TA instruments, New Castle, DE, USA).

The temperature and pH optima of PbChi7A (15 nM) were determined using the HRP assay described above against chitobiose (2 mM). For determination of the pH optimum, a Britton Robinson buffer (100 mM each) was adjusted to pH 4–11 and the oxidation activity was measured after 1 h incubation at 4 °C. The temperature stability of PbChi7A (15 nM) was measured by oxidation of 2 mM chitobiose within 1 h at different temperatures (25 °C–55 °C).

**NMR Spectroscopy of FgCelDH7C, FgChi7B and MoChi7A.** Initially, the analysis of post oxidation reactions was conducted on a reaction mixture containing FgChi7B or MoChi7A (1 µM of each), 25 mM GalNAc or GlcNAc as substrates, buffered with 50 mM NaOAc buffer pH 5.2 (500 µL, 25 °C, 30 min). The reaction mixtures were ice-cooled and immediately analysed at 15 °C. Spectra were collected using an 800 MHz Bruker Avance II spectrometer equipped with an Oxford 18.7 T magnet and a TCI cryoprobe by transferring 500 µL of the post-reaction mixtures to NMR sample tubes (5 mm) and performing ¹H, ¹H-¹H TOCSY, ¹H-¹H COSY, ¹H-¹³C HMBC NMR and multiplicity-edited ¹H-¹³C HSQC NMR spectroscopy analyses. ¹H-¹H TOCSY spectra were acquired as data matrices of 1024 × 256 complex data points sampling 128 and 64 ms in the direct and indirect dimensions, respectively. ¹H-¹H COSY spectra were acquired as data matrices of 2048 × 128 complex data points sampling 320 and 20 ms in the direct and indirect dimensions, respectively, while ¹H-¹³C HMBC NMR spectra were acquired as data matrices of 2048 × 256 complex data points sampling 212 and 5.8 ms in the direct and indirect dimensions, respectively. Multiplicity-edited ¹H-¹³C HSQC spectra were collected by sampling the NMR spectra for 106 and 15.9 ms in the direct (¹H) and indirect (¹³C) dimensions, respectively. Once intermediate and product identification was established (Supplementary Fig. 3), a series of ¹H-¹³C HSQC spectra was acquired at 15 °C to follow the conversion of GlcNAc and GalNAc by MoChi7A (1 µM) in 50 mM NaOAc buffer pH 5.2 after 5, 15, 60 and 180 min experiment time. The ¹H-¹³C HSQC spectra were collected by sampling the NMR spectra for 159 and 8.3 ms in the direct (¹H) and indirect (¹³C) dimensions, respectively, with experiments of 19 min duration to yield the time course shown in Supplementary Fig. 4 upon integration of the C2H2 NMR cross peaks.

For cellobiose, a time series of ¹H-¹³C HSQC spectra for the enzymatic conversion was acquired using FgCelDH7C (0.56 µM) to yield assignments of cellobionolactone and of cellobionic acid signals as displayed in Supplementary Figure 5. Time-resolved ¹H NMR spectra sampling 16384 complex data points of the FID for 1.27 s were subsequently used to track the reaction kinetics for cellobiose (2 mM) conversion by FgCelDH7C (0.56 µM) at 25 °C in the presence or in the absence of 1.3 mM DCIP as an electron acceptor using the same equipment as above. The ¹H singlet signal for DCIP as well as signals for H5 of cellobionolactone and H2 of cellobionic acid were integrated alongside the H2 and
H3 signals for the reducing end of cellobiose. Signal areas were normalized to the number of contributing hydrogen atoms. The same buffer, enzyme concentration and same NMR setup were subsequently used to track the dependence of cellotetraose oxidation (6 mM cellotetraose) by FgCelDH7C in the presence of 0.025, 0.05, 0.1, 0.2, 0.5 or 1 mM DCIP. All NMR spectra were acquired, processed and integrated in Topspin 3.5 pl6 with ample zero filling in all dimensions.

### Crystallization and structure Determination of FgCelDH7C and FgChi7B.

Crystallization screens PACT++ (Jena Bioscience, Jena, Germany), JCSG+ and Morpheus screens (Molecular Dimensions, Sheffield, UK) and optimisation plates of FgCelDH7C (9.0 mg mL⁻¹, deglycosylated with EndoH, New England Biolab, Ipswich, MA, USA) and FgChi7B (10.5 mg mL⁻¹), both in 20 mM NaOAc pH 5.5, were performed with a Crystal Gryphon liquid handling robot (Art Robbins Instruments, Sunnyvale, USA). Drops containing equal volumes of protein and reservoir solutions were mixed. Crystals of FgCelDH7C appeared in the F7 condition (NaSCN, 20% (w/v) PEG 3350), which was optimised by varying PEG content (12.5–25% w/v). The best crystals formed within 4 days at 20 °C, using 150 nL of reservoir and protein solutions (1:1 ratio) in 50 µl reservoir solution. These crystals were cryo-protected with PEG400 in reservoir solution and cryo-cooled in a nylon loop. Crystals of FgChi7B were initially observed in PACT++ condition E10 (PEG 3350, 22% w/v, 200 mM sodium potassium phosphate using a 1:1 protein:reservoir ratio. Crystals were optimised by micro-seeding using the above condition at 22 and 24% PEG (w/v). Drops contained 1:1:0.2 of protein:reservoir: seed stock (100 fold diluted seed stock) ratio, respectively. The crystals were cryo-cooled using a mixture of reservoir and PEG400 (0.7:0.3 ratio) added to the drop prior freezing in liquid nitrogen. Diffraction data for FgCelDH7C and FgChi7B crystals were collected in BioMax beamline at MaxIV (Lund, Sweden) using MXCuBE v3 and id30A3 (MASSIF-3) beamline at ESRF (Grenoble, France) using to 1.64 Å and 2.38 Å resolution, respectively. Both structures were solved using molecular replacement in Phaser⁵⁶ with model coordinates of *TrXylO* and initial automated model building and refinement using PHENIX.autobuild⁵⁷ and PHE-NIX refine⁵⁸, respectively. Manual model rebuilding and map inspection were performed in Coot⁵⁹ and analysed using MolProbity⁶⁰. The final validated structural models of FgCelDH7C and FgChi7B were deposited to the Protein Data Bank under the entry codes 6YJI and 6YJO, respectively (Data collection and refinement statistics are in Supplementary Table 4).

### Cellulose degradation assays based on the interplay between LPMOs and AA7 enzymes.
The assays were performed using suspensions of either 1% PASC⁶¹ and/or 0.5% Avicel (Honeywell Fluka, Morris Plains, NJ, U.S.A) (both w/v) in 250 or 500 µL containing 4.4 µM of PaLPMO9E²⁰, PaLPMO9H²⁰ or LsAA9A³⁵ as well as 0.4 µM FgCelDH7C (0.2 µM in case of LsAA9A) in 50 mM NaOAc, pH 5.2 in the presence of cellotetraose (DP4, 0.8 mM), unless stated otherwise. For the comparison of FgCelDH7C and MoChi7A, 0.2 µM of each was used. Controls were also performed in absence of individual reaction components or in the presence of 0.6 mM DCIP. Moreover, the AA7-mediated LPMO activity was compared to reactions in the presence of ascorbate (1 mM), L-Cys (1 mM), reactions containing horseradish peroxidase (HRP, 0.2–38 nM in the presence of 0.1 mM AAP, 1 mM DCHBS), superoxide dismutase (SOD, recombinantly expressed in *E. coli*, 0.07–0.5 µM, Sigma Aldrich) in 100 µM phosphate, pH 7.0, or CDH from *Podospora anserina* (*PaCDH*, 1.2 µM)³⁴.

The reactions were performed in 2 mL Eppendorf tubes with 850 rpm stirring in thermomixer (Eppendorf, Hamburg, Germany) at 35 °C for 2–24 h and thereafter quenched by boiling (100 °C, 10 min) and/or fourfold dilution in NaOH to 0.1 M followed by centrifugation (15,000 × g, 15 min, 4 °C) to separate the released saccharides in the supernatants from the insoluble fraction prior to analysis. The solubilised oxidised and non-oxidised oligosaccharides, released by the enzymatic treatments, were analysed by high-performance anion-exchange chromatography coupled with pulsed amperometric detection (HPAEC-PAD) (ThermoFischer Scientific) using a CarboPac™ PA1 (ThermoFischer Scientific) and a CarboPac™ PA1 guard column (2 × 50 mm) at 0.25 mL min⁻¹ and an eluent system as described by Westereng et al.⁶². The Chromleon 6.80 SR10 software was used to interface the instrument and analyse the chromatograms. Non-oxidised oligosaccharides were used as standards (Megazyme, Wicklow, Ireland), whereas the corresponding C1-oxidised standards (from DP2-DP6) were provided from native cello-oligosaccharides using PaCDHB treatment³⁴.

### EPR spectroscopy.
EPR experiments were performed on a Bruker ELEXSYS E500 X-band spectrometer equipped with an ER4102ST standard rectangular Bruker EPR cavity fitted to an Oxford Instruments ESR 900 helium flow cryostat. All EPR studies were carried out at 50 K with 4 mW microwave power at 9.479 GHz and 3 mT modulation amplitude at 100 kHz. The presented spectra were averaged from the accumulation of 4 scans to improve the signal-to-noise ratio. The EPR samples were prepared under anaerobic conditions in glove box (Jacomex) and frozen prior to the EPR analysis. In a first experiment, two 160 µL protein samples were prepared in 50 mM NaOAc pH 5.2 and incubated at 25 °C for 5 min: (i) control mixtures of PaLPMO9H (100 µM) and cellotriose (DP3, 1 mM), (ii) same as (i), but supplemented with FgCelDH7C (10 µM).

In the second set of experiments, PaLPMO9H (20 µM) was prepared in the same buffer in the absence or presence of 40 µM dithionite and incubated at 25 °C for 15 min to achieve the full reduction of the catalytic Cu(II) center of PaLPMO9H into Cu(I) by dithionite. Next, the following solutions were prepared and incubated anaerobically for 5 min at 25 °C: (i) FgCelDH7C (20 µM) with sub-stoichiometric amount of sodium dithionite (16 µM) to reduce 80% of the flavo-enzyme, while avoiding excess dithionite in the solution, (ii) PaLPMO9H (20 µM) without additives in the same buffer. Equal volumes of the two solutions were then pooled into a 160 µL EPR sample and incubated for few minutes prior freezing and EPR analysis.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

### Data availability

The atomic coordinates of FgCelDH7C and FgChi7B have been deposited in the Protein Data Bank (https://www.rcsb.org) under the PDB accessions 6YJI and 6YJO, respectively (see also Supplementary Table 4). The GenPept accession IDs of the enzymes characterised in the study are given in Supplementary Table 1. All the data are available from the corresponding authors upon request. Source data are provided with this paper.

Received: 26 October 2020; Accepted: 9 March 2021;

Published online: 09 April 2021

### References

1. Eriksson, K. E., Pettersson, B. & Westermark, U. Oxidation—important enzyme reaction in fungal degradation of cellulose. *FEBS Lett.* **49**, 282–285 (1974).

2. Lo Leggio, L. et al. Structure and boosting activity of a starch-degrading lytic polysaccharide monooxygenase. *Nat. Commun.* **6**, 5961 (2015).

3. Vu, V. V., Beeson, W. T., Span, E. A., Farquhar, E. R., and Marletta, M. A. A family of starch-active polysaccharide monooxygenases. Proceedings of the National Academy of Sciences 111, 13822–13827 https://doi.org/10.1073/pnas.1408090111 (2014).

4. Haddad Momeni, M. et al. Loss of AA13 LPMOs impairs degradation of resistant starch and reduces the growth of *Aspergillus nidulans*. *Biotechnol. Biofuels* **13**, 135 (2020).

5. Vaaja-Kolstad, G. et al. An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. *Science* **330**, 219–222 (2010).

6. Quinlan, R. J. et al. Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. *Proc. Natl Acad. Sci. USA* **108**, 15079–15084 (2011).

7. Beeson, W. T., Phillips, C. M., Cate, J. H. D. & Marletta, M. A. Oxidative cleavage of cellulose by fungal copper-dependent polysaccharide monooxygenases. *J. Am. Chem. Soc.* **134**, 890–892 (2012).

8. Filatriault-Chastel, C. et al. AA16, a new lytic polysaccharide monooxygenase family identified in fungal secretomes. *Biotechnol. Biofuels* **12**, 55 (2019).

9. Nekiunaite, L. et al. FgLPMO9A from *Fusarium graminearum* cleaves xyloglucan independently of the backbone substitution pattern. *FEBS Lett.* **590**, 3346–3356 (2016).

10. Couturier, M. et al. Lytic xylan oxidases from wood-decay fungi unlock biomass degradation. *Nat. Chem. Biol.* **14**, 306–310 (2018).

11. Bissaro, B., Várnaí, A., Røhr, Å. K. & Eijsink, V. G. H. Oxidoreductases and reactive oxygen species in conversion of lignocellulosic biomass. *Microbiol. Mol. Biol. Rev.* **82**, e00029–00018 (2018).

12. Lombard, V., Ramulu, H. G., Drula, E., Coutinho, P. M. & Henrissat, B. The carbohydrate-active enzymes database (CAZy) in 2013. *Nucleic Acids Res.* **42**, D490–D495 (2014).

13. Benedetti, M. et al. Four *Arabidopsis* berberine bridge enzyme-like proteins are specific oxidases that inactivate the elicitor-active oligogalacturonides. *Plant J.* **94**, 260–273 (2018).

14. Locci, F. et al. An *Arabidopsis* berberine bridge enzyme-like protein specifically oxidizes cellulose oligomers and plays a role in immunity. *Plant J.* **98**, 540–554 (2019).

15. Lee, M. H. et al. Structural characterization of glucooligosaccharide oxidase from *Acremonium strictum*. *Appl Environ. Microbiol.* **71**, 8881–8887 (2005).

16. Xu, F. et al. A novel carbohydrate:acceptor oxidoreductase from *Microdochium nivale*. *Eur. J. Biochem.* **268**, 1136–1142 (2001).

17. Ferrari, A. R. et al. Discovery of a xylooligosaccharide oxidase from *Myceliophthora thermophila* C1. *J. Biol. Chem.* **291**, 23709–23718 (2016).

18. Heuts, D. P., Janssen, D. B. & Fraaije, M. W. Changing the substrate specificity of a chitooligosaccharide oxidase from *Fusarium graminearum* by model-inspired site-directed mutagenesis. *FEBS Lett.* **581**, 4905–4909 (2007).
19. Savino, S., Jensen, S., Terwisscha van Scheltinga, A. & Fraaije, M. W. Analysis of the structure and substrate scope of chitooligosaccharide oxidase reveals high affinity for C2-modified glucosamines. *FEBS Lett.* **594**, 2819–2828 (2020).

20. Leferink, N. G., Heuts, D. P., Fraaije, M. W. & van Berkel, W. J. The growing VAO flavoprotein family. *Arch. Biochem. Biophys.* **474**, 292–301 (2008).

21. Bissaro, B. et al. Oxidative cleavage of polysaccharides by monocopper enzymes depends on H₂O₂. *Nat. Chem. Biol.* **13**, 1123–1128 (2017).

22. Kuusk, S. et al. Kinetics of H₂O₂-driven degradation of chitin by a bacterial lytic polysaccharide monooxygenase. *J. Biol. Chem.* **293**, 523–531 (2018).

23. Jones, S. M., Transue, W. J., Meier, K. K., Kelemen, B. & Solomon, E. I. Kinetic analysis of amino acid radicals formed in H₂O₂-driven Cu¹ LPMO reoxidation implicates dominant homolytic reactivity. *Proc. Natl Acad. Sci.* **20192499** (2020).

24. Varnai, A., Umezawa, K., Yoshida, M. & Eijsink, V. G. H. The pyrroloquinoline-quinone-dependent pyranose dehydrogenase from *Coprinopsis cinerea* drives lytic polysaccharide monooxygenase action. *Appl. Environ. Microbiol* **84**, e00156–00118 (2018).

25. Garajova, S. et al. Single-domain flavoenzymes trigger lytic polysaccharide monooxygenases for oxidative degradation of cellulose. *Sci. Rep.* **6**, 28276 (2016).

26. Phillips, C. M., Beeson, W. T., Cate, J. H. & Marletta, M. A. Cellobiose dehydrogenase and a copper-dependent polysaccharide monooxygenase potentiate cellulose degradation by *Neurospora crassa*. *ACS Chem. Biol.* **6**, 1399–1406 (2011).

27. Kracher, D. et al. Extracellular electron transfer systems fuel cellulose oxidative degradation. *Science* **352**, 1098–1101 (2016).

28. Miyauchi, S. et al. Conserved white-rot enzymatic mechanism for wood decay in the Basidiomycota genus *Pycnoporus*. *DNA Res.* **27**, dsaa011 (2020).

29. Kracher, D. et al. Polysaccharide oxidation by lytic polysaccharide monooxygenase is enhanced by engineered cellobiose dehydrogenase. *FEBS J.* **287**, 897–908 (2020).

30. Berrin, J. G., Rosso, M. N. & Abou Hachem, M. Fungal secretomics to probe the biological functions of lytic polysaccharide monooxygenases. *Carbohydr. Res.* **448**, 155–160 (2017).

31. Nekiunaite, L., Arntzen, M. Ø., Svensson, B., Vaaje-Kolstad, G. & Abou Hachem, M. Lytic polysaccharide monooxygenases and other oxidative enzymes are abundantly secreted by *Aspergillus nidulans* grown on different starches. *Biotechnol. Biofuels* **9**, 187 (2016).

32. Winkler, A. et al. Structural and mechanistic studies reveal the functional role of bivalent flavinylation in berberine bridge enzyme. *J. Biol. Chem.* **284**, 19993–20001 (2009).

33. Daniel, B. et al. Oxidation of monolignols by members of the berberine bridge enzyme family suggests a role in plant cell wall metabolism. *J. Biol. Chem.* **290**, 18770–18781 (2015).

34. Bennati-Granier, C. et al. Substrate specificity and regioselectivity of fungal AA9 lytic polysaccharide monooxygenases secreted by *Podospora anserina*. *Biotechnol. Biofuels* **8**, 90 (2015).

35. Frandsen, K. E. H. et al. The molecular basis of polysaccharide cleavage by lytic polysaccharide monooxygenases. *Nat. Chem. Biol.* **12**, 298–303 (2016).

36. Frandsen, K. E. H. et al. Identification of the molecular determinants driving the substrate specificity of fungal lytic polysaccharide monooxygenases (LPMOs). *J. Biol. Chem.* **296**, 100086 (2021).

37. Kracher, D., Andlar, M., Furtmuller, P. G. & Ludwig, R. Active-site copper reduction promotes substrate binding of fungal lytic polysaccharide monooxygenase and reduces stability. *J. Biol. Chem.* **293**, 1676–1687 (2018).

38. Bissaro, B., Kommedal, E., Rohr, Å. K. & Eijsink, V. G. H. Controlled depolymerization of cellulose by light-driven lytic polysaccharide oxygenases. *Nat. Commun.* **11**, 890 (2020).

39. Winkler, A., Kutchan, T. M. & Macheroux, P. 6-S-cysteinylation of bio-covalently attached FAD in berberine bridge enzyme tunes the redox potential for optimal activity. *J. Biol. Chem.* **282**, 24437–24443 (2007).

40. Huang, C.-H. et al. Functional roles of the 6-S-cysteinyl, 8alpha-N1-histidyl FAD in glucooligosaccharide oxidase from *Acremonium strictum*. *J. Biol. Chem.* **283**, 30990–30996 (2008).

41. Dermoun, Z. et al. The NADP-reducing hydrogenase from *Desulfovibrio fructosovorans*: functional interaction between the C-terminal region of HndA and the N-terminal region of HndD subunits. *Biochim. Biophys. Acta* **1556**, 217–225 (2002).

42. Kutchan, T. M. & Dittrich, H. Characterization and mechanism of the berberine bridge enzyme, a covalently flavinylated oxidase of benzophenanthridine alkaloid biosynthesis in plants. *J. Biol. Chem.* **270**, 24475–24481 (1995).

43. Romero, E., Gomez Castellanos, J. R., Gadda, G., Fraaije, M. W. & Mattevi, A. Same substrate, many reactions: oxygen activation in flavoenzymes. *Chem. Rev.* **118**, 1742–1769 (2018).

44. Zafred, D. et al. Rationally engineered flavin-dependent oxidase reveals steric control of dioxygen reduction. *FEBS J.* **282**, 3060–3074 (2015).

45. Hernández-Ortega, A. et al. Role of active site histidines in the two half-reactions of the aryl-alcohol oxidase catalytic cycle. *Biochemistry* **51**, 6595–6608 (2012).

46. Johansen, K. S. Lytic polysaccharide monooxygenases: the microbial power tool for lignocellulose degradation. *Trends Plant Sci.* **21**, 926–936 (2016).

47. Loose, J. S. M. et al. Multipoint precision binding of substrate protects lytic polysaccharide monooxygenases from self-destructive off-pathway processes. *Biochemistry* **57**, 4114–4124 (2018).

48. Massey, V. Activation of molecular oxygen by flavins and flavoproteins. *J. Biol. Chem.* **269**, 22459–22462 (1994).

49. Tan, T.-C. et al. Structural basis for cellobiose dehydrogenase action during oxidative cellulose degradation. *Nat. Commun.* **6**, 7542 (2015).

50. Breslmayr, E. et al. Protein conformational change is essential for reductive activation of lytic polysaccharide monooxygenase by cellobiose dehydrogenase. *ACS Catal.* **10**, 4842–4853 (2020).

51. Katoh, K., Rozewicki, J. & Yamada, K. D. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. *Brief. Bioinform.* **20**, 1160–1166 (2019).

52. Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. *Mol. Biol. Evol.* **17**, 540–552 (2000).

53. Lemoine, F. et al. NGPhylogeny.fr: new generation phylogenetic services for non-specialists. *Nucleic Acids Res.* **47**, W260–W265 (2019).

54. Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. *Nucleic Acids Res.* **47**, W256–W259 (2019).

55. Crooks, G. E., Hon, G., Chandonia, J. M. & Brenner, S. E. WebLogo: a sequence logo generator. *Genome Res.* **14**, 1188–1190 (2004).

56. Bunkoczi, G. et al. MRage: automated molecular replacement. *Acta Crystallogr. Sect. D* **69**, 2276–2286 (2013).

57. Terwilliger, T. C. et al. Iterative model building, structure refinement and density modification with the PHENIX AutoBuild wizard. *Acta Crystallogr. Sect. D* **64**, 61–69 (2008).

58. Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. *Acta Crystallogr. Sect. D. Biol. Crystallogr.* **66**, 213–221 (2010).

59. Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. Features and development of Coot. *Acta Crystallogr. Sect. D. Biol. Crystallogr.* **66**, 486–501 (2010).

60. Williams, C. J. et al. MolProbity: more and better reference data for improved all-atom structure validation. *Protein Sci.* **27**, 293–315 (2018).

61. Wood, T. M. In *Methods in Enzymology* (ed^^(eds)). (Academic Press, 1988).

62. Westereng, B. et al. Efficient separation of oxidized cello-oligosaccharides generated by cellulose degrading lytic polysaccharide monooxygenases. *J. Chromatogr. A* **1271**, 144–152 (2013).

## Acknowledgements

This study was supported by a Novo Nordisk Foundation Post-doc grant within the Biotechnology Based Synthesis and Production Research Programme to MHM (NNF17OC0025642). F.F. and D.H.W. were funded by The Novo Nordisk Foundation Grant No. NNF10CC1016517. NMR spectra were acquired at the NMR center DTU supported by the Villum and Carlsberg Foundations. J.G.B., B.G. (Bruno G) and F.B. (Frédéric B) thank IM2B (Institut de Microbiologie, Bioénergies et Biotechnologie). The DSC instrument is funded by the Carlsberg Foundation Grant (2013-01-0112) to M.A.H. The authors are also grateful to the EPR facilities available at the French EPR network (IR CNRS 3443) and the Aix-Marseille University EPR center. Prof. Peter Westh and Dr. Jeppi Kari are acknowledged for discussion on the redox-activities of AA7.

## Author contributions

M.H.M., J.G.B. and M.A.H. conceived the research and designed the experiments.

M.H.M. and F.F. crystallised the enzymes, F.F. collected the crystallographic data and generated the first models, M.H.M. refined the structures with F.F. M.H.M. performed the sequence and phylogenetic analysis with help from V.L. Recombinant enzyme production was performed by M.H.M. and M.H. M.H.M. generated the initial biochemical characterisation except for PbChi7A. M.H.M. and T.S.N. carried out the kinetic analyses and TSN measured the pH, temperature and DCIP activity profiles. All biochemical data on PbChi7A were generated by O.R. and T.V.V. ERM provided the funding and the facilities for the work on PbChi7A. The NMR analyses were performed and interpreted by S.M. MHM performed the cellulose degradation assays and the HPAEC-PAD with support from S.G. M.H.M. and B.B. designed the SOD/HRP effect experiments, which were carried out by M.H.M. B.B. and M.H.M., F.B. and B.G. performed the EPR analysis.

D.H.W. and B.H. provided funding and interpretation of the crystallography and sequence analysis part. M.A.H. and M.H.M. analysed the data and wrote the first version of the manuscript. J.G.B. and B.B. contributed to the writing of later versions of the manuscript. All authors contributed and approved the final version of the manuscript.

## Competing interests

A provisional patent has been filed based on data from this paper by M.H.M., J.G.B. and M.A.H. The other authors declare that they have no competing interests.
**Additional information**

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-021-22372-0.

**Correspondence** and requests for materials should be addressed to J.-G.B. or M.A.H.

**Peer review information** *Nature Communications* thanks Igor Polikarpov and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021