The crystal structure of the catalytic core domain of endoglucanase I from Trichoderma reesei at 3.6 A resolution, and a comparison with related enzymes

Genomics review of holocellulose deconstruction by aspergilli

SUMMARY

Biomass is constructed of dense recalcitrant polymeric materials: proteins, lignin, and holocellulose, a fraction constituting fibrous cellulose wrapped in hemicellulose-pectin. Bacteria and fungi are abundant in soil and forest floors, actively recycling biomass mainly by extracting sugars from holocellulose degradation. Here we review the genome-wide contents of seven Aspergillus species and unravel hundreds of gene models encoding holocellulose-degrading enzymes. Numerous apparent gene duplications followed functional evolution, grouping similar genes into smaller coherent functional families according to specialized structural features, domain organization, biochemical activity, and genus genome distribution. Aspergilli contain about 37 cellulase gene models, clustered in two mechanistic categories: 27 hydrolyze and 10 oxidize glycosidic bonds. Within the oxidative enzymes, we found two cellobiose dehydrogenases that produce oxygen radicals utilized by eight lytic polysaccharide monooxygenases that oxidize glycosidic linkages, breaking crystalline cellulose chains and making them accessible to hydrolytic enzymes. Among the hydrolases, six cellobiohydrolases with a tunnel-like structural fold embrace single crystalline cellulose chains and cooperate at nonreducing or reducing end termini, splitting off cellobiose. Five endoglucanases group into four structural families and interact randomly and internally with cellulose through an open cleft catalytic domain, and finally, seven extracellular β-glucosidases cleave cellobiose and related oligomers into glucose. Aspergilli contain, on average, 30 hemicellulase and 7 accessory gene models, distributed among 9 distinct functional categories: the backbone-attacking enzymes xylanase, mannosidase, arabinase, and xyloglucanase, the short-side-chain-removing enzymes xylan α-1,2-glucuronidase, arabinofuranosidase, and xylosidase, and the accessory enzymes acetyl xylan and feruloyl esterases.

Characterization of the Dielectric Constant in the Trichoderma reesei Cel7B Active Site

An attempt is made to evaluate the dielectric constant of the Trichoderma reesei Cel7B active site. Through kinetic measurements, the pKa value of the catalytic acid E201 is determined. Mutations (away from E201) with net charge changes are introduced to perturb the E201 pKa. It is shown that the mutation with a +1 charge change (including G225R, G230R, and A335R) decreases the pKa of E201, whereas the mutation with a -1 charge change (including Q149E, A222D, G225D, and G230D) increases the pKa. This effect is consistent with the electrostatic interaction between the changed charge and the E201 side chain. The fitting of the experimental data yields an apparent dielectric constant of 25-80. Molecular dynamics simulations with explicit water molecules indicate that the high solvent accessibility of the active site contributes largely to the high dielectric constant. ONIOM calculations show that high dielectric constant benefits the catalysis through decreasing the energy of the transition state relative to that of the enzyme substrate complex.

Mutagenesis of Trichoderma reesei endoglucanase I: impact of expression host on activity and stability at elevated temperatures

BACKGROUND

Trichoderma reesei is a key cellulase source for economically saccharifying cellulosic biomass for the production of biofuels. Lignocellulose hydrolysis at temperatures above the optimum temperature of T. reesei cellulases (~50°C) could provide many significant advantages, including reduced viscosity at high-solids loadings, lower risk of microbial contamination during saccharification, greater compatibility with high-temperature biomass pretreatment, and faster rates of hydrolysis. These potential advantages motivate efforts to engineer T. reesei cellulases that can hydrolyze lignocellulose at temperatures ranging from 60-70°C.

RESULTS

A B-factor guided approach for improving thermostability was used to engineer variants of endoglucanase I (Cel7B) from T. reesei (TrEGI) that are able to hydrolyze cellulosic substrates more rapidly than the recombinant wild-type TrEGI at temperatures ranging from 50-70°C. When expressed in T. reesei, TrEGI variant G230A/D113S/D115T (G230A/D113S/D115T Tr_TrEGI) had a higher apparent melting temperature (3°C increase in Tm) and improved half-life at 60°C (t1/2 = 161 hr) than the recombinant (T. reesei host) wild-type TrEGI (t1/2 = 74 hr at 60°C, Tr_TrEGI). Furthermore, G230A/D113S/D115T Tr_TrEGI showed 2-fold improved activity compared to Tr_TrEGI at 65°C on solid cellulosic substrates, and was as efficient in hydrolyzing cellulose at 60°C as Tr_TrEGI was at 50°C. The activities and stabilities of the recombinant TrEGI enzymes followed similar trends but differed significantly in magnitude depending on the expression host (Escherichia coli cell-free, Saccharomyces cerevisiae, Neurospora crassa, or T. reesei). Compared to N.crassa-expressed TrEGI, S. cerevisiae-expressed TrEGI showed inferior activity and stability, which was attributed to the lack of cyclization of the N-terminal glutamine in Sc_TrEGI and not to differences in glycosylation. N-terminal pyroglutamate formation in TrEGI expressed in S. cerevisiae was found to be essential in elevating its activity and stability to levels similar to the T. reesei or N. crassa-expressed enzyme, highlighting the importance of this ubiquitous modification in GH7 enzymes.

