Correlation of structure, function and protein dynamics in GH7 cellobiohydrolases from Trichoderma atroviride, T. reesei and T. harzianum

Biochemical Characterization of an Endoglucanase GH7 from Thermophile Thermothielavioides terrestris Expressed on Aspergillus nidulans

Endoglucanases (EC 3.2.1.4) are important enzymes involved in the hydrolysis of cellulose, acting randomly in the β-1,4-glycosidic bonds present in the amorphous regions of the polysaccharide chain. These biocatalysts have been classified into 14 glycosyl hydrolase (GH) families. The GH7 family is of particular interest since it may act on a broad range of substrates, including cellulose, β-glucan, and xylan, an attractive feature for biotechnological applications, especially in the renewable energy field. In the current work, a gene from the thermophilic fungus Thermothielavioides terrestris, encoding an endoglucanase GH7 (TtCel7B), was cloned in the secretion vector pEXPYR and transformed into the high-protein-producing strain Aspergillus nidulans A773. Purified TtCel7B has a molecular weight of approximately 66 kDa, evidenced by SDS-PAGE. Circular dichroism confirmed the high β-strand content consistent with the canonical GH7 family β-jellyroll fold, also observed in the 3D homology model of TtCel7B. Biochemical characterization assays showed that TtCel7B was active over a wide range of pH values (3.5–7.0) and temperatures (45–70 °C), with the highest activity at pH 4.0 and 65 °C. TtCel7B also was stable over a wide range of pH values (3.5–9.0), maintaining more than 80% of its activity after 24 h. The KM and Vmax values in low-viscosity carboxymethylcellulose were 9.3 mg mL−1 and 2.5 × 104 U mg−1, respectively. The results obtained in this work provide a basis for the development of applications of recombinant TtCel7B in the renewable energy field.

TtCel7A: A Native Thermophilic Bifunctional Cellulose/Xylanase Exogluclanase from the Thermophilic Biomass-Degrading Fungus Thielavia terrestris Co3Bag1, and Its Application in Enzymatic Hydrolysis of Agroindustrial Derivatives

The biomass-degrading thermophilic ascomycete fungus Thielavia terrestris Co3Bag1 produces TtCel7A, a native bifunctional cellulase/xylanase GH7 family. The purified TtCel7A, with an estimated molecular weight of 71 kDa, was biochemically characterized. TtCel7A displayed an optimal pH of 5.5 for both activities and an optimal temperature of 60 and 50 °C for cellulolytic and xylanolytic activities, respectively. The half-lives determined for cellulase activity were 140, 106, and 41 min at 50, 60, and 70 °C, respectively, whereas the half-lives observed for xylanase activity were 24, 10, and 1.4 h at 50, 60, and 70 °C, respectively. The K_M and V_max values were 3.12 mg/mL and 50 U/mg for cellulase activity and 0.17 mg/mL and 42.75 U/mg for xylanase activity. Circular dichroism analysis suggests changes in the secondary structure of TtCel7A in the presence of CMC as the substrate, whereas no modifications were observed with beechwood xylan. TtCel7A displayed the excellent capability to hydrolyze CMC, beechwood xylan, and complex substrates such as oat bran, wheat bran, and sugarcane bagasse, with glucose and cellobiose being the main products released; also, slightly less endo cellulase and xylanase activities were observed. Thus, suggesting TtCel7A has an exo- and endomode of action. Based on the characteristics of the enzyme, it might be considered a good candidate for industrial applications.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Protein Crystallization of Two Recombinant Lpt Proteins

The prerequisite for 3D structure determination of macromolecules via X-ray crystallography is well-ordered, diffracting crystals. Here, we report the recombinant production, biophysical/biochemical protein sample characterization, and vapor diffusion sitting drop crystallization protocols for two lipopolysaccharide transport proteins: LptH from Pseudomonas aeruginosa (Pa-LptH) and an inactive LptC mutant (G153R) from Escherichia coli (EcLptC_24-191G153R).

Heterologous production of cellulose- and starch-degrading hydrolases to expand Saccharomyces cerevisiae substrate utilization: Lessons learnt

