Thermoprotection by a cell membrane-localized metacaspase in a green alga.

Caspases are restricted to animals, while other organisms, including plants, possess metacaspases (MCAs), a more ancient and broader class of structurally related yet biochemically distinct proteases. Our current understanding of plant MCAs is derived from studies in streptophytes, and mostly in Arabidopsis (Arabidopsis thaliana) with 9 MCAs with partially redundant activities. In contrast to streptophytes, most chlorophytes contain only 1 or 2 uncharacterized MCAs, providing an excellent platform for MCA research. Here we investigated CrMCA-II, the single type-II MCA from the model chlorophyte Chlamydomonas (Chlamydomonas reinhardtii). Surprisingly, unlike other studied MCAs and similar to caspases, CrMCA-II dimerizes both in vitro and in vivo. Furthermore, activation of CrMCA-II in vivo correlated with its dimerization. Most of CrMCA-II in the cell was present as a proenzyme (zymogen) attached to the plasma membrane (PM). Deletion of CrMCA-II by genome editing compromised thermotolerance, leading to increased cell death under heat stress. Adding back either wild-type or catalytically dead CrMCA-II restored thermoprotection, suggesting that its proteolytic activity is dispensable for this effect. Finally, we connected the non-proteolytic role of CrMCA-II in thermotolerance to the ability to modulate PM fluidity. Our study reveals an ancient, MCA-dependent thermotolerance mechanism retained by Chlamydomonas and probably lost during the evolution of multicellularity.

Structure-function study of a Ca2+-independent metacaspase involved in lateral root emergence.

Metacaspases are part of an evolutionarily broad family of multifunctional cysteine proteases, involved in disease and normal development. As the structure-function relationship of metacaspases remains poorly understood, we solved the X-ray crystal structure of an Arabidopsis thaliana type II metacaspase (AtMCA-IIf) belonging to a particular subgroup not requiring calcium ions for activation. To study metacaspase activity in plants, we developed an in vitro chemical screen to identify small molecule metacaspase inhibitors and found several hits with a minimal thioxodihydropyrimidine-dione structure, of which some are specific AtMCA-IIf inhibitors. We provide mechanistic insight into the basis of inhibition by the TDP-containing compounds through molecular docking onto the AtMCA-IIf crystal structure. Finally, a TDP-containing compound (TDP6) effectively hampered lateral root emergence in vivo, probably through inhibition of metacaspases specifically expressed in the endodermal cells overlying developing lateral root primordia. In the future, the small compound inhibitors and crystal structure of AtMCA-IIf can be used to study metacaspases in other species, such as important human pathogens, including those causing neglected diseases.

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.

The dissociation mechanism of processive cellulases.

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

The hydrolysis mechanism of a GH45 cellulase and its potential relation to lytic transglycosylase and expansin function.

Family 45 glycoside hydrolases (GH45) are endoglucanases that are integral to cellulolytic secretomes, and their ability to break down cellulose has been successfully exploited in textile and detergent industries. In addition to their industrial relevance, understanding the molecular mechanism of GH45-catalyzed hydrolysis is of fundamental importance because of their structural similarity to cell wall-modifying enzymes such as bacterial lytic transglycosylases (LTs) and expansins present in bacteria, plants, and fungi. Our understanding of the catalytic itinerary of GH45s has been incomplete because a crystal structure with substrate spanning the -1 to +1 subsites is currently lacking. Here we constructed and validated a putative Michaelis complex in silico and used it to elucidate the hydrolytic mechanism in a GH45, Cel45A from the fungus Humicola insolens, via unbiased simulation approaches. These molecular simulations revealed that the solvent-exposed active-site architecture results in lack of coordination for the hydroxymethyl group of the substrate at the -1 subsite. This lack of coordination imparted mobility to the hydroxymethyl group and enabled a crucial hydrogen bond with the catalytic acid during and after the reaction. This suggests the possibility of a nonhydrolytic reaction mechanism when the catalytic base aspartic acid is missing, as is the case in some LTs (murein transglycosylase A) and expansins. We calculated reaction free energies and demonstrate the thermodynamic feasibility of the hydrolytic and nonhydrolytic reaction mechanisms. Our results provide molecular insights into the hydrolysis mechanism in HiCel45A, with possible implications for elucidating the elusive catalytic mechanism in LTs and expansins.

Advantages of a distant cellulase catalytic base.