CONCLUSION

Structure-guided evolution of T. reesei EGI was used to engineer enzymes with increased thermal stability and activity on solid cellulosic substrates. Production of TrEGI enzymes in four hosts highlighted the impact of the expression host and the role of N-terminal pyroglutamate formation on the activity and stability of TrEGI enzymes.

Direct and up-close views of plant cell walls show a leading role for lignin-modifying enzymes on ensuing xylanases

BACKGROUND

A key barrier that limits the full potential of biological processes to create new, sustainable materials and fuels from plant fibre is limited enzyme accessibility to polysaccharides and lignin that characterize lignocellulose networks. Moreover, the heterogeneity of lignocellulosic substrates means that different enzyme combinations might be required for efficient transformation of different plant resources. Analytical techniques with high chemical sensitivity and spatial resolution that permit direct characterization of solid samples could help overcome these challenges by allowing direct visualization of enzyme action within plant fibre, thereby identify barriers to enzyme action.

RESULTS

In the current study, the high spatial resolution (about 30 nm) of scanning transmission X-ray microscopy (STXM), and the detection sensitivity (ppm) of time-of-flight secondary ion mass spectrometry (ToF-SIMS), were harnessed for the first time to investigate the progression of laccase, cellulase and xylanase activities through wood samples, and to evaluate complementary action between lignin-modifying and polysaccharide-degrading enzymes. In particular, complementary insights from the STXM and ToF-SIMS analyses revealed the key role of laccase in promoting xylanase activity throughout and between plant cell walls.

CONCLUSIONS

The spatial resolution of STXM clearly revealed time-dependent progression and spatial distribution of laccase and xylanase activities, whereas ToF-SIMS analyses confirmed that laccase promoted protein penetration into fibre samples, leading to an overall increase in polysaccharide degradation. Spectromicroscopic visualizations of plant cell wall chemistry allowed simultaneous tracking of changes to lignin and polysaccharide contents, which provides new possibilities for investigating the complementary roles of lignin-modifying and carbohydrate-active enzymes.

Carbohydrate-protein interactions that drive processive polysaccharide translocation in enzymes revealed from a computational study of cellobiohydrolase processivity

Translocation of carbohydrate polymers through protein tunnels and clefts is a ubiquitous biochemical phenomenon in proteins such as polysaccharide synthases, glycoside hydrolases, and carbohydrate-binding modules. Although static snapshots of carbohydrate polymer binding in proteins have long been studied via crystallography and spectroscopy, the molecular details of polysaccharide chain processivity have not been elucidated. Here, we employ simulation to examine how a cellulose chain translocates by a disaccharide unit during the processive cycle of a glycoside hydrolase family 7 cellobiohydrolase. Our results demonstrate that these biologically and industrially important enzymes employ a two-step mechanism for chain threading to form a Michaelis complex and that the free energy barrier to chain threading is significantly lower than the hydrolysis barrier. Taken with previous studies, our findings suggest that the rate-limiting step in enzymatic cellulose degradation is the glycosylation reaction, not chain processivity. Based on the simulations, we find that strong electrostatic interactions with polar residues that are conserved in GH7 cellobiohydrolases, but not in GH7 endoglucanases, at the leading glucosyl ring provide the thermodynamic driving force for polysaccharide chain translocation. Also, we consider the role of aromatic-carbohydrate interactions, which are widespread in carbohydrate-active enzymes and have long been associated with processivity. Our analysis suggests that the primary role for these aromatic residues is to provide tunnel shape and guide the carbohydrate chain to the active site. More broadly, this work elucidates the role of common protein motifs found in carbohydrate-active enzymes that synthesize or depolymerize polysaccharides by chain translocation mechanisms coupled to catalysis.

The significance of cellulolytic enzymes produced by Trichoderma in opportunistic lifestyle of this fungus

The degradation of native cellulose to glucose monomers is a complex process, which requires the synergistic action of the extracellular enzymes produced by cellulolytic microorganisms. Among fungi, the enzymatic systems that can degrade native cellulose have been extensively studied for species belonging to the genera of Trichoderma. The majority of the cellulolytic enzymes described so far have been examples of Trichoderma reesei, extremely specialized in the efficient degradation of plant cell wall cellulose. Other Trichoderma species, such as T. harzianum, T. koningii, T. longibrachiatum, and T. viride, known for their capacity to produce cellulolytic enzymes, have been isolated from various ecological niches, where they have proved successful in various heterotrophic interactions. As saprotrophs, these species are considered to make a contribution to the degradation of lignocellulosic plant material. Their cellulolytic potential is also used in interactions with plants, especially in plant root colonization. However, the role of cellulolytic enzymes in species forming endophytic associations with plants or in those existing in the substratum for mushroom cultivation remains unknown. The present review discusses the current state of knowledge about cellulolytic enzymes production by Trichoderma species and the encoding genes, as well as the involvement of these proteins in the lifestyle of Trichoderma.

