Molecular modeling suggests induced fit of Family I carbohydrate-binding modules with a broken-chain cellulose surface

Insights into the effect of protein glycosylation on carbohydrate substrate binding

The role of glycosylation in the binding of glycoproteins to carbohydrate substrates has not been well understood. The present study addresses this knowledge gap by elucidating the links between the glycosylation patterns of a model glycoprotein, a Family 1 carbohydrate-binding module (TrCBM1), and the thermodynamic and structural properties of its binding to different carbohydrate substrates using isothermal titration calorimetry and computational simulation. The variations in glycosylation patterns cause a gradual transition of the binding to soluble cellohexaose from an entropy-driven process to an enthalpy-driven one, a trend closely correlated with the glycan-induced shift of the predominant binding force from hydrophobic interactions to hydrogen bonding. However, when binding to a large surface of solid cellulose, glycans on TrCBM1 have a more dispersed distribution and thus have less adverse impact on the hydrophobic interaction forces, leading to overall improved binding. Unexpectedly, our simulation results also suggest an evolutionary role of O-mannosylation in transforming the substrate binding features of TrCBM1 from those of type A CBMs to those of type B CBMs. Taken together, these findings provide new fundamental insights into the molecular basis of the role of glycosylation in protein-carbohydrate interactions and are expected to better facilitate further studies in this area.

Direct Evidence for Aligned Binding of Cellulase Enzymes to Cellulose Surfaces

The conversion of biomass into green fuels and chemicals is of great societal interest. Engineers have been designing new cellulase enzymes for the breakdown of otherwise insoluble cellulose materials. A barrier to the rational design of new enzymes has been our lack of a molecular picture of how cellulase binding occurs. A critical factor is the attachment via the enzyme's carbohydrate binding module (CBM). To elucidate the structural and mechanistic details of cellulase adsorption, we have combined experimental data from sum frequency generation spectroscopy with molecular dynamics simulations to probe the equilibrium structure and surface alignment of a 14-residue peptide mimicking the CBM. The data show that binding is driven by hydrogen bonding and that tyrosine side chains within the CBM align the cellulase with the registry of the cellulose surface. Such an alignment is favorable for the translocation and effective cellulose breakdown and is therefore likely an important parameter for the design of novel enzymes.

Biotechnological Eminence of Chitinases: A Focus on Thermophilic Enzyme Sources, Production Strategies and Prominent Applications

BACKGROUND

Chitin, the second most abundant polysaccharide in nature, is a constantly valuable and renewable raw material after cellulose. Due to advancement in technology, industrial interest has grown to take advantage of the chitin.

OBJECTIVE

Now, biomass is being treated with diverse microbial enzymes or cells for the production of desired products under best industrial conditions. Glycosidic bonds in chitin structure are degraded by chitinase enzymes, which are characterized into number of glycoside hydrolase (GHs) families.

METHODS

Thermophilic microorganisms are remarkable sources of industrially important thermostable enzymes, having ability to survive harsh industrial processing conditions. Thermostable chitinases have an edge over mesophilic chitinases as they can hydrolyse the substrate at relatively high temperatures and exhibit decreased viscosity, significantly reduced contamination risk, thermal and chemical stability and increased solubility. Various methods are employed to purify the enzyme and increase its yield by optimizing various parameters such as temperature, pH, agitation, and by investigating the effect of different chemicals and metal ions etc. Results: Thermostable chitinase enzymes show high specific activity at elevated temperature which distinguish them from mesophiles. Genetic engineering can be used for further improvement of natural chitinases, and unlimited potential for the production of thermophilic chitinases has been highlighted due to advancement in synthetic biological techniques. Thermostable chitinases are then used in different fields such as bioremediation, medicine, agriculture and pharmaceuticals.

CONCLUSION

This review will provide information about chitinases, biotechnological potential of thermostable enzyme and the methods by which they are being produced and optimized for several industrial applications. Some of the applications of thermostable chitinases have also been briefly described.

Interaction of enzymes with lignocellulosic materials: causes, mechanism and influencing factors