The inverting glycoside hydrolase Trichoderma reesei (Hypocrea jecorina) Cel6A is a promising candidate for protein engineering for more economical production of biofuels. Until recently, its catalytic mechanism had been uncertain: The best candidate residue to serve as a catalytic base, Asp-175, is farther from the glycosidic cleavage site than in other glycoside hydrolase enzymes. Recent unbiased transition path sampling simulations revealed the hydrolytic mechanism for this more distant base, employing a water wire; however, it is not clear why the enzyme employs a more distant catalytic base, a highly conserved feature among homologs across different kingdoms. In this work, we describe molecular dynamics simulations designed to uncover how a base with a longer side chain, as in a D175E mutant, affects procession and active site alignment in the Michaelis complex. We show that the hydrogen bond network is tuned to the shorter aspartate side chain, and that a longer glutamate side chain inhibits procession as well as being less likely to adopt a catalytically productive conformation. Furthermore, we draw comparisons between the active site in Trichoderma reesei Cel6A and another inverting, processive cellulase to deduce the contribution of the water wire to the overall enzyme function, revealing that the more distant catalytic base enhances product release. Our results can inform efforts in the study and design of enzymes by demonstrating how counterintuitive sacrifices in chemical reactivity can have worthwhile benefits for other steps in the catalytic cycle.

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

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

Enantioselective Binding of Propranolol and Analogues Thereof to Cellobiohydrolase Cel7A.

At the catalytic site for the hydrolysis of cellulose the enzyme cellobiohydrolase Cel7A binds the enantiomers of the adrenergic beta-blocker propranolol with different selectivity. Methyl-to-hydroxymethyl group modifications of propranolol, which result in higher affinity and improved selectivity, were herein studied by 1 H,1 H and 1 H,13 C scalar spin-spin coupling constants as well as utilizing the nuclear Overhauser effect (NOE) in conjunction with molecular dynamics simulations of the ligands per se, which showed the presence of all-antiperiplanar conformations, except for the one containing a vicinal oxygen-oxygen arrangement governed by the gauche effect. For the ligand-protein complexes investigated by NMR spectroscopy using, inter alia, transferred NOESY and saturation-transfer difference (STD) NMR experiments the S-isomers were shown to bind with a higher affinity and a conformation similar to that preferred in solution, in contrast to the R-isomer. The fact that the S-form of the propranolol enantiomer is pre-arranged for binding to the protein is also observed for a crystal structure of dihydroxy-(S)-propranolol and Cel7A presented herein. Whereas the binding of propranolol is entropy driven, the complexation with the dihydroxy analogue is anticipated to be favored also by an enthalpic term, such as for its enantiomer, that is, dihydroxy-(R)-propranolol, because hydrogen-bond donation replaces the corresponding bonding from hydroxyl groups in glucosyl residues of the natural substrate. In addition to a favorable entropy component, albeit lesser in magnitude, this represents an effect of enthalpy-to-entropy compensation in ligand-protein interactions.

Quantum mechanical calculations suggest that lytic polysaccharide monooxygenases use a copper-oxyl, oxygen-rebound mechanism.

Lytic polysaccharide monooxygenases (LPMOs) exhibit a mononuclear copper-containing active site and use dioxygen and a reducing agent to oxidatively cleave glycosidic linkages in polysaccharides. LPMOs represent a unique paradigm in carbohydrate turnover and exhibit synergy with hydrolytic enzymes in biomass depolymerization. To date, several features of copper binding to LPMOs have been elucidated, but the identity of the reactive oxygen species and the key steps in the oxidative mechanism have not been elucidated. Here, density functional theory calculations are used with an enzyme active site model to identify the reactive oxygen species and compare two hypothesized reaction pathways in LPMOs for hydrogen abstraction and polysaccharide hydroxylation; namely, a mechanism that employs a η(1)-superoxo intermediate, which abstracts a substrate hydrogen and a hydroperoxo species is responsible for substrate hydroxylation, and a mechanism wherein a copper-oxyl radical abstracts a hydrogen and subsequently hydroxylates the substrate via an oxygen-rebound mechanism. The results predict that oxygen binds end-on (η(1)) to copper, and that a copper-oxyl-mediated, oxygen-rebound mechanism is energetically preferred. The N-terminal histidine methylation is also examined, which is thought to modify the structure and reactivity of the enzyme. Density functional theory calculations suggest that this posttranslational modification has only a minor effect on the LPMO active site structure or reactivity for the examined steps. Overall, this study suggests the steps in the LPMO mechanism for oxidative cleavage of glycosidic bonds.

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

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

Glycosylated linkers in multimodular lignocellulose-degrading enzymes dynamically bind to cellulose.

Plant cell-wall polysaccharides represent a vast source of food in nature. To depolymerize polysaccharides to soluble sugars, many organisms use multifunctional enzyme mixtures consisting of glycoside hydrolases, lytic polysaccharide mono-oxygenases, polysaccharide lyases, and carbohydrate esterases, as well as accessory, redox-active enzymes for lignin depolymerization. Many of these enzymes that degrade lignocellulose are multimodular with carbohydrate-binding modules (CBMs) and catalytic domains connected by flexible, glycosylated linkers. These linkers have long been thought to simply serve as a tether between structured domains or to act in an inchworm-like fashion during catalytic action. To examine linker function, we performed molecular dynamics (MD) simulations of the Trichoderma reesei Family 6 and Family 7 cellobiohydrolases (TrCel6A and TrCel7A, respectively) bound to cellulose. During these simulations, the glycosylated linkers bind directly to cellulose, suggesting a previously unknown role in enzyme action. The prediction from the MD simulations was examined experimentally by measuring the binding affinity of the Cel7A CBM and the natively glycosylated Cel7A CBM-linker. On crystalline cellulose, the glycosylated linker enhances the binding affinity over the CBM alone by an order of magnitude. The MD simulations before and after binding of the linker also suggest that the bound linker may affect enzyme action due to significant damping in the enzyme fluctuations. Together, these results suggest that glycosylated linkers in carbohydrate-active enzymes, which are intrinsically disordered proteins in solution, aid in dynamic binding during the enzymatic deconstruction of plant cell walls.