Initial recognition of a cellodextrin chain in the cellulose-binding tunnel may affect cellobiohydrolase directional specificity

Cellobiohydrolases processively hydrolyze glycosidic linkages in individual polymer chains of cellulose microfibrils, and typically exhibit specificity for either the reducing or nonreducing end of cellulose. Here, we conduct molecular dynamics simulations and free energy calculations to examine the initial binding of a cellulose chain into the catalytic tunnel of the reducing-end-specific Family 7 cellobiohydrolase (Cel7A) from Hypocrea jecorina. In unrestrained simulations, the cellulose diffuses into the tunnel from the -7 to the -5 positions, and the associated free energy profiles exhibit no barriers for initial processivity. The comparison of the free energy profiles for different cellulose chain orientations show a thermodynamic preference for the reducing end, suggesting that the preferential initial binding may affect the directional specificity of the enzyme by impeding nonproductive (nonreducing end) binding. Finally, the Trp-40 at the tunnel entrance is shown with free energy calculations to have a significant effect on initial chain complexation in Cel7A.

The synergistic action of accessory enzymes enhances the hydrolytic potential of a "cellulase mixture" but is highly substrate specific

BACKGROUND

Currently, the amount of protein/enzyme required to achieve effective cellulose hydrolysis is still too high. One way to reduce the amount of protein/enzyme required is to formulate a more efficient enzyme cocktail by adding so-called accessory enzymes such as xylanase, lytic polysaccharide monooxygenase (AA9, formerly known as GH61), etc., to the cellulase mixture. Previous work has shown the strong synergism that can occur between cellulase and xylanase mixtures during the hydrolysis of steam pretreated corn stover, requiring lower protein loading to achieve effective hydrolysis. However, relatively high loadings of xylanases were required. When family 10 and 11 endo-xylanases and family 5 xyloglucanase were supplemented to a commercial cellulase mixture varying degrees of improved hydrolysis over a range of pretreated, lignocellulosic substrates were observed.

RESULTS

The potential synergistic interactions between cellulase monocomponents and hemicellulases from family 10 and 11 endo-xylanases (GH10 EX and GH11 EX) and family 5 xyloglucanase (GH5 XG), during hydrolysis of various steam pretreated lignocellulosic substrates, were assessed. It was apparent that the hydrolytic activity of cellulase monocomponents was enhanced by the addition of accessory enzymes although the "boosting" effect was highly substrate specific. The GH10 EX and GH5 XG both exhibited broad substrate specificity and showed strong synergistic interaction with the cellulases when added individually. The GH10 EX was more effective on steam pretreated agriculture residues and hardwood substrates whereas GH5 XG addition was more effective on softwood substrates. The synergistic interaction between GH10 EX and GH5 XG when added together further enhanced the hydrolytic activity of the cellulase enzymes over a range of pretreated lignocellulosic substrates. GH10 EX addition could also stimulate further cellulose hydrolysis when added to the hydrolysis reactions when the rate of hydrolysis had levelled off.

CONCLUSIONS

Endo-xylanases and xyloglucanases interacted synergistically with cellulases to improve the hydrolysis of a range of pretreated lignocellulosic substrates. However, the extent of improved hydrolysis was highly substrate dependent. It appears that those accessory enzymes, such as GH10 EX and GH5 XG, with broader substrate specificities promoted the greatest improvements in the hydrolytic performance of the cellulase mixture on all of the pretreated biomass substrates.

Insights into exo- and endoglucanase activities of family 6 glycoside hydrolases from Podospora anserina

The ascomycete Podospora anserina is a coprophilous fungus that grows at late stages on droppings of herbivores. Its genome encodes a large diversity of carbohydrate-active enzymes. Among them, four genes encode glycoside hydrolases from family 6 (GH6), the members of which comprise putative endoglucanases and exoglucanases, some of them exerting important functions for biomass degradation in fungi. Therefore, this family was selected for functional analysis. Three of the enzymes, P. anserina Cel6A (PaCel6A), PaCel6B, and PaCel6C, were functionally expressed in the yeast Pichia pastoris. All three GH6 enzymes hydrolyzed crystalline and amorphous cellulose but were inactive on hydroxyethyl cellulose, mannan, galactomannan, xyloglucan, arabinoxylan, arabinan, xylan, and pectin. PaCel6A had a catalytic efficiency on cellotetraose comparable to that of Trichoderma reesei Cel6A (TrCel6A), but PaCel6B and PaCel6C were clearly less efficient. PaCel6A was the enzyme with the highest stability at 45°C, while PaCel6C was the least stable enzyme, losing more than 50% of its activity after incubation at temperatures above 30°C for 24 h. In contrast to TrCel6A, all three studied P. anserina GH6 cellulases were stable over a wide range of pHs and conserved high activity at pH values of up to 9. Each enzyme displayed a distinct substrate and product profile, highlighting different modes of action, with PaCel6A being the enzyme most similar to TrCel6A. PaCel6B was the only enzyme with higher specific activity on carboxymethylcellulose (CMC) than on Avicel and showed lower processivity than the others. Structural modeling predicts an open catalytic cleft, suggesting that PaCel6B is an endoglucanase.

