Advertisement
Microbiology| Volume 293, ISSUE 29, P11564-11573, July 20, 2018

Structural and functional characterization of PL28 family ulvan lyase NLR48 from Nonlabens ulvanivorans

  • ThirumalaiSelvi Ulaganathan
    Affiliations
    Department of Biochemistry, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E5, Canada
    Search for articles by this author
  • Ehud Banin
    Affiliations
    Faculty of Life Sciences, Bar Ilan University, Ramat-Gan 5290002, Israel
    Search for articles by this author
  • William Helbert
    Affiliations
    Centre de Recherches sur les Macromolécules Végétales (CERMAV), CNRS and Grenoble Alpes Université, BP53, 38000 Grenoble Cedex 9, France
    Search for articles by this author
  • Miroslaw Cygler
    Correspondence
    To whom correspondence should be addressed: Dept. of Biochemistry, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E5, Canada. Tel.:306-966-4361;.
    Affiliations
    Department of Biochemistry, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E5, Canada

    Department of Biochemistry, McGill University, Montreal, Quebec H3G 0B1, Canada
    Search for articles by this author
  • Author Footnotes
    3 Please note that the JBC is not responsible for the long-term archiving and maintenance of this site or any other third party-hosted site.
    2 The abbreviations used are: GlcAglucuronic acidIdoAiduronic acidPULpolysaccharide utilization lociPLpolysaccharide lyaseRha3S3-O-sulfated rhamnoseMBPmaltose-binding proteinr.m.s.d.root-mean-square deviationΔUA4-deoxy-l-threo-hex-4-enopyranosiduronic acidPDBProtein Data BankTEVtobacco etch virus.
