Expression, crystal structure and cellulase activity of the thermostable cellobiohydrolase Cel7A from the fungus Humicola grisea var. thermoidea

CelS-Catalyzed Processive Cellulose Degradation and Cellobiose Extraction for the Production of Bioethanol

Bacterial cellulase enzymes are potent candidates for the efficient production of bioethanol, a promising alternative to fossil fuels, from cellulosic biomass. These enzymes catalyze the breakdown of cellulose in plant biomass into simple sugars and then to bioethanol. In the absence of the enzyme, the cellulosic biomass is recalcitrant to decomposition due to fermentation-resistant lignin and pectin coatings on the cellulose surface, which make them inaccessible for hydrolysis. Cellobiohydrolase CelS is a microbial enzyme that binds to cellulose fiber and efficiently cleaves it into a simple sugar (cellobiose) by a repeated processive chopping mechanism. The two contributing factors to the catalytic reaction rate and the yield of cellobiose are the efficient product expulsion from the product binding site of CelS and the movement of the substrate or cellulose chain into the active site. Despite progress in understanding product expulsion in other cellulases, much remains to be understood about the molecular mechanism of processive action of these enzymes. Here, nonequilibrium molecular dynamics simulations using suitable reaction coordinates are carried out to investigate the energetics and mechanism of the substrate dynamics and product expulsion in CelS. The calculated free energy barrier for the product expulsion is three times lower than that for the processive action indicating that product removal is relatively easier and faster than the sliding of the substrate to the catalytic active site. The water traffic near the active site in response to the product expulsion and the processive action is also explored.

Metagenomic mining and structure-function studies of a hyper-thermostable cellobiohydrolase from hot spring sediment

Enzymatic breakdown is an attractive cellulose utilisation method with a low environmental load. Its high temperature operation could promote saccharification and lower contamination risk. Here we report a hyper-thermostable cellobiohydrolase (CBH), named HmCel6A and its variant HmCel6A-3SNP that were isolated metagenomically from hot spring sediments and expressed in Escherichia coli. They are classified into glycoside hydrolases family 6 (GH6). HmCel6A-3SNP had three amino acid replacements to HmCel6A (P88S/L230F/F414S) and the optimum temperature at 95 °C, while HmCel6A did it at 75 °C. Crystal structure showed conserved features among GH6, a (β/α)₈-barrel core and catalytic residues, and resembles TfCel6B, a bacterial CBH II of Thermobifida fusca, that had optimum temperature at 60 °C. From structure-function studies, we discuss unique structural features that allow the enzyme to reach its high thermostability level, such as abundance of hydrophobic and charge-charge interactions, characteristic metal bindings and disulphide bonds. Moreover, structure and surface plasmon resonance analysis with oligosaccharides suggested that the contribution of an additional tryptophan located at the tunnel entrance could aid in substrate recognition and thermostability. These results may help to design efficient enzymes and saccharification methods for cellulose working at high temperatures.

Bacteria from hot springs are known for highly thermostable enzymes, which may have industrial potential. Here, a unique thermostable cellobiohydrolase is reported that can breakdown cellulose at temperature up to 95 degrees Celsius.

Insights into the Role of the Fungal Community in Variations of the Antibiotic Resistome in the Soil Collembolan Gut Microbiome

Fertilization is known to affect antibiotic-resistance gene (ARG) patterns in the soil, even in the gut of soil fauna. Here, we conducted a microcosm experiment to investigate differences of effects of different fertilizers on collembolan gut ARG profiles and to further explore the microecological mechanisms that cause the differences. Although fertilization increased the abundance of ARGs, compared with the conventional manure, the application of antibiotic-reduced manure and vermicompost all curbed the enrichment of ARGs in the gut of collembolans. The results of the structural equation model revealed that changes in the microbial community caused by fertilizations have an important contribution to variations in the ARGs. We further found that the fungal community, like bacterial community, is also an important driver of ARG patterns in the collembolan gut. The fungi belonging to Dokmaia and Talaromyces were significantly correlated with the ARGs in the gut of collembolans. In addition, the application of vermicompost significantly increased the abundance of agricultural beneficial microbes in the soil environment. Together, our results provide an insight into the role of the fungal community on ARG patterns in the soil collembolan gut microbiome and highlight environmental friendliness of vermicomposting.

