Small angle neutron scattering reveals pH-dependent conformational changes in Trichoderma reesei cellobiohydrolase I: implications for enzymatic activity

Light-stimulated T. thermophilus two-domain LPMO9H: Low-resolution SAXS model and synergy with cellulases

Lytic polysaccharide monooxygenases (LPMOs), monocopper enzymes that oxidatively cleave recalcitrant polysaccharides, have important biotechnological applications. Thermothelomyces thermophilus is a rich source of biomass-active enzymes, including many members from auxiliary activities family 9 LPMOs. Here, we report biochemical and structural characterization of recombinant TtLPMO9H which oxidizes cellulose at the C1 and C4 positions and shows enhanced activity in light-driven catalysis assays. TtLPMO9H also shows activity against xyloglucan. The addition of TtLPMO9H to endoglucanases from four different glucoside hydrolase families (GH5, GH12, GH45 and GH7) revealed that the product formation was remarkably increased when TtLPMO9H was combined with GH7 endoglucanase. Finally, we determind the first low resolution small-angle X-ray scattering model of the two-domain TtLPMO9H in solution that shows relative positions of its two functional domains and a conformation of the linker peptide, which can be relevant for the catalytic oxidation of cellulose and xyloglucan.

Glycosylation effects on the structure and dynamics of a full-length Cel7A cellulase

Fungi cellulases are used to degrade cellulose-containing biomass for bioethanol production. Industrial cellulases such as Cel7A from Trichoderma reesei (TrCel7A) are critical in this process. Thus, the understanding of structure and dynamics is crucial for engineering variants with improved cellulolytic activity. This cellulase consists of two domains connected by a flexible and highly glycosylated linker. However, the linker flexibility has hindered the determination of Cel7A complete structure. Herein, based on atomic and sparse data, we applied integrative modelling to build a model of the complete enzyme structure. Next, through simulations, we studied the glycosylation effects on the structure and dynamics of a solubilized TrCel7A. Essential dynamics analysis showed that O-glycosylation in the linker led to the stabilization of protein overall dynamics. O-linked glycans seem to restrict protein dihedral angles distribution in this region, selecting more elongated conformations. Besides the reduced flexibility, functional interdomain motions occurred in a more concerted way in the glycosylated system. In contrast, in the absence of glycosylation, we observed vast conformational plasticity with the functional domains frequently collapsing. We report here evidence that targeting Cel7A linker flexibility by point mutations including modification of glycosylation sites could be a promising design strategy to improve cellulase activity.

High-resolution crystal structures of the glycoside hydrolase family 45 endoglucanase EG27II from the snail Ampullaria crossean

Although endogenous animal cellulases have great potential for industrial applications such as bioethanol production, few investigations have focused on these enzymes. In this study, the glycoside hydrolase family 45 (GH45) subfamily B endoglucanase EG27II from the snail Ampullaria crossean was expressed using a Pichia pastoris expression system and the crystal structure of the apo form was determined at 1.00 Å resolution; this is the highest resolution structure of an animal endoglucanase. The results showed that EG27II has a double-ψ six-stranded β-barrel and that the structure of EG27II more closely resembles those of subfamily C enzymes than those of subfamily A enzymes. The structure of EG27II complexed with cellobiose was also determined under cryoconditions and at room temperature at three pH values, pH 4.0, 5.5 and 8.0, and no structural changes were found to be associated with the change in pH. The structural comparison and catalytic activity measurements showed that Asp137 and Asn112 function as the catalytic acid and base, respectively, and that Asp27 is also an important residue for catalysis. These high-resolution structures of EG27II provide a large amount of information for structure-based enzyme modification and cell-surface engineering, which will advance biofuel production using animal-derived cellulases.

Gradual Crossover from Subdiffusion to Normal Diffusion: A Many-Body Effect in Protein Surface Water

