Product binding varies dramatically between processive and nonprocessive cellulase enzymes

Impact of Synergy Partner Cel7B on Cel7A Binding Rates: Insights from Single-Molecule Data

Enzymatic degradation of cellulosic biomass is a well-established route for the sustainable production of biofuels, chemicals, and materials. A strategy employed by nature and industry to achieve an efficient degradation of cellulose is that cellobiohydrolases (or exocellulases), such as Cel7A, work synergistically with endoglucanases, such as Cel7B, to achieve the complete degradation of cellulose. However, a complete mechanistic understanding of this exo-endo synergy is still lacking. Here, we used single-molecule fluorescence microscopy to quantify the binding kinetics of Cel7A on cellulose when it is acting alone on the cellulose fibrils and in the presence of its synergy partner, the endoglucanase Cel7B. To this end, we used a fluorescently tagged Cel7A and studied its binding in the presence of the unlabeled Cel7B. This provided the single-molecule data necessary for the estimation of the rate constants of association kON and dissociation kOFF of Cel7A for the substrate. We show that the presence of Cel7B does not impact the dissociation rate constant, kOFF. But, the association rate of Cel7A decreases by a factor of 2 when Cel7B is present at a molar proportion of 10:1. This ratio has previously been shown to lead to synergy. This decrease in association rate is observed in a wide range of total enzyme concentrations, from sub nM to μM concentrations. This decrease in kON is consistent with the formation of cellulase clusters recently observed by others using atomic force microscopy.

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.

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.

Valorization of outer tunic of the marine filter feeder Ciona intestinalis towards the production of second-generation biofuel and prebiotic oligosaccharides

BACKGROUND

One of the sustainable development goals focuses on the biomass-based production as a replacement for fossil-based commodities. A novel feedstock with vast potentials is tunicate biomass, which can be pretreated and fermented in a similar way to lignocellulose. Ciona intestinalis is a marine filter feeder that is cultivated to produce fish feed. While the inner tissue body is used for feed production, the surrounding tunic remains as a cellulose-rich by-product, which can be further separated into outer and inner tunic. Ethanol production from organosolv-pretreated whole-tunic biomass was recently validated. The aim of the present study was to evaluate the potential of organosolv pretreated outer-tunic biomass for the production of biofuels and cellobiose that is a disaccharide with prebiotic potential.

RESULTS

As a result, 41.4 g/L of ethanol by Saccharomyces cerevisiae, corresponding to a 90.2% theoretical yield, was achieved under the optimal conditions when the tunicate biomass was pretreated at 195 °C for 60 min at a liquid-to-solid ratio of 50. In addition, cellobiose production by enzymatic hydrolysis of the pretreated tunicate biomass was demonstrated with a maximum conversion yield of 49.7 wt. %.

CONCLUSIONS

The utilisation of tunicate biomass offers an eco-friendly and sustainable alternative for value-added biofuels and chemicals. The cultivation of tunicate biomass in shallow coastal sea improves the quality of the water and ensures sustainable production of fish feed. Moreover, there is no competition for arable land, which leaves the latter available for food and feed production.

Computing Cellulase Kinetics with a Two-Domain Linear Interaction Energy Approach

While heterogeneous enzyme reactions play an essential role in both nature and green industries, computational predictions of their catalytic properties remain scarce. Recent experimental work demonstrated the applicability of the Sabatier principle for heterogeneous biocatalysis. This provides a simple relationship between binding strength and the catalytic rate and potentially opens a new way for inexpensive computational determination of kinetic parameters. However, broader implementation of this approach will require fast and reliable prediction of binding free energies of complex two-phase systems, and computational procedures for this are still elusive. Here, we propose a new framework for the assessment of the binding strengths of multidomain proteins, in general, and interfacial enzymes, in particular, based on an extended linear interaction energy (LIE) method. This two-domain LIE (2D-LIE) approach was successfully applied to predict binding and activation free energies of a diverse set of cellulases and resulted in robust models with high accuracy. Overall, our method provides a fast computational screening tool for cellulases that have not been experimentally characterized, and we posit that it may also be applicable to other heterogeneously acting biocatalysts.

Semiprocessive Hyperglycosylation of Adhesin by Bacterial Protein N-Glycosyltransferases

Processivity is an important feature of enzyme families such as DNA polymerases, polysaccharide synthases, and protein kinases, to ensure high fidelity in biopolymer synthesis and modification. Here, we reveal processive character in the family of cytoplasmic protein N-glycosyltransferases (NGTs). Through various activity assays, intact protein mass spectrometry, and proteomics analysis, we established that NGTs from nontypeable Haemophilus influenzae and Actinobacillus pleuropneumoniae modify an adhesin protein fragment in a semiprocessive manner. Molecular modeling studies suggest that the processivity arises from the shallow substrate binding groove in NGT, which promotes the sliding of the adhesin over the surface to allow further glycosylations without temporary dissociation. We hypothesize that the processive character of these bacterial protein glycosyltransferases is the mechanism to ensure multisite glycosylation of adhesins in vivo, thereby creating the densely glycosylated proteins necessary for bacterial self-aggregation and adherence to human cells, as a first step toward infection.

pH profiles of cellulases depend on the substrate and architecture of the binding region