Selected strains of Saccharomyces cerevisiae are used for commercial bioethanol production from cellulose and starch, but the high cost of exogenous enzymes for substrate hydrolysis remains a challenge. This can be addressed through consolidated bioprocessing (CBP) where S. cerevisiae strains are engineered to express recombinant glycoside hydrolases during fermentation. Looking back at numerous strategies undertaken over the past four decades to improve recombinant protein production in S. cerevisiae, it is evident that various steps in the protein production "pipeline" can be manipulated depending on the protein of interest and its anticipated application. In this review, we briefly introduce some of the strategies and highlight lessons learned with regards to improved transcription, translation, post-translational modification and protein secretion of heterologous hydrolases. We examine how host strain selection and modification, as well as enzyme compatibility, are crucial determinants for overall success. Finally, we discuss how lessons from heterologous hydrolase expression can inform modern synthetic biology and genome editing tools to provide process-ready yeast strains in future. However, it is clear that the successful expression of any particular enzyme is still unpredictable and requires a trial-and-error approach.

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.

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

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

Enzymatic processing of lignocellulosic biomass: principles, recent advances and perspectives

Efficient saccharification of lignocellulosic biomass requires concerted development of a pretreatment method, an enzyme cocktail and an enzymatic process, all of which are adapted to the feedstock. Recent years have shown great progress in most aspects of the overall process. In particular, increased insights into the contributions of a wide variety of cellulolytic and hemicellulolytic enzymes have improved the enzymatic processing step and brought down costs. Here, we review major pretreatment technologies and different enzyme process setups and present an in-depth discussion of the various enzyme types that are currently in use. We pay ample attention to the role of the recently discovered lytic polysaccharide monooxygenases (LPMOs), which have led to renewed interest in the role of redox enzyme systems in lignocellulose processing. Better understanding of the interplay between the various enzyme types, as they may occur in a commercial enzyme cocktail, is likely key to further process improvements.

Structure and dynamics of ionic liquid tolerant hyperthermophilic endoglucanase Cel12A from Rhodothermus marinus

Economic deconstruction of lignocellulose remains a challenge due to the complex architecture of cellulose, hemicellulose, and lignin. Advancements in pretreatment processes have introduced ionic liquids (ILs) as promising non-derivatizing solvents for reducing biomass recalcitrance and for promoting enzymatic hydrolysis. However, available commercial cellulases are destabilized or inactivated even in low concentration of residual ILs. Thus, a molecular understanding of IL-enzyme interactions is crucial for developing IL-tolerant enzymes with high catalytic activity. In this study, molecular insight behind the IL tolerance of hyperthermophilic endoglucanase Cel12A from Rhodothermus marinus (RmCel12A) has been investigated in 20%, 40%, and 60% 1-ethyl-3-methylimidazolium acetate (EmimAc) through molecular dynamic simulations at 368 K. Though the enzyme retained its stability in all EmimAc concentrations, the activity was affected due to the loss of essential dynamic motions. A protein structure network was constructed using the snapshots of protein structures from the simulation trajectories and the hub properties of residues R20, Y59, W68, W197, E203, and F220 were found to be lost in 60% EmimAc. Emim cations were observed to intrude the active site tunnel and interact with more number of catalytic residues with higher cumulative fractional occupancy in 60% EmimAc than in 20% or 40% EmimAc. Some non-catalytic residues have also been identified at the active site, which can be probable mutation targets for improving the IL tolerance. Our findings reveal the molecular understanding behind the origin of activity loss of RmCel12A and proposed insights for the further improvement of IL sensitivity.

Understanding the behavior of ionic liquid tolerant hyperthermophilic endoglucanase Cel12A from Rhodothermus marinus in different concentrations of EmimAc.

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

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

Ega3 from the fungal pathogen Aspergillus fumigatus is an endo-α-1,4-galactosaminidase that disrupts microbial biofilms

