Skip to main content

Comparative characterization of all cellulosomal cellulases from Clostridium thermocellum reveals high diversity in endoglucanase product formation essential for complex activity

Abstract

Background

Clostridium thermocellum is a paradigm for efficient cellulose degradation and a promising organism for the production of second generation biofuels. It owes its high degradation rate on cellulosic substrates to the presence of supra-molecular cellulase complexes, cellulosomes, which comprise over 70 different single enzymes assembled on protein-backbone molecules of the scaffold protein CipA.

Results

Although all 24 single-cellulosomal cellulases were described previously, we present the first comparative catalogue of all these enzymes together with a comprehensive analysis under identical experimental conditions, including enzyme activity, binding characteristics, substrate specificity, and product analysis. In the course of our study, we encountered four types of distinct enzymatic hydrolysis modes denoted by substrate specificity and hydrolysis product formation: (i) exo-mode cellobiohydrolases (CBH), (ii) endo-mode cellulases with no specific hydrolysis pattern, endoglucanases (EG), (iii) processive endoglucanases with cellotetraose as intermediate product (pEG4), and (iv) processive endoglucanases with cellobiose as the main product (pEG2). These modes are shown on amorphous cellulose and on model cello-oligosaccharides (with degree of polymerization DP 3 to 6). Artificial mini-cellulosomes carrying combinations of cellulases showed their highest activity when all four endoglucanase-groups were incorporated into a single complex. Such a modeled nonavalent complex (n = 9 enzymes bound to the recombinant scaffolding protein CipA) reached half of the activity of the native cellulosome. Comparative analysis of the protein architecture and structure revealed characteristics that play a role in product formation and enzyme processivity.

Conclusions

The identification of a new endoglucanase type expands the list of known cellulase functions present in the cellulosome. Our study shows that the variety of processivities in the enzyme complex is a key enabler of its high cellulolytic efficiency. The observed synergistic effect may pave the way for a better understanding of the enzymatic interactions and the design of more active lignocellulose-degrading cellulase cocktails in the future.

Background

Due to the complex structure of plant cell walls, biomass-derived polysaccharides embody a rich tapestry of sugars and sugar compositions which are degraded by cellulases and other glycoside-depolymerizing enzymes. These enzymes can be described by three-dimensional structural analysis, sequence-based classification, substrate specificity, hydrolytic reaction mode, kinetic parameters, and product formation. Among carbohydrate-active enzymes, the CAZy database [1] classified 145 different glycoside hydrolase (GH) families as of August 2017, whereas cellulases are represented by 14 different GH families. The ill-defined term “cellulase” is generally taken to describe enzymes that depolymerize β-1,4-glycosidic bonds in β-glucans from cellulosic biomass. However, various cellulase types can be distinguished by their different modes of catalytic action. Exo-acting cellobiohydrolases hydrolyze the polysaccharide chain either from the reducing or non-reducing end, while endoglucanases cleave within the cellulose chain to generate new ends that are susceptible to subsequent hydrolysis by exoglucanase enzymes [2]. Binding of the enzyme to the substrate requires the presence of specific carbohydrate-binding modules (CBM) and sugar-binding residues on the enzyme surface and catalytic cleft.

According to the first description by Koshland, the catalytic reaction is retaining or inverting, depending on the nucleophilic attack at the glycosidic bond of the polysaccharide and the resulting stereochemistry of the anomeric carbon [3]. The measurement and classification of cellulase processivity is a daunting task, due to a variety of available assay techniques and a lack of established standards [2, 4]. Processivity can generally be defined as the average number of cleavages on the cellulose chain, before the enzyme dissociates from the substrate (catalytic rate coefficient k cat divided by dissociation rate coefficient k off) [5]. The key differentiating factors among processivity of cellulases have been studied mainly in fungal cellulases and comprise the following: (i) The presence of loop structures to form a tunnel which covers the active site during the processive movement on the cellulose chain [6], (ii) the presence of certain CBMs linked to the catalytic core of an endoglucanase [7], and (iii) the presence of subsites for sugar binding and affinity [8]. Exo-acting cellulases known to hydrolyze the cellulose chains from the reducing ends are GH7 and GH48 enzymes, while enzymes processively acting from the non-reducing ends are GH9 and GH6 [5]. Other processive endoglucanases have also been reported for certain enzymes from the GH5 and GH9 families, such as Cel5H from Saccharophagus degradans [9] and Cel9I from C. thermocellum [10], respectively. In addition, cellulase actions are dictated by further structure–function–stability relationships, e.g., (N-terminal) extensions for stabilization of the catalytic core [11], the presence of specific ion binding sites for selective thermostabilization [12], or the influence of the quaternary structure on substrate specificity [13]. Instead of measuring the “apparent” processivity of cellulases, computational and structural modeling has been used to explain the “intrinsic” processivity of cellulases on a molecular level, as reviewed by [4, 14].

The cellulosomal complex of Clostridium thermocellum is one of the most efficient cellulase systems discovered to date [15]. This multi-modular enzyme system is based on the immobilization and co-localization of over 70 different proteins on a scaffolding structural protein, whereby different enzyme types act synergistically to efficiently degrade the polysaccharide into soluble sugars [16]. Interestingly, transcriptomic and proteomic analysis revealed that the cellulosome contains redundant sets of different cellulases and that regulation of their expression is a function of the substrate [17,18,19]. Nevertheless, the debate over why C. thermocellum (and other cellulolytic bacteria) express such vast and varied numbers of cellulases remains active.

To our knowledge, a comparative characterization of all β-1,4-glucanases present in the cellulosome has not been reported. In this study we characterize the product formation of 24 cellulases on different soluble and insoluble cellulosic substrates and β-1,4-glucans. Furthermore, a comprehensive comparison of activity profiles and product formation kinetics on model oligosaccharides and PASC (phosphoric acid swollen cellulose) is presented. We were able to differentiate between the apparent product spectra formed by GH5 and GH9 endoglucanases. To this end, a hydrolysis product pattern for Cel9D and four GH5 endoglucanases from sub-family 1 (Cel5O, Cel5B, Cel5G and Cel5L) was identified which distinguishes it from all other endoglucanase or cellobiohydrolase (CBH) hydrolysis patterns. Furthermore, we show that this new type of endoglucanolytic cleavage may have implications on the overall hydrolytic efficiency of synthetic (mini-)cellulosomes towards microcrystalline cellulose. The disparity in apparent processivity and substrate preference between glycoside hydrolases of family 9 (GH9) was supported by molecular docking experiments as well as sequence analysis revealing the presence of carbohydrate-binding modules (CBM) and sugar-binding moieties. Our data contribute to a deeper understanding of the cellulosomal cellulase system and may be of relevance for the design and engineering of more efficient enzyme mixtures for biomass degradation in the future.

Methods

Strains, media, and chemicals

Clostridium thermocellum (in the literature also referred to as “Ruminiclostridium thermocellum” [20]) DSM1237 was grown at 60 °C in prereduced GS-2 medium for liquid cultures containing 0.5% (w/v) cellobiose [21]. Recombinant Escherichia coli strains DH10B and BL21(DE) Star (Invitrogen, Carlsbad, USA) were used for cloning and protein expression, respectively. The cells were grown in Lysogeny broth containing 100 µg/mL ampicillin for pET21a(+) plasmids and 50 µg/mL kanamycin for pET24(+) plasmids. If not stated otherwise, chemical reagents were purchased by Sigma-Aldrich (Taufkirchen, Germany).

DNA manipulation and synthesis

Preparation of chromosomal and plasmid DNA, endonuclease digestion, and ligation was carried out by standard procedures [22]. QIAprep Spin Miniprep Kit and PCR purification kit (Qiagen, Hilden, Germany) were used for purification of plasmids and PCR products. Restriction digests of DNA were done as recommended by the manufacturer (NEB, Ipswich, USA). Chemically competent E. coli DH10B cells were used for transformation with plasmid DNA.