Crystal structure and computational characterization of the lytic polysaccharide monooxygenase GH61D from the Basidiomycota fungus Phanerochaete chrysosporium

Carbohydrate structures are modified and degraded in the biosphere by a myriad of mostly hydrolytic enzymes. Recently, lytic polysaccharide mono-oxygenases (LPMOs) were discovered as a new class of enzymes for cleavage of recalcitrant polysaccharides that instead employ an oxidative mechanism. LPMOs employ copper as the catalytic metal and are dependent on oxygen and reducing agents for activity. LPMOs are found in many fungi and bacteria, but to date no basidiomycete LPMO has been structurally characterized. Here we present the three-dimensional crystal structure of the basidiomycete Phanerochaete chrysosporium GH61D LPMO, and, for the first time, measure the product distribution of LPMO action on a lignocellulosic substrate. The structure reveals a copper-bound active site common to LPMOs, a collection of aromatic and polar residues near the binding surface that may be responsible for regio-selectivity, and substantial differences in loop structures near the binding face compared with other LPMO structures. The activity assays indicate that this LPMO primarily produces aldonic acids. Last, molecular simulations reveal conformational changes, including the binding of several regions to the cellulose surface, leading to alignment of three tyrosine residues on the binding face of the enzyme with individual cellulose chains, similar to what has been observed for family 1 carbohydrate-binding modules. A calculated potential energy surface for surface translation indicates that P. chrysosporium GH61D exhibits energy wells whose spacing seems adapted to the spacing of cellobiose units along a cellulose chain.

The identification of and relief from Fe3+ inhibition for both cellulose and cellulase in cellulose saccharification catalyzed by cellulases from Penicillium decumbens

Lignocellulosic biomass is an underutilized, renewable resource that can be converted to biofuels. The key step in this conversion is cellulose saccharification catalyzed by cellulase. In this work, the effect of metal ions on cellulose hydrolysis by cellulases from Penicillium decumbens was reported for the first time. Fe(3+) and Cu(2+) were shown to be inhibitory. Further studies on Fe(3+) inhibition showed the inhibition takes place on both enzyme and substrate levels. Fe(3+) treatment damages cellulases' capability to degrade cellulose and inhibits all major cellulase activities. Fe(3+) treatment also reduces the digestibility of cellulose, due to its oxidation. Treatment of Fe(3+)-treated cellulose with DTT and supplementation of EDTA to saccharification systems partially relieved Fe(3+) inhibition. It was concluded that Fe(3+) inhibition in cellulose degradation is a complicated process in which multiple inhibition events occur, and that relief from Fe(3+) inhibition can be achieved by the supplementation of reducing or chelating agents.

X-ray structure and molecular dynamics simulations of endoglucanase 3 from Trichoderma harzianum: structural organization and substrate recognition by endoglucanases that lack cellulose binding module

Plant biomass holds a promise for the production of second-generation ethanol via enzymatic hydrolysis, but its utilization as a biofuel resource is currently limited to a large extent by the cost and low efficiency of the cellulolytic enzymes. Considerable efforts have been dedicated to elucidate the mechanisms of the enzymatic process. It is well known that most cellulases possess a catalytic core domain and a carbohydrate binding module (CBM), without which the enzymatic activity can be drastically reduced. However, Cel12A members of the glycosyl hydrolases family 12 (GHF12) do not bear a CBM and yet are able to hydrolyze amorphous cellulose quite efficiently. Here, we use X-ray crystallography and molecular dynamics simulations to unravel the molecular basis underlying the catalytic capability of endoglucanase 3 from Trichoderma harzianum (ThEG3), a member of the GHF12 enzymes that lacks a CBM. A comparative analysis with the Cellulomonas fimi CBM identifies important residues mediating interactions of EG3s with amorphous regions of the cellulose. For instance, three aromatic residues constitute a harboring wall of hydrophobic contacts with the substrate in both ThEG3 and CfCBM structures. Moreover, residues at the entrance of the active site cleft of ThEG3 are identified, which might hydrogen bond to the substrate. We advocate that the ThEG3 residues Asn152 and Glu201 interact with the substrate similarly to the corresponding CfCBM residues Asn81 and Arg75. Altogether, these results show that CBM motifs are incorporated within the ThEG3 catalytic domain and suggest that the enzymatic efficiency is associated with the length and position of the substrate chain, being higher when the substrate interact with the aromatic residues at the entrance of the cleft and the catalytic triad. Our results provide guidelines for rational protein engineering aiming to improve interactions of GHF12 enzymes with cellulosic substrates.

Effects of sulfate groups on the adsorption and activity of cellulases on cellulose substrates