Structural and electronic snapshots during the transition from a Cu(II) to Cu(I) metal center of a lytic polysaccharide monooxygenase by X-ray photoreduction.

Lytic polysaccharide monooxygenases (LPMOs) are a recently discovered class of enzymes that employ a copper-mediated, oxidative mechanism to cleave glycosidic bonds. The LPMO catalytic mechanism likely requires that molecular oxygen first binds to Cu(I), but the oxidation state in many reported LPMO structures is ambiguous, and the changes in the LPMO active site required to accommodate both oxidation states of copper have not been fully elucidated. Here, a diffraction data collection strategy minimizing the deposited x-ray dose was used to solve the crystal structure of a chitin-specific LPMO from Enterococcus faecalis (EfaCBM33A) in the Cu(II)-bound form. Subsequently, the crystalline protein was photoreduced in the x-ray beam, which revealed structural changes associated with the conversion from the initial Cu(II)-oxidized form with two coordinated water molecules, which adopts a trigonal bipyramidal geometry, to a reduced Cu(I) form in a T-shaped geometry with no coordinated water molecules. A comprehensive survey of Cu(II) and Cu(I) structures in the Cambridge Structural Database unambiguously shows that the geometries observed in the least and most reduced structures reflect binding of Cu(II) and Cu(I), respectively. Quantum mechanical calculations of the oxidized and reduced active sites reveal little change in the electronic structure of the active site measured by the active site partial charges. Together with a previous theoretical investigation of a fungal LPMO, this suggests significant functional plasticity in LPMO active sites. Overall, this study provides molecular snapshots along the reduction process to activate the LPMO catalytic machinery and provides a general method for solving LPMO structures in both copper oxidation states.

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.

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.

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

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

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

Glycoside hydrolase family 7 (GH7) cellobiohydrolases (CBHs) play a key role in biomass recycling in nature. They are typically the most abundant enzymes expressed by potent cellulolytic fungi, and are also responsible for the majority of hydrolytic potential in enzyme cocktails for industrial processing of plant biomass. The thermostability of the enzyme is an important parameter for industrial utilization. In this study, Cel7 enzymes from different fungi were expressed in a fungal host and assayed for thermostability, including Hypocrea jecorina Cel7A as a reference. The most stable of the homologues, Humicola grisea var. thermoidea Cel7A, exhibits a 10°C higher melting temperature (T(m) of 72.5°C) and showed a 4-5 times higher initial hydrolysis rate than H. jecorina Cel7A on phosphoric acid-swollen cellulose and showed the best performance of the tested enzymes on pretreated corn stover at elevated temperature (65°C, 24 h). The enzyme shares 57% sequence identity with H. jecorina Cel7A and consists of a GH7 catalytic module connected by a linker to a C-terminal CBM1 carbohydrate-binding module. The crystal structure of the H. grisea var. thermoidea Cel7A catalytic module (1.8 Šresolution; R(work) and R(free) of 0.16 and 0.21, respectively) is similar to those of other GH7 CBHs. The deviations of several loops along the cellulose-binding path between the two molecules in the asymmetric unit indicate higher flexibility than in the less thermostable H. jecorina Cel7A.

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

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

The mechanism of cellulose hydrolysis by a two-step, retaining cellobiohydrolase elucidated by structural and transition path sampling studies.

Glycoside hydrolases (GHs) cleave glycosidic linkages in carbohydrates, typically via inverting or retaining mechanisms, the latter of which proceeds via a two-step mechanism that includes formation of a glycosyl-enzyme intermediate. We present two new structures of the catalytic domain of Hypocrea jecorina GH Family 7 cellobiohydrolase Cel7A, namely a Michaelis complex with a full cellononaose ligand and a glycosyl-enzyme intermediate, that reveal details of the 'static' reaction coordinate. We also employ transition path sampling to determine the 'dynamic' reaction coordinate for the catalytic cycle. The glycosylation reaction coordinate contains components of forming and breaking bonds and a conformational change in the nucleophile. Deglycosylation proceeds via a product-assisted mechanism wherein the glycosylation product, cellobiose, positions a water molecule for nucleophilic attack on the anomeric carbon of the glycosyl-enzyme intermediate. In concert with previous structures, the present results reveal the complete hydrolytic reaction coordinate for this naturally and industrially important enzyme family.