Signal peptides were predicted by SignalP 3.0 server [23]. Genes without the signal sequence were amplified with oligonucleotide primers as listed in Additional file 1 and Phusion DNA Polymerase (NEB, Ipswich, USA) with chromosomal DNA from C. thermocellum DSM1237 as template. The synthesized genes cel124 (cthe_0435), cel9-44J, cel9K, and cel48S were optimized for E. coli codon usage by Eurofins (Ebersberg, Germany). The cellulosomal scaffolding protein CipA was synthesized in optimized E. coli codon usage and optimized DNA sequence, including eight cohesins Coh1-2, the carbohydrate-binding module CBM3, Coh3-8, and the C-terminal X-module from C. thermocellum WP_020458017.1 lacking Coh6 and Dockerin type-II. The resulting construct is referred to as CipA8 (see Additional file 2). The amplicons were digested and ligated in frame into the multiple cloning site of the plasmid pET21a(+). The correct sequence of all constructs was verified by resequencing (MWG, Ebersberg, Germany).

Protein purification

For protein expression, the plasmids were transformed into E. coli BL21(DE) Star. The cells were grown at 37 or 20 °C and protein expression from pET21(+) or pET24(+) plasmids was induced by addition of 1 mM isopropyl-β-d-thiogalactopyranoside (IPTG) to an exponentially growing culture. After further growth at 37 °C for 4 h, the cells were harvested by centrifugation at 3440×g (Sorvall RC 6 +, Thermo Fisher, Waltham, USA) for 10 min at 4 °C.

The cells were resuspended in 20 mL lysis buffer (50 mM MOPS pH 7.3, 100 mM NaCl, 10 mM CaCl2, 20 mM imidazole) with the addition of lysozyme (AppliChem, Darmstadt, Germany) to a final concentration of 10 mg/mL and incubated for 30 min on ice. The cells were sonified twice with Sonifier UP 200S (Hielscher, Teltow, Germany) set at amplitude 60%, interval 0.25 and for 4 min. The supernatant after centrifugation (18,000 rpm, 20 min, 4 °C) was loaded onto an immobilized metal HisTrap affinity column (IMAC) (GE Healthcare, Munich, Germany) and eluted with 0.5 M imidazole, 50 mM MOPS pH 7.3, 100 mM NaCl, 10 mM CaCl2. The proteins were examined by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) and stained with Coomassie brilliant blue R-250. The protein concentration was determined spectrophotometrically by measuring the absorbance at 280 nm in a 5 M urea solution (Additional file 3). All protein preparations contained 20% glycerol (v/v) or sucrose and 0.2% sodium azide (w/v) and were proven to be stable on storage at − 20 °C. Table 1 summarizes all proteins analyzed in this study.

Table 1 Summary of the cellulosomal cellulases from C. thermocellum analyzed in this study

Native cellulosome and SM901 extract preparation

Well-grown cultures of C. thermocellum mutant SM901, also referred to as SM1 [41] were centrifuged twice (13,000 rpm, 20 min). Extracellular proteins were precipitated from the cell-free supernatant using saturated (NH4)2SO4 solution added to a final concentration of 60% (v/v). After overnight incubation at 4 °C the proteins were collected by centrifugation (15,000 rpm, 20 min, 4 °C). Supernatant preparations from mutant SM901 were resuspended in 50 mM MES, 0.1 M NaCl, 5 mM CaCl2, pH 6.0. Cellulosomal preparations from C. thermocellum DSM1237 were obtained by affinity digestion and purification method with modifications [42, 43]. Culture supernatant of 1 L well-grown C. thermocellum culture was spun down and incubated with 100 mg/L phosphoric acid swollen cellulose (PASC) overnight at 4 °C. Cellulosomes bound to amorphous cellulose were collected by centrifugation (13,000 rpm, 15 min, 4 °C) and resuspended in 20 mL dialysis buffer (50 mM Tris, 5 mM CaCl2, 5 mM DTT, pH 7.0). The suspension was incubated at 60 °C and dialyzed in a Slide-A-Lyzer cassette (MW cutoff 10,000 Da) against 2 L of dialysis buffer until the suspension was clear. A pure cellulosome preparation was obtained after spinning down hydrolysis debris. Purified enzymes were concentrated with Vivaspin 500 columns (Sartorius-Stedim, Göttingen, Germany) with a cutoff of 30 to 300 kDa. Sodium azide was added to the protein preparations in a final concentration of 0.02% (w/v).

Substrates

Barley β-glucan was purchased from Megazyme (Wicklow, Ireland), Avicel, and carboxymethylcellulose (CMC) from Sigma-Aldrich (Taufkirchen, Germany). PASC was prepared from Avicel as described by Wood [44]. Substrates were used in enzymatic reactions at final concentrations of 0.5% (Barley β-glucan, CMC, PASC) or 1% (Avicel).

Enzymatic assays

All enzymatic reactions were performed under standard reaction conditions at 60 °C in a total volume of 0.5 mL. The standard reaction buffer contained final concentrations of 0.1 M MOPS, pH 6.5, 50 mM NaCl, 10 mM CaCl2, and 2 mM of Tris(2-carboxyethyl)phosphine (TCEP) as reducing agent. The activity of single cellulases was determined with barley β-glucan, CMC, PASC, or Avicel under standard reaction conditions. The activity of complexed cellulases was determined with Avicel (0.25% final concentration) with a standard enzyme load of 2 µg/mL. The enzyme kinetics were performed with 2.5% Avicel and 2 µg/mL of the enzymes. To avoid inhibition of the complexed cellulases by cellobiose, β-glucosidase (TTP0042) from Thermus thermophilus [45] was added to a final concentration of 6 µg/mL. Reducing sugar ends released from the substrates were quantified in triplicates using 3,5-dinitrosalicylic acid method [46]. One enzymatic unit liberates 1 µmol of glucose equivalent per minute.

Binding affinity studies on CipA8 and gel mobility shift assay (EMSA)

Single cellulases were bound to recombinant CipA8 by titrating different stoichiometric ratios of 1:2, 1:4, 1:6, 1:8, and 1:10 (CipA8:enzyme). The assays were performed in 30 µL reaction volume with 10 mM CaCl2 and 0.05 nmol of scaffolding protein CipA8. After 1 h of incubation at room temperature, the dockerin–cohesin interaction resulted in molecular shifts of the unbound cellulases, as visualized by gel mobility shift assay (EMSA) on 6% native gel. Non-complexed CipA8, single enzymes, and native cellulosome were used as standards.

Complex assembly

Cellulase complexes were assembled in gel filtration buffer (50 mM MOPS pH 7.3, 0.5 M NaCl, 20 mM CaCl2) for 1 h at room temperature. The complexes were assembled with a fixed concentration of the structure protein with 8 cohesins type-I and an equimolar amount of cellulases to the number of cohesins. These complexes were purified from non-complexed proteins by size-exclusion chromatography on a Superdex 200 10/300 GL column (GE Healthcare, Little Chalfont, UK) and equilibrated with gel filtration buffer. Size-exclusion chromatography was carried out on an ÄKTA Purifier (GE Healthcare, Munich, Germany). The column was developed with the same buffer at a flow rate of 0.5 mL/min. Fractions of 1 mL were collected and concentrated with Vivaspin 500 columns with a cutoff of 50 kDa. Protein concentration was determined by the BCA method [47] using bovine serum albumin as a standard.

Product analysis