Pretreatment of lignocellulosic biomass with sulfuric acid may leave sulfate groups on its surface that may hinder its biochemical conversion. This study investigates the effects of sulfate groups on cellulase adsorption onto cellulose substrates and the enzymatic hydrolysis of these substrates. Substrates with different sulfate group densities were prepared from H2SO4- and HCl-hydrolyzed and partially and fully desulfated cellulose nanocrystals. Adsorption onto and hydrolysis of the substrates was analyzed by quartz crystal microbalance with dissipation monitoring (QCM-D). The surface roughness of the substrates, measured by atomic force microscopy, increased with decreasing sulfate group density, but their surface accessibilities, measured by QCM-D H2O/D2O exchange experiments, were similar. The adsorption of cellulose binding domains onto sulfated substrates decreased with increasing sulfate group density, but the adsorption of cellulases increased. The rate of hydrolysis of sulfated substrates decreased with increasing sulfate group density. The results indicated an inhibitory effect of sulfate groups on the enzymatic hydrolysis of cellulose, possibly due to nonproductive binding of the cellulases onto the substrates through electrostatic interactions instead of their cellulose binding domains.

Crystallization and preliminary diffraction analysis of the catalytic domain of major extracellular endoglucanase from Xanthomonas campestris pv. campestris

Cellulases, such as endoglucanases, exoglucanases and β-glucosidases, are important enzymes used in the process of enzymatic hydrolysis of plant biomass. The bacteria Xanthomonas campestris pv. campestris expresses a large number of hydrolases and the major endoglucanase (XccEG), a member of glycoside hydrolase family 5 (GH5), is the most strongly secreted extracellularly. In this work, the native XccEG was purified from the extracellular extract and crystallization assays were performed on its catalytic domain. A complete data set was collected on an in-house X-ray source. The crystal diffracted to 2.7 Å resolution and belonged to space group C2, with unit-cell parameters a = 174.66, b = 141.53, c = 108.00 Å, β = 110.49°. The Matthews coefficient suggests a solvent content of 70.1% and the presence of four protein subunits in the asymmetric unit.

Structural, biochemical, and computational characterization of the glycoside hydrolase family 7 cellobiohydrolase of the tree-killing fungus Heterobasidion irregulare

Root rot fungi of the Heterobasidion annosum complex are the most damaging pathogens in temperate forests, and the recently sequenced Heterobasidion irregulare genome revealed over 280 carbohydrate-active enzymes. Here, H. irregulare was grown on biomass, and the most abundant protein in the culture filtrate was identified as the only family 7 glycoside hydrolase in the genome, which consists of a single catalytic domain, lacking a linker and carbohydrate-binding module. The enzyme, HirCel7A, was characterized biochemically to determine the optimal conditions for activity. HirCel7A was crystallized and the structure, refined at 1.7 Å resolution, confirms that HirCel7A is a cellobiohydrolase rather than an endoglucanase, with a cellulose-binding tunnel that is more closed than Phanerochaete chrysosporium Cel7D and more open than Hypocrea jecorina Cel7A, suggesting intermediate enzyme properties. Molecular simulations were conducted to ascertain differences in enzyme-ligand interactions, ligand solvation, and loop flexibility between the family 7 glycoside hydrolase cellobiohydrolases from H. irregulare, H. jecorina, and P. chrysosporium. The structural comparisons and simulations suggest significant differences in enzyme-ligand interactions at the tunnel entrance in the -7 to -4 binding sites and suggest that a tyrosine residue at the tunnel entrance of HirCel7A may serve as an additional ligand-binding site. Additionally, the loops over the active site in H. jecorina Cel7A are more closed than loops in the other two enzymes, which has implications for the degree of processivity, endo-initiation, and substrate dissociation. Overall, this study highlights molecular level features important to understanding this biologically and industrially important family of glycoside hydrolases.

Dissecting and reconstructing synergism: in situ visualization of cooperativity among cellulases

Cellulose is the most abundant biopolymer and a major reservoir of fixed carbon on earth. Comprehension of the elusive mechanism of its enzymatic degradation represents a fundamental problem at the interface of biology, biotechnology, and materials science. The interdependence of cellulose disintegration and hydrolysis and the synergistic interplay among cellulases is yet poorly understood. Here we report evidence from in situ atomic force microscopy (AFM) that delineates degradation of a polymorphic cellulose substrate as a dynamic cycle of alternating exposure and removal of crystalline fibers. Direct observation shows that chain-end-cleaving cellobiohydrolases (CBH I, CBH II) and an internally chain-cleaving endoglucanase (EG), the major components of cellulase systems, take on distinct roles: EG and CBH II make the cellulose surface accessible for CBH I by removing amorphous-unordered substrate areas, thus exposing otherwise embedded crystalline-ordered nanofibrils of the cellulose. Subsequently, these fibrils are degraded efficiently by CBH I, thereby uncovering new amorphous areas. Without prior action of EG and CBH II, CBH I was poorly active on the cellulosic substrate. This leads to the conclusion that synergism among cellulases is morphology-dependent and governed by the cooperativity between enzymes degrading amorphous regions and those targeting primarily crystalline regions. The surface-disrupting activity of cellulases therefore strongly depends on mesoscopic structural features of the substrate: size and packing of crystalline fibers are key determinants of the overall efficiency of cellulose degradation.

Cloning, molecular characterization, and mRNA expression of the thermostable family 3 β-glucosidase from the rare fungus Stachybotrys microspora