Machine learning reveals sequence-function relationships in family 7 glycoside hydrolases

Family 7 glycoside hydrolases (GH7) are among the principal enzymes for cellulose degradation in nature and industrially. These enzymes are often bimodular, including a catalytic domain and carbohydrate-binding module (CBM) attached via a flexible linker, and exhibit an active site that binds cello-oligomers of up to ten glucosyl moieties. GH7 cellulases consist of two major subtypes: cellobiohydrolases (CBH) and endoglucanases (EG). Despite the critical importance of GH7 enzymes, there remain gaps in our understanding of how GH7 sequence and structure relate to function. Here, we employed machine learning to gain data-driven insights into relationships between sequence, structure, and function across the GH7 family. Machine-learning models, trained only on the number of residues in the active-site loops as features, were able to discriminate GH7 CBHs and EGs with up to 99% accuracy, demonstrating that the lengths of loops A4, B2, B3, and B4 strongly correlate with functional subtype across the GH7 family. Classification rules were derived such that specific residues at 42 different sequence positions each predicted the functional subtype with accuracies surpassing 87%. A random forest model trained on residues at 19 positions in the catalytic domain predicted the presence of a CBM with 89.5% accuracy. Our machine learning results recapitulate, as top-performing features, a substantial number of the sequence positions determined by previous experimental studies to play vital roles in GH7 activity. We surmise that the yet-to-be-explored sequence positions among the top-performing features also contribute to GH7 functional variation and may be exploited to understand and manipulate function.

Identification and Characterization of a Novel Hyperthermostable Bifunctional Cellobiohydrolase- Xylanase Enzyme for Synergistic Effect With Commercial Cellulase on Pretreated Wheat Straw Degradation

The novel cellobiohydrolase gene ctcel7 was identified from Chaetomium thermophilum, and its recombinant protein CtCel7, a member of glycoside hydrolase family 7, was heterologously expressed in Pichia pastoris and biochemically characterized. Compared with commercial hydrolases, purified CtCel7 exhibited superior bifunctional cellobiohydrolase and xylanase activities against microcrystalline cellulose and xylan, respectively, under optimal conditions of 60°C and pH 4.0. Moreover, CtCel7 displayed remarkable thermostability with over 90% residual activity after heat (60°C) treatment for 180 min. CtCel7 was insensitive to most detected cations and reagents and preferentially cleaved the β-1,4-glycosidic bond to generate oligosaccharides through the continuous saccharification of lignocellulosic substrates, which are crucial for various practical applications. Notably, the hydrolysis effect of a commercial cellulase cocktail on pretreated wheat straw was substantively improved by its combination with CtCel7. Taken together, these excellent properties distinguish CtCel7 as a robust candidate for the biotechnological production of biofuels and biobased chemicals.

The dissociation mechanism of processive cellulases

Cellulase enzymes deconstruct recalcitrant cellulose into soluble sugars, making them a biocatalyst of biotechnological interest for use in the nascent lignocellulosic bioeconomy. Cellobiohydrolases (CBHs) are cellulases capable of liberating many sugar molecules in a processive manner without dissociating from the substrate. Within the complete processive cycle of CBHs, dissociation from the cellulose substrate is rate limiting, but the molecular mechanism of this step is unknown. Here, we present a direct comparison of potential molecular mechanisms for dissociation via Hamiltonian replica exchange molecular dynamics of the model fungal CBH, Trichoderma reesei Cel7A. Computational rate estimates indicate that stepwise cellulose dethreading from the binding tunnel is 4 orders of magnitude faster than a clamshell mechanism, in which the substrate-enclosing loops open and release the substrate without reversing. We also present the crystal structure of a disulfide variant that covalently links substrate-enclosing loops on either side of the substrate-binding tunnel, which constitutes a CBH that can only dissociate via stepwise dethreading. Biochemical measurements indicate that this variant has a dissociation rate constant essentially equivalent to the wild type, implying that dethreading is likely the predominant mechanism for dissociation.

