The rotamer of the second-sphere histidine in AA9 lytic polysaccharide monooxygenase is pH dependent.

Lytic polysaccharide monooxygenases (LPMOs) catalyze a reaction that is crucial for the biological decomposition of various biopolymers and for the industrial conversion of plant biomass. Despite the importance of LPMOs, the exact molecular-level nature of the reaction mechanism is still debated today. Here, we investigated the pH-dependent conformation of a second-sphere histidine (His) that we call the stacking histidine, which is conserved in fungal AA9 LPMOs and is speculated to assist catalysis in several of the LPMO reaction pathways. Using constant-pH and accelerated molecular dynamics simulations, we monitored the dynamics of the stacking His in different protonation states for both the resting Cu(II) and active Cu(I) forms of two fungal LPMOs. Consistent with experimental crystallographic and neutron diffraction data, our calculations suggest that the side chain of the protonated and positively charged form is rotated out of the active site toward the solvent. Importantly, only one of the possible neutral states of histidine (HIE state) is observed in the stacking orientation at neutral pH or when bound to cellulose. Our data predict that, in solution, the stacking His may act as a stabilizer (via hydrogen bonding) of the Cu(II)-superoxo complex after the LPMO-Cu(I) has reacted with O2 in solution, which, in fine, leads to H2O2 formation. Also, our data indicate that the HIE-stacking His is a poor acid/base catalyst when bound to the substrate and, in agreement with the literature, may play an important stabilizing role (via hydrogen bonding) during the peroxygenase catalysis. Our study reveals the pH titration midpoint values of the pH-dependent orientation of the stacking His should be considered when modeling and interpreting LPMO reactions, whether it be for classical LPMO kinetics or in industry-oriented enzymatic cocktails, and for understanding LPMO behavior in slightly acidic natural processes such as fungal wood decay.

Characterization and engineering of a two-enzyme system for plastics depolymerization.

Plastics pollution represents a global environmental crisis. In response, microbes are evolving the capacity to utilize synthetic polymers as carbon and energy sources. Recently, Ideonella sakaiensis was reported to secrete a two-enzyme system to deconstruct polyethylene terephthalate (PET) to its constituent monomers. Specifically, the I. sakaiensis PETase depolymerizes PET, liberating soluble products, including mono(2-hydroxyethyl) terephthalate (MHET), which is cleaved to terephthalic acid and ethylene glycol by MHETase. Here, we report a 1.6 Å resolution MHETase structure, illustrating that the MHETase core domain is similar to PETase, capped by a lid domain. Simulations of the catalytic itinerary predict that MHETase follows the canonical two-step serine hydrolase mechanism. Bioinformatics analysis suggests that MHETase evolved from ferulic acid esterases, and two homologous enzymes are shown to exhibit MHET turnover. Analysis of the two homologous enzymes and the MHETase S131G mutant demonstrates the importance of this residue for accommodation of MHET in the active site. We also demonstrate that the MHETase lid is crucial for hydrolysis of MHET and, furthermore, that MHETase does not turnover mono(2-hydroxyethyl)-furanoate or mono(2-hydroxyethyl)-isophthalate. A highly synergistic relationship between PETase and MHETase was observed for the conversion of amorphous PET film to monomers across all nonzero MHETase concentrations tested. Finally, we compare the performance of MHETase:PETase chimeric proteins of varying linker lengths, which all exhibit improved PET and MHET turnover relative to the free enzymes. Together, these results offer insights into the two-enzyme PET depolymerization system and will inform future efforts in the biological deconstruction and upcycling of mixed plastics.

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.

Enabling microbial syringol conversion through structure-guided protein engineering.