Open AccessPublished:June 06, 2018DOI:https://doi.org/10.1074/jbc.RA118.003659
      Ulvan is a complex sulfated polysaccharide present in the cell wall of green algae of the genus Ulva (Chlorophyta). The first ulvan-degrading polysaccharide lyases were identified several years ago, and more were discovered through genome sequencing of marine bacteria. Ulvan lyases are now grouped in three polysaccharide lyase (PL) families in the CAZy database, PL24, PL25, and PL28. The recently determined structures of the representative lyases from families PL24 and PL25 show that they adopt a seven-bladed β-propeller fold and utilize the His/Tyr catalytic mechanism. No structural information is yet available for PL28 ulvan lyases. NLR48 from Nonlabens ulvanivorans belongs to PL28 together with its close paralog, NLR42. Biochemical studies of NLR42 have revealed that it can cleave ulvan next to both uronic acid epimers. We report the crystal structure of ulvan lyase NLR48 at 1.9-Å resolution. It has a β-jelly roll fold with an extended, deep, and positively charged substrate-binding cleft. Putative active-site residues were identified from the sequence conservation pattern, and their role was confirmed by site-directed mutagenesis. The structure of an inactive K162M mutant with a tetrasaccharide substrate showed the substrate occupying the “−” subsites. Comparison with lyases from other PL families with β-jelly roll folds supported assignment of the active site and explained its ability to degrade ulvan next to either epimer of uronic acid. NLR48 contains the His/Tyr catalytic machinery with Lys162 and Tyr281 playing the catalytic base/acid roles.

      Introduction

      Red, brown, and green algae are major members of the marine ecosystem and contribute to the huge aquatic biomass (
      • Shalaby E.A.
      Algae as promising organisms for environment and health.
      ). Polysaccharides present in the algal cell wall play a structural role, providing rigidity, and represent a major component of their biomass. They are also a rich source of carbon for cohabiting microorganisms. Carrageenan is the main cell wall polysaccharide of many red algae, alginate is a cell wall constituent of brown algae, and ulvan is present in the green algae belonging to the genus Ulva (Chlorophyta). Carrageenan and alginate have been studied extensively and reported to have many industrial applications (
      • Sudha P.N.
      • Aisverya S.
      • Nithya R.
      • Vijayalakshmi K.
      Industrial applications of marine carbohydrates.
      ), and the enzymes degrading these polysaccharides have been well-characterized (
      • Hehemann J.-H.
      • Boraston A.B.
      • Czjzek M.
      A sweet new wave: structures and mechanisms of enzymes that digest polysaccharides from marine algae.
      ,
      • Michel G.
      • Czjzek M.
      Polysaccharide-degrading enzymes from marine bacteria.
      ). Ulvan is biosynthesized by marine green algae Ulva spp. and Enteromorpha spp. and contributes to 8–30% of the dry weight of the algae. This sulfated polysaccharide has recently gained attention for a variety of industrial applications in the agriculture, food, pharmaceutical, chemical, and biomaterial industries (
      • Shalaby E.A.
      Algae as promising organisms for environment and health.
      ,
      • Cunha L.
      • Grenha A.
      Sulfated seaweed polysaccharides as multifunctional materials in drug delivery applications.
      ,
      • Popa E.G.
      • Reis R.L.
      • Gomes M.E.
      Seaweed polysaccharide-based hydrogels used for the regeneration of articular cartilage.
      ).
      Major components of the ulvan polysaccharide are 3-O-sulfated α-l-rhamnose, α-l-iduronic acid, β-d-glucuronic acid, and xylose. Ulvan is made of two main repetitive disaccharide units. Type A ulvanobiourinic 3-sulfate has [→4)-β-d-GlcA-(1→4)-α-l-Rha3S-(1→]
      The abbreviations used are: GlcA
      glucuronic acid
      IdoA
      iduronic acid
      PUL
      polysaccharide utilization loci
      PL
      polysaccharide lyase
      Rha3S
      3-O-sulfated rhamnose
      MBP
      maltose-binding protein
      r.m.s.d.
      root-mean-square deviation
      ΔUA
      4-deoxy-l-threo-hex-4-enopyranosiduronic acid
      PDB
      Protein Data Bank
      TEV
      tobacco etch virus.
      and type B ulvanobiouronic 3-sulfate has [→4)-α-l-IdoA-(1→4)-α-l-Rha3S(1→] disaccharide units (Scheme 1) (
      • Lahaye M.
      • Robic A.
      Structure and functional properties of ulvan, a polysaccharide from green seaweeds.
      ). The occurrence of a sulfate group in ulvan resembles the glycosaminoglycans present in vertebrates, and the presence of uronic acid makes it subject to depolymerization by a lytic mechanism.
      Figure thumbnail sc1
      Scheme 1Diagram of a tetrasaccharide fragment of ulvan containing either GlcA or IdoA. The product of the cleavage by lyase is shown on the right.
      Ulvan lyases degrade ulvan through a β-elimination mechanism, cleaving the 1–4-glycosidic linkage on the nonreducing end of a uronic acid residue with the formation of an unsaturated product at the newly formed nonreducing end of the polysaccharide. The first identified ulvan lyase and the first representative of the PL28 family was NLR42 (NCBI accession number AEN28574; 46 kDa with two domains) isolated from Nonlabens ulvanivorans (synonym, Persicivirga ulvanivorans) (
      • Nyvall Collén P.
      • Sassi J.F.
      • Rogniaux H.
      • Marfaing H.
      • Helbert W.
      Ulvan lyases isolated from the flavobacteria Persicivirga ulvanivorans are the first members of a new polysaccharide lyase family.
      ). The enzyme was able to cleave ulvan next to both GlcA and IdoA. In addition to the catalytic domain, NLR42 also harbors carbohydrate-binding and sorting domains (
      • Melcher R.L.
      • Neumann M.
      • Fuenzalida Werner J.P.
      • Gröhn F.
      • Moerschbacher B.M.
      Revised domain structure of ulvan lyase and characterization of the first ulvan binding domain.
      ).
      Genome sequencing of marine microorganisms undertaken in the last few years led to the identification of more genes exhibiting ulvan lyase activity (
      • Kopel M.
      • Helbert W.
      • Henrissat B.
      • Doniger T.
      • Banin E.
      Draft genome sequence of Pseudoalteromonas sp. strain PLSV, an ulvan-degrading bacterium.
      ,
      • Kopel M.
      • Helbert W.
      • Henrissat B.
      • Doniger T.
      • Banin E.
      Draft genome sequences of two ulvan-degrading isolates, strains LTR and LOR, that belong to the Alteromonas genus.
      ,
      • Kopel M.
      • Helbert W.
      • Henrissat B.
      • Doniger T.
      • Banin E.
      Draft genome sequence of Nonlabens ulvanivorans, an ulvan-degrading bacterium.
      ). Polysaccharide-degrading machineries occur often as coregulated gene clusters called polysaccharide utilization loci (PUL). The PUL regions include genes coding for various transport proteins, highly specific carbohydrate-acting enzymes (glycoside hydrolases, polysaccharide lyases, and sulfatases), and recognition proteins. The complete degradation of the polysaccharides requires the coordination of these different components. We had previously characterized the ulvan PUL from Alteromonas and Pseudoalteromonas species (
      • Foran E.
      • Buravenkov V.
      • Kopel M.
      • Mizrahi N.
      • Shoshani S.
      • Helbert W.
      • Banin E.
      Functional characterization of a novel “ulvan utilization loci” found in Alteromonas sp. LOR genome.
      ), which led to addition of the families PL24 (LOR107) and PL25 (PLSV3936) to the CAZy database (http://www.cazy.org)
      Please note that the JBC is not responsible for the long-term archiving and maintenance of this site or any other third party-hosted site.
      (
      • Lombard V.
      • Golaconda Ramulu H.
      • Drula E.
      • Coutinho P.M.
      • Henrissat B.
      The carbohydrate-active enzymes database (CAZy) in 2013.
      ). Both LOR107 and PLSV3936 display <14% sequence identity; nevertheless, they both adopt a β-propeller scaffold. LOR107 is reported to degrade ulvan specifically after the IdoA site (
      • Kopel M.
      • Helbert W.
      • Belnik Y.
      • Buravenkov V.
      • Herman A.
      • Banin E.
      New family of ulvan lyases identified in three isolates from the Alteromonadales order.
      ), whereas PLSV3936 can cleave after both GlcA and IdoA sites (
      • Ulaganathan T.
      • Boniecki M.T.
      • Foran E.
      • Buravenkov V.
      • Mizrachi N.
      • Banin E.
      • Helbert W.
      • Cygler M.
      New ulvan-degrading polysaccharide lyase family: structure and catalytic mechanism suggests convergent evolution of active site architecture.
      ). Recently, a new ulvan polysaccharide family, PL28, was added to the CAZy database.
      Here, we report first crystal structure and catalytic mechanism of ulvan lyase NLR48 from N. ulvanivorans grouped in the PL28 family. NLR48 is a 30-kDa protein that was first identified as an ulvan lyase together with its larger homolog, NLR42 (
      • Nyvall Collén P.
      • Sassi J.F.
      • Rogniaux H.
      • Marfaing H.
      • Helbert W.
      Ulvan lyases isolated from the flavobacteria Persicivirga ulvanivorans are the first members of a new polysaccharide lyase family.
      ). NLR48 shows 37% sequence identity with the catalytic module of NLR42, suggesting structural and functional similarity of these enzymes. NLR42 can cleave the ulvan polysaccharide next to both GlcA and IdoA (
      • Nyvall Collén P.
      • Sassi J.F.
      • Rogniaux H.
      • Marfaing H.
      • Helbert W.
      Ulvan lyases isolated from the flavobacteria Persicivirga ulvanivorans are the first members of a new polysaccharide lyase family.
      ). We have determined the crystal structure of NLR48 at 1.9-Å resolution. Unlike the previously characterized ulvan lyases LOR107 and PLSV3936, NLR48 has a β-jelly roll fold. The structure of the NLR active-site mutant K162M in complex with a tetrasaccharide substrate facilitated the identification of active-site residues and allowed us to propose a catalytic mechanism for ulvan degradation.

      Results and discussion

      We have determined the three-dimensional structure of NLR48, the first representative of the polysaccharide lyase family PL28, which has recently been added to the CAZy database. Several constructs were cloned and tested for expression. Expression yields of His-NLR48 and a shorter construct starting at Asp53, called His-NLR48(53D), were high, but the proteins aggregated in solution and were not suitable for crystallization. Therefore, the MBP and GST fusions were tested. The expression level of His-MBP-NLR48 was low in BL21(DE3), but that of the GST-NLR48 fusion was moderate, and no aggregation was observed. It was essential to purify the latter construct at 4 °C to prevent aggregation. This protein yielded well-diffracting crystals.

      Crystal structure of NLR48

      The three-dimensional crystal structure of NLR48 was solved by a single anomalous dispersion method using selenomethionine (SeMet)-incorporated protein crystals. There are two molecules in the asymmetric unit, which are related by a noncrystallographic 2-fold rotation. Both molecules are nearly identical with a root-mean-square deviation (r.m.s.d.) of 0.2 Å for all Cα atoms. NLR48 displays a β-sandwich jelly roll fold with two concave β-sheets stacking on each other (Fig. 1A). Each β-sheet has seven antiparallel, strongly bent β-strands. The inner β-sheet is formed from strands β1-β4-β13-β6-β7-β8-β9, whereas the outer β-sheet contains strands β2-β3-β14-β5-β10-β11-β12 (Fig. 1B). There are seven crossovers between the two sheets. The N-terminal segment, amino acids 21–50, is devoid of secondary structure elements and is fixed to the loop connecting strands β1 (inner sheet) and β2 (outer sheet) by a disulfide bond between Cys32 and Cys59. The inner β-sheet forms a deep cleft extending across the entire β-sheet in a direction nearly perpendicular to the strands. The bottom of the cleft is formed by the centers of the β-stands; one wall is formed from the loops connecting strands β3 to β4 and β5 to β6 (intersheet), whereas the other wall is made from the loops connecting β6 to β7 and β8 to β9 (intra-inner sheet). The cleft is lined predominantly with basic and polar side chains and displays a positive electrostatic surface potential as calculated with APBS software (
      • Jurrus E.
      • Engel D.
      • Star K.
      • Monson K.
      • Brandi J.
      • Felberg L.E.
      • Brookes D.H.
      • Wilson L.
      • Chen J.
      • Liles K.
      • Chun M.
      • Li P.
      • Gohara D.W.
      • Dolinsky T.
      • Konecny R.
      • et al.
      Improvements to the APBS biomolecular solvation software suite.
      ) (Fig. 1C). Two phosphate molecules acquired from the crystallization solution are found in equivalent positions in both molecules. This cleft is likely the substrate-binding site of NLR48. A large positive difference map peak was identified in each molecule near the disulfide bond described above. Based on the environment and the coordination geometry, these peaks are assigned as Ca2+ ions. Each ion is heptacoordinated in a pentagonal bipyramid arrangement by Oδ1Asp79, Oδ2Asp79, Oδ1Asn85, O=CGly59, and O=CSer81 in the equatorial plane and Oδ1Asn61 and O=CAla84 in axial positions. The Ca–O distances are in the 2.25–2.50-Å range (Fig. 2A). This Ca2+ ion provides additional stability to the N-terminal region of NLR48. This Ca2+-binding site is separated from the positively charged cleft, suggesting that Ca2+ plays a structural rather than catalytic role in NLR48.
      Figure thumbnail gr1
      Figure 1Structural details of NLR48. A, stereoview of the cartoon representation of the structure. The molecule is rainbow-colored from blue at the N terminus to red at the C terminus. The β-strands are numbered sequentially along the protein sequence. This and other figures were made with PyMOL (Schrödinger, LLC). B, topology diagram showing the connectivity of the strands. The strands are similarly rainbow-colored as in the cartoon figure. C, view of the electrostatic potential projected on the protein surface within the cleft in the concave side of the protein. The cleft displays strongly positive potential matching the negatively charged polysaccharide substrate. D, residue conservation among sequence homologs of NLR48 as calculated by ConSurf (
      • Landau M.
      • Mayrose I.
      • Rosenberg Y.
      • Glaser F.
      • Martz E.
      • Pupko T.
      • Ben-Tal N.
      ConSurf 2005: the projection of evolutionary conservation scores of residues on protein structures.
      ). The highly conserved residues are colored burgundy. The surface representation shows their concentration within the canyon.
      Figure thumbnail gr2
      Figure 2Ion-binding sites in the native structure of NLR48. A, the coordination sphere of the Ca2+ ion, which plays a structural role in this enzyme. B, the environment of the PO43− ions.
      Two phosphate (PO43−) ions are located within the positively charged cleft and make many contacts with the surrounding residues (Fig. 2B). PO43− at site 1 binds in the middle of the cleft and forms hydrogen bonds with NζLys162, Nϵ2Gln160, Oϵ1Gln160, Nη1Arg117, and OηTyr281. The PO43− at site 2 makes H-bonds with Nη2Arg216, OγSer171 and backbone NHAla175. The phosphate in site 1 has lower B-factors, suggesting a somewhat stronger binding than in site 2.

      Putative substrate-binding and catalytic sites

      The β-sandwich jelly roll fold is not unique within the polysaccharide lyases. It is also observed in alginate lyase families PL7, PL14, and PL18; in heparin lyase family PL13; and in exoglucuronase family PL20 (
      • Garron M.L.
      • Cygler M.
      Structural and mechanistic classification of uronic acid-containing polysaccharide lyases.
      ). Representative structures from these PL families were compared with NLR48. Representative structures for these families were: PL7, alginate lyase (A1-II′;Aly) complexed with tetrasaccharide, Protein Data Bank (PDB) code 2ZAA (
      • Ogura K.
      • Yamasaki M.
      • Mikami B.
      • Hashimoto W.
      • Murata K.
      Substrate recognition by family 7 alginate lyase from Sphingomonas sp. A1.
      ); PL13, heparin lyase I with octasaccharide, Protein Data Bank code 3INA (
      • Han Y.H.
      • Garron M.L.
      • Kim H.Y.
      • Kim W.S.
      • Zhang Z.
      • Ryu K.S.
      • Shaya D.
      • Xiao Z.
      • Cheong C.
      • Kim Y.S.
      • Linhardt R.J.
      • Jeon Y.H.
      • Cygler M.
      Structural snapshots of heparin depolymerization by heparin lyase I.
      ); PL18, alginate lyase aly-SJ02, PDB code 4Q8K (
      • Dong S.
      • Wei T.-D.
      • Chen X.-L.
      • Li C.-Y.
      • Wang P.
      • Xie B.-B.
      • Qin Q.-L.
      • Zhang X.-Y.
      • Pang X.-H.
      • Zhou B.-C.
      • Zhang Y.-Z.
      Molecular insight into the role of the N-terminal extension in the maturation, substrate recognition, and catalysis of a bacterial alginate lyase from polysaccharide lyase family 18.
      ); PL20, glucuronan lyase, PDB code 2ZZJ (
      • Konno N.
      • Ishida T.
      • Igarashi K.
      • Fushinobu S.
      • Habu N.
      • Samejima M.
      • Isogai A.
      Crystal structure of polysaccharide lyase family 20 endo-β-1,4-glucuronan lyase from the filamentous fungus Trichoderma reesei.
      ); and PL14, Chlorella virus enzyme vAL-1, PDB code 3A0N (
      • Ogura K.
      • Yamasaki M.
      • Yamada T.
      • Mikami B.
      • Hashimoto W.
      • Murata K.
      Crystal structure of family 14 polysaccharide lyase with pH-dependent modes of action.
      ). The number of β-strands in the inner and outer sheets varies between these enzymes, but the core structures (∼100 residues; Fig. 3B) superpose well with that of NLR48 (Table 1). They differ mainly in the loop regions (Fig. 3A). The PL14 lyase is the least similar and is excluded from further analysis. Despite sharing the same fold, the sequence identity between NLR48 and these other lyases is less than 14%. Nevertheless, despite low sequence identity, these enzymes have similar active sites and a shared catalytic mechanism belonging to the His/Tyr class (
      • Garron M.L.
      • Cygler M.
      Structural and mechanistic classification of uronic acid-containing polysaccharide lyases.
      ). The active site is located in the cleft extending along the inner β-sheet.
      Figure thumbnail gr3
      Figure 3Stereoview of the superposition of polysaccharide lyases that display the β-sandwich jelly roll fold. A, representative structures from families PL7, PL13, PL14, PL18, and PL20 were superimposed on NLR48 using the Swiss-PdbViewer (
      • Guex N.
      • Peitsch M.C.
      SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling.
      ). NLR48 is colored green, heparinase I (PL13) is colored magenta, alginate lyase (PL7) is colored cyan, PL18 is colored forest green, and PL20 is colored gray. B, the core ∼100 residues of NLR48 that superimpose well among the PL lyases are colored red, and the less structurally conserved segments are colored green. Residues Arg117, Gln160, Lys162, Arg216, and Tyr281 that were mutagenized are shown explicitly in ball-and-stick representation.
      Table 1Superposition of polysaccharide lyase with a β-sandwich jelly roll fold on NLR48
      PDB codePL familyNo. of overlapping Cαr.m.s.d.
      Å
      4Q8KPL181001.47
      2ZZJPL20971.55
      3INAPL13931.60
      2ZAAPL171251.49
      3A0NPL14831.84
      To identify the putative catalytic site, we first utilized ConSurf software (
      • Landau M.
      • Mayrose I.
      • Rosenberg Y.
      • Glaser F.
      • Martz E.
      • Pupko T.
      • Ben-Tal N.
      ConSurf 2005: the projection of evolutionary conservation scores of residues on protein structures.
      ), which showed the highly conserved residues among ∼30 NLR48 homologs (Fig. 1D). The solvent-accessible and conserved residues located within the cleft are: Glu203, Leu178, Gln160, and Arg277 on the bottom of the cleft; Lys162, Asp173, and Arg216 along one wall; and Tyr157, Tyr281, Arg117, Arg80, and Asp78 along the other wall. Next, we compared the structure of NLR48 with the lyases from the aforementioned PL families (excluding PL14). The structural alignment of these four lyases showed that within the cleft only Glu119, Gln160, and Tyr281 had counterparts in four other lyase families, whereas Arg117 is replaced by a histidine in PL20; Lys162 is a histidine in PL7, PL13, and PL18 but a leucine in PL20; and Arg277 has lysine counterparts in PL7, PL13, and PL18 (Fig. S1). PL14 lyases differ in this respect, and their putative active-site residues are located in a different region of the cleft (
      • Ogura K.
      • Yamasaki M.
      • Yamada T.
      • Mikami B.
      • Hashimoto W.
      • Murata K.
      Crystal structure of family 14 polysaccharide lyase with pH-dependent modes of action.
      ).

      Site-directed mutagenesis of putative active-site and substrate-binding residues

      Based on these observations the cleft region was considered as a substrate-binding site. As such, Arg117, Gln160, Tyr281, and Lys162 were deemed potential candidates for active-site residues with Arg216 likely being involved in substrate binding (H-bonds to phosphates in the native structure). These residues were mutated as follows: R117N, Q160A, K162M, R216N, and Y281F. The single mutants were expressed and purified, and their activity on the ulvan substrate was measured. The enzymatic activity was monitored as a function of time by following the appearance of unsaturated sugars in the reaction product with a characteristic absorbance at 232 nm. Although the native NLR48 showed easily measurable activity, no measurable activity was observed for the mutants (Fig. 4). This confirmed their involvement in either catalysis or substrate binding.
      Figure thumbnail gr4
      Figure 4The activity measurements of the WT NLR48 and several mutants on ulvan polysaccharide. Only the naive NLR48 showed measurable activity.

      Structure of NLR48 mutant with a tetrasaccharide substrate

      Three of the aforementioned mutants, K162M, Q160A, and Y281F, were screened for crystallization. Well-diffracting crystals could only be obtained for K162M mutant, and these crystals required high concentration of phosphate ions. These crystals were soaked for 2 h in a solution containing the tetrasaccharide substrate and diffracted to 2.2-Å resolution. They belong to the space group P212121 and are isomorphous to the native NLR48 crystals. The K162M mutation did not introduce any significant changes into the protein structure; the r.m.s.d. between the native and mutant structures for all Cα atoms is 0.3 Å. The electron density map calculated for the complex structure showed a blob of positive density associated with each molecule. The entire tetrasaccharide could be built into this density (Fig. 5A). The tetrasaccharide does not fully enter the cleft but is docked in the entrance and does not reach into the part where the putative active site is located (residues conserved between the three lyases). We adopt the nomenclature of Davies et al. (
      • Davies G.J.
      • Wilson K.S.
      • Henrissat B.
      Nomenclature for sugar-binding subsites in glycosyl hydrolases.
      ) where the sugars are marked as +1, +2, etc. starting from the first sugar on the reducing end of the cleaved bond and −1, −2, etc. on the nonreducing end. Correspondingly, the sugar-binding sites on the protein are referred to as ± subsites. The tetrasaccharide is oriented with its reducing end directed toward the center of the cleft, indicating that it occupies the “−” subsites in the substrate-binding sites.
      Figure thumbnail gr5
      Figure 5Tetrasaccharide substrate bound to NLR48. A, the final 2FoDFc map contoured at 1.5 σ with the superimposed model of the tetrasaccharide. B, close-up of the tetrasaccharide and its contacts with NLR48.
      In contrast to the native structure, the NLR48 mutant harbors only one phosphate at site 1. The second phosphate, observed in site 2 of the native structure, is replaced by Rha3S of the substrate in the mutant. The phosphate in site 1 was refined to an occupancy of ∼0.7. Nevertheless, a high concentration of phosphate in the crystallization conditions might effectively block the center of the cleft from the incoming tetrasaccharide and prevent it from occupying the +1 subsite. Indeed, a comparison with the lyase–oligosaccharide complexes from PL7 and PL13 confirms that the tetrasaccharide in NLR48 occupies − subsites with the Rha3S on the reducing end occupying the −1 subsite. The first three sugars make multiple contacts with the lyase. Rha3S sugar at the −1 subsite interacts with Tyr281 (at the beginning of the crossover loop following β13), Arg216 (beginning of β9), and Tyr157 (beginning of β6), whereas its sulfate H-bonds with Arg282 (loop following β13). GlcA at the −2 subsite interacts with Lys284 (crossover loop after β13), and its carboxylic group H-bonds with Ser211 (loop connecting β8 and β9). Finally, Rha3S at the −3 subsite interacts through its sulfate group with Ser213 (loop connecting β8 and β9) (Fig. 5B). The last sugar, ΔUA, makes no contacts with this molecule of NLR48 or with symmetry-related molecules, but its rotational freedom is restricted through the constraints imposed by the presence of bulky substituents at C3 and C5 positions of the preceding Rha3S sugar. Consequently, this sugar is more mobile than the others and has higher B-factors. No aromatic stacking interactions are observed in the lyase–tetrasaccharide contacts. Another factor that may restrict passage of the tetrasaccharide into the cleft is the dynamics of walls motion: they may have to move out to allow the substrate to slide in. In particular, the loop 164–174 with the Asp173 protruding into the cleft could play a gating function in solution. However, in the crystal structure, this loop butts against its counterpart in the neighboring molecule related by the noncrystallographic 2-fold axis (Fig. 6) and has no ability to expand the cleft. Because the complex was obtained by soaking the substrate into a preformed crystal, this could be an additional obstacle to the presence of a phosphate in the center of the cleft.
      Figure thumbnail gr6
      Figure 6Packing of a neighboring molecule against the loop forming one wall of the substrate-binding site. This tight packing restricts the movement of this loop and partially blocks the entry into the active site.
      Interestingly, an additional phosphate ion was observed in this structure near Arg124 and Arg277 that define the “+” subsites. Although these side chains have the same conformation in the native structure, there are only water molecules in this area. This new phosphate-binding site might reflect a subtle change of the electrostatic potential when the substrate in occupies the − sites.

      Catalytic mechanism

      The β-elimination catalytic mechanism requires neutralization of the negative charge on the C5 carboxylic group. This reduces the pKa of the C5 proton and makes it vulnerable to abstraction by a base. Concomitant donation of a proton by a protonated amino acid or water molecule to the leaving group results in the cleavage of the O–C4 bond. The end product carries an unsaturated molecule at the newly formed nonreducing end (
      • Gacesa P.
      Alginate-modifying enzymes. A proposed unified mechanism of action for the lyases and epimerases.
      ). Two general mechanisms have been identified among polysaccharide lyases: (i) metal-assisted catalysis and (ii) His/Tyr–dependent catalysis (
      • Garron M.L.
      • Cygler M.
      Structural and mechanistic classification of uronic acid-containing polysaccharide lyases.
      ,
      • Garron M.L.
      • Cygler M.
      Uronic polysaccharide degrading enzymes.
      ). Lyases adopting the β-jelly roll fold utilize the latter mechanism. Based on residue conservation, it is likely that NLR48 utilizes the same mechanism.
      To define the roles of residues within the cleft, we superimposed the three lyase–oligosaccharide complexes, NLR48, heparinase I, and alginate lyase (A1-II′;Aly), based only on the backbone atoms of the five conserved residues representing the local environment of the active site (Fig. 7A). In this superposition, the only sugar that has a similar position and orientation is the +1 uronate (Fig. 7B). The +2 and −1 sugars are oriented differently in these structures, reflecting different local shapes of the clefts. Nevertheless, the O1 of Rha3S in the NLR48 complex is in nearly the same place as the corresponding anomeric oxygen atoms in −1 sugars of the other two complexes (Fig. 7B). It is therefore very likely that the +1 GlcA/IdoA sugar is placed in the same way as in the other two lyases. In such a case, Gln160 would serve to neutralize the acidic group on the +1 sugar, promoting protonation through formation of two hydrogen bonds (Fig. 7A). Additional stabilization of the acidic group of the +1 sugar is provided by Arg117, which forms a hydrogen bond to the acidic group of +1 sugar. The Arg117-Gln160 arrangement is further stabilized by the formation of hydrogen bonds with a bridging Glu119 (Fig. 7A). The same arrangement is present in the other two lyases. The hydroxyl group of Tyr281 would point toward the bridging oxygen. With GlcA in the +1 site, the Tyr281 hydroxyl would also be in close proximity to the C5 proton. This residue could therefore perform the proton abstraction and donate it to the leaving group, the oxygen bridging the −1 and +1 sugars (Scheme 2, upper left). In the case of IdoA, the C5 proton is on the opposite side of the sugar ring to Tyr281 (Scheme 2, bottom left). The only side chain close to this proton is Lys162, which was mutated to methionine in the NLR48–tetrasaccharide complex (Fig. 5A). There is a histidine in an equivalent position in PL13 lyases that break the bond next to IdoA and in PL7 lyases that break the bond next to guluronate with the same configuration of the C5 carboxylic group as IdoA. In PL20 lyases, which break the bond next to GlcA, the C5 proton is facing the tyrosine, and the residue in a position equivalent to Lys162 does not participate in the reaction. Indeed, an isoleucine is found in the PL20 family. Because NRL48 can likely cleave the bond next to either epimer of the uronic acid, like its homolog NLR42 (
      • Nyvall Collén P.
      • Sassi J.F.
      • Rogniaux H.
      • Marfaing H.
      • Helbert W.
      Ulvan lyases isolated from the flavobacteria Persicivirga ulvanivorans are the first members of a new polysaccharide lyase family.
      ), Lys162 must be able to act as a base for proton abstraction. This would require an unusual pKa for this side chain in the context of the enzyme. The pKa could be influenced by the proximity to 3-sulfate of +2 Rha3S. In heparinase I, the equivalent residue is a histidine. The role of a lysine residue as a proton acceptor is not unusual; it fulfills the same role in lyases from several other families with different folds, namely PL3, PL9, PL6, and PL11 (
      • Garron M.L.
      • Cygler M.
      Structural and mechanistic classification of uronic acid-containing polysaccharide lyases.
      ). Overall, both the fold and catalytic machinery are conserved among families PL28, PL7, PL13, PL18, and PL20 (Fig. 7).
      Figure thumbnail gr7
      Figure 7Stereoview of the superposition of enzyme–oligosaccharide complexes based on the backbone of five conserved residues in the cleft (Arg117, Glu119, Gln160, Arg277, and Tyr281 and their equivalences in other lyases). A, the key residues are displayed in a stick representation. NLR48, green carbons; PL7 (PDB code 2ZAA), cyan carbons; PL13 (PDB code 2INA), yellow carbons. Only the +1 sugar in PL7 and PL13 is shown for clarity. The Rha3S in −1 subsite of NLR48 is shown. The hydrogen-bond network among the acidic group of +1 sugar, Gln, Arg, and Glu is shown with cyan dashed lines. B, overlay of the oligosaccharides from the three lyases from the above superposition. Only −2 to +2 sugars are shown. NLR48, green carbons; PL7 (PDB code 2ZAA), pink carbons; PL13 (PDB code 2INA), yellow carbons. Only the +1 sugar has a similar position and orientation as is required for the catalysis. The orientations of the − and + sugars differ to be accommodated well into the clefts of individual enzymes. The O1 of Rha3S in NLR48 is very close to the oxygen bridging −1 and +1 sugars, suggesting that the +1 sugar would take a similar position to that observed in the other two lyases.
      Figure thumbnail sc2
      Scheme 2Proposed mechanism of catalysis. Tyr281 acts as a base and an acid for the cleavage near GlcA, whereas for IdoA the base is Lys162, and the acid is Tyr281.

      Comparison with other ulvan lyases

      Lyases that degrade ulvan polysaccharide have been identified and characterized only in the last several years. They have led to the addition of three new families to the CAZy classification, namely PL24, PL25, and PL28 families, in the CAZy database. We have previously characterized the structures and catalytic mechanisms of LOR107 (PL24) and PLSV3936 (PL25). Both ulvan lyases adopt a seven-bladed β-propeller fold. PLSV3936 utilizes the His/Tyr catalytic machinery with Arg204, Tyr188, and His123 as neutralizer, catalytic base, and acid, respectively (
      • Ulaganathan T.
      • Boniecki M.T.
      • Foran E.
      • Buravenkov V.
      • Mizrachi N.
      • Banin E.
      • Helbert W.
      • Cygler M.
      New ulvan-degrading polysaccharide lyase family: structure and catalytic mechanism suggests convergent evolution of active site architecture.
      ). LOR107 uses Arg259 to neutralize the C5 carboxylic group, and His146 appears to play both catalytic base/acid roles aided by His167 and Tyr243 (
      • Ulaganathan T.
      • Helbert W.
      • Kopel M.
      • Banin E.
      • Cygler M.
      Structure-function analyses of a PL24 family ulvan lyase reveal key features and suggest its catalytic mechanism.
      ). Moreover, LOR107 binds the substrate in two stages with Arg320 initiating substrate binding at the top of the cleft followed by substrate sliding deeper into the cleft toward the active site (
      • Ulaganathan T.
      • Helbert W.
      • Kopel M.
      • Banin E.
      • Cygler M.
      Structure-function analyses of a PL24 family ulvan lyase reveal key features and suggest its catalytic mechanism.
      ).
      Despite its different fold, NLR48 shares similarities with the PL24 and PL25 lyases in their catalytic mechanisms. All of them utilize the His/Tyr mechanism; however, the neutralization of the C5 acidic group is achieved by different means. NLR48 uses glutamine that faces the carboxyl and forms a two H-bonds and is supported by a glutamate and an arginine, whereas PL24/25 utilize an arginine to form two H-bonds with the carboxyl.

      Conclusions

      The marine genome sequencing project resulted in the identification of many ulvan-degrading enzymes. These are classified under PL24, PL25, and PL28 families in the CAZy database. Although PL24 and PL25 ulvan lyases have a seven-bladed β-propeller fold, PL28 lyases have a β-jelly roll fold, adding another family to the list of PL enzymes with this architecture. The catalytic mechanism of the PL28 representative lyase, NLR48 from N. ulvanivorans, is very similar to that of other lyase families with a β-jelly roll fold despite very low sequence identity and is a modification of the His/Tyr catalytic mechanism. The fold similarity of PL28 to several other families of polysaccharide lyases indicates their common evolutionary origin with adaptation to different substrates. Finally, the three families of ulvan lyases utilize a different mode of neutralization of the carboxylic group of the uronic acid.

      Materials and methods

      Cloning

      The gene encoding NLR48 (NCBI accession number KEZ94336), lacking the first 25 residues as they are predicted to be a signal peptide, was cloned into three expression vectors with different tags to improve solubility and to aid in purification. NLR48(26–303) was inserted into pET28 vector using the restriction cloning method, whereas the ligation-independent cloning protocol (
      • Aslanidis C.
      • de Jong P.J.
      Ligation-independent cloning of PCR products (LIC-PCR).
      ) was used to insert NLR48(26–303) into pRL652 (modified version of pGEX with N-terminal GST tag followed by TEV cleavage site). All mutants used in this study were based on GST-NLR48(26–303) as a template by following the QuikChange site-directed mutagenesis protocol. All constructs were verified by DNA sequencing.

      Expression and purification of NLR48

      GST-NLR48(26–303) was expressed in BL21(DE3) strain. Overnight inoculum was subcultured into 6 liters of Terrific Broth medium. The cells were induced at ∼1.5 A600 with 0.5 mm isopropyl β-d-1-thiogalactopyranoside. After 16 h of growth at 18 °C, cells were spun down at 6200 × g for 20 min. The cell pellet was resuspended in lysis buffer (20 mm Tris, pH 8, 400 mm NaCl, 5% glycerol) and disrupted by sonication. The supernatant was separated from cell debris by centrifugation at 39,000 × g for 55 min. The clarified supernatant was incubated for 2 h with pre-equilibrated GST resin. Unbound proteins were removed by washing with high-salt buffer (lysis buffer supplemented with 200 mm NaCl). All purification steps were performed at 4 °C. The GST tag was cleaved on the column by incubation for ∼10 h with His-TEV protease, added at a mass ratio of 1:50. The protein eluted from the column was loaded on a nickel-nitrilotriacetic acid column to remove the TEV protease. The flow-through containing highly pure NLR48 was loaded on a size-exclusion column for final purification.

      Crystallization of NLR48

      NLR48 was concentrated to 16 mg/ml and screened for crystallization by the sitting-drop vapor-diffusion method in 96-well plate format by mixing 300 nl of the protein with 300 nl of well solution using the Gryphon crystallization robot (ArtRobbins Instruments, Sunnyvale, CA). The duplicate plates were maintained at temperatures of 9 and 15 °C. The best crystals grew from solution containing 0.056 m sodium phosphate monobasic, 1.344 m potassium phosphate dibasic, pH 8.2. Crystals grown at 9 °C were somewhat bigger than crystals grown at 15 °C.

      Expression, purification, and crystallization of selenomethionine-labeled NLR48

      Expression plasmid carrying GST-NLR48 fusion construct was transformed into methionine auxotroph cells (B834(DE3), Novagen). Overnight inoculum grown in 50 ml of M9 medium supplemented with methionine was subcultured into 1 liter of M9 medium. For SeMet labeling, the cells were grown at 37 °C until A600 reached 1.0. The culture was spun gently at 3000 × g for 20 min, the cell pellet was resuspended in 1 liter of M9 medium, and the culture continued to grow at 37 °C for another 4 h to deplete the remaining methionine. Subsequently, 50 mg/liter SeMet was added to the cell culture 30 min prior to the addition of isopropyl β-d-1-thiogalactopyranoside. The cells were further grown at 18 °C. After ∼16 h of growth, the cells were collected by centrifugation, and SeMet-incorporated protein was purified following the same protocol as used for the native enzyme purification. SeMet-incorporated NLR48 was concentrated to 16 mg/ml for crystallization experiments using the hanging-drop vapor-diffusion method. The crystals grown from 0.49 m sodium phosphate monobasic monohydrate, 0.91 m potassium phosphate dibasic, pH 6.9, at 9 °C were cryoprotected with 30% ethylene glycol, and diffraction data were collected at the 08BM beamline at the Canadian Light Source. Data were processed with the XDS software package (
      • Kabsch W.
      XDS.
      ) to 1.9-Å resolution, and the data collection statistics are shown in Table 2. These crystals diffracted significantly better and showed better spot shape than the native crystals and were used for structure determination and refinement.
      Table 2Data collection and refinement statistics
      Crystal form
      NLR48(SeMet)NLR48(K162M)–tetrasaccharide
      Data collection
       Space groupP21 21 21P21 21 21
      a, b, c (Å)82.5, 102.1, 103.782.5, 102.7, 103.0
       Wavelength (Å)0.98013
      The data set for structure determination was collected at the wavelength of 0.97864. The data set used for refinement was collected at the wavelength 0.98013 as shown in the table.
      0.97949
       Resolution (Å)50–1.950–2.2
       Observed hkl1,001,093432,911
       No. of unique hkl131,85444,398
       Completeness (%)99.8 (97.7)
      The numbers in parentheses refer to the highest resolution shell.
      99.8 (97)
       Redundancy7.69.8
      Rmeas0.143 (1.25)0.100 (1.16)
      I/(σI)10.01 (1.82)14.86 (2.0)
      Refinement statistics
      Rwork0.1690.189
      Rfree0.2010.220
       Wilson B (Å2)26.2849.57
       Ramachandran plot (%)
       Favored96.5795.35
       Allowed2.673.88
       Outliers0.760.78
       r.m.s.d.
       Bonds (Å)0.0150.008
       Angles (°)1.2141.091
      PDB code6D2C6D3U
      a The data set for structure determination was collected at the wavelength of 0.97864. The data set used for refinement was collected at the wavelength 0.98013 as shown in the table.
      b The numbers in parentheses refer to the highest resolution shell.

      Structure solution and refinement of WT NLR48

      The structure of NLR48 was solved by the single anomalous dispersion method with the SeMet data set using the AutoSol procedure in Phenix software (
      • Adams P.D.
      • Afonine P.V.
      • Bunkóczi G.
      • Chen V.B.
      • Davis I.W.
      • Echols N.
      • Headd J.J.
      • Hung L.W.
      • Kapral G.J.
      • Grosse-Kunstleve R.W.
      • McCoy A.J.
      • Moriarty N.W.
      • Oeffner R.
      • Read R.J.
      • Richardson D.C.
      • et al.
      PHENIX: a comprehensive Python-based system for macromolecular structure solution.
      ). There are two molecules in the asymmetric unit. All six expected selenium sites were identified. The initial model contained ∼90% of the expected residues. The remaining residues were built manually using Coot (
      • Emsley P.
      • Lohkamp B.
      • Scott W.G.
      • Cowtan K.
      Features and development of Coot.
      ). Because the SeMet-containing crystal diffracted better than WT NLR48, this data set was used not only for solving the structure but also for refinement. The refinement was performed with Phenix software interspersed with manual rebuilding using Coot. The final Rwork was 0.169, and Rfree was 0.201. Several N-terminal residues are not visible in the electron density. The model contains residues 44–303 in both chains, two Ca2+ ions, four phosphates, and 710 solvent molecules. The validation by MolProbity (
      • Chen V.B.
      • Arendall 3rd, W.B.
      • Headd J.J.
      • Keedy D.A.
      • Immormino R.M.
      • Kapral G.J.
      • Murray L.W.
      • Richardson J.S.
      • Richardson D.C.
      MolProbity: all-atom structure validation for macromolecular crystallography.
      ) identified Arg282 and Val283 as Ramachandran outliers. However, inspection of electron density map confirms the correct modeling of these residues. Data collection and statistics are shown in Table 2. The coordinates and structure factors have been deposited to the Protein Data Bank with accession code 6D2C.

      Site-directed mutagenesis

      Selected residues were mutated to assess their role in catalysis. The mutants R117A, Q160A, K162M, R217N, and Y281F were made using the GST-NLR48 fusion construct as a template by following the QuikChange site-directed mutagenesis protocol. All the constructs were confirmed by sequencing. The mutants were expressed in BL21(DE3) cells and purified using the same protocol as for the WT protein.

      Enzyme activity

      The activity of NLR48 and its various mutants was measured by monitoring the absorption at 232 nm, indicative of the formation of an unsaturated double bond in the sugar on the nonreducing end of one of the two reaction products. The assay solution contained 0.5 mg/ml ulvan polysaccharide in 20 mm Hepes, pH 7.5, 150 mm NaCl. 1 μg of the enzyme was added directly to the assay solution with constant stirring, and the absorbance was measured as a function of time. The absorbance of the substrate alone was used as a control.

      Expression, purification, and crystallization of NLR48 mutants

      Mutants showing no activity, namely K162M, Y281F, and R117N, were expressed and purified. The GST fusion constructs of the mutants were expressed and purified using the same protocol as for WT NLR48. Of these, only K162M was expressed in high yield and was screened for crystallization at 17 mg/ml. Crystals were obtained in the same conditions as for WT NLR48 (0.49 m sodium phosphate monobasic monohydrate, 0.91 m potassium phosphate dibasic, pH 6.9, at 9 °C).

      Structure determination of NLR48–tetrasaccharide complex

      The K162M mutant crystals were used for soaking experiments with an ulvan tetrasaccharide, which contained a mixture of ΔUA-Rha3S-GlcA-Rha3Sα/β and ΔUA-Rha3S-IdoA-Rha3Sα/β. Briefly, a drop containing 3 μl of substrate solution (∼50 mm tetrasaccharide solution in water) mixed with 3 μl of crystallization well solution on a coverslip was equilibrated against 1 ml of well solution. The tetrasaccharide was prepared as described before (
      • Ulaganathan T.
      • Helbert W.
      • Kopel M.
      • Banin E.
      • Cygler M.
      Structure-function analyses of a PL24 family ulvan lyase reveal key features and suggest its catalytic mechanism.
      ). The crystals were soaked in the substrate solution for ∼2 h, cryoprotected in well solution supplemented with 30% ethylene glycol, and flash frozen in liquid N2. Diffraction data were collected at the 08ID beamline, Canadian Light Source. These crystals diffracted to 2.2-Å resolution and were isomorphous with the WT NLR48 crystals. One round of rigid-body refinement was done on the complex dataset with the Phenix refinement program. The tetrasaccharide was built into the positive difference electron density with Coot. Further refinement was done using Phenix software. The complex structure was evaluated for stereochemistry using MolProbity. The final model has 44–303 residues in each molecule, 212 water molecules, two calcium ions, two tetrasaccharides, and four phosphate ions. Data collection and statistics are shown in Table 2. The coordinates and structure factors have been deposited to the Protein Data Bank with accession code 6D3U.

      Author contributions

      T. U. and M. C. conceptualization; T. U. data curation; T. U. and M. C. validation; T. U. investigation; T. U., E. B., and W. H. methodology; T. U. writing-original draft; E. B. and W. H. resources; W. H. and M. C. writing-review and editing; M. C. supervision; M. C. funding acquisition.

      Acknowledgments

      We thank Kevin Voth for comments. We acknowledge the Protein Characterization and Crystallization Facility, College of Medicine, University of Saskatchewan for access to the crystallization robot. Research described in this paper was performed using beamline 08ID-1 at the Canadian Light Source synchrotron, which is supported by the Canada Foundation for Innovation, Natural Sciences and Engineering Research Council of Canada, the University of Saskatchewan, the Government of Saskatchewan, Western Economic Diversification Canada, the National Research Council Canada, and the Canadian Institutes of Health Research.

      Supplementary Material

      References

        • Shalaby E.A.
        Algae as promising organisms for environment and health.
        Plant Signal. Behav. 2011; 6 (21862867): 1338-1350
        • Sudha P.N.
        • Aisverya S.
        • Nithya R.
        • Vijayalakshmi K.
        Industrial applications of marine carbohydrates.
        Adv. Food Nutr. Res. 2014; 73 (25300546): 145-181
        • Hehemann J.-H.
        • Boraston A.B.
        • Czjzek M.
        A sweet new wave: structures and mechanisms of enzymes that digest polysaccharides from marine algae.
        Curr. Opin. Struct. Biol. 2014; 28 (25136767): 77-86
        • Michel G.
        • Czjzek M.
        Polysaccharide-degrading enzymes from marine bacteria.
        in: Trincone A. Marine Enzymes for Biocatalysis. Woodhead Publishing Ltd., Cambridge, UK2013: 429-464
        • Cunha L.
        • Grenha A.
        Sulfated seaweed polysaccharides as multifunctional materials in drug delivery applications.
        Mar. Drugs. 2016; 14 (26927134): E42
        • Popa E.G.
        • Reis R.L.
        • Gomes M.E.
        Seaweed polysaccharide-based hydrogels used for the regeneration of articular cartilage.
        Crit. Rev. Biotechnol. 2015; 35 (24646368): 410-424
        • Lahaye M.
        • Robic A.
        Structure and functional properties of ulvan, a polysaccharide from green seaweeds.
        Biomacromolecules. 2007; 8 (17458931): 1765-1774
        • Nyvall Collén P.
        • Sassi J.F.
        • Rogniaux H.
        • Marfaing H.
        • Helbert W.
        Ulvan lyases isolated from the flavobacteria Persicivirga ulvanivorans are the first members of a new polysaccharide lyase family.
        J. Biol. Chem. 2011; 286 (22009751): 42063-42071
        • Melcher R.L.
        • Neumann M.
        • Fuenzalida Werner J.P.
        • Gröhn F.
        • Moerschbacher B.M.
        Revised domain structure of ulvan lyase and characterization of the first ulvan binding domain.
        Sci. Rep. 2017; 7 (28327560): 44115
        • Kopel M.
        • Helbert W.
        • Henrissat B.
        • Doniger T.
        • Banin E.
        Draft genome sequence of Pseudoalteromonas sp. strain PLSV, an ulvan-degrading bacterium.
        Genome Announc. 2014; 2 (25502665): e01257-14
        • Kopel M.
        • Helbert W.
        • Henrissat B.
        • Doniger T.
        • Banin E.
        Draft genome sequences of two ulvan-degrading isolates, strains LTR and LOR, that belong to the Alteromonas genus.
        Genome Announc. 2014; 2 (25342689): e01081-14
        • Kopel M.
        • Helbert W.
        • Henrissat B.
        • Doniger T.
        • Banin E.
        Draft genome sequence of Nonlabens ulvanivorans, an ulvan-degrading bacterium.
        Genome Announc. 2014; 2 (25125644): e00793-14
        • Foran E.
        • Buravenkov V.
        • Kopel M.
        • Mizrahi N.
        • Shoshani S.
        • Helbert W.
        • Banin E.
        Functional characterization of a novel “ulvan utilization loci” found in Alteromonas sp. LOR genome.
        Algal Res. 2017; 25: 39-46
        • Lombard V.
        • Golaconda Ramulu H.
        • Drula E.
        • Coutinho P.M.
        • Henrissat B.
        The carbohydrate-active enzymes database (CAZy) in 2013.
        Nucleic Acids Res. 2014; 42 (24270786): D490-D495
        • Kopel M.
        • Helbert W.
        • Belnik Y.
        • Buravenkov V.
        • Herman A.
        • Banin E.
        New family of ulvan lyases identified in three isolates from the Alteromonadales order.
        J. Biol. Chem. 2016; 291 (26763234): 5871-5878
        • Ulaganathan T.
        • Boniecki M.T.
        • Foran E.
        • Buravenkov V.
        • Mizrachi N.
        • Banin E.
        • Helbert W.
        • Cygler M.
        New ulvan-degrading polysaccharide lyase family: structure and catalytic mechanism suggests convergent evolution of active site architecture.
        ACS Chem. Biol. 2017; 12 (28290654): 1269-1280
        • Jurrus E.
        • Engel D.
        • Star K.
        • Monson K.
        • Brandi J.
        • Felberg L.E.
        • Brookes D.H.
        • Wilson L.
        • Chen J.
        • Liles K.
        • Chun M.
        • Li P.
        • Gohara D.W.
        • Dolinsky T.
        • Konecny R.
        • et al.
        Improvements to the APBS biomolecular solvation software suite.
        Protein Sci. 2018; 27 (28836357): 112-128
        • Garron M.L.
        • Cygler M.
        Structural and mechanistic classification of uronic acid-containing polysaccharide lyases.
        Glycobiology. 2010; 20 (20805221): 1547-1573
        • Ogura K.
        • Yamasaki M.
        • Mikami B.
        • Hashimoto W.
        • Murata K.
        Substrate recognition by family 7 alginate lyase from Sphingomonas sp. A1.
        J. Mol. Biol. 2008; 380 (18514736): 373-385
        • Han Y.H.
        • Garron M.L.
        • Kim H.Y.
        • Kim W.S.
        • Zhang Z.
        • Ryu K.S.
        • Shaya D.
        • Xiao Z.
        • Cheong C.
        • Kim Y.S.
        • Linhardt R.J.
        • Jeon Y.H.
        • Cygler M.
        Structural snapshots of heparin depolymerization by heparin lyase I.
        J. Biol. Chem. 2009; 284 (19801541): 34019-34027
        • Dong S.
        • Wei T.-D.
        • Chen X.-L.
        • Li C.-Y.
        • Wang P.
        • Xie B.-B.
        • Qin Q.-L.
        • Zhang X.-Y.
        • Pang X.-H.
        • Zhou B.-C.
        • Zhang Y.-Z.
        Molecular insight into the role of the N-terminal extension in the maturation, substrate recognition, and catalysis of a bacterial alginate lyase from polysaccharide lyase family 18.
        J. Biol. Chem. 2014; 289 (25210041): 29558-29569
        • Konno N.
        • Ishida T.
        • Igarashi K.
        • Fushinobu S.
        • Habu N.
        • Samejima M.
        • Isogai A.
        Crystal structure of polysaccharide lyase family 20 endo-β-1,4-glucuronan lyase from the filamentous fungus Trichoderma reesei.
        FEBS Lett. 2009; 583 (19306878): 1323-1326
        • Ogura K.
        • Yamasaki M.
        • Yamada T.
        • Mikami B.
        • Hashimoto W.
        • Murata K.
        Crystal structure of family 14 polysaccharide lyase with pH-dependent modes of action.
        J. Biol. Chem. 2009; 284 (19846561): 35572-35579
        • Landau M.
        • Mayrose I.
        • Rosenberg Y.
        • Glaser F.
        • Martz E.
        • Pupko T.
        • Ben-Tal N.
        ConSurf 2005: the projection of evolutionary conservation scores of residues on protein structures.
        Nucleic Acids Res. 2005; 33 (15980475): W299-W302
        • Davies G.J.
        • Wilson K.S.
        • Henrissat B.
        Nomenclature for sugar-binding subsites in glycosyl hydrolases.
        Biochem. J. 1997; 321 (9020895): 557-559
        • Gacesa P.
        Alginate-modifying enzymes. A proposed unified mechanism of action for the lyases and epimerases.
        FEBS Lett. 1987; 212: 199-202
        • Garron M.L.
        • Cygler M.
        Uronic polysaccharide degrading enzymes.
        Curr. Opin. Struct. Biol. 2014; 28 (25156747): 87-95
        • Ulaganathan T.
        • Helbert W.
        • Kopel M.
        • Banin E.
        • Cygler M.
        Structure-function analyses of a PL24 family ulvan lyase reveal key features and suggest its catalytic mechanism.
        J. Biol. Chem. 2018; 293 (29382716): 4026-4036
        • Aslanidis C.
        • de Jong P.J.
        Ligation-independent cloning of PCR products (LIC-PCR).
        Nucleic Acids Res. 1990; 18 (2235490): 6069-6074
        • Kabsch W.
        XDS.
        Acta Crystallogr. D Biol. Crystallogr. 2010; 66 (20124692): 125-132
        • Adams P.D.
        • Afonine P.V.
        • Bunkóczi G.
        • Chen V.B.
        • Davis I.W.
        • Echols N.
        • Headd J.J.
        • Hung L.W.
        • Kapral G.J.
        • Grosse-Kunstleve R.W.
        • McCoy A.J.
        • Moriarty N.W.
        • Oeffner R.
        • Read R.J.
        • Richardson D.C.
        • et al.
        PHENIX: a comprehensive Python-based system for macromolecular structure solution.
        Acta Crystallogr. D Biol. Crystallogr. 2010; 66 (20124702): 213-221
        • Emsley P.
        • Lohkamp B.
        • Scott W.G.
        • Cowtan K.
        Features and development of Coot.
        Acta Crystallogr. D Biol. Crystallogr. 2010; 66 (20383002): 486-501
        • Chen V.B.
        • Arendall 3rd, W.B.
        • Headd J.J.
        • Keedy D.A.
        • Immormino R.M.
        • Kapral G.J.
        • Murray L.W.
        • Richardson J.S.
        • Richardson D.C.
        MolProbity: all-atom structure validation for macromolecular crystallography.
        Acta Crystallogr. D Biol. Crystallogr. 2010; 66 (20057044): 12-21
        • Guex N.
        • Peitsch M.C.
        SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling.
        Electrophoresis. 1997; 18 (9504803): 2714-2723