Homology Modeling for Enzyme Design

Homology modeling is a very powerful tool in the absence of atomic structures for understanding the general fold of the enzyme, conserved residues, catalytic tunnel/pocket as well as substrate and product binding sites. This information is useful for structure-assisted enzyme design approach for the development of robust enzymes especially for industrial applications.

Engineering enhanced cellobiohydrolase activity

Glycoside Hydrolase Family 7 cellobiohydrolases (GH7 CBHs) catalyze cellulose depolymerization in cellulolytic eukaryotes, making them key discovery and engineering targets. However, there remains a lack of robust structure–activity relationships for these industrially important cellulases. Here, we compare CBHs from Trichoderma reesei (TrCel7A) and Penicillium funiculosum (PfCel7A), which exhibit a multi-modular architecture consisting of catalytic domain (CD), carbohydrate-binding module, and linker. We show that PfCel7A exhibits 60% greater performance on biomass than TrCel7A. To understand the contribution of each domain to this improvement, we measure enzymatic activity for a library of CBH chimeras with swapped subdomains, demonstrating that the enhancement is mainly caused by PfCel7A CD. We solve the crystal structure of PfCel7A CD and use this information to create a second library of TrCel7A CD mutants, identifying a TrCel7A double mutant with near-equivalent activity to wild-type PfCel7A. Overall, these results reveal CBH regions that enable targeted activity improvements.

Cellobiohydrolases (CBHs) are critical for natural and industrial biomass degradation but their structure–activity relationships are not fully understood. Here, the authors present the biochemical and structural characterization of two CBHs, identifying protein regions that confer enhanced CBH activity.

Biochemical and structural insights into a thermostable cellobiohydrolase from Myceliophthora thermophila

Cellobiohydrolases hydrolyze cellulose, a linear polymer with glucose monomers linked exclusively by β-1,4 glycosidic linkages. The widespread hydrogen bonding network tethers individual cellulose polymers forming crystalline cellulose, which prevent the access of hydrolytic enzymes and water molecules. The most abundant enzyme secreted by Myceliophthora thermophila M77 in response to the presence of biomass is the cellobiohydrolase MtCel7A, which is composed by a GH7-catalytic domain (CD), a linker, and a CBM1-type carbohydrate-binding module. GH7 cellobiohydrolases have been studied before, and structural models have been proposed. However, currently available GH7 crystal structures only define separate catalytic domains and/or cellulose-binding modules and do not include the full-length structures that are involved in shaping the catalytic mode of operation. In this study, we determined the 3D structure of catalytic domain using X-ray crystallography and retrieved the full-length enzyme envelope via small-angle X-ray scattering (SAXS) technique. The SAXS data reveal a tadpole-like molecular shape with a rigid linker connecting the CD and CBM. Our biochemical studies show that MtCel7A has higher catalytic efficiency and thermostability as well as lower processivity when compared to the well-studied TrCel7A from Trichoderma reesei. Based on a comparison of the crystallographic structures of CDs and their molecular dynamic simulations, we demonstrate that MtCel7A has considerably higher flexibility than TrCel7A. In particular, loops that cover the active site are more flexible and undergo higher conformational fluctuations, which might account for decreased processivity and enhanced enzymatic efficiency. Our statistical coupling analysis suggests co-evolution of amino acid clusters comprising the catalytic site of MtCel7A, which correlate with the steps in the catalytic cycle of the enzyme.

DATABASE

