- Open Access
Discovery, activity and characterisation of an AA10 lytic polysaccharide oxygenase from the shipworm symbiont Teredinibacter turnerae
Biotechnology for Biofuels volume 12, Article number: 232 (2019)
The quest for novel enzymes for cellulosic biomass-degradation has recently been focussed on lytic polysaccharide monooxygenases (LPMOs/PMOs), Cu-containing proteins that catalyse the oxidative degradation of otherwise recalcitrant polysaccharides using O2 or H2O2 as a co-substrate.
Although classical saprotrophic fungi and bacteria have been a rich source of lytic polysaccharide monooxygenases (LPMOs), we were interested to see if LPMOs from less evident bio-environments could be discovered and assessed for their cellulolytic activity in a biofuel context. In this regard, the marine shipworm Lyrodus pedicellatus represents an interesting source of new enzymes, since it must digest wood particles ingested during its natural tunnel boring behaviour and plays host to a symbiotic bacterium, Teredinibacter turnerae, the genome of which has revealed a multitude of enzymes dedicated to biomass deconstruction. Here, we show that T. turnerae encodes a cellulose-active AA10 LPMO. The 3D structure, at 1.4 Å resolution, along with its EPR spectrum is distinct from other AA10 polysaccharide monooxygenases insofar as it displays a “histidine-brace” catalytic apparatus with changes to the surrounding coordination sphere of the copper. Furthermore, TtAA10A possesses a second, surface accessible, Cu site 14 Å from the classical catalytic centre. Activity measurements show that the LPMO oxidises cellulose and thereby significantly augments the rate of degradation of cellulosic biomass by classical glycoside hydrolases.
Shipworms are wood-boring marine molluscs that can live on a diet of lignocellulose. Bacterial symbionts of shipworms provide many of the enzymes needed for wood digestion. The shipworm symbiont T. turnerae produces one of the few LPMOs yet described from the marine environment, notably adding to the capability of shipworms to digest recalcitrant polysaccharides.
Harnessing the activity of enzymes for the breakdown of cellulose and related plant-cell wall polysaccharides is important in the quest for environmentally sustainable fuels in the form of second-generation biofuel, from cellulosic wastes and bespoke energy crops. In this context, one of the most significant breakthroughs has been the discovery [1,2,3] of chain-breaking “lytic” polysaccharide monooxygenases (LPMOs, sometimes PMOs) that oxidatively deconstruct recalcitrant polysaccharides, most notably cellulose, chitin, hemicelluloses and retrograded starch (generic LPMOs have been extensively reviewed, examples include [4,5,6,7,8]).
When used as part of enzyme cocktails, LPMOs significantly boost the activity of classical glycoside hydrolases, thereby offering great potential for the sustainable breakdown of recalcitrant biomass within a commercial setting. As such, there is an on-going search for new variants of LPMOs, particularly from organisms which are effective and voracious decomposers of biomass (for example ).
Until now, most cellulose-active LPMOs have been isolated and characterised from wood-decaying fungi and soil bacteria. As a complement to these sources of new enzymes, more complex animals and their microbial symbionts present a different biological context for biomass-enzyme discovery. In the present work, we examined marine xylophagous bivalve molluscs—called shipworms—that burrow through submerged wood from which the organism derives nutrients, mainly sugars [10, 11]. Shipworms are a major cause of damage to submerged timber structures.
Within their gills, shipworms harbour communities of endosymbiotic bacteria in specialised eukaryotic cells (bacteriocytes) . Recent transcriptomics and proteomics analyses of dissected organs revealed that wood digestion in shipworms is accomplished through the combined action of enzymes of both endogenous and symbiotic origin [13, 14]. These studies demonstrated that shipworm gill endosymbionts produce a multitude of carbohydrate active enzymes, in which the bacterial genes coding for predicted (CAZY family, http://www.cazy.org ) AA10 LPMOs and GH6 cellobiohydrolases are among the most highly expressed. Compositional analysis of undigested wood and shipworm faeces (frass) also revealed that, while hemicellulose and lignin were virtually unaffected, over 40% of the cellulose was removed after digestion by the shipworm Lyrodus pedicellatus . Such high levels of activity signal that the constituent enzymes may be particularly fruitful targets within the context of the on-going search for biomass degrading enzymes for biofuel production.
Here, we report the in-depth biochemical and structural characterisation of the first LPMO from a shipworm endosymbiont (Teredinibacter turnerae ), the recombinant protein of which was produced heterologously in Escherichia coli. Through mass spectrometry and HPLC analysis of released products, we show that this enzyme has mixed C1–C4 oxidising activity on cellulose and boosts the breakdown of this recalcitrant polysaccharide by glycoside hydrolases. X-ray diffraction and multi-frequency EPR spectroscopy studies reveal a near classical AA10 LPMO secondary/tertiary structure, yet one that contains a second Cu-binding site. The coordinating residues of both copper sites are conserved in close homologs of Teredinibacter turnerae. The histidine brace site, composed of the N-terminal histidine and a second histidine side chain, is conserved across nearly all known LPMOs. The second Cu site is, however, not conserved across the wider LPMO family and its function within the Teredinibacter turnerae LPMO is unclear, although it may be related to charge transfer pathways through the protein.
Expression and enzymatic characterization of the AA10 LPMO from T. turnerae
The gamma-proteobacterium T. turnerae is the only endosymbiont found within the shipworm gills to have been successfully isolated, cultured and had its genome mapped . Through automated annotation and manual BLAST  searches of the predicted T. turnerae proteome, we identified one gene (NCBI Reference Sequence: WP_019602454.1) coding for an AA10 LPMO (hereafter TtAA10A). The predicted protein sequence features an N-terminal signal peptide, the LPMO domain and a serine-rich linker region followed by a carbohydrate binding module (CBM) 10 domain (Fig. 1a). AA10s have been found with appended CBM2, CBM3, CBM5, CBM10, CBM12, CBM18 and CBM73 domains (Bernard Henrissat personal communication) and are known to be active on cellulose or chitin. CBM10 domains are thought to be cellulose binding and so may provide cellulose recognition that is unlikely to be associated with a catalytic event . Whilst different from the CBMs commonly found attached to AA10 proteins , its presence in the domain structure of the TtAA10A gene gives an indication that this protein may be primarily active on glucose-based polysaccharides.
After multiple attempts to express the gene with various affinity, solubility tags and different secretion signals, sufficient protein for analysis was finally obtained through the production of a C-terminally strep-tagged LPMO catalytic domain (from His25 to Gly228) in E. coli (Fig. 1b). The purified tagged protein was loaded with excess copper, de-salted through size-exclusion chromatography, analysed for purity through SDS-PAGE (Fig. 1c) and mass spectrometry-based protein ID (not shown), and used for subsequent experiments.
Recombinant TtAA10A (catalytic domains, 25–228, only) displays the hallmarks of a correctly folded AA10. Thermal shift analysis (Thermofluor) of purified, Cu-loaded TtAA10A indicates a melting temperature (Tm) of 50.4 °C. Stripping copper with 10 mM EDTA lowers the Tm to 42.5 °C, suggesting a protein stabilising effect by the metal cofactor, as reported in previous literature for other LPMOs (for example [20, 21], Fig. 1d). We also noticed a variability in protein preparations, with some preparations containing a single (active-centre) Cu, whilst others contained two Cu atoms, described below.
Activity assays, on both single and double Cu site samples, were carried out on a range of commercial polysaccharide substrates (Avicel, β-chitin from squid pen, α-chitin from shrimp shell, cellohexaose, corn starch, pachyman, beechwood xylan, glucomannan, xyloglucan, lichenan, galactan, galactomannan and mannan) in the presence of the reducing co-factor, gallic acid. Samples were analysed after 24 h by MALDI-TOF MS and peak masses of the reaction products compared to previously published data, revealing a mixed C1–C4 oxidation pattern, exclusively on cellulose, and dependent on the presence of the electron donor (Fig. 2a, b). Products were not detected in any of the negative controls (Additional file 1: Figure S1). MALDI-TOF MS analysis of crude extract from activity assays carried out with Cu-loaded TtAA10A in the presence of 10 mM EDTA failed to detect the release of products (data not shown), indicating that, as expected, copper is essential for activity.
Synergy experiments were performed by co-incubating TtAA10A and commercial glycoside hydrolases (GH6 and GH9) in the presence of Avicel and gallic acid, and the resulting mono- and oligosaccharides were quantified using high-performance anion-exchange chromatography (HPAEC). While reactions containing either LPMO or GH alone released negligible amounts of free sugars, co-incubation reactions showed a strong synergistic effect, further enhanced by the presence of the electron donor (Fig. 2c, d, Additional file 2: Figure S2). It is worth noting that both commercial GHs (GH6 and GH9) tested during these experiments belong to families that were identified as among the most abundant in the digestive proteome of shipworms , strengthening the biological relevance of the abovementioned activity assays in the context of wood digestion in the shipworm milieu.
Electron paramagnetic resonance spectroscopy
Our first evidence that some protein preparations contained two Cu sites came from EPR analyses. The frozen solution (165 K) X-band CW-EPR spectrum of Cu-saturated TtAA10A (Fig. 3) exhibited two sets of hyperfine peaks in the parallel region of the spectrum, indicating the presence of two distinct copper coordination geometries, arising either from different coordination environments within a single site (e.g. differences in protonation states of ligands) or a distinct second copper-binding site. Indeed, an accurate simulation of the parallel region of the spectrum could be obtained with two different species, each of which afforded a different set of spin Hamiltonian parameters, gz = 2.267 and |Az| = 425 MHz (species 1), and gz = 2.314 and |Az| = 465 MHz (species 2), Table 1, with a ratio between species 1 and 2 of approximately 3:2. The gz value of species 2 is high compared to what one might expect for the typical AA10 LPMO copper coordination in the active site (spectroscopy of LPMOs recently reviewed in Ref. ), on the basis of which we assign species 1 to a copper bound to the canonical histidine brace active site. Its spin Hamiltonian values are typical of an axial Cu coordination geometry which contains a mixture of N and O-donating ligands . (Note that species 2 cannot be from a free copper species in solution as all small molecule species are removed during the protein preparation; hence, all copper signals in the EPR arise from protein-bound copper.)
To determine whether the two signals arose from a single copper-binding site with different coordination geometries, or from two distinct copper sites, an X-band CW-EPR titration experiment was performed. The protein was pre-treated with EDTA (10× the protein concentration) to remove any copper and then buffer-exchanged to remove any EDTA. This copper-free protein sample was tested and, as expected, showed no copper-based signal. Addition of 0.2 equivalents of copper (compared to the protein concentration) showed a single signal in the parallel region, assigned to the copper (II) ion within the histidine brace active site (species 1). Further additions of copper increased this histidine brace copper signal, with concomitant growth of the signal for species 2, already evident after 0.4 equivalents of copper (Additional file 3: Figure S3). These titration experiments were performed at a fixed pH and show that the two species in the EPR spectrum of copper-saturated TtAA10A represent two different Cu-binding sites with slightly different copper-binding affinities, where species 1 is the higher affinity site. Furthermore, the sample of TtAA10A with 0.4 equivalents of Cu was left at 4 °C for 48 h, and its EPR spectrum re-examined. This sample showed no difference in the ratio of copper species, demonstrating that the different binding sites did not arise because of large differences in the kinetics of copper binding.
Noticeably, in several preparations that we produced, a sample of TtAA10A was isolated which exhibited only a single copper signal in the EPR spectrum. The reasons for this difference in copper stoichiometry of the isolated protein are not clear as these samples were ostensibly prepared using identical conditions as those that afforded TtAA10A with two distinct Cu signals in the X-band EPR spectrum (species 1 and species 2). We were not able to detect any differences in activity for these singly copper occupied preparations relative to previous samples, but we could take advantage of these samples to measure both X-band and Q-band CW-EPR spectra for the active-centre Cu of TtAA10A (Fig. 4) with only the histidine brace occupied, as judged by reference to previous spectra. This sample, therefore, allowed us to perform a simultaneous fit of both the X-band and Q-band spectra to yield more accurate spin Hamiltonian parameters for the copper ion in the histidine brace active site (species 1). These values are reported in Table 2. Addition of PASC to TtAA10A did not cause any change of the EPR spectra (data not shown).
3D structure of TtAA10A
To gain further insights into the molecular basis for the biochemical properties of TtAA10A, and to probe this, potentially unusual, dual Cu structure, we determined the crystal structure for the recombinantly expressed protein to 1.4 Å resolution (Additional file 4: Table S1). The overall structure revealed a core immunoglobulin-like fold decorated by loops and a helical bundle as typically observed for enzymes from this family (Fig. 5). Indeed, structural comparisons using the DALI server  reveal closest structural matches to Cellvibrio japonicus AA10A (PDB ID 5fjq)  and Serratia marcescens CBP21 (PDB ID 2bem)  with RMSDs of 2.4 Å and 2.3 Å over 180 and 170 Cα-positions, respectively, representing only 30% identity at the sequence level. Given the high activity of TtAA10A on cellulose, it may be somewhat surprising that the two closest structural matches to this enzyme are chitin-active AA10s. The third closest structural match, however, was the AA10 from Streptomyces coelicolor (ScAA10, PDB ID 4oy7 ) which is a cellulose-specific AA10 giving an RMSD of 2.5 Å over 160 Cα atoms. TtAA10A and ScAA10 share only 26% sequence identity even though they are active on the same substrate, which further highlights the difficulty in relating LPMO substrate specificity based on sequence and overall structure alone (further discussed in the context of AA9 in Ref. ).
As for all LPMOs studied so far (as defined by their activity), the active site of TtAA10A is formed by the “histidine brace”  motif which is located at the centre of an almost flat surface (Fig. 5a, b). A single copper ion was modelled at this position coordinated in a typical T-shaped geometry by the amino terminus and side chain of His 1 and the side chain imidazole of His 107. In the axial positions around the active site copper ion, TtAA10A has Phe 195 and Gly 105. These positions are often occupied by a phenylalanine/tyrosine and an alanine side chain, respectively, in other AA10s. The latter has been implicated in creating a steric environment which drives the formation of the distorted active site coordination geometry observed in chitin-specific AA10s in their Cu(II) oxidation state [21, 26, 28,29,30]. Here, the replacement of Ala by Gly allows the active site copper to adopt a slightly more axial coordination geometry, closer to the one typical of AA9s and cellulose-active AA10s and consistent with the spin Hamiltonian parameters of species 1 in our EPR analysis described earlier, Fig. 4. LPMO crystal structures are often dogged by photo-reduction as a result of radiation damage such that the resting geometry of the active site cannot be directly observed in crystal structures (examples include [21, 28, 30]). Analysis of the sphere surrounding the copper ion in TtAA10A reveals only weak density for a water molecule 2.6 Å from the copper and stronger density for a second water molecule 3.2 Å away hydrogen bonding to Glu 53. These water molecules are too distant from the copper to be considered as directly coordinating. It, therefore, appears that the copper in this enzyme has also undergone photo-reduction to the Cu(I) oxidation state.
Our structure reveals the position of the second copper-binding site that was observed in the EPR spectra. This site is 14.4 Å (Cu…Cu) from the histidine brace copper ion in a large negatively charged patch on the protein surface (Fig. 5a, b; Additional file 5: Figure S4). This second copper ion is directly coordinated by His 165, Glu 5, Asp 101, a water molecule and His 207* which is provided by the Strep-tag from an adjacent molecule in the crystal. Due to the observation of this interaction, we checked the oligomeric state of the protein in solution using SEC–MALLS (Additional file 6: Figure S5). This confirmed that the protein is monomeric, suggesting that the copper–His207* interaction is a crystal artefact (although one that may hint at potential protein–protein interaction). Nonetheless, added to the EPR data, the structure suggests that this second site is occupied in solution, with His 207* likely replaced by a water molecule. Multiple sequence alignment of the top 500 TtAA10A orthologues identified through BlastP searches shows that, while an acidic residue in position 5 is not uncommon among AA10s, residues 101 and 165 are largely conserved only within LPMOs from bacteria that are closely related to Teredinibacter.
There has been considerable debate about possible positions in which electron donors, both small molecule and proteinaceous, may bind to LPMOs to allow catalysis when the enzyme is bound to the solid substrate surface (see, for example [19, 31,32,33]). Indeed, examination of the TtAA10A structure for potential charge transfer pathways using the programme EHPath shows that a clear and rapid hole-hopping pathway with a mean hole residence time of only 20 ms exists between histidine 1 and tyrosine 3 (10 Å separation). Tyrosine 3 is adjacent (5.3 Å) to the second Cu site, thus providing an efficient charge transfer pathway between the two copper sites . Therefore, given the potential charge transfer pathway between the two copper sites, we investigated whether the second metal site (in our case occupied by copper, although we could not displace the Cu with Fe, Ni, Zn and Mn salts) represents a binding site for a proteinaceous redox partner (the binding of another protein to this site is hinted at by the Strep-tag association with a neighbouring molecule in the crystalline lattice), and we attempted to pull down proteins from the T. turnerae predicted secretome that may stably interact with TtAA10A using affinity column (StrepTrap HP) immobilised TtAA10A. These experiments (data not shown) did not lead to the isolation of any candidate protein-based activators for TtAA10A, but it cannot be ruled out that an activating enzyme could bind transiently in this region to allow electron transfer to the LPMO and hence the initiation of catalysis. It should be noted, however, that during the structure refinement, we also identified a sodium binding site on the protein surface (Additional file 7: Figure S6). Whether these additional binding sites are a result of the increased charge that this protein may have to be stable in the saline environment in which T. turnerae resides remains an open question. Nonetheless, these surface features of TtAA10A may be of interest to enzyme engineers if LPMOs are to be stabilised or adapted to specific conditions to be deployed in industrial bioreactors.
The structural and biochemical characterisation of TtAA10A illuminates the key role of symbiont LPMOs in wood digestion by shipworms. Interestingly, the transcriptome of L. pedicellatus contains putative members of two LPMO families (AA10 and AA15, see “Methods” for more details). Predicted AA10 sequences have best matches in the genomes of marine bacteria and are expressed only in the shipworm gills (as reported in previous literature for this and other shipworm species [13, 14]), while putative AA15s have high similarity to sequences from invertebrate genomes (particularly from molluscs) and feature much lower gene expression levels. Shotgun proteomics analysis of the content of the shipworm digestive system has also shown the presence of mature bacterial AA10 LPMOs but no detectable amounts of endogenous AA15s . On the other hand, AA15 LPMOs have high gene expression and protein abundance in the digestive system of the primitive insect Thermobia domestica, and were shown to play a role in its ability to digest plant biomass without microbial assistance, while no bacterial AA10s were detected . Unlike T. domestica, shipworms appear to have co-opted bacterial AA10 towards cellulose digestion, while the role of endogenous AA15s remains to be determined. Interestingly, the study on T. domestica showed that AA15 genes are widespread among invertebrates, including those that do not feed on plant biomass. As such, their ancestral role is most likely to remodel endogenous chitin , a structural polysaccharide required for the formation of insect exoskeletons and mollusc shells.
The in vitro characterisation of TtAA10A revealed the high specificity of this enzyme towards cellulose, a result that we expected based on the presence of a putative cellulose binding domain (CBM10) at the C-terminus of the protein sequence in T. turnerae. The products generated by TtAA10A, its EPR spectrum and its X-ray structure, however, were unexpected. While most peaks detected through MALDI-TOF MS are consistent with a mixed C1–C4 oxidation on the cellulose backbone, we also detected a product that could be compatible with higher levels of oxidation/prolonged incubation. The significance of this side reaction in TtAA10A is beyond the scope of the current manuscript. However, both X-ray structural features and EPR spectra support the existence of a second copper-binding site (distinct from the conserved histidine brace) that could potentially influence LPMO mechanism of action, electron donation, and thereby the nature of the generated products. The identification of a sodium-binding site on the surface of TtAA10A is also of interest, as it might play a role in stabilising the enzyme and help it cope with the marine environment conditions, thus opening up new potential opportunities in industrial applications requiring high salt. Shipworm symbionts, therefore, offer an exciting environmental niche in which to hunt for LPMOs with unusual properties.
Heterologous gene expression and protein purification
The LPMO domain of TtAA10A (without its CBM) was codon optimised and cloned using the In-Fusion HD cloning kit (Takara, Saint-Germain-en-Laye, France) into a modified pET22b vector containing at the N-terminus the pelB leader sequence to direct protein production to the periplasm, and a C-terminal Strep-tag. The construct was transformed into Tig Chaperone E. coli cells (Pgro7 Chaperone set Takara). LB liquid medium containing ampicillin (100 μg/mL) and chloramphenicol (35 μg/mL) was inoculated with starter culture, and l-arabinose (0.5 g/L final concentration) was added to induce expression of the Tig chaperone. Cultures were grown at 37 °C until the OD was approximately 0.6. The cultures were left to cool slightly before IPTG was added to a final concentration of 1 mM, and they were incubated overnight at 16 °C with shaking.
Cultures were harvested by centrifugation, 5000g for 30 min. For each 100 mL of original culture the pellet was gently resuspended in 5 mL of ice cold 50 mM Tris HCl buffer pH 8 with 20% v/v sucrose and left on ice for 30 min with occasional mixing. The cell suspension was centrifuged again at 8000 rpm for 10 min, and the cells subjected to osmotic shock; pellet was resuspended in ice-cold 1 mM MgSO4 (5 mL per 100 mL of initial culture) plus AEBSF protease inhibitor, and left on ice for 30 min with occasional mixing. The suspension was centrifuged again, and the supernatant collected, filtered and diluted to make up a 1× PBS solution (using stock 10× PBS, pH 7.4).
The sample was loaded onto a strep column (GE Healthcare) pre-equilibrated in 1× PBS pH 7.4. The column was washed with 1× PBS buffer for 5 column volumes, followed by protein elution using 5 column volumes of 1× PBS pH 7.4, 2.5 mM desthiobiotin.
The dilute TtAA10A-strep protein sample was copper loaded by incubation with excess CuSO4 (1 mM final concentration equivalent to 10× the protein concentration) at 4 °C overnight. To remove excess desthiobiotin and unbound copper, the protein sample was filtered, concentrated and passed through a HiLoadTM 16/60 Superdex 75 gel filtration column (Ge Healthcare) equilibrated with PBS pH 7.4.
In vitro activity assays
Activity assays were carried out on microcrystalline cellulose (Avicel), squid pen chitin, shrimp chitin, glucomannan (low viscosity from konjac, Megazyme), pachyman (Megazyme), xyloglucan (from tamarind, Megazyme), lichenan (from icelandic moss, Megazyme), galactan (from lupin, Megazyme), galactomannan (from carob, Megazyme), mannan (borohydride reduced, Megazyme), cellohexaose (Megazyme, used in final concentration of 40 µM), corn starch (Sigma) and beechwood xylan (Serva).
Typical reactions for LPMO characterization through mass spectrometry were carried out by mixing 4 mg/mL substrate with purified TtAA10A (2 μM) and 4 mM gallic acid in a total volume of 100 μL in 2 mL plastic reaction tubes for 24 h. All reactions analysed via MALDI-TOF MS were carried out in 50 mM ammonium acetate buffer pH 6 and incubated at 28 °C shaking at 600 rpm and the supernatant used for analysis. Samples were analysed by MALDI-TOF MS as described in .
Reactions used for product quantification and boosting experiments with TtAA10A were typically carried out in 50 mM sodium phosphate buffer pH 6 in triplicates of 100 μL each for 3 h at 600 rpm at 28 °C. Each reaction contained 2 μM purified LPMO, 4 mg/mL substrate, and 1 mM gallic acid. Commercial GH6 (cat. number E-CBHIIM, Megazyme, 0.8 mU) and GH9 (cat. number CZ03921, NZYTech, 10 µg) were added to 100 μL reactions. After 3-h incubation, 400 μL of ethanol was added to stop the reaction, spun down and 400 μL of supernatant was transferred to new plastic tubes, dried down and re-suspended in 80 μL of pure water, filtered and analysed via HPAEC as previously described .
Electron paramagnetic resonance (EPR) spectroscopy
Frozen solution CW X- and Q-band EPR spectra of TtAA10A were collected at 165 K on a Bruker micro EMX spectrometer operating at ~ 9.3 GHz, modulation amplitude of 4 G and 10.02 mW microwave power or at 77 K on a Jeol JES-X320 spectrometer operating at ~ 34.7 GHz with microwave power of 1 mW, respectively. Protein samples used during single EPR analysis were all in the concentration range 100–200 μM for X-band and 1 mM for Q-band, in 1X PBS buffer at pH 7.4. Simulations of the collected spectra were carried out in Easy Spin 5.2.6  integrated into MatLab 2016a software to determine the g and A-tensor parameters. Accurate determination of the spin Hamiltonian parameters for species 1 was obtained by simultaneous fitting of both X- and Q-band spectra. The superhyperfine coupling values for the nitrogen atoms could not be determined accurately, although it was noted that the fitting was highly improved by addition of two nitrogen atoms with coupling in the range 35–40 MHz. Accurate determination of the gx, gy, |Ax| and |Ay| values for species 2 was not possible due to the overlap of the two species in the perpendicular region. EPR of different batches of protein showed that some samples contained two copper species, whilst other samples contained only one. EPR Copper titrations (CuSO4 1 M) were carried out on protein free of copper (10 mM EDTA treated, followed by extensive buffer exchange), with spectra taken before the titration was started confirming the lack of coordinated copper in the protein. Copper was added to the protein solution whilst contained within the EPR tubes, in additions of 0.2 equivalents (to the concentration of the protein in the sample). Spectra were measured after each addition of copper solution. Raw EPR data is available on request through Research Data York with DOI: https://doi.org/10.15124/4d239806-7529-45fc-95a1-c14bcadc24d2.
TtAA10A crystallisation, structure solution and refinement
TtAA10A was screened for crystallisation at 7 mg/mL using the Hampton HT (Hampton Research) and PEG/Ion (Qiagen) screens in 96 well sitting drops using Mosquito robotics (TTP Labtech). Initial hits were obtained in the PEG/Ion screen condition B8 (0.2 M magnesium formate dehydrate and 20% w/v PEG 3350, pH 7), which were subsequently optimised in hanging drops, screening across a narrow range of concentrations for both magnesium formate dehydrate and PEG 3350. A single crystal was cryo-cooled for data collection by plunging in liquid nitrogen without the addition of cryo-protectant.
X-ray diffraction data were collected at Diamond Light Source, beamline I04 and were processed in the CCP4i2 pipeline  to 1.4 Å resolution (see Additional file 2: Table S1, for processing statistics). The structure was determined by single wavelength anomalous dispersion (SAD) in CRANK2  using the anomalous signal from protein-bound copper ions. Briefly, SHELX  was used to solve the phase problem and to generate an initial structural model which was further elaborated upon in BUCCANEER . The final model was generated using iterative cycles of rebuilding and refinement in COOT  and REFMAC , respectively (Additional file 2: Table S1). The final model and accompanying structure factors have been deposited in the Protein Data Bank with accession code 6RW7.
Availability of data and materials
Coordinates and observed data have been deposited in the PDB with accession code 6RW7.
degree of polymerisation
electron paramagnetic resonance
high-pressure anion-exchange chromatography
high-pressure liquid chromatography
(lytic) polysaccharide monooxygenase
- MALDI (TOF):
matrix-assisted laser desorption/ionisation (time of flight)
phosphoric acid swollen cellulose
single-wavelength anomalous dispersion
Vaaje-Kolstad G, Westereng B, Horn SJ, Liu Z, Zhai H, Sørlie M, Eijsink VGH. An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Science. 2010;330(6001):219–22.
Quinlan RJ, Sweeney MD, Lo Leggio L, Otten H, Poulsen J-CN, Johansen KS, Krogh KBRM, Jørgensen CI, Tovborg M, Anthonsen A, et al. Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. Proc Natl Acad Sci USA. 2011;108(37):15079–84.
Beeson WT, Phillips CM, Cate JH, Marletta MA. Oxidative cleavage of cellulose by fungal copper-dependent polysaccharide monooxygenases. J Am Chem Soc. 2012;134(2):890–2.
Hemsworth GR, Déjean G, Davies GJ, Brumer H. Learning from microbial strategies for polysaccharide degradation. Biochem Soc Trans. 2016;44(1):94–108.
Hemsworth GR, Johnston EM, Davies GJ, Walton PH. Lytic polysaccharide monooxygenases in biomass conversion. Trends Biotechnol. 2015;33(12):747–61.
Ciano L, Davies GJ, Tolman WB, Walton PH. Bracing copper for the catalytic oxidation of C–H bonds. Nat Catal. 2018;1(8):571–7.
Beeson WT, Vu VV, Span EA, Phillips CM, Marletta MA. Cellulose degradation by polysaccharide monooxygenases. Annu Rev Biochem. 2015;84:923–46.
Horn SJ, Vaaje-Kolstad G, Westereng B, Eijsink VG. Novel enzymes for the degradation of cellulose. Biotechnol Biofuels. 2012;5(1):45.
Harris PV, Welner D, McFarland KC, Re E, Navarro Poulsen JC, Brown K, Salbo R, Ding H, Vlasenko E, Merino S, et al. Stimulation of lignocellulosic biomass hydrolysis by proteins of glycoside hydrolase family 61: structure and function of a large, enigmatic family. Biochemistry. 2010;49(15):3305–16.
Cragg SM, Beckham GT, Bruce NC, Bugg TDH, Distel DL, Dupree P, Etxabe AG, Goodell BS, Jellison J, McGeehan JE, et al. Lignocellulose degradation mechanisms across the Tree of Life. Curr Opin Chem Biol. 2015;29:108–19.
Distel DL. The biology of marine wood boring bivalves and their bacterial endosymbionts. Acs Sym Ser. 2003;845:253–71.
Distel DL, Beaudoin DJ, Morrill W. Coexistence of multiple proteobacterial endosymbionts in the gills of the wood-boring bivalve Lyrodus pedicellatus (Bivalvia: Teredinidae). Appl Environ Microb. 2002;68(12):6292–9.
Sabbadin F, Pesante G, Elias L, Besser K, Li Y, Steele-King C, Stark M, Rathbone DA, Dowle AA, Bates R, et al. Uncovering the molecular mechanisms of lignocellulose digestion in shipworms. Biotechnol Biofuels. 2018;11:59.
O’Connor RM, Fung JM, Sharp KH, Benner JS, McClung C, Cushing S, Lamkin ER, Fomenkov AI, Henrissat B, Londer YY, et al. Gill bacteria enable a novel digestive strategy in a wood-feeding mollusk. Proc Natl Acad Sci USA. 2014;111(47):E5096–104.
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42(1):D490–5.
Distel DL, Morrill W, MacLaren-Toussaint N, Franks D, Waterbury J. Teredinibacter turnerae gen. nov., sp. nov., a dinitrogen-fixing, cellulolytic, endosymbiotic gamma-proteobacterium isolated from the gills of wood-boring molluscs (Bivalvia: Teredinidae). Int J Syst Evol Microbiol. 2002;52(Pt 6):2261–9.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Millward-Sadler SJ, Davidson K, Hazlewood GP, Black GW, Gilbert HJ, Clarke JH. Novel cellulose-binding domains, NodB homologues and conserved modular architecture in xylanases from the aerobic soil bacteria Pseudomonas fluorescens subsp. cellulosa and Cellvibrio mixtus. Biochem J. 1995;312(Pt 1):39–48.
Book AJ, Yennamalli RM, Takasuka TE, Currie CR, Phillips GN Jr, Fox BG. Evolution of substrate specificity in bacterial AA10 lytic polysaccharide monooxygenases. Biotechnol Biofuels. 2014;7:109.
Sabbadin F, Hemsworth GR, Ciano L, Henrissat B, Dupree P, Tryfona T, Marques RDS, Sweeney ST, Besser K, Elias L, et al. An ancient family of lytic polysaccharide monooxygenases with roles in arthropod development and biomass digestion. Nat Commun. 2018;9(1):756.
Glyn RH, Robbert QK, Rebecca CG, Sally JL, Johan PT, Alison JP, Gideon JD, Paul HW. The copper active site of CBM33 polysaccharide oxygenases. J Am Chem Soc. 2013;135(16):6069–77.
Peisach J, Blumberg WE. Structural implications derived from the analysis of electron paramagnetic resonance spectra of natural and artificial copper proteins. Arch Biochem Biophys. 1974;165(2):691–708.
Holm L, Laakso LM. Dali server update. Nucleic Acids Res. 2016;44(W1):W351–5.
Forsberg Z, Nelson CE, Dalhus B, Mekasha S, Loose JS, Crouch LI, Rohr AK, Gardner JG, Eijsink VG, Vaaje-Kolstad G. Structural and functional analysis of a lytic polysaccharide monooxygenase important for efficient utilization of chitin in Cellvibrio japonicus. J Biol Chem. 2016;291(14):7300–12.
Vaaje-Kolstad G, Houston DR, Riemen AHK, Eijsink VGH, van Aalten DMF. Crystal structure and binding properties of the Serratia marcescens chitin-binding protein CBP21. J Biol Chem. 2005;280(12):11313–9.
Forsberg Z, Mackenzie AK, Sorlie M, Rohr AK, Helland R, Arvai AS, Vaaje-Kolstad G, Eijsink VG. Structural and functional characterization of a conserved pair of bacterial cellulose-oxidizing lytic polysaccharide monooxygenases. Proc Natl Acad Sci USA. 2014;111(23):8446–51.
Lenfant N, Hainaut M, Terrapon N, Drula E, Lombard V, Henrissat B. A bioinformatics analysis of 3400 lytic polysaccharide oxidases from family AA9. Carbohydr Res. 2017;448:166–74.
Gudmundsson M, Kim S, Wu M, Ishida T, Momeni MH, Vaaje-Kolstad G, Lundberg D, Royant A, Ståhlberg J, Eijsink VGH, et al. Structural and electronic snapshots during the transition from a Cu(II) to Cu(I) metal center of a lytic polysaccharide monooxygenase by X-ray photoreduction. J Biol Chem. 2014;289(27):18782–92.
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(10):1647–56.
Gregory RC, Hemsworth GR, Turkenburg JP, Hart SJ, Walton PH, Davies GJ. Activity, stability and 3-D structure of the Cu(ii) form of a chitin-active lytic polysaccharide monooxygenase from Bacillus amyloliquefaciens. Dalton Trans. 2016;45(42):16904–12.
Li X, Beeson William T, Phillips Christopher M, Marletta Michael A, Cate Jamie HD. Structural basis for substrate targeting and catalysis by fungal polysaccharide monooxygenases. Structure. 2012;20(6):1051–61.
Tan TC, Kracher D, Gandini R, Sygmund C, Kittl R, Haltrich D, Hallberg BM, Ludwig R, Divne C. Structural basis for cellobiose dehydrogenase action during oxidative cellulose degradation. Nat Commun. 2015;6:7542.
Courtade G, Wimmer R, Rohr AK, Preims M, Felice AK, Dimarogona M, Vaaje-Kolstad G, Sorlie M, Sandgren M, Ludwig R, et al. Interactions of a fungal lytic polysaccharide monooxygenase with beta-glucan substrates and cellobiose dehydrogenase. Proc Natl Acad Sci USA. 2016;113(21):5922–7.
Teo RD, Wang R, Smithwick ER, Migliore A, Therien MJ, Beratan DN. Mapping hole hopping escape routes in proteins. Proc Natl Acad Sci USA. 2019;116(32):15811–6.
Stoll S, Schweiger A. EasySpin, a comprehensive software package for spectral simulation and analysis in EPR. J Mag Res (San Diego, Calif: 1997). 2006;178(1):42–55.
Potterton L, Agirre J, Ballard C, Cowtan K, Dodson E, Evans PR, Jenkins HT, Keegan R, Krissinel E, Stevenson K, et al. CCP4i2: the new graphical user interface to the CCP4 program suite. Acta Crystallogr D Struct Biol. 2018;74(Pt 2):68–84.
Ness SR, de Graaff RA, Abrahams JP, Pannu NS. CRANK: new methods for automated macromolecular crystal structure solution. Structure. 2004;12(10):1753–61.
Sheldrick GM. A short history of SHELX. Acta Crystallogr A. 2008;64(Pt 1):112–22.
Cowtan K. The Buccaneer software for automated model building. 1. Tracing protein chains. Acta Crystallogr D Biol Crystallogr. 2006;62(Pt 9):1002–11.
Emsley P, Cowtan K. Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr. 2004;60(Pt 12 Pt 1):2126–32.
Murshudov GN, Vagin AA, Dodson EJ. Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D Biol Crystallogr. 1997;53(Pt 3):240–55.
Zhang H, Yohe T, Huang L, Entwistle S, Wu P, Yang Z, Busk PK, Xu Y, Yin Y. dbCAN2: a meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2018;46(W1):W95–101.
Lo Leggio L, Simmons TJ, Poulsen JC, Frandsen KE, Hemsworth GR, Stringer MA, von Freiesleben P, Tovborg M, Johansen KS, De Maria L, et al. Structure and boosting activity of a starch-degrading lytic polysaccharide monooxygenase. Nat Commun. 2015;6:5961.
We thank Diamond Light Source for access to beamline i04 (proposal number mx-13587) that contributed to the results presented here.
This work was supported by the Biotechnology and Biological Sciences Research Council (BBSRC through Grants BB/R007705/1 and BB/L001926/1). GJD is supported by the Royal Society through the Ken Murray Research Professorship.
Ethics approval and consent to participate
Consent for publication
All authors have seen the manuscript and consent to publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1. MALDI-TOF MS analysis of in vitro negative control activity assays with purified TtAA10A, under the same experimental conditions as in Fig. 2a. The panels show spectra of products obtained after incubation of 4 mg/mL Avicel (a), Avicel with 4 mM gallic acid (b) and Avicel with 2 µM TtAA10A (c). The spectra show no detectable amounts of native or oxidised cello-oligosaccharides. Relative intensity represents 1.23 × 103.
Additional file 2: Figure S2. HPAEC chromatograms showing the release of cellobiose from Avicel during boosting experiments with TtAA10A, commercial GH6 and gallic acid. The identity and concentration of cellobiose (retention time approx. 4.5 mins) was determined by analysis of a commercial standard (Std). The small amount of additional cellobiose released with LPMO in the absence of external reducing substrate is assumed to be derived from low quantities of unknown reductant, as has been observed on other systems such as AA13  and AA15 .
Additional file 3: Figure S3. EPR Cu titration reveals the presence of two distinct Cu sites. EPR copper titration experiment in which the protein was pre-treated with EDTA to remove the copper (as shown by the lack of EPR signal). Addition of 0.2 equivalents of Cu (to the concentration of protein) produced a single Cu species indicative of the histidine brace. Addition of a further 0.2 equivalents caused a change in the spectra with appearance of a second set of hyperfine peaks in the parallel region. The clearest change is seen after 0.6 equivalents of Cu have been added, where multiple peaks in the spectra indicate two Cu species. Loss of resolution in the last two spectra may be due to signal dilution.
Additional file 5: Figure S4. Electrostatic Surface Potential for TtAA10A. The electrostatic charge distribution for TtAA10A has been mapped onto the protein surface using the APBS plug-in for PyMol at ± 10 KBT/e. The histidine brace, secondary copper binding site and sodium binding site are outlined with black lines showing the negatively charged areas to which the ions bind.
Additional file 6: Figure S5. SEC-MALLS analysis of TtAA10A, where the solid red line indicates the refractive index, the dashed line is the light scattering (these are essentially identical so indistinguishable), dotted line the UV response at 280 nm. The central line inside the peak is representative of the molar mass (2.44 × 104 Da), indicating that TtAA10A forms a monomer in solution.
Additional file 7: Figure S6. The Na+ Site on the TtAA10A Surface. The modelled sodium ion is shown as a grey sphere with the groups that coordinate it shown as sticks coloured by atom type. The 2Fobs–Fcalc map is shown as a blue wire mesh contoured at 1σ. The sodium ion is octahedrally coordinated with three main chain carbonyl groups and three water molecules. Sodium was assigned at this site based on the coordinating bond lengths and this was the ion that gave a B-factor following refinement closest to the surrounding protein atoms.
About this article
Cite this article
Fowler, C.A., Sabbadin, F., Ciano, L. et al. Discovery, activity and characterisation of an AA10 lytic polysaccharide oxygenase from the shipworm symbiont Teredinibacter turnerae. Biotechnol Biofuels 12, 232 (2019). https://doi.org/10.1186/s13068-019-1573-x
- 3D structure