The kinetics of product formation were studied on PASC and β-1,4-gluco-oligosaccharides (cello-oligosaccharides) from DP2 (cellobiose) to DP6 (cellohexaose) by thin-layer chromatography. Aliquots were taken at different time points during an enzymatic reaction, and the enzyme was inactivated by incubation at 95 °C for 15 min and subsequently stored at − 20 °C for further analysis. One to 5 µL of the aliquots was spotted on TLC silica gel 60 aluminum plates (Merck, Darmstadt, Germany) using acetonitrile/water (80:20, v/v) as the mobile phase. A mixture of DP1—DP6 cello-oligosaccharides was used as standard. Detection was performed according to De Stefanis and Ponte [48], documentation and density plot calculation was performed with ImageJ (http://imagej.net/). Glucose tetramer type B (G4G3G4G) and type C (G4G4G3G) were analyzed using a high-performance anion-exchange chromatography with pulsed amperometric detection (HPAEC-PAD) on an ICS 3000 Dionex chromatography system with a CarboPac PA1 column (4 × 250 mm) and a PA1-precolumn (4 × 50 mm). The column temperature was set to 30 °C and the injection volume was 25 µL at a flow rate of 1 mL/min. The eluent gradient for analyte separation was 7.5 mM sodium acetate with 100 mM NaOH at 0 min and increased linearly up to 100 mM sodium acetate with 100 mM NaOH at 67.5 min. After each run, the washing step consisted of 650 mM sodium acetate during 4 min and equilibration with 100 mM NaOH for 16.3 min. Carbohydrate detection based on the waveform “standard carbohydrate quad” was set to 1 Hz. Samples were diluted by factor 10 with Milli-Q water before analyzing the polysaccharide hydrolysates by HPAEC-PAD. All oligosaccharides were purchased at Megazyme, Bray, Ireland.

Structural sequence alignments and molecular docking

Multiple sequence alignments were performed with T-Coffee (http://tcoffee.crg.cat/) [49] and ESPript 3 (http://espript.ibcp.fr/ESPript/ESPript/) [50]. The sequence similarity tree was visualized with Mega 5.2 [51]. Structure prediction was performed using RaptorX (http://raptorx.uchicago.edu/), and models obtained were visualized as surface plots and amino acid overlay with the Visual Molecular Dynamics program. In silico docking experiments with cellohexaose and selected cellulases were performed with AutoDock Vina (version 1.1.2) [52] using the following procedure: Water molecules and ligands were deleted manually and structural alignments were performed using MultiSeq in the Visual Molecular Dynamics program, resulting in aligned pdb files. Aligned molecules were rotated with PyMOL (x55, y20, z-24) and saved separately. Polar hydrogens were added using AutoDockTools (version 1.5.6) [53], macromolecule was chosen under flexible, residues were selected, and rotational bonds were defined. AutoDock was performed with flexible residues (exhaustiveness 24), and the results were loaded in chimera and saved from ViewDock and converted with OpenBabel (automated bonding disabled). All input molecules were joined into a singular output molecule. Proteins and the cellohexaose sugar substrate were visualized as surface model representation.

Results

Characterizing the cellulosomal cellulases

According to genome sequence analysis and to proteomics data of the extracellular cellulosomal complex of C. thermocellum [17, 18, 33, 54], in total 24 cellulase-encoding genes were selected for subsequent enzyme characterization (Table 1). The ORFs encoding putative cellulolytic proteins were subjected to PCR-cloning or gene synthesis. Enzyme preparations were obtained by heterologous expression without the predicted N-terminal signal peptide sequence and subsequent His-tag purification (purified enzymes are summarized in Additional file 3). As the proteins were expressed with an intact type-I dockerin binding module, the binding capacity of each protein was tested on the recombinant scaffolding protein CipA8 with eight single cohesin modules. All tested proteins assembled with CipA8 via cohesin–dockerin interaction. However, the molar ratio of full stoichiometric binding varied for each enzyme (see Additional file 4).

In order to identify true β-1,4-glucanases, the degradation capability of glucose tetramer type B (G4G3G4G) and type C (G4G4G3G) was determined by HPAEC-PAD. Only Cel5-26H was specifically cleaving the β-1,3-glycosidic bond, whereas the other enzymes had no detectable activity on this type of glycosidic bond (see Additional file 5). Concomitantly, the products formed from model cello-oligosaccharides (cellotriose to cellohexaose), and activity on various cellulosic substrates were assessed (Figs. 1, 2). Unmodified substrate preparations were amorphous PASC and insoluble Avicel. In order to distinguish exo- from endo-acting cellulases, various β-glucan backbones were tested either with mixed-linkage β-1,3/1,4-glucan (barley) or side chain-modified CMC. Cleavage of these substrates is an indication for endo-acting cellulases which hydrolyze randomly at the β-1,4-linkages of the polysaccharide chain. In contrast, exo-acting cellobiohydrolases thread the cellulose molecule from its free cellulose chain end through a tunnel built by loop structures around the active site. Modified and mixed-linkage β-glucans block the enzymes’ processive activity by steric hindrance. Hence, significant activity is only observed on unmodified cellulose. The specific enzyme activities (µmol of reducing sugar ends per minute and per nmol of protein) were obtained under the optimal conditions for cellulosome activity (at 60 °C and pH 5.8; see “Methods” section and Additional file 6).

Fig. 1
figure 1

Comparison of all β-1,4-glucanases from the C. thermocellum cellulosome (listed in first column). Second column: Intermediate and final product analysis of different cellulases on various cello-oligosaccharides and PASC as substrate. The arrays show the oligosaccharide products (degree of polymerization ranging from glucose DP 1 to cellohexaose DP 6) on the Y-axis and the kinetic product shift over time (X-axis). The sugar amount detected by thin-layer chromatography is depicted as heat map representation, with relative intensities of the sugar products ranging from 1% (light gray) to 100% (black). Explanations of time points: 0.5 = 0.5 min; 2 = 2 min; 5 = 5 min; ¼ h = 15 min; 1 h = 60 min; 2 h = 120 min; on = overnight incubation. Empty fields (white) indicate that no products were formed, or products were below the detection limit of thin-layer chromatography. The pattern of protein CtCel124 is not shown due to its low activity. Third column: Activity of recombinant cellulases on various substrates (average values from triplicate measurements) at optimal cellulosome activity parameters (60 °C, pH 5.8; see Additional file 6). Fourth column: presence of glycoside hydrolase (GH) families, carbohydrate-binding modules (CBM), and Ig-like modules (Ig). For continuation of this figure, please see Fig. 2

Fig. 2
figure 2

Comparison of all β-1,4-glucanases from the C. thermocellum cellulosome (continuation of Fig. 1)

The substrate preference and sugar product spectrum of the cellulosomal cellulases vary substantially, regardless of enzyme family and module architecture. As expected, for the CBHs Cbh9A, Cel48S, and Cel9K no or very weak activities on modified substrates were observed, whereas on PASC substantial product formation was found. In contrast, Cel9D, Cel9-44J, Cel8A, and Cel5E were most active on CMC. Other proteins like Cel5L and Cel5G released the highest amount of reducing sugar ends on microcrystalline cellulose.

The apparent hydrolysis pattern of these enzymes was further studied on various cello-oligosaccharide standards and PASC with TLC over time (Figs. 1, 2). A suitable enzyme dilution was chosen to visualize the presence of all intermediate products formed during the hydrolysis reaction. To this end, we were able to identify four different product patterns. As expected, CBHs (exo-acting from the sugar ends) released specifically cellobiose as the only product over time (Cel48S, Cel9K, and Cbh9A). In contrast, endo-acting β-1,4-glucanases showed a more diverse product pattern. On PASC, the apparent random cleavage mode of non-processive endoglucanases (EG) is indicated by the formation of diverse cello-oligosaccharides and longer chain dextrins like cellopentaose (DP ≥ 5) at the beginning of the hydrolysis reaction with no preferred product at any time. This pattern is found with different GH family proteins such as Cel8A, Cel5E, Cel5-26H, Cel9-44J, and Cel9T. After prolonged incubation times (overnight), the final products are mainly cellobiose and cellotriose (DP 2 to DP 3). In contrast, processively acting endoglucanases are characterized by specifically cleaving off short-chain oligosaccharides of defined length (DP 2 or 4) at the beginning of the hydrolysis on PASC. This can be interpreted as an internal cut into the cellulose chain followed by a processive cleavage of even-numbered short cello-oligosaccharides before the enzyme falls off. Two different groups of processive endoglucanases can be distinguished, depending on the main product formed during hydrolysis: pEG4 and pEG2.

The cellotetraose-type processive endoglucanase (pEG4) group demonstrates cleavage and release of defined cello-oligosaccharides with DP 4 as intermediate product at the beginning of the hydrolysis reaction on the tested substrates. All the members of this group belong to glycoside hydrolase family 9. In cellobiose-type endoglucanases (pEG2), only cellobiose and small amounts of cellotriose as intermediate and final products were observed, e.g., all members of GH5 sub-family 1 (Cel5B, Cel5G, Cel5L, and Cel5O). Interestingly, the pEG2 hydrolysis pattern is also demonstrated by Cel9D, which resulted in cellobiose and a small amount of glucose as the only and final degradation products. This result was confirmed by the hydrolysis products from cello-oligosaccharides as substrate, which also produced cellobiose as major degradation product, whereas glucose was released to a lesser extent (Figs. 1, 2).

Role of endoglucanase processivity in synthetic protein complexes

The presence and selective attachment of single enzymatic functions to the scaffolding protein has been discussed to be the key factor for effective cellulose degradation by the native cellulosome and synthetic multi-enzyme complexes [43, 55]. The discovery of different processivity groups of cellulases (Figs. 1, 2) prompted us to construct di-, tri-, and tetravalent mini-cellulosomal complexes to test their efficiency and synergism on microcrystalline cellulose: Different combinations of endo- and exo-active cellulases were bound to the scaffolding carrier protein CipA8 in equal stoichiometric loadings via the specific dockerin–cohesin protein–protein interaction. Upon loading of the scaffolding protein to saturation (all binding positions are bound by single cellulases), the high-molecular weight fractions were separated from unbound single cellulases by size-exclusion chromatography and pooled. The complex activity resulted in the release of soluble reducing sugar products from the insoluble substrate Avicel (Fig. 3a). As a result, after 2 days of incubation at 60 °C, divalent cellulase combinations of endo/exo as well as exo/exo components (basic complex with cellulases SK, KA, and SA, meaning Cel48S/Cel9K, Cel9K/Cbh9A, and Cel48S/Cbh9A complexes, respectively) showed the lowest activities with less than 500 µM reducing sugar end products per reaction, as compared with trivalent complexes comprising two CBH enzymes (one from the reducing sugar end and one from the non-reducing end-type, Cel48S and Cel9K, respectively) and one endoglucanase. To further analyze the impact of the type of endoglucanase incorporated in the complex, we further compared the presence of non-processive endoglucanases (complex SK with Cel5-26H) with processive ones (complex SK with Cel9R, Cel9L, and Cel9D, respectively). Interestingly the complexes containing the processive endoglucanases Cel5L (complex SKL) gave the best result (up to 736.6 µM) from all trivalent mini-cellulosomal complexes. Even a complex of four different enzymes (SKAR) including two different endoglucanase functions (non-processive and cellotetraose-releasing endoglucanase, whereas a pEG2-type was missing) did not result in higher productivity (566.2 µM).

Fig. 3
figure 3

Hydrolytic efficiency of multi-enzyme complexes on Avicel as substrate. a End-point activity of 1 µg of the enzymes on 0.25% Avicel after 42 h in dependence of the endoglucanase functions present in the complex. The bars represent amount of reducing sugar ends (glucose equivalents) as average values from biological replicates (at least duplicate measurements with standard deviations represented as × 1 SD). The endoglucanase product pattern present (+) or absent (−) in the complex was non-processive (EG) or processive with cellobiose (pEG2) or cellotetraose (pEG4) as intermediate or main product. The complex “all EG types” consists of 9 different enzymes, whereas each cellulase function is present in this complex (cellobiohydrolases, non-processive endoglucanases, members of pEG2 and pEG4, respectively). Each of the complexed mixtures comprised equal stoichiometric loading and statistical distribution of eight single enzymes on CipA8 by cohesin–dockerin protein interaction. The enzyme complexes were purified by gel filtration to exclude the impact of unbound single cellulases. As controls, the activity of complexed and non-complexed enzyme extracts of C. thermocellum mutant SM901 [41] is shown together with the native cellulosome. Abbreviations: Cel8A (A), Cel9D (D), Cel5-26H (H), Cel9-44J (J), Cel9K (K), Cel5L (L), Cel9R (R), and Cel48S (S). b Enzyme kinetics of cellulosomal complexes on 2.5% Avicel. c Electrophoretic mobility shift showing the binding capacity of recombinant scaffolding protein CipA8 made possible by its eight cohesin binding modules. Complex formation by cohesin–dockerin interaction is visible by up-shifted protein bands in the native gel. 10 µM of CipA8 was titrated with 80 µM of a nonavalent cellulase mixture (all EG types + CipA8) for statistically binding all free cohesin modules. As another control, the SM901 enzyme extract was also completely bound (SM901 + CipA8). The 6% native PAGE gel was stained with Coomassie R-250

In order to analyze the influence of more endoglucanase functions on a complex, we designed a fully synthetic cellulosomal nonavalent complex (“all EG types”) containing 25% of Cel48S and Cel9K, Cbh9A, Cel8A, Cel9Q, Cel9T (each 12.5%, corresponding a stoichiometric binding to one cohesin module) and a mixture of Cel5G, Cel9R, and Cel9-44J (each 4.2%) which most closely resembles the cellulase composition of the native cellulosome complex. The fully recombinant enzyme mixture (termed “all EG types”) contains all different classes of endoglucanase functions and showed on average 52.6 ± 1.4% of the activity of the native cellulosome enzyme preparation from C. thermocellum on 2.5% microcrystalline cellulose (Fig. 3b). The single enzyme components as well as the native enzyme mixture from C. thermocellum mutant SM901 assembled with recombinant scaffolding protein CipA8 to form enzyme complexes, whereas the stoichiometric binding capacity equals 1:8 (CipA8: single enzyme ratio) (Fig. 3c).

Comparative sequence analysis and structural modeling

In order to predict certain sequence signatures that trigger the processive status of the endoglucanases, the module architecture, the presence of carbohydrate binding and other modules as well as tertiary/secondary structure prediction and sugar-binding moieties was compared. The multiple sequence alignment analysis of all 24 full-length protein sequences (including catalytic core and adjacent modules like CBMs, immunoglobulin-like modules, and others) could not differentiate between the apparent processivity status and the product specificity between the cellulosomal endoglucanases (data not shown). Noteworthy, this is also the case for the subset of cellulases belonging to GH9 which represent the majority of all cellulosomal cellulases (13 out of 24 cellulases in total).

Structure-based multiple sequence alignments and molecular modeling analysis of representative GH9 catalytic modules with different product spectra were performed: cellobiohydrolase Cbh9A [56], non-processive endoglucanase Cel9T [57], and the processive endoglucanase Cel9D [58]. The catalytic module of Cel9A (formerly called E4) from Thermobifida fusca (formerly known as Thermomonospora fusca) was chosen as it has been intensively characterized and as it shares relatively high sequence identity to Cel9F (57.2%) and Cel9T (35.9%), respectively [7, 59] (Fig. 4). The comparative analysis revealed 12 α-helices forming the (α/α)6-barrel fold typical for GH9 catalytic modules and amino acid residues that may be involved in substrate-binding, according to available structural data [56,57,58,59] and molecular docking simulations (Fig. 5). The active site comprises the conserved catalytic triad of the nucleophile/base (two aspartic acid residues in the DAGD-motif) and glutamic acid as catalytic proton donor. Sugar-binding moieties that are conserved in the sequence alignment share aromatic properties (tyrosine Y, tryptophan W) or are amino acids with electrically charged side chains (arginine R, histidine H, aspartic acid D, and glutamic acid E). The number of predicted substrate-binding residues varies between the Cbh9A with 14 residues, followed by Cel9F and Cel9T (12 residues each) and Cel9D comprising 10 interaction partners. Subsites G553, Y555, W616, W678, H737, and R739 of cellobiohydrolase Cbh9A are conserved among the compared structures covering the interactions of carbohydrate-binding positions + 2 to − 3 relative to the glycosidic linkage cleaved, while W473, L476, G546, S547, and T797 are unique sites with binding to position + 2 to − 2 cello-oligosaccharides. One of two loop regions that confer exo-activity in Cbh9A comprises E606 as another binding residue. In contrast, aromatic residues needed for interaction with larger sugars at position − 3 and − 4 were found to be present in endo-mode acting enzymes only, but are absent in Cbh9A and Cel9D. As putative binding residues we identified the residues W281, Y343 for Cel9F and W314, Y395 for Cel9T, respectively. Both aromatic amino acids are strictly conserved in this particular position among all other cellulosomal endoglucanases of family 9 (data not shown). In similarity to Cbh9A, the cellobiose-releasing processive endo-acting cellulase Cel9D lacks these aromatic residues binding to cello-oligosaccharides at positions − 3 and − 4, whereas unique aromatic sugar-binding residues are predicted, e.g., F276 instead of histidine at subsite + 2 and W560 instead of tyrosine at subsite − 2. Again, all other endoglucanases including Cel9T and Cel9F share conserved histidines or tyrosines at these particular positions as a common feature.

Fig. 4
figure 4

Structure-based multiple sequence alignment of GH9 family catalytic modules of four C. thermocellum cellulases: Cel9D, Cbh9A, Cel9T, and Cel9F. α-Helices (α- and η-helices), β-sheets, and loops in Cbh9A are indicated and numbered above the sequences as squiggles and arrows, respectively. Strict α-turns are indicated with TTT, strict β-turns with TT. The catalytic triad in the active sites is indicated with asterisks. Amino acids of the endoglucanase TfCel9A from Thermobifida fusca known to be involved in substrate-binding [59, 60] are shown as black triangles, those identified from cellobiohydrolase Cbh9A [56] are marked as gray triangles. The numbers below indicate the corresponding cello-oligosaccharide positions reported to interact/bind. Carbohydrate positions + 1 and + 2 are the expected product sites. Loop regions conferring exo-activity of Cbh9A are highlighted in light blue [56]

Fig. 5
figure 5

Structural comparison of four different catalytic clefts from glycoside hydrolase GH9 cellulases. The figure depicts cellobiohydrolase Cbh9A (PDB structure 1RQ5), processive endoglucanase pEG2 (Cel9D, PDB 1CLC), pEG4 (Cel9F, predicted structure), and non-processive endoglucanase Cel9T (PDB accession 2YIK) as gray surface plots and their corresponding sugar-binding moieties (red sticks) and catalytic triad (in blue). Cellohexaose (Glc6) was taken from PDB 7CEL. Numbers in black depict cello-oligosaccharide positions (+ 2 to − 2) according to the nomenclature for sugar-binding subsites [61] in the catalytic cleft from protein–ligand interaction data for Cbh9A according to [56], Cel9T [57], Cel9A from T. fusca [59] and the structural sequence alignment (Fig. 4)

Discussion

The recalcitrant nature and heterogeneous physical structure of cellulose selects for a varied arsenal of enzymatic machinery to efficiently degrade this kind of biomass. The native cellulosome of C. thermocellum is a model for co-localization of single enzymes on carrier proteins for synergistic activity on crystalline cellulose. The steric proximity of different enzyme classes seems to be the key feature of the cellulosomal system, inspiring researchers to systematically study and develop modified in vitro cellulase complexes [16, 43, 55, 62, 63]. From all of the more than 70 identified enzyme components identified via genome, transcriptome, and proteome analysis, 24 different enzymes are associated with the scission of cellulosic β-1,4-glycosidic bonds by exhibiting β-1,4-glucanohydrolase activity [17,18,19, 54]. Despite hydrolyzing an identical chemical bond, these cellulases are generally distinguished by their protein fold, mode of hydrolysis, and substrate specificity, as documented in the CAZy online database of glycoside hydrolase family proteins (http://www.cazy.org/). The cellulases present in the cellulosome of C. thermocellum are found in five different GH families (Table 1), namely families 5, 8, 9, 48, and the recently identified family 124 [1, 33]. Although all single enzymes have been reported before (Table 1), the lack of standardized experimental conditions (enzyme and substrate loading) has hindered any meaningful inter-laboratory comparisons of the available biochemical data. In this study activity, parameters like temperature, buffer, and pH were chosen in accordance with the optimum reaction conditions of the native cellulosome (Additional file 6).

Analysis of intermediate product kinetics and product ratios was employed to distinguish different processivity groups with the aid of thin-layer chromatograms of all 24 cellulosomal cellulases. This approach allowed for qualitative and semiquantitative discrimination of distinct product patterns [4]. Four such pattern types were obtained: (i) cellobiohydrolases (CBHs), (ii) non-processive endoglucanases without predominant hydrolysis products (EGs), apparent processive endoglucanases with (iii) cellotetraose as the intermediate product (pEG4), and (iv) cellobiose as the major product during substrate hydrolysis (pEG2).

Cellulosomal GH9 proteins were shown to produce all four types of cellulase product patterns and seem to be the most diverse enzyme family with regard to composition of the module architecture, product spectrum, and activity mode (Figs. 1, 2). Cbh9A, Cel9K, and Cel48S are CBHs, specifically releasing cellobiose from unmodified cellulose and cellodextrins, whereas they do not efficiently hydrolyze CMC and mixed-linkage β-glucan from barley (Figs. 1, 2). The processive action of CBHs, leading to the release of cellobiose, is favored by the 180° rotation of the glucose moieties within the cellulose chain [4]. Non-processive endo-acting β-1,4-glucanases (EGs) are characterized by their indiscriminate scission of cello-oligosaccharides and an acceptance of substrates with side chain modifications or mixed-linkage substrates. Thin-layer chromatographic product analysis revealed that this endoglucanase group generates cellodextrins with no preferential hydrolysis pattern when tested on the different types of substrates. Three GH families were found to show this type of endoglucanase activity, with the highest activities seen on CMC and barley β-glucan, namely GH5 (Cel5E and Cel5H), GH8 (Cel8A), and GH9 proteins (Cel9-44J and Cel9T). The results from TLC analysis support this finding, as long-chain products (e.g., cellopentaose or larger, DP ≥ 5) which are characteristic to non-processive endoglucanases were observed in their digestion patterns.

In contrast, processively acting endoglucanases regularly show low activity on CMC and barley β-glucan. This can be explained by steric hindrance inhibiting further substrate cleavage, or by immobilization of the enzymes as carbohydrate-binding modules inhibit dissociation from the tightly bound substrate. Interestingly, about half of the cellulosomal endoglucanases produce cellotetraose as the intermediate product (i.e., pEG4-type cellulases: Cel9F, Cel9N, Cel9P, Cel9Q, Cel9R, Cel9T, Cel9U, Cel9V, Cel9W, Lec9A, and Lec9B). With the exception of Cel9P, they all share identical module architecture with a GH9 catalytic module connected to a CBM3c. A major functional role of the CBM is to decrease the enzyme dissociation constant k off by interaction of the polysaccharide chain with a diverse set of binding residues on the CBM surface [5]. In processive endoglucanases, the catalytic module is joined to a family 3c carbohydrate-binding module that is aligned with the active site cleft. The endoglucanase Cel9A from T. fusca was shown to be processive upon the presence of the CBM3c module, whereby the truncation of the binding module converted the enzyme into a non-processive endoglucanase [5, 59]. In terms of bioenergetics it seems reasonable to infer that C. thermocellum expresses a redundant and large set of processive endoglucanases, as cellotetraose was shown to be preferably assimilated during growth on cellulose [64].

The most interesting observation of this study was the detection of cellobiose as main product of the pEG2-type cellulases, which was found in endoglucanase GH5 sub-family 1 proteins (Cel5B, Cel5G, Cel5L, Cel5O) and one representative of GH9 (Cel9D). Cel5O is the only representative of cellobiose-producing endoglucanases of type pEG2 that comprises a CBM3b module. In this study, Cel5O shows characteristics of a processively active endoglucanase rather than the suggested cellobiohydrolase function that has been reported previously [36].

Of particular note is that a mixture of non-processive and processive-type endoglucanases within a nonavalent complex (n = 9 different enzymes, currently named “all EG types”), which reconstitutes the intricate cellulosome, achieved the most efficient degradation of cellulose with a recombinant enzyme complex in this study (Fig. 3a, b). A native enzyme mixture from the cipA-deficient C. thermocellum mutant SM901 [41] complexed with recombinant CipA8 reached almost the same activity as the native cellulosomal complex. These data are in accordance with previously published results, where a higher cellulolytic efficiency was observed with a more diverse complex composition [43, 63, 65, 66]. The observed diversity of the hydrolysis pattern and substrate specificity of the cellulosomal cellulases may be an adaptation of the cellulosome complex to avoid stalling (also referred as jamming) of cellulases during substrate degradation [67]. Our results therefore indicate that different endoglucanase types present in the cellulosome complex may contribute to its high efficiency in lignocellulosic biomass degradation.

Sequence and structural comparison of cellulosomal GH9 cellulases allow identifying binding residues that may interact with cello-oligosaccharide sugar moieties entering the catalytic cleft upon hydrolysis (Fig. 5; Additional file 7). The (α/α)6-barrel fold of T. fusca cellulase Cel9A, a cellobiose-producing enzyme, contains an open active site cleft and at least 9 sugar-binding subsites to bind positions + 4 to − 2 [59]. The lack of substrate-binding residues from subsites − 1 to − 4 results in weaker binding. The dissociation of the sugar chain bound to the enzyme rather than entering the empty subsites after cleavage results in decreased cellulase processivity [5, 33]. In the cellulosomal GH9 cellulases, comprising most of the pEG4 enzymes with cellotetraose as the intermediate product, conserved aromatic and electrically charged residues were identified that may correlate with the observed product formation pattern: non-processive endoglucanases and pEG4 comprise additional tryptophan or tyrosine residues that were shown to bind the − 3 and − 4 sugar moieties, and are absent in the CBHs Cel9K and Cbh9A and in the pEG2 enzyme Cel9D. These additional binding subsites may explain the production of longer oligosaccharide products during hydrolysis (such as DP 4) by binding a larger portion of the cellulose chain. In turn, the presence of more binding residues at the − 2 to + 2 subsites may result in an increased processivity via higher affinity to the sugar chain after cleavage. From molecular docking models, this stronger binding capacity causes conformational changes to the cello-oligosaccharide (see Additional file 7). Indeed, Cbh9A and Cel9D share two additional amino acid positions for a tighter binding of the + 1/+ 2 subsites, specifically F276 (binding + 1) and F279 (+ 2) in Cel9D and W473 (+ 1) and L476 (+ 2) in Cbh9A, respectively. These amino acid residues are absent in the other cellulase types (pEG4 and EG) and may trigger the release of cellobiose as the main hydrolysis product (+ positions are the subsites of an enzymatically bound sugar chain that are released as products after hydrolysis).

Strikingly, structural similarities were also found. Cbh9A and Cel9D both share an immunoglobulin-like module that was shown to stabilize the catalytic module in Cbh9A [11]. In another study, the effect of a N-terminal extremity of Cel5F from S. degradans was shown to protrude into the active site of the neighboring enzyme within a trimeric quaternary structure [13], thereby influencing the substrate specificity of the cellulase. Although Cbh9A exhibits a higher sequence similarity with Cel9D than to the other cellulases (29% amino acid identity) and a similar product spectrum, Cel9D lacks the characteristic loop structure from Cbh9A, which blocks the active site after the − 2 subsite [56], thus allowing the initial endo-attack of Cel9D. Cel9D comprises less binding residues than Cbh9A which leads to a lower binding affinity for the substrate as shown by molecular docking analysis (Additional file 7). This could be due to the structure of the catalytic cleft which is flatter and broader in Cel9D than for the other glycoside family GH9 proteins.

Conclusions

From a comparative analysis of all 24 cellulosomal β-1,4-glucanases from C. thermocellum, four different product formation patterns are observed that coincide with the apparent processivity of these enzymes. The data suggest that the presence of each processivity type is necessary for peak complex activity and therefore contributes to the high efficiency of the cellulosome. Our study paves the way for the future optimization of cellulosomal complexes by supporting a deeper understanding of the synergistic action of cellulases of different processivity types. These results may help to target efficient enzyme mixtures for industrial degradation of lignocelluloses as a basis for second generation biofuels.

Abbreviations

CBH:

cellobiohydrolase

CBM:

carbohydrate-binding module

CMC:

carboxymethylcellulose

Coh:

cohesin

DP:

degree of polymerization

EG:

endoglucanase (endo-mode cellulases with no specific hydrolysis pattern)

EMSA:

gel mobility shift assay

GH:

glycoside hydrolase

HPAEC-PAD:

high-performance anion-exchange chromatography with pulsed amperometric detection

PASC:

phosphoric acid swollen cellulose

pEG2:

processive endoglucanases with cellobiose as the main product

pEG4:

processive endoglucanases with cellotetraose as intermediate product

SDS-PAGE:

sodium dodecyl sulfate–polyacrylamide gel electrophoresis

TCEP:

tris (2-carboxyethyl) phosphine

TLC:

thin-layer chromatography

References

  1. Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42(D1):D490–5.

    Article  CAS  Google Scholar 

  2. Payne CM, Knott BC, Mayes HB, Hansson H, Himmel ME, Sandgren M, et al. Fungal cellulases. Chem Rev. 2015;115(3):1308–448.

    Article  CAS  Google Scholar 

  3. Koshland DE. Stereochemistry and the mechanism of enzymatic reactions. Biol Rev. 1953;28(4):416–36.

    Article  CAS  Google Scholar 

  4. Horn SJ, Sorlie M, Vårum KM, Väljamäe P, Eijsink VGH. Chapter five—measuring processivity. In: Harry JG, editor. Methods in Enzymology, vol. 510. London: Academic Press; 2012. p. 69–95.

    Google Scholar 

  5. Wilson DB, Kostylev M. Cellulase processivity. In: Himmel EM, editor. Biomass conversion: methods and protocols. Totowa: Humana Press; 2012. p. 93–9.

    Chapter  Google Scholar 

  6. Rouvinen J, Bergfors T, Teeri T, Knowles J, Jones T. Three-dimensional structure of cellobiohydrolase II from Trichoderma reesei. Science. 1990;249(4967):380–6.

    Article  CAS  Google Scholar 

  7. Irwin D, Shin D-H, Zhang S, Barr BK, Sakon J, Karplus PA, et al. Roles of the catalytic domain and two cellulose binding domains of Thermomonospora fusca E4 in cellulose hydrolysis. J Bacteriol. 1998;180(7):1709–14.

    CAS  Google Scholar 

  8. Parsiegla G, Belaïch A, Belaïch JP, Haser R. Crystal structure of the cellulase Cel9M enlightens structure/function relationships of the variable catalytic modules in glycoside hydrolases. Biochemistry. 2002;41(37):11134–42.

    Article  CAS  Google Scholar 

  9. Watson BJ, Zhang H, Longmire AG, Moon YH, Hutcheson SW. Processive endoglucanases mediate degradation of cellulose by Saccharophagus degradans. J Bacteriol. 2009;191(18):5697–705.

    Article  CAS  Google Scholar 

  10. Zverlov VV, Velikodvorskaya GA, Schwarz WH. Two new cellulosome components encoded downstream of cell in the genome of Clostridium thermocellum: the non-processive endoglucanase CeIN and the possibly structural protein CseP. Microbiology (UK). 2003;149:515–24.

    Article  CAS  Google Scholar 

  11. Kataeva IA, Uversky VN, Brewer JM, Schubot F, Rose JP, Wang BC, et al. Interactions between immunoglobulin-like and catalytic modules in Clostridium thermocellum cellulosomal cellobiohydrolase CbhA. Protein Eng Des Sel. 2004;17(11):759–69.

    Article  CAS  Google Scholar 

  12. Santos CR, Paiva JH, Sforça ML, Neves JL, Navarro RZ, Cota J, Akao PK, Hoffmam ZB, Meza AN, Smetana JH, Nogueira ML, Polikarpov I, Xavier-Neto J, Squina FM, Ward RJ, Ruller R, Zeri AC, Murakami MT. Dissecting structure-function-stability relationships of a thermostable GH5-CBM3 cellulase from Bacillus subtilis 168. Biochem J. 2012;441(1):95–104. doi:10.1042/BJ20110869.

    Article  CAS  Google Scholar 

  13. Lafond M, Sulzenbacher G, Freyd T, Henrissat B, Berrin JG, Garron ML. The quaternary structure of a glycoside hydrolase dictates specificity toward β-glucans. J Biol Chem. 2016;291(13):7183–94. doi:10.1074/jbc.M115.695999 (Epub 2016 Jan 11).

    Article  CAS  Google Scholar 

  14. Beckham GT, Stahlberg J, Knott BC, Himmel ME, Crowley MF, Sandgren M, et al. Towards a molecular-level theory of carbohydrate processivity in glycoside hydrolases. Curr Opin Biotechnol. 2014;27:96–106.

    Article  CAS  Google Scholar 

  15. Bayer EA, Kenig R, Lamed R. Adherence of Clostridium thermocellum to cellulose. J Bacteriol. 1983;156(2):818–27.

    CAS  Google Scholar 

  16. Artzi L, Bayer EA, Morais S. Cellulosomes: bacterial nanomachines for dismantling plant polysaccharides. Nat Rev Micro. 2017;15(2):83–95.

    Article  CAS  Google Scholar 

  17. Zverlov VV, Kellermann J, Schwarz WH. Functional subgenomics of Clostridium thermocellum cellulosomal genes: identification of the major catalytic components in the extracellular complex and detection of three new enzymes. Proteomics. 2005;5(14):3646–53.

    Article  CAS  Google Scholar 

  18. Raman B, Pan C, Hurst GB, Rodriguez M Jr, McKeown CK, Lankford PK, et al. Impact of pretreated switchgrass and biomass carbohydrates on Clostridium thermocellum ATCC 27405 cellulosome composition: a quantitative proteomic analysis. PLoS ONE. 2009;4(4):e5271.

    Article  Google Scholar 

  19. Wei H, Fu Y, Magnusson L, Baker JO, Maness P-C, Xu Q, et al. Comparison of transcriptional profiles of Clostridium thermocellum grown on cellobiose and pretreated yellow poplar using RNA-Seq. Front Microbiol. 2014;5:142.

    Article  Google Scholar 

  20. Yutin N, Galperin MY. A genomic update on clostridial phylogeny: Gram-negative spore-formers and other misplaced clostridia. Environ Microbiol. 2013;15(10):2631–41.

    CAS  Google Scholar 

  21. Johnson EA, Sakajoh M, Halliwell G, Madia A, Demain AL. Saccharification of complex cellulosic substrates by the cellulase system from Clostridium thermocellum. Appl Environ Microbiol. 1982;43(5):1125–32.

    CAS  Google Scholar 

  22. Sambrook J, Russell DW. Molecular cloning: a laboratory manual. 3rd ed. Coldspring: Harbour Laboratory Press; 2001.

    Google Scholar 

  23. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Meth. 2011;8(10):785–6.

    Article  CAS  Google Scholar 

  24. Kruus K, Wang WK, Ching JT, Wu JHD. Exoglucanase activities of the recombinant Clostridium thermocellum CelS, a major cellulosome component. J Bacteriol. 1995;177(6):1641–4.

    Article  CAS  Google Scholar 

  25. Kataeva I, Li XL, Chen HZ, Choi SK, Ljungdahl LG. Cloning and sequence analysis of a new cellulase gene encoding CelK, a major cellulosome component of Clostridium thermocellum: evidence for gene duplication and recombination. J Bacteriol. 1999;181(17):5288–95.

    CAS  Google Scholar 

  26. Zverlov VV, Velikodvorskaya GA, Schwarz WH, Kellermann J, Staudenbauer WL. Duplicated Clostridium thermocellum cellobiohydrolase gene encoding cellulosomal subunits S3 and S5. Appl Microbiol Biotechnol. 1999;51(6):852–9.

    Article  CAS  Google Scholar 

  27. Zverlov VV, Velikodvorskaya GV, Schwarz WH, Bronnenmeier K, Kellermann J, Staudenbauer WL. Multidomain structure and cellulosomal localization of the Clostridium thermocellum cellobiohydrolase CbhA. J Bacteriol. 1998;180(12):3091–9.

    CAS  Google Scholar 

  28. Schwarz WH, Grabnitz F, Staudenbauer WL. Properties of a Clostridium thermocellum endoglucanase produced in Escherichia coli. Appl Environ Microbiol. 1986;51(6):1293–9.

    CAS  Google Scholar 

  29. Hall J, Hazlewood GP, Barker PJ, Gilbert HJ. Conserved reiterated domains in Clostridium thermocellum endoglucanases are not essential for catalytic activity. Gene. 1988;69(1):29–38.

    Article  CAS  Google Scholar 

  30. Yague E, Beguin P, Aubert JP. Nucleotide sequence and deletion analysis of the cellulase-encoding gene celH of Clostridium thermocellum. Gene. 1990;89(1):61–7.

    Article  CAS  Google Scholar 

  31. Ahsan MM, Kimura T, Karita S, Sakka K, Ohmiya K. Cloning, DNA sequencing, and expression of the gene encoding Clostridium thermocellum cellulase CelJ, the largest catalytic component of the cellulosome. J Bacteriol. 1996;178(19):5732–40.

    Article  CAS  Google Scholar 

  32. Kurokawa J, Hemjinda E, Arai T, Kimura T, Sakka K, Ohmiya K. Clostridium thermocellum cellulase CelT, a family 9 endoglucanase without an Ig-like domain or family 3c carbohydrate-binding module. Appl Microbiol Biotechnol. 2002;59(4–5):455–61.

    CAS  Google Scholar 

  33. Bras JLA, Cartmell A, Carvalho ALM, Verze G, Bayer EA, Vazana Y, et al. Structural insights into a unique cellulase fold and mechanism of cellulose hydrolysis. Proc Natl Acad Sci USA. 2011;108(13):5237–42.

    Article  CAS  Google Scholar 

  34. Grepinet O, Beguin P. Sequence of the cellulase gene of Clostridium thermocellum coding for endoglucanase-b. Nucleic Acids Res. 1986;14(4):1791–9.

    Article  CAS  Google Scholar 

  35. Lemaire M, Beguin P. Nucleotide sequence of the celG gene of Clostridium thermocellum and characterization of its product, endoglucanase CelG. J Bacteriol. 1993;175(11):3353–60.

    Article  CAS  Google Scholar 

  36. Zverlov VV, Velikodvorskaya GA, Schwarz WH. A newly described cellulosomal cellobiohydrolase, CelO, from Clostridium thermocellum: investigation of the exo-mode of hydrolysis, and binding capacity to crystalline cellulose. Microbiology (UK). 2002;148:247–55.

    Article  CAS  Google Scholar 

  37. Joliff G, Beguin P, Aubert JP. Nucleotide sequence of the cellulase gene celD encoding endoglucanase-D of Clostridium thermocellum. Nucleic Acids Res. 1986;14(21):8605–13.

    Article  CAS  Google Scholar 

  38. Navarro A, Chebrou MC, Beguin P, Aubert JP. Nucleotide sequence of the cellulase gene celF of Clostridium thermocellum. Res Microbiol. 1991;142(9):927–36.

    Article  CAS  Google Scholar 

  39. Arai T, Ohara H, Karita S, Kimura T, Sakka K, Ohmiya K. Sequence of celQ and properties of CelQ, a component of the Clostridium thermocellum cellulosome. Appl Microbiol Biotechnol. 2001;57(5–6):660–6.

    Article  CAS  Google Scholar 

  40. Zverlov VV, Schantz N, Schwarz WH. A major new component in the cellulosome of Clostridium thermocellum is a processive endo-beta-1,4-glucanase producing cellotetraose. FEMS Microbiol Lett. 2005;249(2):353–8.

    Article  CAS  Google Scholar 

  41. Zverlov VV, Klupp M, Krauss J, Schwarz WH. Mutations in the scaffoldin gene, cipA, of Clostridium thermocellum with impaired cellulosome formation and cellulose hydrolysis: insertions of a new transposable element, IS1447, and implications for cellulase synergism on crystalline cellulose. J Bacteriol. 2008;190(12):4321–7.

    Article  CAS  Google Scholar 

  42. Morag E, Bayer EA, Lamed R. Affinity digestion for the near-total recovery of purified cellulosome from Clostridium thermocellum. Enzyme Microb Technol. 1992;14(4):289–92.

    Article  CAS  Google Scholar 

  43. Krauss J, Zverlov VV, Schwarz WH. In vitro reconstitution of the complete Clostridium thermocellum cellulosome and synergistic activity on crystalline cellulose. Appl Environ Microbiol. 2012;78(12):4301–7.

    Article  CAS  Google Scholar 

  44. Wood TM. Preparation of crystalline, amorphous, and dyed cellulase substrates. Methods in enzymology, vol. 160. London: Academic Press; 1988. p. 19–25.

    Google Scholar 

  45. Ohta T, Tokishita S-I, Imazuka R, Mori I, Okamura J, Yamagata H. β-Glucosidase as a reporter for the gene expression studies in Thermus thermophilus and constitutive expression of DNA repair genes. Mutagenesis. 2006;21(4):255–60.

    Article  CAS  Google Scholar 

  46. Wood TM, Bhat KM. Methods for measuring cellulase activities. Methods in enzymology, vol. 160. London: Academic Press; 1988. p. 87–112.

    Google Scholar 

  47. Smith PK, Krohn RI, Hermanson GT, Mallia AK, Gartner FH, Provenzano MD, et al. Measurement of protein using bicinchoninic acid. Anal Biochem. 1985;150(1):76–85.

    Article  CAS  Google Scholar 

  48. Destefanis VA, Ponte JG. Separation of sugars by thin-layer chromatography. J Chromatogr. 1968;34(1):116–20.

    Article  Google Scholar 

  49. Notredame C, Higgins DG, Heringa J. T-coffee: a novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;302(1):205–17.

    Article  CAS  Google Scholar 

  50. Robert X, Gouet P. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 2014;42(W1):W320–4.

    Article  CAS  Google Scholar 

  51. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.

    Article  CAS  Google Scholar 

  52. Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31(2):455–61.

    CAS  Google Scholar 

  53. Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J Comput Chem. 2009;30(16):2785–91.

    Article  CAS  Google Scholar 

  54. Gold ND, Martin VJJ. Global view of the Clostridium thermocellum cellulosome revealed by quantitative proteomic analysis. J Bacteriol. 2007;189(19):6787–95.

    Article  CAS  Google Scholar 

  55. Stern J, Kahn A, Vazana Y, Shamshoum M, Moraïs S, Lamed R, et al. Significance of relative position of cellulases in designer cellulosomes for optimized cellulolysis. PLoS ONE. 2015;10(5):e0127326.

    Article  Google Scholar 

  56. Schubot FD, Kataeva IA, Chang J, Shah AK, Ljungdahl LG, Rose JP, et al. Structural basis for the exocellulase activity of the cellobiohydrolase CbhA from Clostridium thermocellum. Biochemistry. 2004;43(5):1163–70.

    Article  CAS  Google Scholar 

  57. Kesavulu MM, Tsai JY, Lee HL, Liang PH, Hsiao CD. Structure of the catalytic domain of the Clostridium thermocellum cellulase CelT. Acta Crystallogr Sect D-Biol Crystallogr. 2012;68:310–20.

    Article  CAS  Google Scholar 

  58. Juy M, Amit AG, Alzari PM, Poljak RJ, Claeyssens M, Beguin P, et al. 3-Dimensional structure of a thermostable bacterial cellulase. Nature. 1992;357(6373):89–91.

    Article  CAS  Google Scholar 

  59. Sakon J, Irwin D, Wilson DB, Karplus PA. Structure and mechanism of endo/exocellulase E4 from Thermomonospora fusca. Nat Struct Biol. 1997;4(10):810–8.

    Article  CAS  Google Scholar 

  60. Sakon J, Karplus PA. Structure and mechanism of endo/exocellulase E4 from Thermomonospora fusca. Biophys J. 1998;74(2):A255.

    Google Scholar 

  61. Davies GJ, Wilson KS, Henrissat B. Nomenclature for sugar-binding subsites in glycosyl hydrolases. Biochem J. 1997;321(2):557–9.

    Article  CAS  Google Scholar 

  62. Fierobe H-P, Bayer EA, Tardif C, Czjzek M, Mechaly A, Bélaı̈ch A, et al. Degradation of cellulose substrates by cellulosome chimeras: substrate targeting versus proximity of enzyme components. J Biol Chem. 2002;277(51):49621–30.

    Article  CAS  Google Scholar 

  63. Hirano K, Kurosaki M, Nihei S, Hasegawa H, Shinoda S, Haruki M, et al. Enzymatic diversity of the Clostridium thermocellum cellulosome is crucial for the degradation of crystalline cellulose and plant biomass. Sci Rep. 2016;6:35709.

    Article  CAS  Google Scholar 

  64. Zhang Y-HP, Lynd LR. Cellulose utilization by Clostridium thermocellum: bioenergetics and hydrolysis product assimilation. Proc Natl Acad Sci USA. 2005;102(20):7321–5.

    Article  CAS  Google Scholar 

  65. Hirano K, Nihei S, Hasegawa H, Haruki M, Hirano N. Stoichiometric assembly of cellulosome generates maximum synergy for the degradation of crystalline cellulose, as revealed by in vitro reconstitution of the Clostridium thermocellum cellulosome. Appl Environ Microbiol. 2015;81(14):4756–66.

    Article  CAS  Google Scholar 

  66. Fendri I, Tardif C, Fierobe H-P, Lignon S, Valette O, Pagès S, et al. The cellulosomes from Clostridium cellulolyticum. FEBS J. 2009;276(11):3076–86.

    Article  CAS  Google Scholar 

  67. Igarashi K, Uchihashi T, Koivula A, Wada M, Kimura S, Okamoto T, et al. Traffic jams reduce hydrolytic efficiency of cellulase on cellulose surface. Science. 2011;333(6047):1279–82.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

VVZ, WHS, and BL planned and designed the research. CH, BL, FB, KD, RS, AR, and MMe performed the experiments. MMo and CH performed the structural modeling. BL, CH, FB, and VVZ analyzed the data. BL, CH, VVZ, WHS, WL, and SG wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors acknowledge English proofreading of the manuscript by Arman Schwarz.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its Additional files.

Consent for publication

Not applicable.

Ethical approval and consent to participate

Not applicable.

Funding

Publication of this work was supported by the German Research Foundation (DFG) and the Technische Universität München within the funding program Open Access Publishing. This work was supported by the German Federal Ministry of Education and Research (BMBF) within the National Research Strategy Bioeconomy 2030 and the GO-Bio funding measure (research Grant Numbers 0316147 and 031A383).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Benedikt Leis or Vladimir V. Zverlov.

Additional files

Additional file 1. Primer sequences used in this study.

13068_2017_928_MOESM2_ESM.docx

Additional file 2. Schematic representation of functional modules and primary amino acid sequence of synthesized scaffolding protein A from C. thermocellum.

Additional file 3. SDS-PAGE gel documentation of all proteins used in this study after the purification process.

Additional file 4. Binding properties of cellulosomal cellulases on the recombinant scaffolding protein CipA8.

13068_2017_928_MOESM5_ESM.docx

Additional file 5. Product degradation pattern of glucose tetramers with selected cellulosomal cellulases using HPAED-PAD.

Additional file 6. Purification and characterization of native cellulosome complex from C. thermocellum.

Additional file 7. Molecular docking of cellohexaose in the catalytic cleft of selected cellulosomal cellulases.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Leis, B., Held, C., Bergkemper, F. et al. Comparative characterization of all cellulosomal cellulases from Clostridium thermocellum reveals high diversity in endoglucanase product formation essential for complex activity. Biotechnol Biofuels 10, 240 (2017). https://doi.org/10.1186/s13068-017-0928-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-017-0928-4

Keywords