The atomic coordinates and structural factors of MtCel7A have been deposited in the Protein Data Bank with accession number 5W11.

Correlation of structure, function and protein dynamics in GH7 cellobiohydrolases from Trichoderma atroviride, T. reesei and T. harzianum

BACKGROUND

The ascomycete fungus Trichoderma reesei is the predominant source of enzymes for industrial conversion of lignocellulose. Its glycoside hydrolase family 7 cellobiohydrolase (GH7 CBH) TreCel7A constitutes nearly half of the enzyme cocktail by weight and is the major workhorse in the cellulose hydrolysis process. The orthologs from Trichoderma atroviride (TatCel7A) and Trichoderma harzianum (ThaCel7A) show high sequence identity with TreCel7A, ~ 80%, and represent naturally evolved combinations of cellulose-binding tunnel-enclosing loop motifs, which have been suggested to influence intrinsic cellobiohydrolase properties, such as endo-initiation, processivity, and off-rate.

RESULTS

The TatCel7A, ThaCel7A, and TreCel7A enzymes were characterized for comparison of function. The catalytic domain of TatCel7A was crystallized, and two structures were determined: without ligand and with thio-cellotriose in the active site. Initial hydrolysis of bacterial cellulose was faster with TatCel7A than either ThaCel7A or TreCel7A. In synergistic saccharification of pretreated corn stover, both TatCel7A and ThaCel7A were more efficient than TreCel7A, although TatCel7A was more sensitive to thermal inactivation. Structural analyses and molecular dynamics (MD) simulations were performed to elucidate important structure/function correlations. Moreover, reverse conservation analysis (RCA) of sequence diversity revealed divergent regions of interest located outside the cellulose-binding tunnel of Trichoderma spp. GH7 CBHs.

CONCLUSIONS

We hypothesize that the combination of loop motifs is the main determinant for the observed differences in Cel7A activity on cellulosic substrates. Fine-tuning of the loop flexibility appears to be an important evolutionary target in Trichoderma spp., a conclusion supported by the RCA data. Our results indicate that, for industrial use, it would be beneficial to combine loop motifs from TatCel7A with the thermostability features of TreCel7A. Furthermore, one region implicated in thermal unfolding is suggested as a primary target for protein engineering.

Improving the thermal stability of cellobiohydrolase Cel7A from Hypocrea jecorina by directed evolution

Secreted mixtures of Hypocrea jecorina cellulases are able to efficiently degrade cellulosic biomass to fermentable sugars at large, commercially relevant scales. H. jecorina Cel7A, cellobiohydrolase I, from glycoside hydrolase family 7, is the workhorse enzyme of the process. However, the thermal stability of Cel7A limits its use to processes where temperatures are no higher than 50 °C. Enhanced thermal stability is desirable to enable the use of higher processing temperatures and to improve the economic feasibility of industrial biomass conversion. Here, we enhanced the thermal stability of Cel7A through directed evolution. Sites with increased thermal stability properties were combined, and a Cel7A variant (FCA398) was obtained, which exhibited a 10.4 °C increase in T_m and a 44-fold greater half-life compared with the wild-type enzyme. This Cel7A variant contains 18 mutated sites and is active under application conditions up to at least 75 °C. The X-ray crystal structure of the catalytic domain was determined at 2.1 Å resolution and showed that the effects of the mutations are local and do not introduce major backbone conformational changes. Molecular dynamics simulations revealed that the catalytic domain of wild-type Cel7A and the FCA398 variant exhibit similar behavior at 300 K, whereas at elevated temperature (475 and 525 K), the FCA398 variant fluctuates less and maintains more native contacts over time. Combining the structural and dynamic investigations, rationales were developed for the stabilizing effect at many of the mutated sites.

Comparative insights into the saccharification potentials of a relatively unexplored but robust Penicillium funiculosum glycoside hydrolase 7 cellobiohydrolase

BACKGROUND