Dynamics of hydration water is essential for the function of biomacromolecules. Previous studies have demonstrated that water molecules exhibit subdiffusion on the surface of biomacromolecules; yet the microscopic mechanism remains vague. Here, by performing neutron scattering, molecular dynamics simulations, and analytic modeling on hydrated perdeuterated protein powders, we found water molecules jump randomly between trapping sites on protein surfaces, whose waiting times obey a broad distribution, resulting in subdiffusion. Moreover, the subdiffusive exponent gradually increases with observation time towards normal diffusion due to a many-body volume-exclusion effect.

Concerted motions and large-scale structural fluctuations of Trichoderma reesei Cel7A cellobiohydrolase

Cellobiohydrolases (CBHs) are key enzymes for the saccharification of cellulose and play major roles in industrial settings for biofuel production. The catalytic core domain of these enzymes exhibits a long and narrow binding tunnel capable of binding glucan chains from crystalline cellulose and processively hydrolyze them. The binding cleft is topped by a set of loops, which are believed to play key roles in substrate binding and cleavage processivity. Here, we present an analysis of the loop motions of the Trichoderma reesei Cel7A catalytic core domain (TrCel7A) using conventional and accelerated molecular dynamics simulations. We observe that the loops exhibit highly coupled fluctuations and cannot move independently of each other. In the absence of a substrate, the characteristic large amplitude dynamics of TrCel7A consists of breathing motions, where the loops undergo open-and-close fluctuations. Upon substrate binding, the open-close fluctuations of the loops are quenched and one of the loops moves parallel to the binding site, possibly to allow processive motion along the glucan chain. Using microsecond accelerated molecular dynamics, we observe large-scale fluctuations of the loops (up to 37 Å) and the entire exposure of the TrCel7A binding site in the absence of the substrate, resembling an endoglucanase. These results suggest that the initial CBH-substrate contact and substrate recognition by the enzyme are similar to that of endoglucanases and, once bound to the substrate, the loops remain closed for proper enzymatic activity.

Modification of the nanostructure of lignocellulose cell walls via a non-enzymatic lignocellulose deconstruction system in brown rot wood-decay fungi

Wood decayed by brown rot fungi and wood treated with the chelator-mediated Fenton (CMF) reaction, either alone or together with a cellulose enzyme cocktail, was analyzed by small angle neutron scattering (SANS), sum frequency generation (SFG) spectroscopy, Fourier transform infrared (FTIR) analysis, X-ray diffraction (XRD), atomic force microscopy (AFM), and transmission electron microscopy (TEM). Results showed that the CMF mechanism mimicked brown rot fungal attack for both holocellulose and lignin components of the wood. Crystalline cellulose and lignin were both depolymerized by the CMF reaction. Porosity of the softwood cell wall did not increase during CMF treatment, enzymes secreted by the fungi did not penetrate the decayed wood. The enzymes in the cellulose cocktail also did not appear to alter the effects of the CMF-treated wood relative to enhancing cell wall deconstruction. This suggests a rethinking of current brown rot decay models and supports a model where monomeric sugars and oligosaccharides diffuse from the softwood cell walls during non-enzymatic action. In this regard, the CMF mechanism should not be thought of as a "pretreatment" used to permit enzymatic penetration into softwood cell walls, but instead it enhances polysaccharide components diffusing to fungal enzymes located in wood cell lumen environments during decay. SANS and other data are consistent with a model for repolymerization and aggregation of at least some portion of the lignin within the cell wall, and this is supported by AFM and TEM data. The data suggest that new approaches for conversion of wood substrates to platform chemicals in biorefineries could be achieved using the CMF mechanism with >75% solubilization of lignocellulose, but that a more selective suite of enzymes and other downstream treatments may be required to work when using CMF deconstruction technology. Strategies to enhance polysaccharide release from lignocellulose substrates for enhanced enzymatic action and fermentation of the released fraction would also aid in the efficient recovery of the more uniform modified lignin fraction that the CMF reaction generates to enhance biorefinery profitability.

Deleting the Ig-Like Domain of Alicyclobacillus acidocaldarius Endoglucanase Cel9A Causes a Simultaneous Increase in the Activity and Stability