Understanding the pH effect of cellulolytic enzymes is of great technological importance. In this study, we have examined the influence of pH on activity and stability for central cellulases (Cel7A, Cel7B, Cel6A from Trichoderma reesei, and Cel7A from Rasamsonia emersonii). We systematically changed pH from 2 to 7, temperature from 20°C to 70°C, and used both soluble (4-nitrophenyl β- d-lactopyranoside [pNPL]) and insoluble (Avicel) substrates at different concentrations. Collective interpretation of these data provided new insights. An unusual tolerance to acidic conditions was observed for both investigated Cel7As, but only on real insoluble cellulose. In contrast, pH profiles on pNPL were bell-shaped with a strong loss of activity both above and below the optimal pH for all four enzymes. On a practical level, these observations call for the caution of the common practice of using soluble substrates for the general characterization of pH effects on cellulase activity. Kinetic modeling of the experimental data suggested that the nucleophile of Cel7A experiences a strong downward shift in pK_a upon complexation with an insoluble substrate. This shift was less pronounced for Cel7B, Cel6A, and for Cel7A acting on the soluble substrate, and we hypothesize that these differences are related to the accessibility of water to the binding region of the Michaelis complex.

Binding Process and Free Energy Characteristics of Cellulose Chain into the Catalytic Domain of Cellobiohydrolase TrCel7A

It was observed in experiments that the catalytic domain (CD) of Trichoderma reesei Cel7A (TrCel7A) hydrolyzes crystalline cellulose in a processive manner, but the underlying binding mechanism is still unknown. Here, through replica-exchange molecular dynamics simulations, we find that the loading and sucking-in process of the cellulose chain into CD is entropy-driven and enthalpy-unfavorable, which firmly relate to the desolvation of the binding channel of CD. During the loading process, hydrophobic interactions play a dominant role because several aromatic residues have been identified to guide the cellulose chain processing. At the active site, a transition from enthalpy- to entropy-driven is detected for the driving force. Such a finding reveals the indispensability of the catalytic reaction of the glycosidic bond to provide the energy to drive the movements of the cellulose chain. Our study reveals the interaction pictures between the cellulose chain and TrCel7A at the atomic level, which helps better understand the catalytic mechanism of TrCel7A.

Thermodynamic Signatures of Substrate Binding for Three Thermobifida fusca Cellulases with Different Modes of Action

The enzymatic breakdown of recalcitrant polysaccharides is achieved by synergistic enzyme cocktails of glycoside hydrolases (GHs) and accessory enzymes. Many GHs are processive, meaning that they stay bound to the substrate between subsequent catalytic interactions. Cellulases are GHs that catalyze the hydrolysis of cellulose [β-1,4-linked glucose (Glc)]. Here, we have determined the relative subsite binding affinity for a glucose moiety as well as the thermodynamic signatures for (Glc)₆ binding to three of the seven cellulases produced by the bacterium Thermobifida fusca. TfCel48A is exo-processive, TfCel9A endo-processive, and TfCel5A endo-nonprocessive. Initial hydrolysis of (Glc)₅ and (Glc)₆ was performed in H₂¹⁸O enabling the incorporation of an ¹⁸O atom at the new reducing end anomeric carbon. A matrix-assisted laser desorption ionization time-of-flight mass spectrometry analysis of the products reveals the intensity ratios of otherwise identical ¹⁸O- and ¹⁶O-containing products to provide insight into how the substrate is placed during productive binding. The two processive cellulases have significant binding affinity in subsites where products dissociate during processive hydrolysis, aligned with a need to have a pushing potential to remove obstacles on the substrate. Moreover, we observed a correlation between processive ability and favorable binding free energy, as previously postulated. Upon ligand binding, the largest contribution to the binding free energy is desolvation for all three cellulases as determined by isothermal titration calorimetry. The two endo-active cellulases show a more favorable solvation entropy change compared to the exo-active cellulase, while the two processive cellulases have less favorable changes in binding enthalpy compared to the nonprocessive TfCel5A.

Systematic deletions in the cellobiohydrolase (CBH) Cel7A from the fungus Trichoderma reesei reveal flexible loops critical for CBH activity

Glycoside hydrolase family 7 (GH7) cellulases are some of the most efficient degraders of cellulose, making them particularly relevant for industries seeking to produce renewable fuels from lignocellulosic biomass. The secretome of the cellulolytic model fungus Trichoderma reesei contains two GH7s, termed TrCel7A and TrCel7B. Despite having high structural and sequence similarities, the two enzymes are functionally quite different. TrCel7A is an exolytic, processive cellobiohydrolase (CBH), with high activity on crystalline cellulose, whereas TrCel7B is an endoglucanase (EG) with a preference for more amorphous cellulose. At the structural level, these functional differences are usually ascribed to the flexible loops that cover the substrate-binding areas. TrCel7A has an extensive tunnel created by eight peripheral loops, and the absence of four of these loops in TrCel7B makes its catalytic domain a more open cleft. To investigate the structure-function relationships of these loops, here we produced and kinetically characterized several variants in which four loops unique to TrCel7A were individually deleted to resemble the arrangement in the TrCel7B structure. Analysis of a range of kinetic parameters consistently indicated that the B2 loop, covering the substrate-binding subsites -3 and -4 in TrCel7A, was a key determinant for the difference in CBH- or EG-like behavior between TrCel7A and TrCel7B. Conversely, the B3 and B4 loops, located closer to the catalytic site in TrCel7A, were less important for these activities. We surmise that these results could be useful both in further mechanistic investigations and for guiding engineering efforts of this industrially important enzyme family.