For the production of biofuel (bioethanol), enzymatic adsorption onto a lignocellulosic biomass surface is a prior condition for the enzymatic hydrolysis process to occur. Lignocellulosic substances are mainly composed of cellulose, hemicellulose and lignin. The polysaccharide matrix (cellulose and hemicellulose) is capable of producing bioethanol. Therefore, lignin is removed or its concentration is reduced from the adsorption substrates by pretreatments. Selected enzymes are used for the production of reducing sugars from cellulosic materials, which in turn are converted to bioethanol. Adsorption of enzymes onto the substrate surface is a complicated process. A large number of research have been performed on the adsorption process, but little has been done to understand the mechanism of adsorption process. This article reviews the mechanisms of adsorption of enzymes onto the biomass surfaces. A conceptual adsorption mechanism is presented which will fill the gaps in literature and help researchers and industry to use adsorption more efficiently. The process of enzymatic adsorption starts with the reciprocal interplay of enzymes and substrates and ends with the establishment of molecular and cellular binding. The kinetics of an enzymatic reaction is almost the same as that of a characteristic chemical catalytic reaction. The influencing factors discussed in detail are: surface characteristics of the participating materials, the environmental factors, such as the associated flow conditions, temperature, concentration, etc. Pretreatment of lignocellulosic materials and optimum range of shear force and temperature for getting better results of adsorption are recommended.

Molecular dynamics investigations of oligosaccharides recognized by family 16 and 22 carbohydrate binding modules

As a non-catalytic domain, carbohydrate binding modules (CBMs) are often considered to play some key roles in the degradation and recognition of polysaccharides catalyzed by cellulases. In this work, we investigated the recognition dynamics of cello- or xylo-saccharides by two typical CBMs (CBM16-1 and CBM22-2), which are grouped into Type B CBMs. By combining extensive molecular dynamics, principle component analysis, and binding free energy calculations, we constructed several complex models of the two CBMs in both complex cello- and xylo-oligosaccharides. The corresponding substrate recognition affinity and critical residues having significant contributions were systematically investigated. The residues containing aromatic side chain groups were shown to contribute significantly to substrate binding. The calculated binding free energies were in fairly good agreement with the experimental measurements with the absolute mean error of 0.69 kcal mol^-1. The overall electrostatic interactions were shown to have negative effects on substrate recognition. Further metadynamics simulations revealed the substrate dissociation process.

Interaction of carbohydrate-binding modules with poly(ethylene terephthalate)

Poly(ethylene terephthalate) (PET) is one of the most widely applied synthetic polymers, but its hydrophobicity is challenging for many industrial applications. Biotechnological modification of PET surface can be achieved by PET hydrolyzing cutinases. In order to increase the adsorption towards their unnatural substrate, the enzymes are fused to carbohydrate-binding modules (CBMs) leading to enhanced activity. In this study, we identified novel PET binding CBMs and characterized the CBM-PET interplay. We developed a semi-quantitative method to detect CBMs bound to PET films. Screening of eight CBMs from diverse families for PET binding revealed one CBM that possesses a high affinity towards PET. Molecular dynamics (MD) simulations of the CBM-PET interface revealed tryptophan residues forming an aromatic triad on the peptide surface. Their interaction with phenyl rings of PET is stabilized by additional hydrogen bonds formed between amino acids close to the aromatic triad. Furthermore, the ratio of hydrophobic to polar contacts at the interface was identified as an important feature determining the strength of PET binding of CBMs. The interaction of CBM tryptophan residues with PET was confirmed experimentally by tryptophan quenching measurements after addition of PET nanoparticles to CBM. Our findings are useful for engineering PET hydrolyzing enzymes and may also find applications in functionalization of PET.

Molecular dynamics investigations of cello-oligosaccharide recognition by Cel9G-CBM3c from Clostridium cellulovorans

The processive mechanism of cellulases against cellulose represents one of the key mechanisms in the conversion of biomass. A reliable model of substrate binding in a multidomain cellulase is a prerequisite for fully understanding this mechanism. In this study, the specificity of the recognition of the polysaccharide by the multidomain endoglucanase Cel9G from Clostridium cellulovorans was investigated by molecular dynamics simulations. Aromatic ring-containing residues were found to be critical for stabilizing the substrate. The calculated subtotal contributions of polar residues close to the active site, e.g., D58, E244, R315 and D420, also have some critical functions in substrate binding. Unlike other members of the carbohydrate-binding module family, CBM3c alone is shown not to bind cellulose very well, which is also consistent with experimental conclusions.

Comparison of catalytically non-productive adsorption of fungal proteins to lignins and pseudo-lignin using isobaric mass tagging