The filamentous fungus Stachybotrys microspora possess a rich β-glucosidase system composed of five β-glucosidases. Three of them were already purified to homogeneity and characterized. In order to isolate the β-glucosidase genes from S. microspora and study their regulation, a PCR strategy using consensus primers was used as a first step. This approach enabled the isolation of three different fragments of family 3 β-glucosidase gene. A representative genomic library was constructed and probed with one amplified fragment gene belonging to family 3 of β-glucosidase. After two rounds of hybridization, seven clones were obtained and the analysis of DNA plasmids leads to the isolation of one clone (CF3) with the largest insert of 7 kb. The regulatory region shows multiple TC-rich elements characteristic of constitutive promoter, explaining the expression of this gene under glucose condition, as shown by zymogram and RT-PCR analysis. The tertiary structure of the deduced amino acid sequence of Smbgl3 was predicted and has shown three conserved domains: an (α/β)8 triose phosphate isomerase (TIM) barrel, (α/β)5 sandwich, and fibronectin type III domain involved in protein thermostability. Zymogram analysis highlighted such thermostable character of this novel β-glucosidase.

Hallmarks of processivity in glycoside hydrolases from crystallographic and computational studies of the Serratia marcescens chitinases

Degradation of recalcitrant polysaccharides in nature is typically accomplished by mixtures of processive and nonprocessive glycoside hydrolases (GHs), which exhibit synergistic activity wherein nonprocessive enzymes provide new sites for productive attachment of processive enzymes. GH processivity is typically attributed to active site geometry, but previous work has demonstrated that processivity can be tuned by point mutations or removal of single loops. To gain additional insights into the differences between processive and nonprocessive enzymes that give rise to their synergistic activities, this study reports the crystal structure of the catalytic domain of the GH family 18 nonprocessive endochitinase, ChiC, from Serratia marcescens. This completes the structural characterization of the co-evolved chitinolytic enzymes from this bacterium and enables structural analysis of their complementary functions. The ChiC catalytic module reveals a shallow substrate-binding cleft that lacks aromatic residues vital for processivity, a calcium-binding site not previously seen in GH18 chitinases, and, importantly, a displaced catalytic acid (Glu-141), suggesting flexibility in the catalytic center. Molecular dynamics simulations of two processive chitinases (ChiA and ChiB), the ChiC catalytic module, and an endochitinase from Lactococcus lactis show that the nonprocessive enzymes have more flexible catalytic machineries and that their bound ligands are more solvated and flexible. These three features, which relate to the more dynamic on-off ligand binding processes associated with nonprocessive action, correlate to experimentally measured differences in processivity of the S. marcescens chitinases. These newly defined hallmarks thus appear to be key dynamic metrics in determining processivity in GH enzymes complementing structural insights.

A coarse-grained model for synergistic action of multiple enzymes on cellulose

BACKGROUND

Degradation of cellulose to glucose requires the cooperative action of three classes of enzymes, collectively known as cellulases. Endoglucanases randomly bind to cellulose surfaces and generate new chain ends by hydrolyzing β-1,4-D-glycosidic bonds. Exoglucanases bind to free chain ends and hydrolyze glycosidic bonds in a processive manner releasing cellobiose units. Then, β-glucosidases hydrolyze soluble cellobiose to glucose. Optimal synergistic action of these enzymes is essential for efficient digestion of cellulose. Experiments show that as hydrolysis proceeds and the cellulose substrate becomes more heterogeneous, the overall degradation slows down. As catalysis occurs on the surface of crystalline cellulose, several factors affect the overall hydrolysis. Therefore, spatial models of cellulose degradation must capture effects such as enzyme crowding and surface heterogeneity, which have been shown to lead to a reduction in hydrolysis rates.

RESULTS

We present a coarse-grained stochastic model for capturing the key events associated with the enzymatic degradation of cellulose at the mesoscopic level. This functional model accounts for the mobility and action of a single cellulase enzyme as well as the synergy of multiple endo- and exo-cellulases on a cellulose surface. The quantitative description of cellulose degradation is calculated on a spatial model by including free and bound states of both endo- and exo-cellulases with explicit reactive surface terms (e.g., hydrogen bond breaking, covalent bond cleavages) and corresponding reaction rates. The dynamical evolution of the system is simulated by including physical interactions between cellulases and cellulose.

CONCLUSIONS

Our coarse-grained model reproduces the qualitative behavior of endoglucanases and exoglucanases by accounting for the spatial heterogeneity of the cellulose surface as well as other spatial factors such as enzyme crowding. Importantly, it captures the endo-exo synergism of cellulase enzyme cocktails. This model constitutes a critical step towards testing hypotheses and understanding approaches for maximizing synergy and substrate properties with a goal of cost effective enzymatic hydrolysis.

Real-time observation of the swelling and hydrolysis of a single crystalline cellulose fiber catalyzed by cellulase 7B from Trichoderma reesei