GH7 cellobiohydrolases (CBH1) are vital for the breakdown of cellulose. We had previously observed the enzyme as the most dominant protein in the active cellulose-hydrolyzing secretome of the hypercellulolytic ascomycete-Penicillium funiculosum (NCIM1228). To understand its contributions to cellulosic biomass saccharification in comparison with GH7 cellobiohydrolase from the industrial workhorse-Trichoderma reesei, we natively purified and functionally characterized the only GH7 cellobiohydrolase identified and present in the genome of the fungus.

RESULTS

There were marginal differences observed in the stability of both enzymes, with P. funiculosum (PfCBH1) showing an optimal thermal midpoint (Tm) of 68 °C at pH 4.4 as against an optimal Tm of 65 °C at pH 4.7 for T. reesei (TrCBH1). Nevertheless, PfCBH1 had an approximate threefold lower binding affinity (Km), an 18-fold higher turnover rate (kcat), a sixfold higher catalytic efficiency as well as a 26-fold higher enzyme-inhibitor complex equilibrium dissociation constant (Ki) than TrCBH1 on p-nitrophenyl-β-d-lactopyranoside (pNPL). Although both enzymes hydrolyzed cellooligomers (G2-G6) and microcrystalline cellulose, releasing cellobiose and glucose as the major products, the propensity was more with PfCBH1. We equally observed this trend during the hydrolysis of pretreated wheat straws in tandem with other core cellulases under the same conditions. Molecular dynamic simulations conducted on a homology model built using the TrCBH1 structure (PDB ID: 8CEL) as a template enabled us to directly examine the effects of substrate and products on the protein dynamics. While the catalytic triads-EXDXXE motifs-were conserved between the two enzymes, subtle variations in regions enclosing the catalytic path were observed, and relations to functionality highlighted.

CONCLUSION

To the best of our knowledge, this is the first report about a comprehensive and comparative description of CBH1 from hypercellulolytic ascomycete-P. funiculosum NCIM1228, against the backdrop of the same enzyme from the industrial workhorse-T. reesei. Our study reveals PfCBH1 as a viable alternative for CBH1 from T. reesei in industrial cellulase cocktails.

Identification of Glutaminyl Cyclase Genes Involved in Pyroglutamate Modification of Fungal Lignocellulolytic Enzymes

The breakdown of plant biomass to simple sugars is essential for the production of second-generation biofuels and high-value bioproducts. Currently, enzymes produced from filamentous fungi are used for deconstructing plant cell wall polysaccharides into fermentable sugars for biorefinery applications. A post-translational N-terminal pyroglutamate modification observed in some of these enzymes occurs when N-terminal glutamine or glutamate is cyclized to form a five-membered ring. This modification has been shown to confer resistance to thermal denaturation for CBH-1 and EG-1 cellulases. In mammalian cells, the formation of pyroglutamate is catalyzed by glutaminyl cyclases. Using the model filamentous fungus Neurospora crassa, we identified two genes (qc-1 and qc-2) that encode proteins homologous to mammalian glutaminyl cyclases. We show that qc-1 and qc-2 are essential for catalyzing the formation of an N-terminal pyroglutamate on CBH-1 and GH5-1. CBH-1 and GH5-1 produced in a Δqc-1 Δqc-2 mutant, and thus lacking the N-terminal pyroglutamate modification, showed greater sensitivity to thermal denaturation, and for GH5-1, susceptibility to proteolytic cleavage. QC-1 and QC-2 are endoplasmic reticulum (ER)-localized proteins. The pyroglutamate modification is predicted to occur in a number of additional fungal proteins that have diverse functions. The identification of glutaminyl cyclases in fungi may have implications for production of lignocellulolytic enzymes, heterologous expression, and biotechnological applications revolving around protein stability.