Catalytically non-productive adsorption of fungal enzymes to pseudo-lignin (PL) was compared to adsorption to lignin preparations derived from different sources (SL, spruce; BL, birch; OL, beech) using different methods [steam pretreatment/enzymatic saccharification (SL, BL) and organosolv processing (OL)]. The protein adsorption to the SL was more extensive than the adsorption to the hardwood lignins, which was relatively similar to the adsorption to the PL. The adsorption patterns of 13 individual proteins were studied using isobaric mass tagging with TMTsixplex reagent and LC-MS/MS analysis. The results suggest that, on an average, adsorption of proteins equipped with carbohydrate-binding modules, such as the cellulases CBHI, EGII, and EGIV, was less dependent on the quality of the lignin/PL than adsorption of other proteins, such as β-Xyl, Xyn-1, and Xyn-2, which are involved in xylan degradation.

Approximate quantum chemical methods for modelling carbohydrate conformation and aromatic interactions: β-cyclodextrin and its adsorption on a single-layer graphene sheet

Adsorption of carbohydrates on graphene has the potential to improve graphene dispersibility in water. Here we assess the ability of DFTB-based and NDDO-based quantum chemical methods to model β-cyclodextrin conformations and interactions with graphene.

Carbohydrate-binding modules of fungal cellulases: occurrence in nature, function, and relevance in industrial biomass conversion

In this review, the present knowledge on the occurrence of cellulases, with a special emphasis on the presence of carbohydrate-binding modules (CBMs) in various fungal strains, has been summarized. The importance of efficient fungal cellulases is growing due to their potential uses in biorefinery processes where lignocellulosic biomasses are converted to platform sugars and further to biofuels and chemicals. Most secreted cellulases studied in detail have a bimodular structure containing an active core domain attached to a CBM. CBMs are traditionally been considered as essential parts in cellulases, especially in cellobiohydrolases. However, presently available genome data indicate that many cellulases lack the binding domains in cellulose-degrading organisms. Recent data also demonstrate that CBMs are not necessary for the action of cellulases and they solely increase the concentration of enzymes on the substrate surfaces. On the other hand, in practical industrial processes where high substrate concentrations with low amounts of water are employed, the enzymes have been shown to act equally efficiently with and without CBM. Furthermore, available kinetic data show that enzymes without CBMs can desorb more readily from the often lignaceous substrates, that is, they are not stuck on the substrates and are thus available for new actions. In this review, the available data on the natural habitats of different wood-degrading organisms (with emphasis on the amount of water present during wood degradation) and occurrence of cellulose-binding domains in their genome have been assessed in order to identify evolutionary advantages for the development of CBM-less cellulases in nature.

Measurements of single molecular affinity interactions between carbohydrate-binding modules and crystalline cellulose fibrils

Combining atomic force microscopy (AFM) recognition imaging and single molecule dynamic force spectroscopy (SMDFS), we studied the single molecule affinity interactions between the carbohydrate-binding module (CBM) and plant cell wall cellulose using the CBM3a (from Clostridium thermocellum) and CBM2a (from Cellvibrio japonicus) functionalized AFM tips. The binding efficiencies of the CBMs to the cellulose were determined by the binding areas on the crystalline cellulose fibrils surface using the recognition imaging. Several dynamic and kinetic parameters, such as the reconstructed free energy change, energy barrier and bond lifetime constant, were also obtained based on the measured single molecule unbinding forces, which are used to illuminate the affinity of the CBMs binding to the natural and single cellulose surface from a totally different aspect. It was found that CBM3a has a little higher binding efficiency and affinity than CBM2a to both natural and extracted cellulose surfaces and both the CBMs have higher affinities to the natural cell wall cellulose compared to the extracted single cellulose. The in-depth understanding of the binding mechanisms of the CBM-cellulose interactions of this study may pave the way for more efficient plant cell wall degradation and eventually facilitate biofuel production.

Two family 11 xylanases from Achaetomium sp. Xz-8 with high catalytic efficiency and application potentials in the brewing industry

This study identified two family-11 xylanase genes (xynC81 and xynC83) in Achaetomium sp. Xz-8, a thermophilic strain from a desert area with substantial xylanase activity, and successfully expressed them in Pichia pastoris . Their deduced amino acid sequences showed the highest identity of ≤90% to known fungal xylanases and of ≤62% with each other. The purified recombinant xylanases showed optimal activities at pH 5.5 and 60-65 °C and exhibited stability over pH 5.0-10.0 and temperatures at 55 °C and below. XynC81 had high catalytic efficiency (6082 mL/s/mg), and XynC83 was favorable for xylooligosaccharide production. Under simulated mashing conditions, combination of XynC83 and a commercial β-glucanase improved the filtration rate by 34.76%, which is much better than that of Novozymes Ultraflo (20.71%). XynC81 and XynC83 had a synergistic effect on viscosity reduction (7.08%), which is comparable with that of Ultraflo (8.47%). Thus, XynC81 and XynC83 represent good candidates for application in the brewing industry.