Microbial conversion of aromatic compounds is an emerging and promising strategy for valorization of the plant biopolymer lignin. A critical and often rate-limiting reaction in aromatic catabolism is O-aryl-demethylation of the abundant aromatic methoxy groups in lignin to form diols, which enables subsequent oxidative aromatic ring-opening. Recently, a cytochrome P450 system, GcoAB, was discovered to demethylate guaiacol (2-methoxyphenol), which can be produced from coniferyl alcohol-derived lignin, to form catechol. However, native GcoAB has minimal ability to demethylate syringol (2,6-dimethoxyphenol), the analogous compound that can be produced from sinapyl alcohol-derived lignin. Despite the abundance of sinapyl alcohol-based lignin in plants, no pathway for syringol catabolism has been reported to date. Here we used structure-guided protein engineering to enable microbial syringol utilization with GcoAB. Specifically, a phenylalanine residue (GcoA-F169) interferes with the binding of syringol in the active site, and on mutation to smaller amino acids, efficient syringol O-demethylation is achieved. Crystallography indicates that syringol adopts a productive binding pose in the variant, which molecular dynamics simulations trace to the elimination of steric clash between the highly flexible side chain of GcoA-F169 and the additional methoxy group of syringol. Finally, we demonstrate in vivo syringol turnover in Pseudomonas putida KT2440 with the GcoA-F169A variant. Taken together, our findings highlight the significant potential and plasticity of cytochrome P450 aromatic O-demethylases in the biological conversion of lignin-derived aromatic compounds.

Kinetic and molecular dynamics study of inhibition and transglycosylation in Hypocrea jecorina family 3 ß-glucosidases.

ß-Glucosidases enhance enzymatic biomass conversion by relieving cellobiose inhibition of endoglucanases and cellobiohydrolases. However, the susceptibility of these enzymes to inhibition and transglycosylation at high glucose or cellobiose concentrations severely limits their activity and, consequently, the overall efficiency of enzyme mixtures. We determined the impact of these two processes on the hydrolytic activity of the industrially relevant family 3 ß-glucosidases from Hypocrea jecorina, HjCel3A and HjCel3B, and investigated the underlying molecular mechanisms through kinetic studies, binding free energy calculations, and molecular dynamics (MD) simulations. HjCel3B had a 7-fold higher specificity for cellobiose than HjCel3A but greater tendency for glucose inhibition. Energy decomposition analysis indicated that cellobiose has relatively weak electrostatic interactions with binding site residues, allowing it to be easily displaced by glucose and free to inhibit other hydrolytic enzymes. HjCel3A is, thus, preferable as an industrial ß-glucosidase despite its lower activity caused by transglycosylation. This competing pathway to hydrolysis arises from binding of glucose or cellobiose at the product site after formation of the glycosyl-enzyme intermediate. MD simulations revealed that binding is facilitated by hydrophobic interactions with Trp-37, Phe-260, and Tyr-443. Targeting these aromatic residues for mutation to reduce substrate affinity at the product site would therefore potentially mitigate transglycosidic activity. Engineering improved variants of HjCel3A and other structurally similar ß-glucosidases would have a significant economic effect on enzymatic biomass conversion in terms of yield and production cost as the process can be consequently conducted at higher substrate loadings.

Improving Enzyme Optimum Temperature Prediction with Resampling Strategies and Ensemble Learning.

Accurate prediction of the optimal catalytic temperature (Topt) of enzymes is vital in biotechnology, as enzymes with high Topt values are desired for enhanced reaction rates. Recently, a machine learning method (temperature optima for microorganisms and enzymes, TOME) for predicting Topt was developed. TOME was trained on a normally distributed data set with a median Topt of 37 °C and less than 5% of Topt values above 85 °C, limiting the method's predictive capabilities for thermostable enzymes. Due to the distribution of the training data, the mean squared error on Topt values greater than 85 °C is nearly an order of magnitude higher than the error on values between 30 and 50 °C. In this study, we apply ensemble learning and resampling strategies that tackle the data imbalance to significantly decrease the error on high Topt values (>85 °C) by 60% and increase the overall R2 value from 0.527 to 0.632. The revised method, temperature optima for enzymes with resampling (TOMER), and the resampling strategies applied in this work are freely available to other researchers as Python packages on GitHub.

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)6 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)5 and (Glc)6 was performed in H218O enabling the incorporation of an 18O 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 18O- and 16O-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.

Evaluation of the Immunomodulatory Effects of All-Trans Retinoic Acid Solid Lipid Nanoparticles and Human Mesenchymal Stem Cells in an A549 Epithelial Cell Line Model.