Predominant Nonproductive Substrate Binding by Fungal Cellobiohydrolase I and Implications for Activity Improvement

Enzymatic conversion of the most abundant renewable source of organic compounds, cellulose to fermentable sugars is attractive for production of green fuels and chemicals. The major component of industrial enzyme systems, cellobiohydrolase I from Hypocrea jecorina (Trichoderma reesei) (HjCel7A) processively splits disaccharide units from the reducing ends of tightly packed cellulose chains. HjCel7A consists of a catalytic domain (CD) and a carbohydrate-binding module (CBM) separated by a linker peptide. A tunnel-shaped substrate-binding site in the CD includes nine subsites for β-d-glucose units, seven of which (-7 to -1) precede the catalytic center. Low catalytic activity of Cel7A is the bottleneck and the primary target for improvement. Here it is shown for the first time that, in spite of much lower apparent k_cat of HjCel7A at the hydrolysis of β-1,4-glucosidic linkages in the fluorogenic cellotetra- and -pentaose compared to the structurally related endoglucanase I (HjCel7B), the specificity constants (catalytic efficiency) k_cat /K_m for both enzymes are almost equal in these reactions. The observed activity difference appears from strong nonproductive substrate binding by HjCel7A, particularly significant for MU-β-cellotetraose (MUG₄ ). Interaction of substrates with the subsites -6 and -5 proximal to the nonconserved Gln101 residue in HjCel7A decreases K_m,ap by >1500 times. HjCel7A can be nonproductively bound onto cellulose surface with K_d ≈2-9 nM via CBM and CD that captures six terminal glucose units of cellulose chain. Decomposition of this nonproductive complex can determine the rate of cellulose conversion. MUG₄ is a promising substrate to select active cellobiohydrolase I variants with reduced nonproductive substrate binding.

The inhibition of hemicellulosic sugars on cellulose hydrolysis are highly dependant on the cellulase productive binding, processivity, and substrate surface charges

In this study, the influence of major hemicellulosic sugars (mannose and xylose) on cellulose hydrolysis and major enzyme activities were evaluated by using both commercial enzyme cocktail and purified cellulase monocomponents over a "library" of cellulosic substrates. Surprisingly, the results showed that unlike glucose, mannose/xylose did not inhibit individual cellulase activities but significantly decreased their hydrolytic performance on cellulose substrates. When various enzyme-substrate interactions (e.g. adsorption/desorption, productive binding, and processive moving) were evaluated, it appeared that these hemicellulosic sugars significantly reduced the productive binding and processivity of Cel7A, which in turn limited cellulase hydrolytic efficacy. Among a range of major cellulose characteristics (e.g. crystallinity, degree of polymerization, accessibility, and surface charges), the acid group content of the cellulosic substrates seemed to be the main driver that determined the extent of hemicellulosic sugar inhibition. Our results provided new insights for better understanding the sugar inhibition mechanisms of cellulose hydrolysis.

Insight into the process of product expulsion in cellobiohydrolase Cel6A from Trichoderma reesei by computational modeling

Glycoside hydrolase cellulase family 6 from Trichoderma reesei (TrCel6A) is an important cellobiohydrolase to hydrolyze cellooligosaccharide into cellobiose. The knowledge of enzymatic mechanisms is critical for improving the conversion efficiency of cellulose into ethanol or other chemicals. However, the process of product expulsion, a key component of enzymatic depolymerization, from TrCel6A has not yet been described in detail. Here, conventional molecular dynamics and steered molecular dynamics (SMD) were applied to study product expulsion from TrCel6A. Tyr103 may be a crucial residue in product expulsion given that it exhibits two different posthydrolytic conformations. In one conformation, Tyr103 rotates to open the -3 subsite. However, Tyr103 does not rotate in the other conformation. Three different routes for product expulsion were proposed on the basis of the two different conformations. The total energy barriers of the three routes were calculated through SMD simulations. The total energy barrier of product expulsion through Route 1, in which Tyr103 does not rotate, was 22.2 kcal·mol^-1. The total energy barriers of product expulsion through Routes 2 and 3, in which Tyr103 rotates to open the -3 subsite, were 10.3 and 14.4 kcal·mol^-1, respectively. Therefore, Routes 2 and 3 have lower energy barriers than Route 1, and Route 2 is the thermodynamically optimal route for product expulsion. Consequently, the rotation of Tyr103 may be crucial for product release from TrCel6A. Results of this work have potential applications in cellulase engineering.

Mechanistic kinetic models of enzymatic cellulose hydrolysis-A review