Introduction of a unique tryptophan residue into various positions of Bacillus licheniformis DnaK

Site-directed mutagenesis together with biochemical and biophysical techniques were used to probe effects of single-tryptophan-incorporated mutations on a bacterial molecular chaperone, Bacillus licheniformis DnaK (BlDnaK). Specifically, five phenylalanine residues (Phe(120), Phe(174), Phe(186), Phe(378) and Phe(396)) of BlDnaK were individually replaced by single tryptophans, thus creating site-specific probes for the fluorescence analysis of the protein. The steady-state ATPase activity for BlDnaK, F120W, F174W, F186W, F378W, and F396W was determined to be 76.01, 52.82, 25.32, 53.31, 58.84, and 47.53 nmol Pi/min/mg, respectively. Complementation test revealed that the single mutation at codons 120, 186, and 378 of the dnaK gene still allowed an Escherichia coli dnaK756-Ts strain to grow at a stringent temperature of 44°C. Simultaneous addition of co-chaperones and NR-peptide did not synergistically stimulate the ATPase activity of F174W and F396W, and these two proteins were unable to assist the refolding of GdnHCl-denatured luciferase. The heat-induced denaturation of all variants could be fitted adequately to a three-state model, in agreement with the observation for the wild-type protein. By CD spectral analysis, GdnHCl-induced unfolding transition for BlDnaK was 1.51 M corresponding to ΔG(N-U) of 1.69 kcal/mol; however, the transitions for mutant proteins were 1.07-1.55 M equivalent to ΔG(N-U) of 0.94-2.93 kcal/mol. The emission maximum of single-tryptophan-incorporated variants was in the range of 333.2-335.8 nm. Acrylamide quenching analysis showed that the mutant proteins had a dynamic quenching constant of 3.0-4.2 M(-1). Taken together, these results suggest that the molecular properties of BlDnaK have been significantly changed upon the individual replacement of selected phenylalanine residues by tryptophan.

Experimental evidence for the involvement of amino acid residue Glu398 in the autocatalytic processing of Bacillus licheniformis γ-glutamyltranspeptidase

The role of glutamate 398 in the autocatalytic processing of Bacillus licheniformis γ-glutamyltranspeptidase (BlGGT) was explored by site-directed mutagenesis. This glutamate was substituted by either alanine, aspartate, arginine or glutamine and the expressed mutant enzymes were purified to apparent homogeneity with metal-affinity chromatography. SDS–PAGE analysis showed that E398A, E398D and E398K were unable to process themselves into a large and a small subunit. However, E398Q was not only able to process itself, but also had a catalytic activity comparable to that of BlGGT. As compared with the wild-type enzyme, no significant change in circular dichroism spectra was observed for the mutant proteins. Thermal unfolding of BlGGT, E398A, E398D, E398K and E398Q followed the two-state unfolding process with a transition point (T_m) of 47.7–69.4 °C. Tryptophan fluorescence spectra of the mutant enzymes were different from the wild-type protein in terms of fluorescence intensity. Native BlGGT started to unfold beyond ∼1.92 M guanidine hydrochloride (GdnHCl) and reached an unfolded intermediate, [GdnHCl]_{0.5, N–U}, at 3.07 M equivalent to free energy change (Δ⁢GN−UH2⁢O) of 14.53 kcal/mol for the N → U process, whereas the denaturation midpoints for the mutant enzymes were 1.31–2.99 M equivalent to Δ⁢GN−UH2⁢O of 3.29–12.05 kcal/mol. Taken together, these results strongly suggest that the explored glutamate residue is indeed important for the autocatalytic processing of BlGGT.

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.

Exploring the conformational space of amorphous cellulose using NMR chemical shifts

(13)C-labeled amorphous cellulose and (13)C NMR chemical shifts by 2D (13)C-(13)C correlation spectroscopy were obtained in the regenerated solid-state from ionic liquids. On the basis of the assigned chemical shifts, combined with information from molecular dynamics and quantum chemistry computer simulations a twisted structure for amorphous cellulose is proposed exposing more hydrophilic surface than that of extended crystalline cellulose.