Endoglucanase Cel9A from Alicyclobacillus acidocaldarius (AaCel9A) is a monomeric enzyme with 537 residues. This enzyme has an Ig-like domain in the N-terminus of the catalytic domain. In this study, the role of the Ig-like domain on the activity, stability, and structural rigidity of AaCel9A and the effect of calcium on enzyme activity and stability were examined by comparing a truncated enzyme with deletion of the Ig-like domain (AaCel9AΔN) to the wild-type enzyme. Our results showed that the deletion of the Ig-like domain increased the catalytic efficiency of the truncated enzyme up to threefold without any significant changes in the K m of the enzyme. Furthermore, pH and temperature optimum for activity were shifted from 6.5 to 7.5 and from 65 to 60 °C, respectively, by deletion of the Ig-like domain. The thermal stability and fluorescence quenching results indicated that the stability and rigidity of the truncated enzyme have been more than that of the wild-type enzyme. Calcium similarly increased the catalytic efficiency of the enzymes (up to 40 %) and remarkably raised the stability of the AaCel9A compared to the AaCel9AΔN. This shows that Ig-like domain has a role in the increase of the enzyme stability by calcium in the wild-type enzyme.

Biophysical and biochemical studies of a major endoglucanase secreted by Xanthomonas campestris pv. campestris

Endoglucanases are the main cellulolytic enzymes secreted by the bacterium Xanthomonas campestris pv. campestris (Xcc). The major endoglucanase exported by this bacterium into an external milieu is an enzyme XccCel5A, which belongs to GH5 family subfamily 1 and is encoded by the gene engXCA. We purified XccCel5A using ammonium sulfate precipitation followed by size exclusion chromatography and identified it by zymogram analysis. Circular dichroism and fluorescence spectroscopy studies showed that XccCel5A is stable in a wide pH range and up to about 55°C and denatures at the higher temperatures. The optimal conditions for enzyme activity were identified as T=45°C and pH=7.0. Under the optimum conditions the catalytic efficiency (kcat/KM) of the enzyme was determined as 5.16×10(4)s(-1)M(-1) using carboxymethylcellulose (CMC) as a substrate. Our SAXS studies revealed extended tadpole-shape molecular assembly, typical for cellulases, and allowed to determine an overall shape of the enzyme and a relative position of the catalytic and cellulose binding domains.

Chemical stability of a cold-active cellulase with high tolerance toward surfactants and chaotropic agent

CelE1 is a cold-active endo-acting glucanase with high activity at a broad temperature range and under alkaline conditions. Here, we examined the effects of pH on the secondary and tertiary structures, net charge, and activity of CelE1. Although variation in pH showed a small effect in the enzyme structure, the activity was highly influenced at acidic conditions, while reached the optimum activity at pH 8. Furthermore, to estimate whether CelE1 could be used as detergent additives, CelE1 activity was evaluated in the presence of surfactants. Ionic and nonionic surfactants were not able to reduce CelE1 activity significantly. Therefore, CelE1 was found to be promising candidate for use as detergent additives. Finally, we reported a thermodynamic analysis based on the structural stability and the chemical unfolding/refolding process of CelE1. The results indicated that the chemical unfolding proceeds as a reversible two-state process. These data can be useful for biotechnological applications.

Production of bacterial cellulose with controlled deuterium-hydrogen substitution for neutron scattering studies

Isotopic enrichment of biomacromolecules is a widely used technique that enables the investigation of the structural and dynamic properties to provide information not accessible with natural abundance isotopic composition. This study reports an approach for deuterium incorporation into bacterial cellulose. A media formulation for growth of Acetobacter xylinus subsp. sucrofermentans and Gluconacetobacter hansenii was formulated that supports cellulose production in deuterium (D) oxide. The level of D incorporation can be varied by altering the ratio of deuterated and protiated glycerol used during cell growth in the D2O-based growth medium. Spectroscopic analysis and mass spectrometry show that the level of deuterium incorporation is high (>90%) for the perdeuterated form of bacterial cellulose. The small-angle neutron scattering profiles of the cellulose with different amounts of D incorporation are all similar indicating that there are no structural changes in the cellulose due to substitution of deuterium for hydrogen. In addition, by varying the amount of deuterated glycerol in the media it was possible to vary the scattering length density of the deuterated cellulose. The ability to control deuterium content of cellulose extends the range of experiments using techniques such as neutron scattering to reveal information about the structure and dynamics of cellulose, and its interactions with other biomacromolecules as well as synthetic polymers used for development of composite materials.