The biodegradation of cellulose involves the enzymatic action of cellulases (endoglucanases), cellobiohydrolases (exoglucanases), and β-glucosidases that act synergistically. The rate and efficiency of enzymatic hydrolysis of crystalline cellulose in vitro decline markedly with time, limiting the large-scale, cost-effective production of cellulosic biofuels. Several factors have been suggested to contribute to this phenomenon, but there is considerable disagreement regarding the relative importance of each. These earlier investigations were hampered by the inability to observe the disruption of crystalline cellulose and its subsequent hydrolysis directly. Here, we show the application of high-resolution atomic force microscopy to observe the swelling of a single crystalline cellulose fiber and its-hydrolysis in real time directly as catalyzed by a single cellulase, the industrially important cellulase 7B from Trichoderma reesei. Volume changes, the root-mean-square roughness, and rates of hydrolysis of the surfaces of single fibers were determined directly from the images acquired over time. Hydrolysis dominated the early stage of the experiment, and swelling dominated the later stage. The high-resolution images revealed that the combined action of initial hydrolysis followed by swelling exposed individual microfibrils and bundles of microfibrils, resulting in the loosening of the fiber structure and the exposure of microfibrils at the fiber surface. Both the hydrolysis and swelling were catalyzed by the native cellulase; under the same conditions, its isolated carbohydrate-binding module did not cause changes to crystalline cellulose. We anticipate that the application of our AFM-based analysis on other cellulolytic enzymes, alone and in combination, will provide significant insight into the process of cellulose biodegradation and greatly facilitate its application for the efficient and economical production of cellulosic ethanol.

Binding preferences, surface attachment, diffusivity, and orientation of a family 1 carbohydrate-binding module on cellulose

Cellulase enzymes often contain carbohydrate-binding modules (CBMs) for binding to cellulose. The mechanisms by which CBMs recognize specific surfaces of cellulose and aid in deconstruction are essential to understand cellulase action. The Family 1 CBM from the Trichoderma reesei Family 7 cellobiohydrolase, Cel7A, is known to selectively bind to hydrophobic surfaces of native cellulose. It is most commonly suggested that three aromatic residues identify the planar binding face of this CBM, but several recent studies have challenged this hypothesis. Here, we use molecular simulation to study the CBM binding orientation and affinity on hydrophilic and hydrophobic cellulose surfaces. Roughly 43 μs of molecular dynamics simulations were conducted, which enables statistically significant observations. We quantify the fractions of the CBMs that detach from crystal surfaces or diffuse to other surfaces, the diffusivity along the hydrophobic surface, and the overall orientation of the CBM on both hydrophobic and hydrophilic faces. The simulations demonstrate that there is a thermodynamic driving force for the Cel7A CBM to bind preferentially to the hydrophobic surface of cellulose relative to hydrophilic surfaces. In addition, the simulations demonstrate that the CBM can diffuse from hydrophilic surfaces to the hydrophobic surface, whereas the reverse transition is not observed. Lastly, our simulations suggest that the flat faces of Family 1 CBMs are the preferred binding surfaces. These results enhance our understanding of how Family 1 CBMs interact with and recognize specific cellulose surfaces and provide insights into the initial events of cellulase adsorption and diffusion on cellulose.

Enzymatic synthesis of cellulose II-like substance via cellulolytic enzyme-mediated transglycosylation in an aqueous medium

The enzymatic synthesis of cellulose-like substance via a non-biosynthetic pathway has been achieved by transglycosylation in an aqueous system of the corresponding substrate, cellotriose for cellulolytic enzyme endo-acting endoglucanase I (EG I) from Hypocrea jecorina. A significant amount of water-insoluble product precipitated out from the reaction system. MALDI-TOF mass analysis showed that the resulting precipitate had a degree of polymerization (DP) of up to 16 from cellotriose. Solid-state (13)C NMR spectrum of the resulting water-insoluble product revealed that all carbon resonance lines were assigned to two kinds of anhydroglucose residues in the corresponding structure of cellulose II. X-ray diffraction (XRD) measurement as well as (13)C NMR analysis showed that the crystal structure corresponds to cellulose II with a high degree of crystallinity. We propose the multiple oligomers form highly crystalline cellulose II as a result of self-assembly via oligomer-oligomer interaction when they precipitate.

Cellulases dig deep: in situ observation of the mesoscopic structural dynamics of enzymatic cellulose degradation

Enzymatic hydrolysis of cellulose is key for the production of second generation biofuels, which represent a long-standing leading area in the field of sustainable energy. Despite the wealth of knowledge about cellulase structure and function, the elusive mechanism by which these enzymes disintegrate the complex structure of their insoluble substrate, which is the gist of cellulose saccharification, is still unclear. We herein present a time-resolved structural characterization of the action of cellulases on a nano-flat cellulose preparation, which enabled us to overcome previous limitations, using atomic force microscopy (AFM). As a first step in substrate disintegration, elongated fissures emerge which develop into coniform cracks as disintegration continues. Detailed data analysis allowed tracing the surface evolution back to the dynamics of crack morphology. This, in turn, reflects the interplay between surface degradation inside and outside of the crack. We observed how small cracks evolved and initially increased in size. At a certain point, the crack diameter stagnated and then started decreasing again. Stagnation corresponds with a decrease in the total amount of surface which is fissured and thus leads to the conclusion that the surface hydrolysis “around” the cracks is proceeding more rapidly than inside the cracks. The mesoscopic view presented here is in good agreement with various mechanistic proposals from the past and allows a novel insight into the structural dynamics occurring on the cellulosic substrate through cellulase action.