Aspergillus fumigatus is an opportunistic fungal pathogen that causes both chronic and acute invasive infections. Galactosaminogalactan (GAG) is an integral component of the A. fumigatus biofilm matrix and a key virulence factor. GAG is a heterogeneous linear α-1,4-linked exopolysaccharide of galactose and GalNAc that is partially deacetylated after secretion. A cluster of five co-expressed genes has been linked to GAG biosynthesis and modification. One gene in this cluster, ega3, is annotated as encoding a putative α-1,4-galactosaminidase belonging to glycoside hydrolase family 114 (GH114). Herein, we show that recombinant Ega3 is an active glycoside hydrolase that disrupts GAG-dependent A. fumigatus and Pel polysaccharide-dependent Pseudomonas aeruginosa biofilms at nanomolar concentrations. Using MS and functional assays, we demonstrate that Ega3 is an endo-acting α-1,4-galactosaminidase whose activity depends on the conserved acidic residues, Asp-189 and Glu-247. X-ray crystallographic structural analysis of the apo Ega3 and an Ega3-galactosamine complex, at 1.76 and 2.09 Å resolutions, revealed a modified (β/α)₈-fold with a deep electronegative cleft, which upon ligand binding is capped to form a tunnel. Our structural analysis coupled with in silico docking studies also uncovered the molecular determinants for galactosamine specificity and substrate binding at the -2 to +1 binding subsites. The findings in this study increase the structural and mechanistic understanding of the GH114 family, which has >600 members encoded by plant and opportunistic human pathogens, as well as in industrially used bacteria and fungi.

Synthesis of Human Milk Oligosaccharides: Protein Engineering Strategies for Improved Enzymatic Transglycosylation

Human milk oligosaccharides (HMOs) signify a unique group of oligosaccharides in breast milk, which is of major importance for infant health and development. The functional benefits of HMOs create an enormous impetus for biosynthetic production of HMOs for use as additives in infant formula and other products. HMO molecules can be synthesized chemically, via fermentation, and by enzymatic synthesis. This treatise discusses these different techniques, with particular focus on harnessing enzymes for controlled enzymatic synthesis of HMO molecules. In order to foster precise and high-yield enzymatic synthesis, several novel protein engineering approaches have been reported, mainly concerning changing glycoside hydrolases to catalyze relevant transglycosylations. The protein engineering strategies for these enzymes range from rationally modifying specific catalytic residues, over targeted subsite -1 mutations, to unique and novel transplantations of designed peptide sequences near the active site, so-called loop engineering. These strategies have proven useful to foster enhanced transglycosylation to promote different types of HMO synthesis reactions. The rationale of subsite -1 modification, acceptor binding site matching, and loop engineering, including changes that may alter the spatial arrangement of water in the enzyme active site region, may prove useful for novel enzyme-catalyzed carbohydrate design in general.

Gene Co-expression Network Reveals Potential New Genes Related to Sugarcane Bagasse Degradation in Trichoderma reesei RUT-30

The biomass-degrading fungus Trichoderma reesei has been considered a model for cellulose degradation, and it is the primary source of the industrial enzymatic cocktails used in second-generation (2G) ethanol production. However, although various studies and advances have been conducted to understand the cellulolytic system and the transcriptional regulation of T. reesei, the whole set of genes related to lignocellulose degradation has not been completely elucidated. In this study, we inferred a weighted gene co-expression network analysis based on the transcriptome dataset of the T. reesei RUT-C30 strain aiming to identify new target genes involved in sugarcane bagasse breakdown. In total, ~70% of all the differentially expressed genes were found in 28 highly connected gene modules. Several cellulases, sugar transporters, and hypothetical proteins coding genes upregulated in bagasse were grouped into the same modules. Among them, a single module contained the most representative core of cellulolytic enzymes (cellobiohydrolase, endoglucanase, β-glucosidase, and lytic polysaccharide monooxygenase). In addition, functional analysis using Gene Ontology (GO) revealed various classes of hydrolytic activity, cellulase activity, carbohydrate binding and cation:sugar symporter activity enriched in these modules. Several modules also showed GO enrichment for transcription factor activity, indicating the presence of transcriptional regulators along with the genes involved in cellulose breakdown and sugar transport as well as other genes encoding proteins with unknown functions. Highly connected genes (hubs) were also identified within each module, such as predicted transcription factors and genes encoding hypothetical proteins. In addition, various hubs contained at least one DNA binding site for the master activator Xyr1 according to our in silico analysis. The prediction of Xyr1 binding sites and the co-expression with genes encoding carbohydrate active enzymes and sugar transporters suggest a putative role of these hubs in bagasse cell wall deconstruction. Our results demonstrate a vast range of new promising targets that merit additional studies to improve the cellulolytic potential of T. reesei strains and to decrease the production costs of 2G ethanol.