Pyroglutamate modification is the post-translational conversion of N-terminal glutamine or glutamate into a cyclized amino acid derivative. This modification is well studied in animal systems but poorly explored in fungal systems. In Neurospora crassa, we show that this modification takes place in the ER and is catalyzed by two well-conserved enzymes, ubiquitously conserved throughout the fungal kingdom. We demonstrate that the modification is important for the structural stability and aminopeptidase resistance of CBH-1 and GH5-1, two important cellulase enzymes utilized in industrial plant cell wall deconstruction. Many additional fungal proteins predicted in the genome of N. crassa and other filamentous fungi are predicted to carry an N-terminal pyroglutamate modification. Pyroglutamate addition may also be a useful way to stabilize secreted proteins and peptides, which can be easily produced in fungal production systems.

Biochemical and Structural Characterizations of Two Dictyostelium Cellobiohydrolases from the Amoebozoa Kingdom Reveal a High Level of Conservation between Distant Phylogenetic Trees of Life

Glycoside hydrolase family 7 (GH7) cellobiohydrolases (CBHs) are enzymes commonly employed in plant cell wall degradation across eukaryotic kingdoms of life, as they provide significant hydrolytic potential in cellulose turnover. To date, many fungal GH7 CBHs have been examined, yet many questions regarding structure-activity relationships in these important natural and commercial enzymes remain. Here, we present the crystal structures and a biochemical analysis of two GH7 CBHs from social amoeba: Dictyostelium discoideum Cel7A (DdiCel7A) and Dictyostelium purpureum Cel7A (DpuCel7A). DdiCel7A and DpuCel7A natively consist of a catalytic domain and do not exhibit a carbohydrate-binding module (CBM). The structures of DdiCel7A and DpuCel7A, resolved to 2.1 Å and 2.7 Å, respectively, are homologous to those of other GH7 CBHs with an enclosed active-site tunnel. Two primary differences between the Dictyostelium CBHs and the archetypal model GH7 CBH, Trichoderma reesei Cel7A (TreCel7A), occur near the hydrolytic active site and the product-binding sites. To compare the activities of these enzymes with the activity of TreCel7A, the family 1 TreCel7A CBM and linker were added to the C terminus of each of the Dictyostelium enzymes, creating DdiCel7ACBM and DpuCel7ACBM, which were recombinantly expressed in T. reesei DdiCel7ACBM and DpuCel7ACBM hydrolyzed Avicel, pretreated corn stover, and phosphoric acid-swollen cellulose as efficiently as TreCel7A when hydrolysis was compared at their temperature optima. The Ki of cellobiose was significantly higher for DdiCel7ACBM and DpuCel7ACBM than for TreCel7A: 205, 130, and 29 μM, respectively. Taken together, the present study highlights the remarkable degree of conservation of the activity of these key natural and industrial enzymes across quite distant phylogenetic trees of life.

IMPORTANCE

GH7 CBHs are among the most important cellulolytic enzymes both in nature and for emerging industrial applications for cellulose breakdown. Understanding the diversity of these key industrial enzymes is critical to engineering them for higher levels of activity and greater stability. The present work demonstrates that two GH7 CBHs from social amoeba are surprisingly quite similar in structure and activity to the canonical GH7 CBH from the model biomass-degrading fungus T. reesei when tested under equivalent conditions (with added CBM-linker domains) on an industrially relevant substrate.

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.

The three-dimensional structure of the cellobiohydrolase Cel7A from Aspergillus fumigatus at 1.5 Å resolution

The enzymatic degradation of plant cell-wall cellulose is central to many industrial processes, including second-generation biofuel production. Key players in this deconstruction are the fungal cellobiohydrolases (CBHs), notably those from family GH7 of the carbohydrate-active enzymes (CAZY) database, which are generally known as CBHI enzymes. Here, three-dimensional structures are reported of the Aspergillus fumigatus CBHI Cel7A solved in uncomplexed and disaccharide-bound forms at resolutions of 1.8 and 1.5 Å, respectively. The product complex with a disaccharide in the +1 and +2 subsites adds to the growing three-dimensional insight into this family of industrially relevant biocatalysts.