Two structurally discrete GH7-cellobiohydrolases compete for the same cellulosic substrate fiber

BACKGROUND

Cellulose consisting of arrays of linear beta-1,4 linked glucans, is the most abundant carbon-containing polymer present in biomass. Recalcitrance of crystalline cellulose towards enzymatic degradation is widely reported and is the result of intra- and inter-molecular hydrogen bonds within and among the linear glucans. Cellobiohydrolases are enzymes that attack crystalline cellulose. Here we report on two forms of glycosyl hydrolase family 7 cellobiohydrolases common to all Aspergillii that attack Avicel, cotton cellulose and other forms of crystalline cellulose.

RESULTS

Cellobiohydrolases Cbh1 and CelD have similar catalytic domains but only Cbh1 contains a carbohydrate-binding domain (CBD) that binds to cellulose. Structural superpositioning of Cbh1 and CelD on the Talaromyces emersonii Cel7A 3-dimensional structure, identifies the typical tunnel-like catalytic active site while Cbh1 shows an additional loop that partially obstructs the substrate-fitting channel. CelD does not have a CBD and shows a four amino acid residue deletion on the tunnel-obstructing loop providing a continuous opening in the absence of a CBD. Cbh1 and CelD are catalytically functional and while specific activity against Avicel is 7.7 and 0.5 U.mg prot-1, respectively specific activity on pNPC is virtually identical. Cbh1 is slightly more stable to thermal inactivation compared to CelD and is much less sensitive to glucose inhibition suggesting that an open tunnel configuration, or absence of a CBD, alters the way the catalytic domain interacts with the substrate. Cbh1 and CelD enzyme mixtures on crystalline cellulosic substrates show a strong combinatorial effort response for mixtures where Cbh1 is present in 2:1 or 4:1 molar excess. When CelD was overrepresented the combinatorial effort could only be partially overcome. CelD appears to bind and hydrolyze only loose cellulosic chains while Cbh1 is capable of opening new cellulosic substrate molecules away from the cellulosic fiber.

CONCLUSION

Cellobiohydrolases both with and without a CBD occur in most fungal genomes where both enzymes are secreted, and likely participate in cellulose degradation. The fact that only Cbh1 binds to the substrate and in combination with CelD exhibits strong synergy only when Cbh1 is present in excess, suggests that Cbh1 unties enough chains from cellulose fibers, thus enabling processive access of CelD.

Overview of computer modeling of cellulose

Although it has a deceptively simple primary structure, the collective organization of bulk cellulose, particularly as it exists in cellulose fibers in the cell walls of living plants and other organisms, is quite diverse and complex. While some experimental techniques, such as vibrational spectroscopy and diffraction from partially crystalline samples, are able to provide insights into the organization of bulk cellulose, its intrinsic complexity has left many questions still unanswered. For this reason, additional probes of cellulose structure would be highly desirable. With the continuing advances in computer power through massive parallelization, and the steady progress in computer codes and force fields for modeling carbohydrate systems, molecular mechanics simulations have become an attractive means of studying cellulosic systems at the atomic and molecular level. The coming decade will almost certainly see remarkable advances in the understanding of cellulose using such simulations.

Contribution of conserved Glu255 and Cys289 residues to catalytic activity of recombinant aldehyde dehydrogenase from Bacillus licheniformis

Based on the sequence homology, we have modeled the three-dimensional structure of Bacillus licheniformis aldehyde dehydrogenase (BlALDH) and identified two different residues, Glu255 and Cys289, that might be responsible for the catalytic function of the enzyme. The role of these residues was further investigated by site-directed mutagenesis and biophysical analysis. The expressed parental and mutant proteins were purified by nickel-chelate chromatography, and their molecular masses were determined to be approximately 53 kDa by SDS-PAGE. As compared with the parental BlALDH, a dramatic decrease or even complete loss of the dehydrogenase activity was observed for the mutant enzymes. Structural analysis showed that the intrinsic fluorescence and circular dichroism spectra of the mutant proteins were similar to the parental enzyme, but most of the variants exhibited a different sensitivity towards thermal- and guanidine hydrochloride-induced denaturation. These observations indicate that residues Glu255 and Cys289 play an important role in the dehydrogenase activity of BlALDH, and the rigidity of the enzyme has been changed as a consequence of the mutations.

A cellular automaton model of crystalline cellulose hydrolysis by cellulases

BACKGROUND