Conformational changes in a hyperthermostable glycoside hydrolase: enzymatic activity is a consequence of the loop dynamics and protonation balance

Endo-β-1, 4-mannanase from Thermotoga petrophila (TpMan) is a modular hyperthermostable enzyme involved in the degradation of mannan-containing polysaccharides. The degradation of these polysaccharides represents a key step for several industrial applications. Here, as part of a continuing investigation of TpMan, the region corresponding to the GH5 domain (TpManGH5) was characterized as a function of pH and temperature. The results indicated that the enzymatic activity of the TpManGH5 is pH-dependent, with its optimum activity occurring at pH 6. At pH 8, the studies demonstrated that TpManGH5 is a molecule with a nearly spherical tightly packed core displaying negligible flexibility in solution, and with size and shape very similar to crystal structure. However, TpManGH5 experiences an increase in radius of gyration in acidic conditions suggesting expansion of the molecule. Furthermore, at acidic pH values, TpManGH5 showed a less globular shape, probably due to a loop region slightly more expanded and flexible in solution (residues Y88 to A105). In addition, molecular dynamics simulations indicated that conformational changes caused by pH variation did not change the core of the TpManGH5, which means that only the above mentioned loop region presents high degree of fluctuations. The results also suggested that conformational changes of the loop region may facilitate polysaccharide and enzyme interaction. Finally, at pH 6 the results indicated that TpManGH5 is slightly more flexible at 65°C when compared to the same enzyme at 20°C. The biophysical characterization presented here is well correlated with the enzymatic activity and provide new insight into the structural basis for the temperature and pH-dependent activity of the TpManGH5. Also, the data suggest a loop region that provides a starting point for a rational design of biotechnological desired features.

Mechanism of lignin inhibition of enzymatic biomass deconstruction

BACKGROUND

The conversion of plant biomass to ethanol via enzymatic cellulose hydrolysis offers a potentially sustainable route to biofuel production. However, the inhibition of enzymatic activity in pretreated biomass by lignin severely limits the efficiency of this process.

RESULTS

By performing atomic-detail molecular dynamics simulation of a biomass model containing cellulose, lignin, and cellulases (TrCel7A), we elucidate detailed lignin inhibition mechanisms. We find that lignin binds preferentially both to the elements of cellulose to which the cellulases also preferentially bind (the hydrophobic faces) and also to the specific residues on the cellulose-binding module of the cellulase that are critical for cellulose binding of TrCel7A (Y466, Y492, and Y493).

CONCLUSIONS

Lignin thus binds exactly where for industrial purposes it is least desired, providing a simple explanation of why hydrolysis yields increase with lignin removal.

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.

Crystallization and preliminary X-ray diffraction analysis of Hypocrea jecorina Cel7A in two new crystal forms

Cel7A (previously known as cellobiohydrolase I) from Hypocrea jecorina was crystallized in two crystalline forms, neither of which have been previously reported. Both forms co-crystallize under the same crystallization conditions. The first crystal form belonged to space group C2, with unit-cell parameters a=152.5, b=44.9, c=57.6 Å, β=101.2°, and diffracted X-rays to 1.5 Å resolution. The second crystal form belonged to space group P6₃22, with unit-cell parameters a=b≃155, c≃138 Å, and diffracted X-rays to 2.5 Å resolution. The crystals were obtained using full-length Cel7A, which consists of a large 434-residue N-terminal catalytic domain capable of cleaving cellulose, a 27-residue flexible linker and a small 36-residue C-terminal carbohydrate-binding module (CBM). However, a preliminary analysis of the electron-density maps suggests that the linker and CBM are disordered in both crystal forms. Complete refinement and structure analysis are currently in progress.