PURPOSE: To investigate two potential strategies aimed at targeting the inflammatory pathogenesis of COPD: a small molecule, all trans retinoic acid (atRA) and human mesenchymal stem cells (hMSCs). METHODS: atRA was formulated into solid lipid nanoparticles (SLNs) via the emulsification-ultrasonication method, and these SLNs were characterised physicochemically. Assessment of the immunomodulatory effects of atRA-SLNs on A549 cells in vitro was determined using ELISA. hMSCs were suspended in a previously developed methylcellulose, collagen and beta-glycerophosphate hydrogel prior to investigating their immunomodulatory effects in vitro. RESULTS: SLNs provided significant encapsulation of atRA and also sustained its release over 72 h. A549 cells were viable following the addition of atRA SLNs and showed a reduction in IL-6 and IL-8 levels. A549 cells also remained viable following addition of the hMSC/hydrogel formulation - however, this formulation resulted in increased levels of IL-6 and IL-8, indicating a potentially pro-inflammatory effect. CONCLUSION: Both atRA SLNs and hMSCs show potential for modulating the environment in inflammatory disease, though through different mechanisms and leading to different outcomes - despite both being explored as strategies for use in inflammatory disease. atRA shows promise by acting in a directly anti-inflammatory manner, whereas further research into the exact mechanisms and behaviours of hMSCs in inflammatory diseases is required.

Molecular Determinants of Substrate Affinity and Enzyme Activity of a Cytochrome P450BM3 Variant.

Cytochrome P450BM3 catalyzes the hydroxylation and/or epoxidation of fatty acids, fatty amides, and alcohols. Protein engineering has produced P450BM3 variants capable of accepting drug molecules normally metabolized by human P450 enzymes. The enhanced substrate promiscuity has been attributed to the greater flexibility of the lid of the substrate channel. However, it is not well understood how structurally different and highly polar drug molecules can stably bind in the active site nor how the activity and coupling efficiency of the enzyme may be affected by the lack of enzyme-substrate complementarity. To address these important aspects of non-native small molecule binding, this study investigated the binding of drug molecules with different size, charge, polar surface area, and human P450 affinity on the promiscuous R47L/F87V/L188Q/E267V/F81I pentuple mutant of P450BM3. Binding free energy data and energy decomposition analysis showed that pentuple mutant P450BM3 stably binds (i.e., negative ΔGb°) a broad range of substrate and inhibitor types because dispersion interactions with active site residues overcome unfavorable repulsive and electrostatic effects. Molecular dynamics simulations revealed that 1) acidic substrates tend to disrupt the heme propionate A-K69 salt bridge, which may reduce heme oxidizing ability, and 2) the lack of complementarity leads to high substrate mobility and water density in the active site, which may lead to uncoupling. These factors must be considered in future developments of P450BM3 as a biocatalyst in the large-scale production of drug metabolites.

Inhibition of Mammalian Glycoprotein YKL-40: IDENTIFICATION OF THE PHYSIOLOGICAL LIGAND.

YKL-40 is a mammalian glycoprotein associated with progression, severity, and prognosis of chronic inflammatory diseases and a multitude of cancers. Despite this well documented association, identification of the lectin's physiological ligand and, accordingly, biological function has proven experimentally difficult. YKL-40 has been shown to bind chito-oligosaccharides; however, the production of chitin by the human body has not yet been documented. Possible alternative ligands include proteoglycans, polysaccharides, and fibers like collagen, all of which makeup the extracellular matrix. It is likely that YKL-40 is interacting with these alternative polysaccharides or proteins within the body, extending its function to cell biological roles such as mediating cellular receptors and cell adhesion and migration. Here, we consider the feasibility of polysaccharides, including cello-oligosaccharides, hyaluronan, heparan sulfate, heparin, and chondroitin sulfate, and collagen-like peptides as physiological ligands for YKL-40. We use molecular dynamics simulations to resolve the molecular level recognition mechanisms and calculate the free energy of binding the hypothesized ligands to YKL-40, addressing thermodynamic preference relative to chito-oligosaccharides. Our results suggest that chitohexaose and hyaluronan preferentially bind to YKL-40 over collagen, and hyaluronan is likely the preferred physiological ligand, because the negatively charged hyaluronan shows enhanced affinity for YKL-40 over neutral chitohexaose. Collagen binds in two locations at the YKL-40 surface, potentially related to a role in fibrillar formation. Finally, heparin non-specifically binds at the YKL-40 surface, as predicted from structural studies. Overall, YKL-40 likely binds many natural ligands in vivo, but its concurrence with physical maladies may be related to associated increases in hyaluronan.

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.

