Towards a molecular-level theory of carbohydrate processivity in glycoside hydrolases

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.

Nanoscale Engineering of Designer Cellulosomes

Biocatalysts showcase the upper limit obtainable for high-speed molecular processing and transformation. Efforts to engineer functionality in synthetic nanostructured materials are guided by the increasing knowledge of evolving architectures, which enable controlled molecular motion and precise molecular recognition. The cellulosome is a biological nanomachine, which, as a fundamental component of the plant-digestion machinery from bacterial cells, has a key potential role in the successful development of environmentally-friendly processes to produce biofuels and fine chemicals from the breakdown of biomass waste. Here, the progress toward so-called "designer cellulosomes", which provide an elegant alternative to enzyme cocktails for lignocellulose breakdown, is reviewed. Particular attention is paid to rational design via computational modeling coupled with nanoscale characterization and engineering tools. Remaining challenges and potential routes to industrial application are put forward.

Aromatic-Mediated Carbohydrate Recognition in Processive Serratia marcescens Chitinases

Microorganisms use a host of enzymes, including processive glycoside hydrolases, to deconstruct recalcitrant polysaccharides to sugars. Processive glycoside hydrolases closely associate with polymer chains and repeatedly cleave glycosidic linkages without dissociating from the crystalline surface after each hydrolytic step; they are typically the most abundant enzymes in both natural secretomes and industrial cocktails by virtue of their significant hydrolytic potential. The ubiquity of aromatic residues lining the enzyme catalytic tunnels and clefts is a notable feature of processive glycoside hydrolases. We hypothesized that these aromatic residues have uniquely defined roles, such as substrate chain acquisition and binding in the catalytic tunnel, that are defined by their local environment and position relative to the substrate and the catalytic center. Here, we investigated this hypothesis with variants of Serratia marcescens family 18 processive chitinases ChiA and ChiB. We applied molecular simulation and free energy calculations to assess active site dynamics and ligand binding free energies. Isothermal titration calorimetry provided further insight into enthalpic and entropic contributions to ligand binding free energy. Thus, the roles of six aromatic residues, Trp-167, Trp-275, and Phe-396 in ChiA, and Trp-97, Trp-220, and Phe-190 in ChiB, have been examined. We observed that point mutation of the tryptophan residues to alanine results in unfavorable changes in the free energy of binding relative to wild-type. The most drastic effects were observed for residues positioned at the "entrances" of the deep substrate-binding clefts and known to be important for processivity. Interestingly, phenylalanine mutations in ChiA and ChiB had little to no effect on chito-oligomer binding, in accordance with the limited effects of their removal on chitinase functionality.

Processivity, Substrate Positioning, and Binding: The Role of Polar Residues in a Family 18 Glycoside Hydrolase

The enzymatic degradation of recalcitrant polysaccharides such as cellulose (β-1,4-linked glucose) and chitin (β-1,4-linked N-acetylglucosamine) by glycoside hydrolases (GHs) is of significant biological and economical importance. In nature, depolymerization is primarily accomplished by processive GHs, which remain attached to the substrate between subsequent hydrolytic reactions. Recent computational efforts have suggested that the processive ability of a GH is directly linked to the ligand binding free energy. The contribution of individual aromatic residues in the active site of these enzymes has been extensively studied. In this study, we offer the first experimental evidence confirming correlation of binding free energy and degree of processivity and evidence that polar residues are essential for maintaining processive ability. Exchanging Thr(276) with Ala in substrate binding subsite -2 in the processive ChiA of Serratia marcescens results in a decrease in both the enthalpy (2.6 and 3.8 kcal/mol) and free energy (0.5 and 2.2 kcal/mol) for the binding to the substrate (GlcNAc)6 and the inhibitor allosamidin, respectively, compared to that of the wild type. Moreover, the initial apparent processivity as measured by [(GlcNAc)2]/[GlcNAc] ratios (17.1 ± 0.4) and chitin degradation efficiency (20%) are greatly reduced for ChiA-T276A versus those of the wild type (30.1 ± 1.5 and 75%, respectively). Mutation of Arg(172) to Ala reduces the level of recognition and positioning of the substrate into the active site. Molecular dynamics simulations indicate ChiA-R172A behaves like the wild type, but the dynamics of ChiA-T276A are greatly influenced by mutation, which is reflective of their influence on processivity.