Predicting enzyme adsorption to lignin films by calculating enzyme surface hydrophobicity

The inhibitory action of lignin on cellulase cocktails is a major challenge to the biological saccharification of plant cell wall polysaccharides. Although the mechanism remains unclear, hydrophobic interactions between enzymes and lignin are hypothesized to drive adsorption. Here we evaluate the role of hydrophobic interactions in enzyme-lignin binding. The hydrophobicity of the enzyme surface was quantified using an estimation of the clustering of nonpolar atoms, identifying potential interaction sites. The adsorption of enzymes to lignin surfaces, measured using the quartz crystal microbalance, correlates to the hydrophobic cluster scores. Further, these results suggest a minimum hydrophobic cluster size for a protein to preferentially adsorb to lignin. The impact of electrostatic contribution was ruled out by comparing the isoelectric point (pI) values to the adsorption of proteins to lignin surfaces. These results demonstrate the ability to predict enzyme-lignin adsorption and could potentially be used to design improved cellulase cocktails, thus lowering the overall cost of biofuel production.

Deciphering the effect of the different N-glycosylation sites on the secretion, activity, and stability of cellobiohydrolase I from Trichoderma reesei

N-linked glycosylation modulates and diversifies the structures and functions of the eukaryotic proteome through both intrinsic and extrinsic effects on proteins. We investigated the significance of the three N-linked glycans on the catalytic domain of cellobiohydrolase I (CBH1) from the filamentous fungus Trichoderma reesei in its secretion and activity. While the removal of one or two N-glycosylation sites hardly affected the extracellular secretion of CBH1, eliminating all of the glycosylation sites did induce expression of the unfolded protein response (UPR) target genes, and secretion of this CBH1 variant was severely compromised in a calnexin gene deletion strain. Further characterization of the purified CBH1 variants showed that, compared to Asn270, the thermal reactivity of CBH1 was significantly decreased by removal of either Asn45 or Asn384 glycosylation site during the catalyzed hydrolysis of soluble substrate. Combinatorial loss of these two N-linked glycans further exacerbated the temperature-dependent inactivation. In contrast, this thermal labile property was less severe when hydrolyzing insoluble cellulose. Analysis of the structural integrity of CBH1 variants revealed that removal of N-glycosylation at Asn384 had a more pronounced effect on the integrity of regular secondary structure compared to the loss of Asn45 or Asn270. These data implicate differential roles of N-glycosylation modifications in contributing to the stability of specific functional regions of CBH1 and highlight the potential of improving the thermostability of CBH1 by tuning proper interactions between glycans and functional residues.

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

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

Structural insight into the role of Streptococcus parasanguinis Fap1 within oral biofilm formation

The fimbriae-associated protein 1 (Fap1) is a major adhesin of Streptococcus parasanguinis, a primary colonizer of the oral cavity that plays an important role in the formation of dental plaque. Fap1 is an extracellular adhesive surface fibre belonging to the serine-rich repeat protein (SRRP) family, which plays a central role in the pathogenesis of streptococci and staphylococci. The N-terminal adhesive region of Fap1 (Fap1-NR) is composed of two domains (Fap1-NR(α) and Fap1-NR(β)) and is projected away from the bacterial surface via the extensive serine-rich repeat region, for adhesion to the salivary pellicle. The adhesive properties of Fap1 are modulated through a pH switch in which a reduction in pH results in a rearrangement between the Fap1-NR(α) and Fap1-NR(β) domains, which assists in the survival of S. parasanguinis in acidic environments. We have solved the structure of Fap1-NR(α) at pH 5.0 at 3.0Ǻ resolution and reveal how subtle rearrangements of the 3-helix bundle combined with a change in electrostatic potential mediates 'opening' and activation of the adhesive region. Further, we show that pH-dependent changes are critical for biofilm formation and present an atomic model for the inter-Fap1-NR interactions which have been assigned an important role in the biofilm formation.