Effect of pH and temperature on the global compactness, structure, and activity of cellobiohydrolase Cel7A from Trichoderma harzianum

Due to its elevated cellulolytic activity, the filamentous fungus Trichoderma harzianum (T. harzianum) has considerable potential in biomass hydrolysis application. Cellulases from Trichoderma reesei have been widely used in studies of cellulose breakdown. However, cellulases from T. harzianum are less-studied enzymes that have not been characterized biophysically and biochemically as yet. Here, we examined the effects of pH and temperature on the secondary and tertiary structures, compactness, and enzymatic activity of cellobiohydrolase Cel7A from T. harzianum (Th Cel7A) using a number of biophysical and biochemical techniques. Our results show that pH and temperature perturbations affect Th Cel7A stability by two different mechanisms. Variations in pH modify protonation of the enzyme residues, directly affecting its activity, while leading to structural destabilization only at extreme pH limits. Temperature, on the other hand, has direct influence on mobility, fold, and compactness of the enzyme, causing unfolding of Th Cel7A just above the optimum temperature limit. Finally, we demonstrated that incubation with cellobiose, the product of the reaction and a competitive inhibitor, significantly increased the thermal stability of Th Cel7A. Our studies might provide insights into understanding, at a molecular level, the interplay between structure and activity of Th Cel7A at different pH and temperature conditions.

A structural study of Hypocrea jecorina Cel5A

Interest in generating lignocellulosic biofuels through enzymatic hydrolysis continues to rise as nonrenewable fossil fuels are depleted. The high cost of producing cellulases, hydrolytic enzymes that cleave cellulose into fermentable sugars, currently hinders economically viable biofuel production. Here, we report the crystal structure of a prevalent endoglucanase in the biofuels industry, Cel5A from the filamentous fungus Hypocrea jecorina. The structure reveals a general fold resembling that of the closest homolog with a high-resolution structure, Cel5A from Thermoascus aurantiacus. Consistent with previously described endoglucanase structures, the H. jecorina Cel5A active site contains a primarily hydrophobic substrate binding groove and a series of hydrogen bond networks surrounding two catalytic glutamates. The reported structure, however, demonstrates stark differences between side-chain identity, loop regions, and the number of disulfides. Such structural information may aid efforts to improve the stability of this protein for industrial use while maintaining enzymatic activity through revealing nonessential and immutable regions.

"Fluctuograms" reveal the intermittent intra-protein communication in subtilisin Carlsberg and correlate mechanical coupling with co-evolution

The mechanism of intra-protein communication and allosteric coupling is key to understanding the structure-property relationship of protein function. For subtilisin Carlsberg, the Ca²⁺-binding loop is distal to substrate-binding and active sites, yet the serine protease function depends on Ca²⁺ binding. The atomic molecular dynamics (MD) simulations of apo and Ca²⁺-bound subtilisin show similar structures and there is no direct evidence that subtilisin has alternative conformations. To model the intra-protein communication due to Ca²⁺ binding, we transform the sequential segments of an atomic MD trajectory into separate elastic network models to represent anharmonicity and nonlinearity effectively as the temporal and spatial variation of the mechanical coupling network. In analogy to the spectrogram of sound waves, this transformation is termed the “fluctuogram” of protein dynamics. We illustrate that the Ca²⁺-bound and apo states of subtilisin have different fluctuograms and that intra-protein communication proceeds intermittently both in space and in time. We found that residues with large mechanical coupling variation due to Ca²⁺ binding correlate with the reported mutation sites selected by directed evolution for improving the stability of subtilisin and its activity in a non-aqueous environment. Furthermore, we utilize the fluctuograms calculated from MD to capture the highly correlated residues in a multiple sequence alignment. We show that in addition to the magnitude, the variance of coupling strength is also an indicative property for the sequence correlation observed in a statistical coupling analysis. The results of this work illustrate that the mechanical coupling networks calculated from atomic details can be used to correlate with functionally important mutation sites and co-evolution.

A hallmark of protein molecules is their machine-like behaviors while carrying out biological functions. At the molecular level, molecular signals such as binding a metal ion at an action site can cause long-range effects and alter protein function. Such phenomena are often referred to as intra-protein communication or allosteric coupling. Elucidating the underlying mechanisms could lead to novel discovery of molecular modulators to regulate protein function in a more specific and effective manner. A long-standing puzzle is the roles of the anharmonicity and nonlinearity in protein dynamics. To incorporate these characters in modeling intra-protein communication, we devise a “fluctuogram” analysis to record the choreography of allosteric coupling in an atomic molecular dynamics simulation. We show that fluctuogram analysis can bridge the results of physics-based simulation and sequence alignment in bioinformatics by capturing the residues that exhibit high correlation in a multiple sequence alignment. We also show that the fluctuograms calculated from atomic details have the potential to be applied as a tool to select mutation sites for modulating protein function.