Bioconversion of lignocellulose forms the basis for renewable, advanced biofuels, and bioproducts. Mechanisms of hydrolysis of cellulose by cellulases have been actively studied for nearly 70 years with significant gains in understanding of the cellulolytic enzymes. Yet, a full mechanistic understanding of the hydrolysis reaction has been elusive. We present a review to highlight new insights gained since the most recent comprehensive review of cellulose hydrolysis kinetic models by Bansal et al. (2009) Biotechnol Adv 27:833-848. Recent models have taken a two-pronged approach to tackle the challenge of modeling the complex heterogeneous reaction-an enzyme-centric modeling approach centered on the molecularity of the cellulase-cellulose interactions to examine rate limiting elementary steps and a substrate-centric modeling approach aimed at capturing the limiting property of the insoluble cellulose substrate. Collectively, modeling results suggest that at the molecular-scale, how rapidly cellulases can bind productively (complexation) and release from cellulose (decomplexation) is limiting, while the overall hydrolysis rate is largely insensitive to the catalytic rate constant. The surface area of the insoluble substrate and the degrees of polymerization of the cellulose molecules in the reaction both limit initial hydrolysis rates only. Neither enzyme-centric models nor substrate-centric models can consistently capture hydrolysis time course at extended reaction times. Thus, questions of the true reaction limiting factors at extended reaction times and the role of complexation and decomplexation in rate limitation remain unresolved. Biotechnol. Bioeng. 2017;114: 1369-1385. © 2017 Wiley Periodicals, Inc.

Human Chitotriosidase Is an Endo-Processive Enzyme

Human chitotriosidase (HCHT) is involved in immune response to chitin-containing pathogens in humans. The enzyme is able to degrade chitooligosaccharides as well as crystalline chitin. The catalytic domain of HCHT is connected to the carbohydrate binding module (CBM) through a flexible hinge region. In humans, two active isoforms of HCHT are found-the full length enzyme and its truncated version lacking CBM and the hinge region. The active site architecture of HCHT is reminiscent to that of the reducing-end exo-acting processive chitinase ChiA from bacterium Serratia marcescens (SmChiA). However, the presence of flexible hinge region and occurrence of two active isoforms are reminiscent to that of non-processive endo-chitinase from S. marcescens, SmChiC. Although the studies on soluble chitin derivatives suggest the endo-character of HCHT, the mode of action of the enzyme on crystalline chitin is not known. Here, we made a thorough characterization of HCHT in terms of the mode of action, processivity, binding, and rate constants for the catalysis and dissociation using α-chitin as substrate. HCHT efficiently released the end-label from reducing-end labelled chitin and had also high probability (95%) of endo-mode initiation of processive run. These results qualify HCHT as an endo-processive enzyme. Processivity and the rate constant of dissociation of HCHT were found to be in-between those, characteristic to processive exo-enzymes, like SmChiA and randomly acting non-processive endo-enzymes, like SmChiC. Apart from increasing the affinity for chitin, CBM had no major effect on kinetic properties of HCHT.

Engineering of d-fructose-6-phosphate aldolase A for improved activity towards cinnamaldehyde

d-Fructose-6-phosphate aldolase A (FSAA) from Escherichia coli was engineered for enhanced catalytic efficiency towards cinnamaldehyde.

Anomeric Selectivity and Product Profile of a Processive Cellulase

Cellobiohydrolases (CBHs) make up an important group of enzymes for both natural carbon cycling and industrial deconstruction of lignocellulosic biomass. The consecutive hydrolysis of one cellulose strand relies on an intricate pattern of enzyme-substrate interactions in the long, tunnel-shaped binding site of the CBH. In this work, we have investigated the initial complexation mode with cellulose of the most thoroughly studied CBH, Cel7A from Hypocrea jecorina (HjCel7A). We found that HjCel7A predominantly produces glucose when it initiates a processive run on insoluble microcrystalline cellulose, confirming the validity of an even and odd product ratio as an estimate of processivity. Moreover, the glucose released from cellulose was predominantly α-glucose. A link between the initial binding mode of the enzyme and the reducing end configuration was investigated by inhibition studies with the two anomers of cellobiose. A clear preference for β-cellobiose in product binding site +2 was observed for HjCel7A, but not the homologous endoglucanase, HjCe7B. Possible relationships between this anomeric preference in the product site and the prevalence of odd-numbered initial-cut products are discussed, and a correlation between processivity and anomer selectivity is proposed.

Who's on base? Revealing the catalytic mechanism of inverting family 6 glycoside hydrolases

In several important classes of inverting carbohydrate-active enzymes, the identity of the catalytic base remains elusive, including in family 6 Glycoside Hydrolase (GH6) enzymes, which are key components of cellulase cocktails for cellulose depolymerization. Despite many structural and kinetic studies with both wild-type and mutant enzymes, especially on the Trichoderma reesei (Hypocrea jecorina) GH6 cellulase (TrCel6A), the catalytic base in the single displacement inverting mechanism has not been definitively identified in the GH6 family. Here, we employ transition path sampling to gain insight into the catalytic mechanism, which provides unbiased atomic-level understanding of key order parameters involved in cleaving the strong glycosidic bond. Our hybrid quantum mechanics and molecular mechanics (QM/MM) simulations reveal a network of hydrogen bonding that aligns two active site water molecules that play key roles in hydrolysis: one water molecule drives the reaction by nucleophilic attack on the substrate and a second shuttles a proton to the putative base (D175) via a short water wire. We also investigated the case where the putative base is mutated to an alanine, an enzyme that is experimentally still partially active. The simulations predict that proton hopping along a water wire via a Grotthuss mechanism provides a mechanism of catalytic rescue. Further simulations reveal that substrate processive motion is 'driven' by strong electrostatic interactions with the protein at the product sites and that the -1 sugar adopts a 2SO ring configuration as it reaches its binding site. This work thus elucidates previously elusive steps in the processive catalytic mechanism of this important class of enzymes.