Cellulose from plant biomass is an abundant, renewable material which could be a major feedstock for low emissions transport fuels such as cellulosic ethanol. Cellulase enzymes that break down cellulose into fermentable sugars are composed of different types - cellobiohydrolases I and II, endoglucanase and β-glucosidase - with separate functions. They form a complex interacting network between themselves, soluble hydrolysis product molecules, solution and solid phase substrates and inhibitors. There have been many models proposed for enzymatic saccharification however none have yet employed a cellular automaton approach, which allows important phenomena, such as enzyme crowding on the surface of solid substrates, denaturation and substrate inhibition, to be considered in the model.

RESULTS

The Cellulase 4D model was developed de novo taking into account the size and composition of the substrate and surface-acting enzymes were ascribed behaviors based on their movements, catalytic activities and rates, affinity for, and potential for crowding of, the cellulose surface, substrates and inhibitors, and denaturation rates. A basic case modeled on literature-derived parameters obtained from Trichoderma reesei cellulases resulted in cellulose hydrolysis curves that closely matched curves obtained from published experimental data. Scenarios were tested in the model, which included variation of enzyme loadings, adsorption strengths of surface acting enzymes and reaction periods, and the effect on saccharide production over time was assessed. The model simulations indicated an optimal enzyme loading of between 0.5 and 2 of the base case concentrations where a balance was obtained between enzyme crowding on the cellulose crystal, and that the affinities of enzymes for the cellulose surface had a large effect on cellulose hydrolysis. In addition, improvements to the cellobiohydrolase I activity period substantially improved overall glucose production.

CONCLUSIONS

Cellulase 4D simulates the enzymatic hydrolysis of cellulose to glucose by surface and solution phase-acting enzymes and accounts for complex phenomena that have previously not been included in cellulose hydrolysis models. The model is intended as a tool for industry, researchers and educators alike to explore options for enzyme engineering and process development and to test hypotheses regarding cellulase mechanisms.

Tissue-specific biomass recalcitrance in corn stover pretreated with liquid hot-water: SEM imaging (part 2)

In the first part of our work, we combined compositional analysis, pretreatment and enzyme hydrolysis for fractionated pith, rind, and leaf tissues from a hybrid stay-green corn, in order to identify the role of structural characteristics on enzyme hydrolysis of cell walls. Hydrolysis experiments coupled with chemical analysis of the different fractions of corn stover showed significant differences in cell wall structure before and after liquid hot water pretreatment. The extent of enzyme hydrolysis followed the sequence rind < leaves < pith with 90% conversion of cellulose to glucose in 24 h in the best cases. Since similar lignin contents remained after liquid hot water pretreatment of leaves, rind, and pith, our results indicated that the amount of lignin alone is not sufficient to explain the different enzymatic hydrolysis characteristics of the fractions. While the role of structural characteristics on enzyme hydrolysis of cell walls is measured as described in part I, the SEM images presented in this part II of our work show that sugar yields from enzymatic hydrolysis of corn fractions correlate with changes in plant cell wall structure both before and after liquid hot water pretreatment.

High-temperature enzymatic breakdown of cellulose

Cellulose is an abundant and renewable biopolymer that can be used for biofuel generation; however, structural entrapment with other cell wall components hinders enzyme-substrate interactions, a key bottleneck for ethanol production. Biomass is routinely subjected to treatments that facilitate cellulase-cellulose contacts. Cellulases and glucosidases act by hydrolyzing glycosidic bonds of linear glucose β-1,4-linked polymers, producing glucose. Here we describe eight high-temperature-operating cellulases (TCel enzymes) identified from a survey of thermobacterial and archaeal genomes. Three TCel enzymes preferentially hydrolyzed soluble cellulose, while two preferred insoluble cellulose such as cotton linters and filter paper. TCel enzymes had temperature optima ranging from 85°C to 102°C. TCel enzymes were stable, retaining 80% of initial activity after 120 h at 85°C. Two modes of cellulose breakdown, i.e., with endo- and exo-acting glucanases, were detected, and with two-enzyme combinations at 85°C, synergistic cellulase activity was observed for some enzyme combinations.

Probing carbohydrate product expulsion from a processive cellulase with multiple absolute binding free energy methods

