Research | Open | Published:
Evolution of substrate specificity in bacterial AA10 lytic polysaccharide monooxygenases
Biotechnology for Biofuelsvolume 7, Article number: 109 (2014)
Understanding the diversity of lignocellulose-degrading enzymes in nature will provide insights for the improvement of cellulolytic enzyme cocktails used in the biofuels industry. Two families of enzymes, fungal AA9 and bacterial AA10, have recently been characterized as crystalline cellulose or chitin-cleaving lytic polysaccharide monooxygenases (LPMOs). Here we analyze the sequences, structures, and evolution of LPMOs to understand the factors that may influence substrate specificity both within and between these enzyme families.
Comparative analysis of sequences, solved structures, and homology models from AA9 and AA10 LPMO families demonstrated that, although these two LPMO families are highly conserved, structurally they have minimal sequence similarity outside the active site residues. Phylogenetic analysis of the AA10 family identified clades with putative chitinolytic and cellulolytic activities. Estimation of the rate of synonymous versus non-synonymous substitutions (dN/dS) within two major AA10 subclades showed distinct selective pressures between putative cellulolytic genes (subclade A) and CBP21-like chitinolytic genes (subclade D). Estimation of site-specific selection demonstrated that changes in the active sites were strongly negatively selected in all subclades. Furthermore, all codons in the subclade D had dN/dS values of less than 0.7, whereas codons in the cellulolytic subclade had dN/dS values of greater than 1.5. Positively selected codons were enriched at sites localized on the surface of the protein adjacent to the active site.
The structural similarity but absence of significant sequence similarity between AA9 and AA10 families suggests that these enzyme families share an ancient ancestral protein. Combined analysis of amino acid sites under Darwinian selection and structural homology modeling identified a subclade of AA10 with diversifying selection at different surfaces, potentially used for cellulose-binding and protein-protein interactions. Together, these data indicate that AA10 LPMOs are under selection to change their function, which may optimize cellulolytic activity. This work provides a phylogenetic basis for identifying and classifying additional cellulolytic or chitinolytic LPMOs.
The two most abundant polysaccharides in nature are cellulose and chitin . Plants, insects, crustaceans, molluscs, and fungi all utilize these two highly stable polymers as primary components of their cell walls. Deconstruction of polysaccharides is essential for ecosystem-level carbon and nitrogen cycling. Moreover, polysaccharides are potential energy sources that could help supplement the current massive demand for fossil fuels . Intensive efforts worldwide focus on conversion of these energy-rich biomolecules into free sugars that can be fermented into biofuels or other value-added bioproducts. However, hydrolysis of these polymers is difficult due to their crystalline structure, the stability of the β-glucosidic bond, and their close association with hemicellulose, lignin, and other modifying molecules [1, 3]. Cellulolytic and chitinolytic enzymes capable of this have been identified in a myriad of organisms, but most often in bacteria and fungi . While the biochemical activities and mechanisms of hydrolytic enzymes have been known for decades, oxygenolytic pathways for deconstruction of chitin and cellulose have only recently been identified [5–8].
CBH1, one of the first representatives of what are now recognized to be lytic polysaccharides monooxygenases (LPMOs), was secreted by Streptomyces olivaceoviridis and interacted with α-chitin, but since it lacked classical hydrolytic activity, it was thus considered to be a non-hydrolytic carbohydrate binding module (CBM) . An ortholog of CBH1, chitin-binding protein 21 (CBP21) was identified in Serratia marcescens  and initially classified as carbohydrate binding module 33 (CBM33, now systematically called Auxiliary Activity 10, AA10).a The function of CBP21 was first demonstrated by Vaaje-Kolstad et al. , who showed cleavage of crystalline chitin in an O2-dependent reaction. Soon after this report, others showed that the eukaryotic counterpart, fungal glycoside hydrolase 61 (GH61, now systematically called Auxiliary Activity 9, AA9) was a Cu2+-dependent enzyme [11–14]. An oxidative function has also been demonstrated for CelS2, an AA10 from Streptomyces coelicolor , which reacts synergistically with hydrolytic cellobiohydrolases and endoglucanases , and more recently for BlAA10A from Bacillus licheniformis and E8 from Thermobifida fusca  which react with chitin and cellulose, respectively, giving four AA10 enzymes whose function has been determined.
AA9 and AA10 incorporate a single 18O from 18O2 into polysaccharide cleavage products, and so are now classified as LPMOs . To date, structures of six AA9 and five AA10 enzymes have been solved, including one nuclear magnetic resonance (NMR) structure [11, 17–19]. Overall, the two LPMO families share a conserved β-sandwich fold , and many residues on the substrate-binding surface are conserved. Moreover, Cu2+ has been identified in the active sites [8, 17, 18]. Although recent computational studies support the involvement of a copper-oxyl radical intermediate , the catalytic mechanism of this reaction is still largely unexplored.
Oxidative polysaccharide cleavage results in the formation of an aldonic acid from C1 oxidation  or a ketoaldose from C4 oxidation [21, 22]. The monooxygenase reaction stoichiometry requires the addition of 2e- from an oxidoreductase or other external electron donor. The presence of oxidoreductases has been reported in various cellulolytic fungi , though an actual, physiological electron partner for LPMOs has not been unambiguously determined.
In this study, we compared amino acid sequences and protein structures in order to explore the evolutionary relatedness of AA9 and AA10. Conserved sequence and structural features were correlated with potential substrate interactions and surfaces potentially used by electron donors. Phylogenetic analysis suggests that cellulose- and chitin-specific enzymes are distributed into different subclades within bacterial AA10, as has been recently reported for the fungal AA9 [18, 21]. Potential evolutionarily pressures within the AA10 family were examined in order to understand how Darwinian selection might have influenced substrate specificity.
Structural comparison of LPMO families AA9 and AA10
Figure 1a shows five crystal structures from the AA9 family. These are from Hypocrea jecorina (Trichoderma reesei, Protein Data Bank (pdb) id: 2VTC) , Thielavia terrestris (pdb id: 3EII) , Thermoascus aurantiacus (pdb id: 2YET) , Neurospora crassa PMO-2 (pdb id: 4EIR) , and N. crassa PMO-3 (pdb id: 4EIS) . Structures of four AA10 enzymes are also shown in Figure 1b. These are from S. marcescens (pdb id: 2BEM) , Vibrio cholerae O1 biovar EI Tor (pdb id: 2XWX) , Burkholderia pseudomallei (pdb id: 3UAM), and Enterococcus faecalis (pdb id: 4A02) . Both AA9 and AA10 have a conserved β-sandwich fold with three to four β-sheet strands (Figure 1c and 1d). The average root mean square (RMS) deviation of the aligned structures is approximately 3 Å (Table 1) In addition to the fold-level similarity between AA9 and AA10, two key histidine (His) residues that coordinate a Cu2+ ion at their active sites are also highly conserved in both families (Table 1 and Figure 1c). The structural superposition of the metal ligands suggests that this configuration is essential for activity (inset in Figure 1c and 1d). A notable difference between AA9 and AA10 is the third, non-coordinating active site residue; being primarily tyrosine in the former and primarily phenylalanine in the latter, with a relatively few exceptions presently also identified.
Six AA10 structures from E. faecalis released in the pdb show copper in the active site, and a recently published structure of AA9 from Phanerochaete chrysosporium (pdb id: 4B5Q) also shows copper bound in the active site . Copper binds with nanomolar affinity to AA10 [8, 17]; its presence is consistent with O2 activation required for the LPMO reaction.
Surface electrostatic potential on the binding surfaces of AA9 and AA10
To explore factors that may contribute to substrate specificity in the AA9 and AA10 families, we characterized the electrostatic potential present at the substrate-binding surface. In both families, the metal-binding histidine residues are part of a planar surface that constitutes the polysaccharide-binding surface . Figure 1a and 1b show the surface electrostatic potential of representatives from both AA9 and AA10 families. For the AA9 proteins, which are biochemically characterized as cellulose monooxygenases, negatively charged residues (shown in red) prominently surround the active site (Figure 1a, yellow circle). In contrast, the AA10 chitin monooxygenases contain both positively charged (shown in blue) and negatively charged residues (shown in red) surrounding the active site (Figure 1b, yellow circle). Aachmann et al.  used NMR to identify residues from the chitinolytic AA10 enzyme from S. marcescens (pdb id: 2BEM) that are involved in chitin binding. These residues are Q53, Y54, E55, Q67, S58, L110, T111, A112, H114, and T116 . The positions of the corresponding residues from the other AA10 enzymes that align with 2BEM are shown as yellow on a grey surface in the lower parts of Figure 1a and 1b. In the other members of the AA10 family, most of these structurally conserved residues are also surface-exposed (Figure 1b, bottom). However, in the AA9 family, only a few are exposed at the polysaccharide-binding surface (Figure 1a, bottom), indicating that different residues from the folded structures will be involved in substrate binding in the AA9 and AA10 families.
Diversity of domain structures in AA9 and AA10 proteins
Another possible determinant of substrate specificity within the AA9 and AA10 families is the domain architecture. LPMO enzymes have a diverse composition of domains: they can be single catalytic domains, associated with various CBMs, or even associated with other catalytic domains (such as glycoside hydrolase (GH) domains). Figure 2 shows a Cytoscape (The Cytoscape Consortium, San Diego, CA) protein sequence homology network accounting for the variations in domains in the AA9 and AA10 families, where nodes represent enzymes or functional classes, and edges represent sequence similarity (bit score >200, evalue <1e-50). In order to prepare this network, sequences were collected from CAZy, compared via pairwise BLAST analysis, and then annotated with secondary CAZy domains. Nodes are colored according to their phylum-level taxonomic identification. The network contains 184 AA9 sequences and 495 AA10 sequences. All AA9 proteins were from eukaryotes, with a vast majority (99%) from the fungal phyla Ascomycota (135 sequences) and Basidiomycota (34 sequences). Of the protein sequence in the AA9 family, 31% include a secondary carbohydrate binding module 1 (CBM1), which has been reported to bind cellulose . Seven AA9 sequences are associated with CBM0, an unclassified CBM family [30, 31].
The AA10 family is exclusively from prokaryotes, with 226 sequences from Proteobacteria, 145 from Actinobacteria, and 132 from Firmicutes (Figure 2). There were no edges linking members of the AA9 and AA10 families at the similarity threshold of evalue <1e-50. Furthermore, when the similarity threshold was relaxed to 1e-5 there were still no connections between the AA9 and AA10 families. While Figure 2 shows that the AA9 network contains interspersed sequences from Ascomycota and Basidiomycota, the AA10 family shows clear taxonomic groupings assembled from different bacterial phyla. These results also show that while the active site residues of the AA9 and AA10 families are mostly conserved (Table 1), these two families do not share any other significant sequence similarities or consistent linkages to other domains.
Figure 2 also shows that the AA10 family is combined with a variety of secondary CBM domains, with 31% of the total sequences including cellulose-binding domains CBM2 and CBM3 [32, 33] or chitin-binding domains CBM5 and CBM12 . Further phylogenetic binning of AA10 showed expansion within the genera of Streptomyces, Bacillus, and Vibrio (Additional file 1: Figure S1). Interestingly, 94% of the AA10 sequences that included a cellulose-binding CBM were from the phylum Actinobacteria, whereas 95% of sequences including a chitin-binding CBM were from the phyla Firmicutes and Proteobacteria. Finally, two genes were identified that also encoded a glycoside hydrolase domain, suggesting a rare but possibly synergistic pairing of glycoside hydrolase and LPMO catalytic activities in a single enzyme.
Phylogenic analysis of LPMO families
To gain further insight into the evolutionary relationship and possible functional roles of the distinct LPMO families, we created phylogenetic trees representing the AA9 and AA10 families (Figures 3 and 4, respectively). Briefly, sequences were collected, curated to remove redundant sequences with 100% identity, aligned, trimmed to the conserved catalytic domain, and then the tree was constructed by MrBayes phylogenetic analysis . The resulting consensus tree was midpoint rooted and annotated with associated carbohydrate-binding modules in addition to the AA9 or AA10 catalytic domains. The five crystal structures determined for AA9, 2YET, 2VTC, 4EIS, 4EIR, and 3EII, were mapped onto the phylogenetic tree. In Figure 3, the surfaces of these structures have been colored to identify highly conserved residues shared across the AA9 family. The tree was also annotated to indicate whenever a putative cellobiose dehydrogenase (AA3 family enzymes) was present in the host genome using a cutoff criterion of 35% identity to N. crassa CDH1. The ability of CDH to act as the proximal electron donor for LPMO in cellulose oxidative cleavage has been demonstrated in this organism [18, 36–38].
The AA9 LPMOs have been classified into four functional types based on their reaction products . These are shown in Figure 3 as red boxes. LPMO1 enzymes hydroxylate the C1 position of pyranose rings and produce an aldonolactone [18, 21], while LPMO2 enzymes hydroxylate the C4 position of pyranose rings and produce a 4-ketoaldose [21, 22]. LPMO3 enzymes are less specific [13, 19, 21, 39], and produce both aldonolactone and non-reducing end oxidized products, while LPMO3* produce only aldonic acids .Mapping of the four LPMO subgroups onto the global AA9 phylogeny showed that the LPMO2, LPMO3, and LPMO3* subgroups are monophyletic, with each having a single phylogenetic clade that corresponds to distinct functional classes (red boxes). In contrast, LPMO1 enzymes span a major evolutionary division as two branches cross into this functional class, indicating more sequence diversity in the LPMO1 family. Examples where all four LPMO functional types were fused to additional CBM domains are identified in Figure 3. Moreover, Figure 3 also shows that the majority of AA9 proteins come from organisms that also contain a cellobiose dehydrogenase homolog.
The AA10 phylogenetic tree was generated in a similar manner using the catalytic domains of all non-redundant sequences present in the CAZy database. The AA10 tree shown in Figure 4 represents 374 non-redundant sequences that are entirely bacterial in origin. The tree was annotated with secondary CBM domains (central column), and divided into two major clades (clade I and clade II) that could be subdivided into four additional subclades (A through D). The biochemically characterized cellulose-oxidizing LPMOs from S. coelicolor (A3) and T. fusca are present in subclade A,b while all other LPMOs with experimental confirmation of their reaction with chitin are present in subclades C and D .The tree was also annotated with microarray-based gene expression data for the six variants of AA10 present in Streptomyces sp. SirexAA-E (SirexAA-E) . Clade I contains a delineated mixture of phyla, with subclade C containing sequences only from Actinobacteria and with subclade D containing sequences from Firmicutes and Proteobacteria. Clade II is primarily composed of Actinobacteria and separates into subclades A and B. Subclades A and B contain only cellulose-binding CBMs (CBM2 and CBM3) associated with the catalytic AA10 domain, whereas subclades C and D contain only chitin-binding CBMs (CBM5 and 12). Furthermore, expression data from SirexAA-E shows that genes from subclades A and B were selectively upregulated only when cells were grown in medium containing cellulose as the sole carbon source, while genes from subclade C were upregulated only during growth on chitin .
The cellulose-oxidizing LPMOs from AA10 are primarily present in Actinobacteria, an aerobic filamentous bacterial phyla found in soil, but also associated with insects and other animals . In Figure 4, the structures of four AA10 enzymes are mapped to the phylogenetic tree: 3UAM, 4A02, 2BEM, and 2XWX. Additionally, predicted protein structures for expressed AA10 from SirexAA-E are mapped onto the tree. There is high amino-acid sequence identity among the AA10 proteins whose structures have been determined, with the highest sequence conservation observed at the active site (magenta color). Interestingly, homology models consistently predict an additional surface exposed loop region on the same side of the protein as the active site in clade II proteins (chitin oxidation), but not in clade I (cellulose oxidation). The position of this loop can be recognized in pdb id: 4GBO, the E7 enzyme from T. fusca . Recently, Vu et al. have identified a role for these extra loops in substrate recognition and control of specificity of reaction in the AA9 family .
Homology modeling of AA10 proteins and conserved sequence motif in LPMOs
Several LPMOs from within the AA10 family have been experimentally verified to be either chitin or cellulose monooxygenases (such as CBP21 and BlAA10A) which react with chitin, and CelS2 and E8, which react with cellulose [6, 25].b To further explore structural determinants that control substrate specificity, we compared homology models for 43 proteins that spanned the AA10 family (Figure 5) across the clade I and clade II sequences shown in Figure 4. Homology modeling using I_TASSER , followed by superposition of the modeled structures showed that the most significant structural differences were located in the substrate binding region (Figure 5, Additional file 2: Table S1). Specifically, the positions of loops (shown for illustration purposes only) on the substrate-binding side of the protein had more variations than other parts of the modeled structures. Correspondingly, the insertion observed in the sequence alignments mapped to loops on the substrate-binding side of the AA10 family. Given the structural variability of clades I and II, and differences in measured catalytic functions, it is likely that these structural differences help to modulate substrate selectivity between chitin and cellulose in AA10, as now predicted for AA9.
To improve our mapping of potential functional determinants onto the modeled structures, Multiple EM for Motif Elicitation (MEME)  was used. This approach identified three sequence motifs among the 43 AA10 proteins (Figure 5). These motifs were mapped back onto the structures and homology models. Simultaneously, MEME was used to determine whether there were significant motifs observed in the published structures of AA9 (Figure 5).In the homology-modeled proteins (shown in Figure 5 corresponding to the four AA10 clades shown in Figure 4), the three MEME motifs ranged from about 25 to 41 residues in length. Motif1 was present in both AA9 and AA10, and contained the variable insertion regions that possibly yield substrate selectivity in the AA10 family. Motif2 and motif3 were observed only in AA10 (Figure 5b). It is interesting to note that the difference in the number of motifs identified in AA9 as compared to AA10 provides an additional line of evidence supporting the possibility of evolutionary selection in these two families. Although experimental evidence as to what these motifs (Figure 5c-e) contribute is currently lacking, it is clear from the superposition of the homology-modeled structures (Figure 5, cyan sequences) that motif1 is well-positioned to play a role in substrate binding and discrimination between binding to chitin or cellulose. Interestingly, motif2 and motif3 span the breadth of the protein and connect the substrate-binding surface to the opposite side of the protein where potential electron donor proteins might interact.
Evolution of chitinolytic and cellulolytic subclades within the AA10 family
To evaluate the selective pressure on these functionally defined clades, the rates of non-synonymous and synonymous codon substitutions (dN and dS, respectively) of the catalytic domain were estimated (Figure 6a). Pairwise comparisons were performed against all genes within either subclade A or subclade D. Estimations with dN values greater than 0.01 and dS values less than 1.5 were reported to allow sufficient mutational signal and to avoid the effect of back mutations that would artificially increase dS and reduce dN with increased sequence divergence . Pairwise comparisons indicate that a group of the chitinolytic genes from subclade D are primarily under negative selection (Figure 6, dN/dS <0.2), while a second group is under more neutral selection (1 > dN/dS >0.2). However, genes from subclade A have a significantly different distribution than from subclade D. Very few genes from cellulolytic subclade A show negative selection, while a significant proportion show increased positive selection (dN/dS >1). To confirm these results, site-specific dN/dS values were estimated for subclades A and D (Figure 6b). The results show that a significant number of residues in subclade A were indeed positively selected, while residues in subclade D were all negatively or neutrally selected. When plotted on the protein structures, the negatively selected sites in both subclades A and D are primarily located around the active-site residues. In contrast, most of the positively selected residues in subclade A are surface exposed, including regions on the putative substrate-binding surface and along the interior of the protein traversing from the substrate-binding surface to the opposite surface of the protein. Interestingly, this latter region may provide a surface for interaction with accessory redox proteins such as cellobiose dehydrogenase (AA3 enzymes).
In this study, we analyzed the AA9 and AA10 families using available protein structures and sequence information to evaluate differences between and within the families, to explore features that influence substrate specificity, and to characterize selective pressures that may have led to functional diversification.
LPMOs share a common structural fold and a spatial conservation of active site residues, as seen by their low root mean square deviation (RMSD) values (ranging up to 3.3 Å, Table 1). While the core structural folds and the active site geometry of these two LPMO families are similar, there is low homology at the amino-acid sequence level, and the surface electrostatic potentials at the substrate-binding surface show considerable differences in charge distributions. Indeed, comparison of all AA9 and AA10 proteins available in the CAZy database failed to identify any sequences from across these two families that have significant homology (evalue <1e-5). Our results indicate that although AA9 and AA10 families share structural similarities, they have so significantly diverged from a common ancestor that the only residue-level homology that remains is in the active site residues.
Due to the low sequence similarity between AA9 and AA10 families we analyzed their phylogenetic relationships separately. The AA9 phylogenetic tree is separated into three major evolutionarily related groups which partially correspond to the four types of enzyme activity observed for LPMOs . LPMO2, LPMO3, and LPMO3* enzyme activities correspond to monophyletic clades, which suggests vertical inheritance and conserved enzyme functions within each clade. In contrast, LPMO1 enzymes are present in a polyphyletic clade, indicating a more diverse sequence space and potentially varied enzyme function. Sequences from Ascomycetes and Basidiomycetes are scattered throughout the three major evolutionarily related groups in AA9, suggesting an ancestral sequence that was shared before these two phyla separated.
The AA10 phylogenetic tree was separated into two major phylogenetic groups. When annotated with known activities, the two clades appear to separate enzymes with different substrate specificities. Clade I contains all biochemically defined chitin monooxygenases, while clade II contains subclades that are either cellulose or chitin monooxygenases. Gene expression data from SirexAA-E grown on either chitin or cellulose as the sole carbon source further corroborates this assessment . We also observed that CBM domain composition varies between clade I and II. Clade I is dominated by CBM5 and 12 domains, which are primarily chitin binding, but possibly can have a lignin-binding function as well . Clade II is enriched in CBM2 domains, which are primarily associated with cellulose binding. Most recently, Forsberg et al. showed the binding specificity of CelS2 either with or without the associated CBM2 domain . Interestingly, although this CBM2 domain was tightly bound to either α or β-chitin, the corresponding AA10 domain (CelS2) only reacted with cellulose. Further biochemical verification will be necessary to extend these observations more broadly into phylogenetic space.
To identify sequence and structural features that may contribute to clade II activity against cellulose, we generated homology-modeled structures for 43 sequences that span the phylogeny in AA10. Using MEME, these homology-modeled structures were identified to have three highly significant motifs, where motif1 shows the largest structural variability. Specifically, this variable motif is contained in a loop of un-modeled sequence at the substrate-binding surface and is only found in subclade A. Subclade A of AA10 contains biochemically characterized cellulose monooxygenases, and also contains the most highly upregulated AA10 enzyme when SirexAA-E is grown on cellulose . We hypothesize that this additional sequence at the binding surface is a defining feature of cellulose-active AA10 enzymes, paralleling the identification of a loop-modulating reaction specificity in the AA9 enzymes . Motif2 and motif3, which span the breadth of the protein, connect the substrate-binding surface to the opposite side of the protein. This suggests a possibility for modulation of electron donor interactions.
Finally, we explored the selective pressures within two clades of the AA10 family to understand how diversification may be distributed in this enzyme family. The results show that chitinolytic enzymes in subclade D (chitinolytic enzymes) have mostly negative selection at both the whole gene and site-specific levels. In contrast, subclade A (both chitinolytic and cellulolytic enzymes) contains more genes with diversifying selection at both the whole gene and site-specific levels. This result indicates that subclade A may have undergone a change in substrate specificity and that genes within this clade are potentially being selected for increased activity. Together, these data suggest that the ancestral form of AA10 may have been a chitin monooxygenase, and that clade II has apparently further specialized for cellulose oxidation. Selection may be towards more favorable substrate binding, better interactions with accessory redox proteins, such as cellobiose dehydrogenase enzymes, or perhaps both.
In summary, this study provides a better understanding of the evolution of functional diversity within the recently discovered AA9 and AA10 LPMO families. Together, these data suggest that AA9 and AA10 families share a distant common ancestor. Furthermore, clades within the AA10 family are specialized for different substrates and subclade A has undergone diversifying selection at surface-exposed regions of the protein.
Materials and methods
Sequence similarity network
AA9 and AA10 protein-coding sequences were identified on the Carbohydrate-Active Enzyme (CAZy) database , and harvested from the National Center for Biotechnology Information (NCBI) protein database. All AA9 and AA10 sequences were compared against each other using BLAST  to identify similar proteins. All sequences were also re-annotated with CAZy families to identify the domain structure of each protein. This data was then used to build a similarity network using Cytoscape 2.8.0 , and visualized as an organic layout. Nodes in the network represent unique protein sequences and CAZy families. Edges represent a BLAST bit score of ≥200 (evalue ≥1 × e-50), or an annotation to a CAZy category. Nodes were annotated with taxonomic information at the phylum level.
Phylogenetic tree construction
AA9 and AA10 phylogenetic trees were constructed by first identifying proteins from the CAZy database, and the harvesting sequence from NCBI. Sequences from either AA9 or AA10 families were aligned using Multiple Sequence Comparison by Log-Expectation (MUSCLE) on the Cyberinfrastructure for Phylogenetic Research (CIPRES, https://www.phylo.org/portal2/login!input.action) Science Gateway . Aligned sequences were then trimmed to retain only the AA9 or AA10 domain; sequences lacking the conserved active site His residues were removed from the alignment. Phylogenetic trees were generated using MrBayes code with a calculated standard deviation of ≤0.05. Non-default parameters were set to mcmc, ngen = 10,000,000, temp = 0.200, burninfrac = 0.25, stoprule = No, sump burnin = 4000, and sumt burnin = 4000. Resulting trees were annotated with pfam, phyla, solved structures, and cellobiose dehydrogenase homolog information.
Evolutionary rate estimation
Coding sequences for subclades A-D from the AA10 family were collected and codon alignments were generated with MUSCLE. Sequences were trimmed to retain only the AA10 domain. Codon alignments were masked with Zorro (http://phylogenomics.wordpress.com/software/zorro/) to generate quality scores for codon positions , and then a phylogenetic tree was generated with RAxML (http://sco.h-its.org/exelixis/web/software/raxml/index.html) using the masking scores . Pairwise codon substitution models (dN/dS values) were estimated using the CODEML program in the PAML package [51, 52]. Variables were set at CodonFreq = 0 and model = 0. Only pairwise dN values with values ≥0.01 and dS values ≤1.5 were reported so as to allow for sufficient mutational signal and to avoid the effects of back mutations. Site-specific codon substitution models were generated using the CODEML program in PAML, with model = 0, NSsites = 3, ncatG = 3 fix_kappa = 0, fix_omega = 0, cleandata = 1, and fix_blength = 2.
Protein three-dimensional structure comparison
The Dali protein structure alignment database (http://ekhidna.biocenter.helsinki.fi/dali_server/) was used to calculate %RMSD and %ID of LPMO enzymes whose structures are known using 2BEM as a query . Structures with the ten best%RMSD are shown in Table 1.
The 43 sequences highlighted in Figure 4 were considered for prediction of three-dimensional structures using iterative threading assembly refinement algorithm (I-TASSER) (http://zhanglab.ccmb.med.umich.edu/I-TASSER/) . For each sequence, the signal peptides and other domains besides the Cu2+-binding catalytic domain were included in the homology modeling. Alignments used for modeling are tabulated in Additional file 2: Table S1. Models obtained with the highest C-score were retained for further analysis. The homology models have been deposited at Model Archive (doi:10.5452/ma-asp8e) .
Structural comparisons were done using the Combinatorial Extension algorithm  implemented in PyMOL (Schrödinger, Portland, OR). Protein surface electrostatics calculations were carried using Adaptive Poisson-Boltzmann Solver (APBS) , where an externally generated pdb (P) file with per-atom charge (Q) and radius (R) (PQR file) file was used to calculate the electrostatics. The parameters used were solvent and protein dielectrics of 78.0 and 2.0 respectively, solvent radius of 1.4, and a monovalent ion concentration of 0.15 M. The visualization was depicted in PyMOL with positive and negative molecular surface ranging from -2kT/e to 2kT/e.
Sequence-based motifs were identified using Multiple Em for Motif Elicitation (MEME) (http://meme.nbcr.net/meme/) . The occurrence of motifs in the sequence was assumed to be distributed either zero or one per sequence. Three motifs were identified for each set of sequences given. The phylogenetic tree of AA10 was divided into four subclades (A to D) based on major phylogenetic clades. For each clade, the motifs were identified using MEME.
aCBM33 has recently been renamed as AA10; likewise GH61 has been renamed as AA9 . These names will be used throughout.
bSACTE_3159 from the highly cellulolytic Streptomyces sp. SirexAA-E, has also been confirmed to contain Cu2+ and have O2-dependent cellulose oxidation activity (M. Mbughuni and BG Fox, unpublished data).
Lytic polysaccharide monooxygenases
Himmel ME, Ding SY, Johnson DK, Adney WS, Nimlos MR, Brady JW, Foust TD: Biomass recalcitrance: engineering plants and enzymes for biofuels production. Science 2007, 315: 804-807.
Faaij APC: Bio-energy in Europe: changing technology choices. Energ Policy 2006, 34: 322-342.
Wilson DB: Microbial diversity of cellulose hydrolysis. Curr Opin Microbiol 2011, 14: 259-263.
Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS: Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev 2002, 66: 506-577. Table of contents
Culpepper MA, Rosenzweig AC: Architecture and active site of particulate methane monooxygenase. Crit Rev Biochem Mol Biol 2012, 47: 483-492.
Forsberg Z, Vaaje-Kolstad G, Westereng B, Bunaes AC, Stenstrom Y, MacKenzie A, Sorlie M, Horn SJ, Eijsink VG: Cleavage of cellulose by a CBM33 protein. Protein Sci 2011, 20: 1479-1483.
Harris PV, Welner D, McFarland KC, Re E, Navarro Poulsen JC, Brown K, Salbo R, Ding H, Vlasenko E, Merino S, Xu F, Cherry J, Larsen S, Lo Leggio L: Stimulation of lignocellulosic biomass hydrolysis by proteins of glycoside hydrolase family 61: structure and function of a large, enigmatic family. Biochemistry 2010, 49: 3305-3316.
Hemsworth GR, Taylor EJ, Kim RQ, Gregory RC, Lewis SJ, Turkenburg JP, Parkin A, Davies GJ, Walton PH: The copper active site of CBM33 polysaccharide oxygenases. J Am Chem Soc 2013, 135: 6069-6077.
Schnellmann J, Zeltins A, Blaak H, Schrempf H: The novel lectin-like protein CHB1 is encoded by a chitin-inducible Streptomyces olivaceoviridis gene and binds specifically to crystalline alpha-chitin of fungi and other organisms. Mol Microbiol 1994, 13: 807-819.
Suzuki K, Suzuki M, Taiyoji M, Nikaidou N, Watanabe T: Chitin binding protein (CBP21) in the culture supernatant of Serratia marcescens 2170. Biosci Biotechnol Biochem 1998, 62: 128-135.
Vaaje-Kolstad G, Westereng B, Horn SJ, Liu Z, Zhai H, Sorlie M, Eijsink VG: An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Science 2010, 330: 219-222.
Levasseur A, Drula E, Lombard V, Coutinho PM, Henrissat B: Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. Biotechnol Biofuels 2013, 6: 41.
Quinlan RJ, Sweeney MD, Lo Leggio L, Otten H, Poulsen JC, Johansen KS, Krogh KB, Jorgensen CI, Tovborg M, Anthonsen A, Tryfona T, Walter CP, Dupree P, Xu F, Davies GJ, Walton PH: Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. Proc Natl Acad Sci U S A 2011, 108: 15079-15084.
Westereng B, Ishida T, Vaaje-Kolstad G, Wu M, Eijsink VG, Igarashi K, Samejima M, Stahlberg J, Horn SJ, Sandgren M: The putative endoglucanase PcGH61D from Phanerochaete chrysosporium is a metal-dependent oxidative enzyme that cleaves cellulose. PLoS One 2011, 6: e27807.
Moser F, Irwin D, Chen S, Wilson DB: Regulation and characterization of Thermobifida fusca carbohydrate-binding module proteins E7 and E8. Biotechnol Bioeng 2008, 100: 1066-1077.
Forsberg Z, Rohr AK, Mekasha S, Andersson KK, Eijsink VG, Vaaje-Kolstad G, Sorlie M: Comparative study of two chitin-active and two cellulose-active AA10-type lytic polysaccharide monooxygenases. Biochemistry 2014, 53: 1647-1656.
Aachmann FL, Sorlie M, Skjak-Braek G, Eijsink VG, Vaaje-Kolstad G: NMR structure of a lytic polysaccharide monooxygenase provides insight into copper binding, protein dynamics, and substrate interactions. Proc Natl Acad Sci U S A 2012, 109: 18779-18784.
Li X, Beeson WT, Phillips CM, Marletta MA, Cate JH: Structural basis for substrate targeting and catalysis by fungal polysaccharide monooxygenases. Structure 2012, 20: 1051-1061.
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: 1399-1406.
Kim S, Stahlberg J, Sandgren M, Paton RS, Beckham GT: Quantum mechanical calculations suggest that lytic polysaccharide monooxygenases use a copper-oxyl, oxygen-rebound mechanism. Proc Natl Acad Sci U S A 2014, 111: 149-154.
Vu VV, Beeson WT, Phillips CM, Cate JHD, Marletta MA: Determinants of regioselective hydroxylation in the fungal polysaccharide monooxygenases. J Am Chem Soc 2014, 136: 562-565.
Isaksen T, Westereng B, Aachmann FL, Agger JW, Kracher D, Kittl R, Ludwig R, Haltrich D, Eijsink VG, Horn SJ: A C4-oxidizing lytic polysaccharide monooxygenase cleaving both cellulose and cello-oligosaccharides. J Biol Chem 2014, 289: 2632-2642.
Hori C, Gaskell J, Igarashi K, Samejima M, Hibbett D, Henrissat B, Cullen D: Genomewide analysis of polysaccharides degrading enzymes in 11 white- and brown-rot Polyporales provides insight into mechanisms of wood decay. Mycologia 2013, 105: 1412-1427.
Karkehabadi S, Hansson H, Kim S, Piens K, Mitchinson C, Sandgren M: The first structure of a glycoside hydrolase family 61 member, Cel61B from Hypocrea jecorina, at 1.6 A resolution. J Mol Biol 2008, 383: 144-154.
Vaaje-Kolstad G, Houston DR, Riemen AH, Eijsink VG, van Aalten DM: Crystal structure and binding properties of the Serratia marcescens chitin-binding protein CBP21. J Biol Chem 2005, 280: 11313-11319.
Wong E, Vaaje-Kolstad G, Ghosh A, Hurtado-Guerrero R, Konarev PV, Ibrahim AF, Svergun DI, Eijsink VG, Chatterjee NS, van Aalten DM: The Vibrio cholerae colonization factor GbpA possesses a modular structure that governs binding to different host surfaces. PLoS Pathog 2012, 8: e1002373.
Vaaje-Kolstad G, Bohle LA, Gaseidnes S, Dalhus B, Bjoras M, Mathiesen G, Eijsink VG: Characterization of the chitinolytic machinery of Enterococcus faecalis V583 and high-resolution structure of its oxidative CBM33 enzyme. J Mol Biol 2012, 416: 239-254.
Wu M, Beckham GT, Larsson AM, Ishida T, Kim S, Payne CM, Himmel ME, Crowley MF, Horn SJ, Westereng B, Igarashi K, Samejima M, Ståhlberg J, Eijsink VG, Sandgren M: Crystal structure and computational characterization of the lytic polysaccharide monooxygenase GH61D from the Basidiomycota fungus Phanerochaete chrysosporium. J Biol Chem 2013, 288: 12828-12839.
Blake AW, McCartney L, Flint JE, Bolam DN, Boraston AB, Gilbert HJ, Knox JP: Understanding the biological rationale for the diversity of cellulose-directed carbohydrate-binding modules in prokaryotic enzymes. J Biol Chem 2006, 281: 29321-29329.
Hogg D, Pell G, Dupree P, Goubet F, Martin-Orue SM, Armand S, Gilbert HJ: The modular architecture of Cellvibrio japonicus mannanases in glycoside hydrolase families 5 and 26 points to differences in their role in mannan degradation. Biochem J 2003, 371: 1027-1043.
Wang N, Zhang Y, Wang Q, Liu J, Wang H, Xue Y, Ma Y: Gene cloning and characterization of a novel alpha-amylase from alkaliphilic Alkalimonas amylolytica. Biotechnol J 2006, 1: 1258-1265.
Lin ES, Wilson DB: Identification of a celE-binding protein and its potential role in induction of the celE gene in Thermomonospora fusca. J Bacteriol 1988, 170: 3843-3846.
Zhang S, Lao G, Wilson DB: Characterization of a Thermomonospora fusca exocellulase. Biochemistry 1995, 34: 3386-3395.
Huang L, Garbulewska E, Sato K, Kato Y, Nogawa M, Taguchi G, Shimosaka M: Isolation of genes coding for chitin-degrading enzymes in the novel chitinolytic bacterium, Chitiniphilus shinanonensis, and characterization of a gene coding for a family 19 chitinase. J Biosci Bioeng 2012, 113: 293-299.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP: MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol 2012, 61: 539-542.
Bey M, Zhou S, Poidevin L, Henrissat B, Coutinho PM, Berrin JG, Sigoillot JC: Cello-oligosaccharide oxidation reveals differences between two lytic polysaccharide monooxygenases (family GH61) from Podospora anserina. Appl Environ Microbiol 2013, 79: 488-496.
Kittl R, Kracher D, Burgstaller D, Haltrich D, Ludwig R: Production of four Neurospora crassa lytic polysaccharide monooxygenases in Pichia pastoris monitored by a fluorimetric assay. Biotechnol Biofuels 2012, 5: 79.
Langston JA, Shaghasi T, Abbate E, Xu F, Vlasenko E, Sweeney MD: Oxidoreductive cellulose depolymerization by the enzymes cellobiose dehydrogenase and glycoside hydrolase 61. Appl Environ Microbiol 2011, 77: 7007-7015.
Beeson WT, Phillips CM, Cate JH, Marletta MA: Oxidative cleavage of cellulose by fungal copper-dependent polysaccharide monooxygenases. J Am Chem Soc 2012, 134: 890-892.
Takasuka TE, Book AJ, Lewin GR, Currie CR, Fox BG: Aerobic deconstruction of cellulosic biomass by an insect-associated Streptomyces. Sci Rep 2013, 3: 1030.
Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc 2010, 5: 725-738.
Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res 2006, 34: W369-W373.
Hurst LD: The Ka/Ks ratio: diagnosing the form of sequence evolution. Trends Genet 2002, 18: 486.
Bianchetti CM, Harmann CH, Takasuka TE, Hura GL, Dyer K, Fox BG: Fusion of dioxygenase and lignin-binding domains in a novel secreted enzyme from cellulolytic streptomyces sp SirexAA-E. J Biol Chem 2013, 288: 18574-18587.
Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B: The carbohydrate-active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res 2009, 37: D233-D238.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25: 3389-3402.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 13: 2498-2504.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004, 32: 1792-1797.
Wu M, Chatterji S, Eisen JA: Accounting for alignment uncertainty in phylogenomics. PLoS One 2012, 7: e30288.
Stamatakis A: RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30: 1312-1313.
Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci 1997, 13: 555-556.
Yang Z, Rannala B: Bayesian phylogenetic inference using DNA sequences: a Markov Chain Monte Carlo Method. Mol Biol Evol 1997, 14: 717-724.
Holm L, Rosenstrom P: Dali server: conservation mapping in 3D. Nucleic Acids Res 2010, 38: W545-W549.
The Model Archive http://www.modelarchive.org/doi/10.5452/ma-asp8e
Shindyalov IN, Bourne PE: Protein structure alignment by incremental combinatorial extension (CE) of the optimal path. Protein Eng 1998, 11: 739-747.
Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA: Electrostatics of nanosystems: application to microtubules and the ribosome. Proc Natl Acad Sci U S A 2001, 98: 10037-10041.
This work was funded by the Department of Energy Great Lakes Bioenergy Research Center (DOE BER Office of Science DE-FC-02-07ER64494), the Department of Energy Bringing Advanced Computational Techniques to Energy Research program (DOE BACTER DE-FG02-04ER25627), the National Institute of Health Protein Structure Initiative (U01 GM098248), and the National Science Foundation GRAPE (CMMI-0941013).
The authors declare that they have no competing interests.
AJB, RMY, and TET, designed and performed the experiments. All authors reviewed, discussed, and interpreted the results. AJB, RMY, TET, GNP, CRC, and BGF wrote and reviewed the manuscript. All authors read and approved the final manuscript.
Adam J Book, Ragothaman M Yennamalli, Taichi E Takasuka contributed equally to this work.