Strategies to reduce end-product inhibition in family 48 glycoside hydrolases

Family 48 cellobiohydrolases are some of the most abundant glycoside hydrolases in nature. They are able to degrade cellulosic biomass and therefore serve as good enzyme candidates for biofuel production. Family 48 cellulases hydrolyze cellulose chains via a processive mechanism, and produce end products composed primarily of cellobiose as well as other cellooligomers (dp ≤ 4). The challenge of utilizing cellulases in biofuel production lies in their extremely slow turnover rate. A factor contributing to the low enzyme activity is suggested to be product binding to enzyme and the resulting performance inhibition. In this study, we quantitatively evaluated the product inhibitory effect of four family 48 glycoside hydrolases using molecular dynamics simulations and product expulsion free-energy calculations. We also suggested a series of single mutants of the four family 48 glycoside hydrolases with theoretically reduced level of product inhibition. The theoretical calculations provide a guide for future experimental studies designed to produce mutant cellulases with enhanced activity.

Reassembly and co-crystallization of a family 9 processive endoglucanase from its component parts: structural and functional significance of the intermodular linker

Non-cellulosomal processive endoglucanase 9I (Cel9I) from Clostridium thermocellum is a modular protein, consisting of a family-9 glycoside hydrolase (GH9) catalytic module and two family-3 carbohydrate-binding modules (CBM3c and CBM3b), separated by linker regions. GH9 does not show cellulase activity when expressed without CBM3c and CBM3b and the presence of the CBM3c was previously shown to be essential for endoglucanase activity. Physical reassociation of independently expressed GH9 and CBM3c modules (containing linker sequences) restored 60-70% of the intact Cel9I endocellulase activity. However, the mechanism responsible for recovery of activity remained unclear. In this work we independently expressed recombinant GH9 and CBM3c with and without their interconnecting linker in Escherichia coli. We crystallized and determined the molecular structure of the GH9/linker-CBM3c heterodimer at a resolution of 1.68 Å to understand the functional and structural importance of the mutual spatial orientation of the modules and the role of the interconnecting linker during their re-association. Enzyme activity assays and isothermal titration calorimetry were performed to study and compare the effect of the linker on the re-association. The results indicated that reassembly of the modules could also occur without the linker, albeit with only very low recovery of endoglucanase activity. We propose that the linker regions in the GH9/CBM3c endoglucanases are important for spatial organization and fixation of the modules into functional enzymes.

Alleviating product inhibition in cellulase enzyme Cel7A

Enzymes that degrade cellulose into glucose are one of the most expensive components of processes for converting cellulosic biomass to fuels and chemicals. Cellulase enzyme Cel7A is the most abundant enzyme naturally employed by fungi to depolymerize cellulose, and like other cellulases is inhibited by its product, cellobiose. There is thus great economic incentive for minimizing the detrimental effects of product inhibition on Cel7A. In this work, we experimentally generated 10 previously proposed site-directed mutant Cel7A enzymes expected to have reduced cellobiose binding energies (the majority of mutations were to alanine). We then tested their resilience to cellobiose as well as their hydrolytic activities on microcrystalline cellulose. Although every mutation tested conferred reduced product inhibition (and abolished it for some), our results confirm a trade-off between Cel7A tolerance to cellobiose and enzymatic activity: Reduced product inhibition was accompanied by lower overall enzymatic activity on crystalline cellulose for the mutants tested. The tempering effect of mutations on inhibition was nearly constant despite relatively large differences in activities of the mutants. Our work identifies an amino acid in the Cel7A product binding site of interest for further mutational studies, and highlights both the challenge and the opportunity of enzyme engineering toward improving product tolerance in Cel7A.

Effects of lytic polysaccharide monooxygenase oxidation on cellulose structure and binding of oxidized cellulose oligomers to cellulases