Understanding the enzymatic mechanism that cellulases employ to degrade cellulose is critical to efforts to efficiently utilize plant biomass as a sustainable energy resource. A key component of cellulase action on cellulose is product inhibition from monosaccharide and disaccharides in the product site of cellulase tunnel. The absolute binding free energy of cellobiose and glucose to the product site of the catalytic tunnel of the Family 7 cellobiohydrolase (Cel7A) of Trichoderma reesei (Hypocrea jecorina) was calculated using two different approaches: steered molecular dynamics (SMD) simulations and alchemical free energy perturbation molecular dynamics (FEP/MD) simulations. For the SMD approach, three methods based on Jarzynski's equality were used to construct the potential of mean force from multiple pulling trajectories. The calculated binding free energies, -14.4 kcal/mol using SMD and -11.2 kcal/mol using FEP/MD, are in good qualitative agreement. Analysis of the SMD pulling trajectories suggests that several protein residues (Arg-251, Asp-259, Asp-262, Trp-376, and Tyr-381) play key roles in cellobiose and glucose binding to the catalytic tunnel. Five mutations (R251A, D259A, D262A, W376A, and Y381A) were made computationally to measure the changes in free energy during the product expulsion process. The absolute binding free energies of cellobiose to the catalytic tunnel of these five mutants are -13.1, -6.0, -11.5, -7.5, and -8.8 kcal/mol, respectively. The results demonstrated that all of the mutants tested can lower the binding free energy of cellobiose, which provides potential applications in engineering the enzyme to accelerate the product expulsion process and improve the efficiency of biomass conversion.

The O-glycosylated linker from the Trichoderma reesei Family 7 cellulase is a flexible, disordered protein

Fungi and bacteria secrete glycoprotein cocktails to deconstruct cellulose. Cellulose-degrading enzymes (cellulases) are often modular, with catalytic domains for cellulose hydrolysis and carbohydrate-binding modules connected by linkers rich in serine and threonine with O-glycosylation. Few studies have probed the role that the linker and O-glycans play in catalysis. Since different expression and growth conditions produce different glycosylation patterns that affect enzyme activity, the structure-function relationships that glycosylation imparts to linkers are relevant for understanding cellulase mechanisms. Here, the linker of the Trichoderma reesei Family 7 cellobiohydrolase (Cel7A) is examined by simulation. Our results suggest that the Cel7A linker is an intrinsically disordered protein with and without glycosylation. Contrary to the predominant view, the O-glycosylation does not change the stiffness of the linker, as measured by the relative fluctuations in the end-to-end distance; rather, it provides a 16 Å extension, thus expanding the operating range of Cel7A. We explain observations from previous biochemical experiments in the light of results obtained here, and compare the Cel7A linker with linkers from other cellulases with sequence-based tools to predict disorder. This preliminary screen indicates that linkers from Family 7 enzymes from other genera and other cellulases within T. reesei may not be as disordered, warranting further study.

Applications of computational science for understanding enzymatic deconstruction of cellulose

Understanding the molecular-level mechanisms that enzymes employ to deconstruct plant cell walls is a fundamental scientific challenge with significant ramifications for renewable fuel production from biomass. In nature, bacteria and fungi use enzyme cocktails that include processive and non-processive cellulases and hemicellulases to convert cellulose and hemicellulose to soluble sugars. Catalyzed by an accelerated biofuels R&D portfolio, there is now a wealth of new structural and experimental insights related to cellulases and the structure of plant cell walls. From this background, computational approaches commonly used in other fields are now poised to offer insights complementary to experiments designed to probe mechanisms of plant cell wall deconstruction. Here we outline the current status of computational approaches for a collection of critical problems in cellulose deconstruction. We discuss path sampling methods to measure rates of elementary steps of enzyme action, coarse-grained modeling for understanding macromolecular, cellulosomal complexes, methods to screen for enzyme improvements, and studies of cellulose at the molecular level. Overall, simulation is a complementary tool to understand carbohydrate-active enzymes and plant cell walls, which will enable industrial processes for the production of advanced, renewable fuels.

Bioenergy research: a new paradigm in multidisciplinary research

The field of biology is becoming increasingly interdisciplinary and cross-cutting. This changing research atmosphere is creating the way for a new kind of enquiry that while building upon the traditional research establishment is providing a new multidisciplinary framework to more effectively address scientific grand challenges. Using the US Department of Energy sponsored BioEnergy Science Center as an example, we highlight how impactful breakthroughs in biofuel science can be achieved within a large cross-disciplinary team environment. Such transformational insights are key to furthering our understanding and in generating models, theories and processes that can be used to overcome recalcitrance of biomass for sustainable biofuel production. Multidisciplinary approaches have an increasingly greater role to play in meeting rising demands for food, fibre, energy, clean environment and good health. Discoveries achieved by diverse minds and cross-applications of tools and analytical approaches have tremendous potential to fill existing knowledge gaps, clear roadblocks and facilitate translation of basic sciences discoveries as solutions towards addressing some of the most pressing global issues.