The role of catalytic residue pKa on the hydrolysis/transglycosylation partition in family 3 ß-glucosidases.

ß-Glucosidases (ßgls) primarily catalyze the hydrolysis of the terminal glycosidic bond at the non-reducing end of ß-glucosides, although glycosidic bond synthesis (called transglycosylation) can also occur in the presence of another acceptor. In the final reaction step, the glucose product or another substrate competes with water for transfer to the glycosyl-enzyme intermediate. The factors governing the balance between the two pathways are not fully known; however, the involvement of ionizable residues in binding and catalysis suggests that their pKa may play a role. Through constant pH molecular dynamics simulations of a glycoside hydrolase Family 3 (GH3) ßgl, we showed that the pKa of the catalytic acid/base residue, E441, is low (∼2) during either reaction due to E441-R125-E128 and E441-R125-E166 hydrogen bond networks. The low basicity of E441 would reduce its ability to deprotonate the acceptor. This may be less critical for transglycosylation because sugars have a lower deprotonation enthalpy than water. Moreover, their acidity would be increased by hydrogen bonding with R169 at the acceptor binding site. In contrast, no such interaction was observed for catalytic water. The results are likely applicable to other GH3 ßgls because R125, E128, E166, and R169 are conserved. As these enzymes are commonly used in biomass degradation, there is interest in developing variants with enhanced hydrolytic activity. This may be accomplished by elevating the acid/base residue pKa by disrupting its hydrogen bond networks and reducing the affinity and reactivity of a sugar acceptor by mutating R169.

Effect of Mutation and Substrate Binding on the Stability of Cytochrome P450BM3 Variants.

Cytochrome P450BM3 is a heme-containing enzyme from Bacillus megaterium that exhibits high monooxygenase activity and has a self-sufficient electron transfer system in the full-length enzyme. Its potential synthetic applications drive protein engineering efforts to produce variants capable of oxidizing nonnative substrates such as pharmaceuticals and aromatic pollutants. However, promiscuous P450BM3 mutants often exhibit lower stability, thereby hindering their industrial application. This study demonstrated that the heme domain R47L/F87V/L188Q/E267V/F81I pentuple mutant (PM) is destabilized because of the disruption of hydrophobic contacts and salt bridge interactions. This was directly observed from crystal structures of PM in the presence and absence of ligands (palmitic acid and metyrapone). The instability of the tertiary structure and heme environment of substrate-free PM was confirmed by pulse proteolysis and circular dichroism, respectively. Binding of the inhibitor, metyrapone, significantly stabilized PM, but the presence of the native substrate, palmitic acid, had no effect. On the basis of high-temperature molecular dynamics simulations, the lid domain, ß-sheet 1, and Cys ligand loop (a ß-bulge segment connected to the heme) are the most labile regions and, thus, potential sites for stabilizing mutations. Possible approaches to stabilization include improvement of hydrophobic packing interactions in the lid domain and introduction of new salt bridges into ß-sheet 1 and the heme region. An understanding of the molecular factors behind the loss of stability of P450BM3 variants therefore expedites site-directed mutagenesis studies aimed at developing thermostability.

Structural and Functional Characterization of a Lytic Polysaccharide Monooxygenase with Broad Substrate Specificity.