In nature, polysaccharide glycosidic bonds are cleaved by hydrolytic enzymes for a vast array of biological functions. Recently, a new class of enzymes that utilize an oxidative mechanism to cleave glycosidic linkages was discovered; these enzymes are called lytic polysaccharide monooxygenases (LPMO). These oxidative enzymes are synergistic with cocktails of hydrolytic enzymes and are thought to act primarily on crystalline regions, in turn providing new sites of productive attachment and detachment for processive hydrolytic enzymes. In the case of cellulose, the homopolymer of β-1,4-d-glucose, enzymatic oxidation occurs at either the reducing end or the nonreducing end of glucose, depending on enzymatic specificity, and results in the generation of oxidized chemical substituents at polymer chain ends. LPMO oxidation of cellulose is thought to produce either a lactone at the reducing end of glucose that can spontaneously or enzymatically convert to aldonic acid or 4-keto-aldose at the nonreducing end that may further oxidize to a geminal diol. Here, we use molecular simulation to examine the effect of oxidation on the structure of crystalline cellulose. The simulations highlight variations in behaviors depending on the chemical identity of the oxidized species and its location within the cellulose fibril, as different oxidized species introduce steric effects that disrupt local crystallinity and in some cases reduce the work needed for polymer decrystallization. Reducing-end oxidations are easiest to decrystallize when located at the end of the fibril, whereas nonreducing end oxidations readily decrystallize from internal cleavage sites despite their lower solvent accessibility. The differential in decrystallization free energy suggests a molecular mechanism consistent with experimentally observed LPMO/cellobiohydrolase synergy. Additionally, the soluble oxidized cellobiose products released by hydrolytic cellulases may bind to the active sites of cellulases with different affinities relative to cellobiose itself, which potentially affects hydrolytic turnover through product inhibition. To examine the effect of oxidation on cello-oligomer binding, we use thermodynamic integration to compute the relative change in binding free energy between the hydrolyzed and oxidized products in the active site of Family 7 and Family 6 processive glycoside hydrolases, Trichoderma reesei Cel7A and Cel6A, which are key industrial cellulases and commonly used model systems for fungal cellulases. Our results suggest that the equilibrium between the two reducing end oxidized products, favoring the linear aldonic acid, may increase product inhibition, which would in turn reduce processive substrate turnover. In the case of LMPO action at the nonreducing end, oxidation appears to lower affinity with the nonreducing end specific cellulase, reducing product inhibition and potentially promoting processive cellulose turnover. Overall, this suggests that oxidation of recalcitrant polysaccharides by LPMOs accelerates degradation not only by increasing the concentration of chain termini but also by reducing decrystallization work, and that product inhibition may be somewhat reduced as a result.

Cellobiohydrolase and endoglucanase respond differently to surfactants during the hydrolysis of cellulose

BACKGROUND

Non-ionic surfactants such as polyethylene glycol (PEG) can increase the glucose yield obtained from enzymatic saccharification of lignocellulosic substrates. Various explanations behind this effect include the ability of PEG to increase the stability of the cellulases, decrease non-productive cellulase adsorption to the substrate, and increase the desorption of enzymes from the substrate. Here, using lignin-free model substrates, we propose that PEG also alters the solvent properties, for example, water, leading the cellulases to increase hydrolysis yields.

RESULTS

The effect of PEG differs for the individual cellulases. During hydrolysis of Avicel and PASC with a processive monocomponent exo-cellulase cellobiohydrolase (CBH) I, the presence of PEG leads to an increase in the final glucose concentration, while PEG caused no change in glucose production with a non-processive endoglucanase (EG). Also, no effect of PEG was seen on the activity of β-glucosidases. While PEG has a small effect on the thermostability of both cellulases, only the activity of CBH I increases with PEG. Using commercial enzyme mixtures, the hydrolysis yields increased with the addition of PEG. In parallel, we observed that the relaxation time of the hydrolysis liquid phase, as measured by LF-NMR, directly correlated with the final glucose yield. PEG was able to boost the glucose production even in highly concentrated solutions of up to 150 g/L of glucose.

CONCLUSIONS

The hydrolysis boosting effect of PEG appears to be specific for CBH I. The mechanism could be due to an increase in the apparent activity of the enzyme on the substrate surface. The addition of PEG increases the relaxation time of the liquid-phase water, which from the data presented points towards a mechanism related to PEG-water interactions rather than PEG-protein or PEG-substrate interactions.

The predominant molecular state of bound enzyme determines the strength and type of product inhibition in the hydrolysis of recalcitrant polysaccharides by processive enzymes

Processive enzymes are major components of the efficient enzyme systems that are responsible for the degradation of the recalcitrant polysaccharides cellulose and chitin. Despite intensive research, there is no consensus on which step is rate-limiting for these enzymes. Here, we performed a comparative study of two well characterized enzymes, the cellobiohydrolase Cel7A from Hypocrea jecorina and the chitinase ChiA from Serratia marcescens. Both enzymes were inhibited by their disaccharide product, namely chitobiose for ChiA and cellobiose for Cel7A. The products behaved as noncompetitive inhibitors according to studies using the (14)C-labeled crystalline polymeric substrates (14)C chitin nanowhiskers and (14)C-labeled bacterial microcrystalline cellulose for ChiA and Cel7A, respectively. The resulting observed Ki (obs) values were 0.45 ± 0.08 mm for ChiA and 0.17 ± 0.02 mm for Cel7A. However, in contrast to ChiA, the Ki (obs) of Cel7A was an order of magnitude higher than the true Ki value governed by the thermodynamic stability of the enzyme-inhibitor complex. Theoretical analysis of product inhibition suggested that the inhibition strength and pattern can be accounted for by assuming different rate-limiting steps for ChiA and Cel7A. Measuring the population of enzymes whose active site was occupied by a polymer chain revealed that Cel7A was bound predominantly via its active site. Conversely, the active-site-mediated binding of ChiA was slow, and most ChiA exhibited a free active site, even when the substrate concentration was saturating for the activity. Collectively, our data suggest that complexation with the polymer chain is rate-limiting for ChiA, whereas Cel7A is limited by dissociation.