Slow Off-rates and Strong Product Binding Are Required for Processivity and Efficient Degradation of Recalcitrant Chitin by Family 18 Chitinases

Processive glycoside hydrolases are the key components of enzymatic machineries that decompose recalcitrant polysaccharides, such as chitin and cellulose. The intrinsic processivity (P(Intr)) of cellulases has been shown to be governed by the rate constant of dissociation from polymer chain (koff). However, the reported koff values of cellulases are strongly dependent on the method used for their measurement. Here, we developed a new method for determining koff, based on measuring the exchange rate of the enzyme between a non-labeled and a (14)C-labeled polymeric substrate. The method was applied to the study of the processive chitinase ChiA from Serratia marcescens. In parallel, ChiA variants with weaker binding of the N-acetylglucosamine unit either in substrate-binding site -3 (ChiA-W167A) or the product-binding site +1 (ChiA-W275A) were studied. Both ChiA variants showed increased off-rates and lower apparent processivity on α-chitin. The rate of the production of insoluble reducing groups on the reduced α-chitin was an order of magnitude higher than koff, suggesting that the enzyme can initiate several processive runs without leaving the substrate. On crystalline chitin, the general activity of the wild type enzyme was higher, and the difference was magnifying with hydrolysis time. On amorphous chitin, the variants clearly outperformed the wild type. A model is proposed whereby strong interactions with polymer in the substrate-binding sites (low off-rates) and strong binding of the product in the product-binding sites (high pushing potential) are required for the removal of obstacles, like disintegration of chitin microfibrils.

O-glycosylation effects on family 1 carbohydrate-binding module solution structures

Family 1 carbohydrate-binding modules (CBMs) are ubiquitous components of multimodular fungal enzymes that degrade plant cell wall polysaccharides and bind specifically to cellulose. Native glycosylation of family 1 CBMs has been shown to substantially impact multiple physical properties, including thermal and proteolytic stability and cellulose binding affinity. To gain molecular insights into the changes in CBM properties upon glycosylation, solution structures of two glycoforms of a Trichoderma reesei family 1 CBM were studied by NMR spectroscopy: a glycosylated family 1 CBM with a mannose group attached to both Thr1 and Ser3 and a second family 1 CBM with single mannose groups attached to Thr1, Ser3 and Ser14. The structures clearly reveal that monosaccharides at both Ser3 and Ser14 on family 1 CBMs present additional cellulose binding platforms, similar to well-characterized aromatic residues at the binding interface, which align to the cellulose surface. These results are in agreement with previous experimental work demonstrating that glycans at Ser3 and Ser14 impart significant improvements in binding affinity. Additionally, detailed analysis of the NMR structures and molecular simulations indicates that the protein backbone of the CBM is not significantly altered by attachment of monosaccharides, and that the mannose attached to Ser14 may be more flexible than the mannose at Ser3. Overall, the present study reveals how family 1 CBM structures are affected by covalent attachment of monosaccharides, which are likely important post-translational modifications of these common subdomains of fungal plant cell wall degrading enzymes.

DATABASE

Structural data have been deposited in the RCSB Protein Data Bank (PDB codes: 2MWJ and 2MWK) and the BioMagRes Bank (BMRB codes: 25331 and 25332) for CBM_M2 and CBM_M3, respectively.

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.

Characterization of Wall Teichoic Acid Degradation by the Bacteriophage ϕ29 Appendage Protein GP12 Using Synthetic Substrate Analogs

The genetics and enzymology of the biosynthesis of wall teichoic acid have been the extensively studied, however, comparatively little is known regarding the enzymatic degradation of this biological polymer. The GP12 protein from the Bacillus subtilis bacteriophage ϕ29 has been implicated as a wall teichoic acid hydrolase. We have studied the wall teichoic acid hydrolase activity of pure, recombinant GP12 using chemically defined wall teichoic acid analogs. The GP12 protein had potent wall teichoic acid hydrolytic activity in vitro and demonstrated ∼13-fold kinetic preference for glycosylated poly(glycerol phosphate) teichoic acid compared with non-glycosylated. Product distribution patterns suggested that the degradation of glycosylated polymers proceeded from the hydroxyl terminus of the polymer, whereas hydrolysis occurred at random sites in the non-glycosylated polymer. In addition, we present evidence that the GP12 protein possesses both phosphodiesterase and phosphomonoesterase activities.

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.

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.