Systematic docking study of the carbohydrate binding module protein of Cel7A with the cellulose Ialpha crystal model

A computer docking study has been carried out on the crystal surfaces of cellulose Ialpha crystal models for the carbohydrate binding module (CBM) protein of the cellobiohydrolase Cel7A produced by Trichoderma reesei. Binding free energy maps between the CBM and the crystal surface were obtained by calculating the noncovalent interactions and the solvation free energy at grid points covering the area of the unit cell dimensions at the crystal surface. The potential maps obtained from grid searches of the hydrophobic (110) crystal surface exhibited two distinct potential wells. These reflected the 2-fold helical symmetry of the cellulose chain and had lower binding energies at the minimum positions than those for the hydrophilic (100) and (010) crystal surfaces. The CBM-cellulose crystal complex models derived from the minimum positions were then subjected to molecular dynamics (MD) simulation under an explicit solvent system. The (110) complex models exhibited larger affinities at the interface than the (100) and (010) ones. The CBM was more stably bound to the (110) surface when it was placed in an antiparallel orientation with respect to the cellulose fiber axis. In the solvated dynamics state, the curved (110) surface resulting from the fiber twist somewhat assisted a complementary fit with the CBM at the interface. In addition to the conventional Generalized Born (GB) method, the three-dimensional reference interaction site model (3D-RISM) theory was adopted to assess a solvent effect for the solvated MD trajectories. Large exothermic values for the noncovalent interactions appeared correlated to and were mostly compensated by endothermic values for the solvation free energy. These gave total binding free energies of -13 to -28 kcal/mol. Results also suggested that the hydrogen bonding scheme was not essential for substrate specificity.

The energy landscape for the interaction of the family 1 carbohydrate-binding module and the cellulose surface is altered by hydrolyzed glycosidic bonds

A multiscale simulation model is used to construct potential and free energy surfaces for the carbohydrate-binding module [CBM] from an industrially important cellulase, Trichoderma reesei cellobiohydrolase I, on the hydrophobic face of a coarse-grained cellulose Ibeta polymorph. We predict from computation that the CBM alone exhibits regions of stability on the hydrophobic face of cellulose every 5 and 10 A, corresponding to a glucose unit and a cellobiose unit, respectively. In addition, we predict a new role for the CBM: specifically, that in the presence of hydrolyzed cellulose chain ends, the CBM exerts a thermodynamic driving force to translate away from the free cellulose chain ends. This suggests that the CBM is not only required for binding to cellulose, as has been known for two decades, but also that it has evolved to both assist the enzyme in recognizing a cellulose chain end and exert a driving force on the enzyme during processive hydrolysis of cellulose.

Computational simulations of the Trichoderma reesei cellobiohydrolase I acting on microcrystalline cellulose Ibeta: the enzyme-substrate complex

Cellobiohydrolases are the dominant components of the commercially relevant Trichoderma reesei cellulase system. Although natural cellulases can totally hydrolyze crystalline cellulose to soluble sugars, the current enzyme loadings and long digestion times required render these enzymes less than cost effective for biomass conversion processes. It is clear that cellobiohydrolases must be improved via protein engineering to reduce processing costs. To better understand cellobiohydrolase function, new simulations have been conducted using charmm of cellobiohydrolase I (CBH I) from T.reesei interacting with a model segment (cellodextrin) of a cellulose microfibril in which one chain from the substrate has been placed into the active site tunnel mimicking the hypothesized configuration prior to final substrate docking (i.e., the +1 and +2 sites are unoccupied), which is also the structure following a catalytic bond scission. No tendency was found for the protein to dissociate from or translate along the substrate surface during this initial simulation, nor to align with the direction of the cellulose chains. However, a tendency for the decrystallized cellodextrin to partially re-anneal into the cellulose surface hints that the arbitrary starting configuration selected was not ideal.

Cellulosic biofuels

The development of sustainable, low-carbon, liquid fuels from cellulosic biomass will require advances in many areas of science and engineering. This review describes the major topics of enquiry concerning cellulosic biofuels with an emphasis on those areas of research and development that include research problems of interest to plant biologists.