Glycosylation of Cellulases: Engineering Better Enzymes for Biofuels

Cellulose in plant cell walls is the largest reservoir of renewable carbon on Earth. The saccharification of cellulose from plant biomass into soluble sugars can be achieved using fungal and bacterial cellulolytic enzymes, cellulases, and further converted into fuels and chemicals. Most fungal cellulases are both N- and O-glycosylated in their native form, yet the consequences of glycosylation on activity and structure are not fully understood. Studying protein glycosylation is challenging as glycans are extremely heterogeneous, stereochemically complex, and glycosylation is not under direct genetic control. Despite these limitations, many studies have begun to unveil the role of cellulase glycosylation, especially in the industrially relevant cellobiohydrolase from Trichoderma reesei, Cel7A. Glycosylation confers many beneficial properties to cellulases including enhanced activity, thermal and proteolytic stability, and structural stabilization. However, glycosylation must be controlled carefully as such positive effects can be dampened or reversed. Encouragingly, methods for the manipulation of glycan structures have been recently reported that employ genetic tuning of glycan-active enzymes expressed from homogeneous and heterologous fungal hosts. Taken together, these studies have enabled new strategies for the exploitation of protein glycosylation for the production of enhanced cellulases for biofuel production.

Enzymatic degradation of (ligno)cellulose

Glycoside-degrading enzymes play a dominant role in the biochemical conversion of cellulosic biomass into low-price biofuels and high-value-added chemicals. New insight into protein functions and substrate structures, the kinetics of recognition, and degradation events has resulted in a substantial improvement of our understanding of cellulose degradation.

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.

Loop motions important to product expulsion in the Thermobifida fusca glycoside hydrolase family 6 cellobiohydrolase from structural and computational studies

Cellobiohydrolases (CBHs) are typically major components of natural enzyme cocktails for biomass degradation. Their active sites are enclosed in a tunnel, enabling processive hydrolysis of cellulose chains. Glycoside hydrolase Family 6 (GH6) CBHs act from nonreducing ends by an inverting mechanism and are present in many cellulolytic fungi and bacteria. The bacterial Thermobifida fusca Cel6B (TfuCel6B) exhibits a longer and more enclosed active site tunnel than its fungal counterparts. Here, we determine the structures of two TfuCel6B mutants co-crystallized with cellobiose, D274A (catalytic acid), and the double mutant D226A/S232A, which targets the putative catalytic base and a conserved serine that binds the nucleophilic water. The ligand binding and the structure of the active site are retained when compared with the wild type structure, supporting the hypothesis that these residues are directly involved in catalysis. One structure exhibits crystallographic waters that enable construction of a model of the α-anomer product after hydrolysis. Interestingly, the product sites of TfuCel6B are completely enclosed by an "exit loop" not present in fungal GH6 CBHs and by an extended "bottom loop". From the structures, we hypothesize that either of the loops enclosing the product subsites in the TfuCel6B active site tunnel must open substantially for product release. With simulation, we demonstrate that both loops can readily open to allow product release with equal probability in solution or when the enzyme is engaged on cellulose. Overall, this study reveals new structural details of GH6 CBHs likely important for functional differences among enzymes from this important family.

Adsorption and hydrolytic activity of the polycatalytic cellulase nanocomplex on cellulose

The formation of polycatalytic enzyme complexes may enhance the effectiveness of enzymes due to improved substrate interaction and synergistic actions of multiple enzymes in proximity. Much effort has been made to develop highly efficient polycatalytic cellulase complexes by immobilizing cellulases on low-cost polymer or nanoparticle scaffolds, aiming at their potential applications in biomass conversion to fuels. However, some key cellulases carry out the hydrolytic reaction on crystalline cellulose in a directional, processive manner. A large, artificial polycatalytic complex is unlikely to undergo a highly coordinated motion to slide on the cellulose surface as a whole unit. The mechanism underlying the activity enhancements observed in some artificial cellulase complexes and the limit of this approach remain elusive. Herein, we report the synthesis of polycatalytic cellulase complexes bound to colloidal polymer nanoparticles with a magnetic core and describe their unique adsorption, hydrolytic activities, and motions on cellulose. The polycatalytic clusters of cellulases on colloidal polymers show an increased rate of hydrolytic reactions on cellulose, but this was observed mainly at relatively low cellulase-to-cellulose ratios. Enhanced efficiency is mainly attributed to increased local concentrations of cellulases on the scaffolds and their polyvalent interactions with cellulose. However, once bound, the polycatalytic complexes can only carry out reactions locally and are not capable of relocating to new sites rapidly due to their lack of long-range surface mobility and their extremely tight binding. The development of highly optimized polycatalytic complexes may arise by developing novel nanoscaffolds that induce concerted motion of the complex as a whole.

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.

Binding and movement of individual Cel7A cellobiohydrolases on crystalline cellulose surfaces revealed by single-molecule fluorescence imaging