Determination of the action modes of cellulases from hydrolytic profiles over a time course using fluorescence-assisted carbohydrate electrophoresis

Fluorescence-assisted carbohydrate electrophoresis (FACE) is a sensitive and simple method for the separation of oligosaccharides. It relies on labeling the reducing ends of oligosaccharides with a fluorophore, followed by PAGE. Concentration changes of oligosaccharides following hydrolysis of a carbohydrate polymer could be quantitatively measured continuously over time using the FACE method. Based on the quantitative analysis, we suggested that FACE was a relatively high-throughput, repeatable, and suitable method for the analysis of the action modes of cellulases. On account of the time courses of their hydrolytic profiles, the apparent processivity was used to show the different action modes of cellulases. Cellulases could be easily differentiated as exoglucanases, β-glucosidases, or endoglucanases. Moreover, endoglucanases from the same glycoside hydrolases family had a variety of apparent processivity, indicating the different modes of action. Endoglucanases with the same binding capacities and hydrolytic activities had similar oligosaccharide profiles, which aided in their classification. The hydrolytic profile of Trichoderma reesei Cel12A, an endoglucanases from T. reesei, contained glucose, cellobiose, and cellotriose, which revealed that it may have a new glucosidase activity, corresponding to that of EC 3.2.1.74. A hydrolysate study of a T. reesei Cel12A-N20A mutant demonstrated that the FACE method was sufficiently sensitive to detect the influence of a single-site mutation on enzymatic activity.

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.

Probing substrate interactions in the active tunnel of a catalytically deficient cellobiohydrolase (Cel7)

Cellobiohydrolases break down cellulose sequentially by sliding along the crystal surface with a single cellulose strand threaded through the catalytic tunnel of the enzyme. This so-called processive mechanism relies on a complex pattern of enzyme-substrate interactions, which need to be addressed in molecular descriptions of processivity and its driving forces. Here, we have used titration calorimetry to study interactions of cellooligosaccharides (COS) and a catalytically deficient variant (E212Q) of the enzyme Cel7A from Trichoderma reesei. This enzyme has ∼10 glucopyranose subsites in the catalytic tunnel, and using COS ligands with a degree of polymerization (DP) from 2 to 8, different regions of the tunnel could be probed. For COS ligands with a DP of 2-3 the binding constants were around 10(5) m(-1), and for longer ligands (DP 5-8) this value was ∼10(7) m(-1). Within each of these groups we did not find increased affinity as the ligands got longer and potentially filled more subsites. On the contrary, we found a small but consistent affinity loss as DP rose from 6 to 8, particularly at the higher investigated temperatures. Other thermodynamic functions (ΔH, ΔS, and ΔCp) decreased monotonously with both temperature and DP. Combined interpretation of these thermodynamic results and previously published structural data allowed assessment of an affinity profile along the length axis of the active tunnel.

Kinetics of cellobiohydrolase (Cel7A) variants with lowered substrate affinity

Cellobiohydrolases are exo-active glycosyl hydrolases that processively convert cellulose to soluble sugars, typically cellobiose. They effectively break down crystalline cellulose and make up a major component in industrial enzyme mixtures used for deconstruction of lignocellulosic biomass. Identification of the rate-limiting step for cellobiohydrolases remains controversial, and recent reports have alternately suggested either association (on-rate) or dissociation (off-rate) as the overall bottleneck. Obviously, this uncertainty hampers both fundamental mechanistic understanding and rational design of enzymes with improved industrial applicability. To elucidate the role of on- and off-rates, respectively, on the overall kinetics, we have expressed a variant in which a tryptophan residue (Trp-38) in the middle of the active tunnel has been replaced with an alanine. This mutation weakens complex formation, and the population of substrate-bound W38A was only about half of the wild type. Nevertheless, the maximal, steady-state rate was twice as high for the variant enzyme. It is argued that these opposite effects on binding and activity can be reconciled if the rate-limiting step is after the catalysis (i.e. in the dissociation process).

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.