The recently discovered lytic polysaccharide monooxygenases (LPMOs) carry out oxidative cleavage of polysaccharides and are of major importance for efficient processing of biomass. NcLPMO9C from Neurospora crassa acts both on cellulose and on non-cellulose ß-glucans, including cellodextrins and xyloglucan. The crystal structure of the catalytic domain of NcLPMO9C revealed an extended, highly polar substrate-binding surface well suited to interact with a variety of sugar substrates. The ability of NcLPMO9C to act on soluble substrates was exploited to study enzyme-substrate interactions. EPR studies demonstrated that the Cu(2+) center environment is altered upon substrate binding, whereas isothermal titration calorimetry studies revealed binding affinities in the low micromolar range for polymeric substrates that are due in part to the presence of a carbohydrate-binding module (CBM1). Importantly, the novel structure of NcLPMO9C enabled a comparative study, revealing that the oxidative regioselectivity of LPMO9s (C1, C4, or both) correlates with distinct structural features of the copper coordination sphere. In strictly C1-oxidizing LPMO9s, access to the solvent-facing axial coordination position is restricted by a conserved tyrosine residue, whereas access to this same position seems unrestricted in C4-oxidizing LPMO9s. LPMO9s known to produce a mixture of C1- and C4-oxidized products show an intermediate situation.

Structural characterization of a unique marine animal family 7 cellobiohydrolase suggests a mechanism of cellulase salt tolerance.

Nature uses a diversity of glycoside hydrolase (GH) enzymes to convert polysaccharides to sugars. As lignocellulosic biomass deconstruction for biofuel production remains costly, natural GH diversity offers a starting point for developing industrial enzymes, and fungal GH family 7 (GH7) cellobiohydrolases, in particular, provide significant hydrolytic potential in industrial mixtures. Recently, GH7 enzymes have been found in other kingdoms of life besides fungi, including in animals and protists. Here, we describe the in vivo spatial expression distribution, properties, and structure of a unique endogenous GH7 cellulase from an animal, the marine wood borer Limnoria quadripunctata (LqCel7B). RT-quantitative PCR and Western blot studies show that LqCel7B is expressed in the hepatopancreas and secreted into the gut for wood degradation. We produced recombinant LqCel7B, with which we demonstrate that LqCel7B is a cellobiohydrolase and obtained four high-resolution crystal structures. Based on a crystallographic and computational comparison of LqCel7B to the well-characterized Hypocrea jecorina GH7 cellobiohydrolase, LqCel7B exhibits an extended substrate-binding motif at the tunnel entrance, which may aid in substrate acquisition and processivity. Interestingly, LqCel7B exhibits striking surface charges relative to fungal GH7 enzymes, which likely results from evolution in marine environments. We demonstrate that LqCel7B stability and activity remain unchanged, or increase at high salt concentration, and that the L. quadripunctata GH mixture generally contains cellulolytic enzymes with highly acidic surface charge compared with enzymes derived from terrestrial microbes. Overall, this study suggests that marine cellulases offer significant potential for utilization in high-solids industrial biomass conversion processes.

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.

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.

Cello-oligomer-binding dynamics and directionality in family 4 carbohydrate-binding modules.

Carbohydrate-binding modules (CBMs) play significant roles in modulating the function of cellulases, and understanding the protein-carbohydrate recognition mechanisms by which CBMs selectively bind substrate is critical to development of enhanced biomass conversion technology. CBMs exhibit a limited range of specificity and appear to bind polysaccharides in a directional fashion dictated by the position of the ring oxygen relative to the protein fold. The two family 4 CBMs of Cellulomonas fimi Cel9B (CfCBM4) are reported to preferentially bind cellulosic substrates. However, experimental evidence suggests that these CBMs may not exhibit a thermodynamic preference for a particular orientation. We use molecular dynamics (MD) and free energy calculations to investigate protein-carbohydrate recognition mechanisms in CfCBM4-1 and CfCBM4-2 and to elucidate preferential ligand-binding orientation. We evaluate four cellopentaose orientations including that of the crystal structure and three others suggested by nuclear magnetic resonance (NMR). These four orientations differ based on position of the ligand reducing end (RE) and pyranose ring orientations relative to the protein core. MD simulations indicate that the plausible orientations reduce to two conformations. Calculated ligand-binding free energy discerns each of the orientations is equally favorable. The calculated free energies are in excellent agreement with isothermal titration calorimetry measurements from the literature. MD simulations further reveal the approximate structural symmetry of the oligosaccharides relative to the amino acids along the binding cleft plays a role in the promiscuity of ligand binding. A survey of ligand-bound structures suggests this phenomenon may be characteristic of the broader class of proteins belonging to the ß-sandwich fold.

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.

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.