The efficient catalytic conversion of biomass to bioenergy would meet a large portion of energy requirements in the near future. A crucial step in this process is the enzyme-catalyzed hydrolysis of cellulose to glucose that is then converted into fuel such as ethanol by fermentation. Here we use single-molecule fluorescence imaging to directly monitor the movement of individual Cel7A cellobiohydrolases from Trichoderma reesei (TrCel7A) on the surface of insoluble cellulose fibrils to elucidate molecular level details of cellulase activity. The motion of multiple, individual TrCel7A cellobiohydrolases was simultaneously recorded with ∼15-nm spatial resolution. Time-resolved localization microscopy provides insights on the activity of TrCel7A on cellulose and informs on nonproductive binding and diffusion. We measured single-molecule residency time distributions of TrCel7A bound to cellulose both in the presence of and absence of cellobiose the major product and a potent inhibitor of Cel7A activity. Combining these results with a kinetic model of TrCel7A binding provides microscopic insight into interactions between TrCel7A and the cellulose substrate.

Computational investigation of the pH dependence of loop flexibility and catalytic function in glycoside hydrolases

Cellulase enzymes cleave glycosidic bonds in cellulose to produce cellobiose via either retaining or inverting hydrolysis mechanisms, which are significantly pH-dependent. Many fungal cellulases function optimally at pH ~5, and their activities decrease dramatically at higher or lower pH. To understand the molecular-level implications of pH in cellulase structure, we use a hybrid, solvent-based, constant pH molecular dynamics method combined with pH-based replica exchange to determine the pK(a) values of titratable residues of a glycoside hydrolase (GH) family 6 cellobiohydrolase (Cel6A) and a GH family 7 cellobiohydrolase (Cel7A) from the fungus Hypocrea jecorina. For both enzymes, we demonstrate that a bound substrate significantly affects the pKa values of the acid residues at the catalytic center. The calculated pK(a) values of catalytic residues confirm their proposed roles from structural studies and are consistent with the experimentally measured apparent pKa values. Additionally, GHs are known to impart a strained pucker conformation in carbohydrate substrates in active sites for catalysis, and results from free energy calculations combined with constant pH molecular dynamics suggest that the correct ring pucker is stable near the optimal pH for both Cel6A and Cel7A. Much longer molecular dynamics simulations of Cel6A and Cel7A with fixed protonation states based on the calculated pK(a) values suggest that pH affects the flexibility of tunnel loops, which likely affects processivity and substrate complexation. Taken together, this work demonstrates several molecular-level effects of pH on GH enzymes important for cellulose turnover in the biosphere and relevant to biomass conversion processes.

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.

Improving the performance of enzymes in hydrolysis of high solids paper pulp derived from MSW

BACKGROUND

The research aimed to improve the overall conversion efficiency of the CTec® family of enzymes by identifying factors that lead to inhibition and seeking methods to overcome these through process modification and manipulation. The starting material was pulp derived from municipal solid waste and processed in an industrial-scale washing plant.

RESULTS

Analysis of the pulp by acid hydrolysis showed a ratio of 55 : 12 : 6 : 24 : 3 of glucan : xylan : araban/galactan/mannan : lignin : ash. At high total solids content (>18.5% TS) single-stage enzyme hydrolysis gave a maximum glucan conversion of 68%. It was found that two-stage hydrolysis could give higher conversion if sugar inhibition was removed by an intermediate fermentation step between hydrolysis stages. This, however, was not as effective as direct removal of the sugar products, including xylose, by washing of the residual pulp at pH 5. This improved the water availability and allowed reactivation of the pulp-bound enzymes. Inhibition of enzyme activity could further be alleviated by replenishment of β-glucosidase which was shown to be removed during the wash step.

CONCLUSIONS

The two-stage hydrolysis process developed could give an overall glucan conversion of 88%, with an average glucose concentration close to 8% in 4 days, thus providing an ideal starting point for ethanol fermentation with a likely yield of 4 wt%. This is a significant improvement over a single-step process. This hydrolysis configuration also provides the potential to recover the sugars associated with residual solids which are diluted when washing hydrolysed pulp.

Endo-exo synergism in cellulose hydrolysis revisited

Synergistic cooperation of different enzymes is a prerequisite for efficient degradation of cellulose. The conventional mechanistic interpretation of the synergism between randomly acting endoglucanases (EGs) and chain end-specific processive cellobiohydrolases (CBHs) is that EG-generated new chain ends on cellulose surface serve as starting points for CBHs. Here we studied the hydrolysis of bacterial cellulose (BC) by CBH TrCel7A and EG TrCel5A from Trichoderma reesei under both single-turnover and "steady state" conditions. Unaccountable by conventional interpretation, the presence of EG increased the rate constant of TrCel7A-catalyzed hydrolysis of BC in steady state. At optimal enzyme/substrate ratios, the "steady state" rate of synergistic hydrolysis became limited by the velocity of processive movement of TrCel7A on BC. A processivity value of 66 ± 7 cellobiose units measured for TrCel7A on (14)C-labeled BC was close to the leveling off degree of polymerization of BC, suggesting that TrCel7A cannot pass through the amorphous regions on BC and stalls. We propose a mechanism of endo-exo synergism whereby the degradation of amorphous regions by EG avoids the stalling of TrCel7A and leads to its accelerated recruitment. Hydrolysis of pretreated wheat straw suggested that this mechanism of synergism is operative also in the degradation of lignocellulose. Although both mechanisms of synergism are used in parallel, the contribution of conventional mechanism is significant only at high enzyme/substrate ratios.