A Recombinant Thermophilic and Glucose-Tolerant GH1 β-Glucosidase Derived from Hehua Hot Spring

As a crucial enzyme for cellulose degradation, β-glucosidase finds extensive applications in food, feed, and bioethanol production; however, its potential is often limited by inadequate thermal stability and glucose tolerance. In this study, a functional gene (lq-bg5) for a GH1 family β-glucosidase was obtained from the metagenomic DNA of a hot spring sediment sample and heterologously expressed in E. coli and the recombinant enzyme was purified and characterized. The optimal temperature and pH of LQ-BG5 were 55 °C and 4.6, respectively. The relative residual activity of LQ-BG5 exceeded 90% at 55 °C for 9 h and 60 °C for 6 h and remained above 100% after incubation at pH 5.0-10.0 for 12 h. More importantly, LQ-BG5 demonstrated exceptional glucose tolerance with more than 40% activity remaining even at high glucose concentrations of 3000 mM. Thus, LQ-BG5 represents a thermophilic β-glucosidase exhibiting excellent thermal stability and remarkable glucose tolerance, making it highly promising for lignocellulose development and utilization.

Thermostabilizing mechanisms of canonical single amino acid substitutions at a GH1 β-glucosidase probed by multiple MD and computational approaches

β-glucosidases play a pivotal role in second-generation biofuel (2G-biofuel) production. For this application, thermostable enzymes are essential due to the denaturing conditions on the bioreactors. Random amino acid substitutions have originated new thermostable β-glucosidases, but without a clear understanding of their molecular mechanisms. Here, we probe by different molecular dynamics simulation approaches with distinct force fields and submitting the results to various computational analyses, the molecular bases of the thermostabilization of the Paenibacillus polymyxa GH1 β-glucosidase by two-point mutations E96K (TR1) and M416I (TR2). Equilibrium molecular dynamic simulations (eMD) at different temperatures, principal component analysis (PCA), virtual docking, metadynamics (MetaDy), accelerated molecular dynamics (aMD), Poisson-Boltzmann surface analysis, grid inhomogeneous solvation theory and colony method estimation of conformational entropy allow to converge to the idea that the stabilization carried by both substitutions depend on different contributions of three classic mechanisms: (i) electrostatic surface stabilization; (ii) efficient isolation of the hydrophobic core from the solvent, with energetic advantages at the solvation cap; (iii) higher distribution of the protein dynamics at the mobile active site loops than at the protein core, with functional and entropic advantages. Mechanisms i and ii predominate for TR1, while in TR2, mechanism iii is dominant. Loop A integrity and loops A, C, D, and E dynamics play critical roles in such mechanisms. Comparison of the dynamic and topological changes observed between the thermostable mutants and the wildtype protein with amino acid co-evolutive networks and thermostabilizing hotspots from the literature allow inferring that the mechanisms here recovered can be related to the thermostability obtained by different substitutions along the whole family GH1. We hope the results and insights discussed here can be helpful for future rational approaches to the engineering of optimized β-glucosidases for 2G-biofuel production for industry, biotechnology, and science.

Role of metagenomics in prospecting novel endoglucanases, accentuating functional metagenomics approach in second-generation biofuel production: a review

As the fossil fuel reserves are depleting rapidly, there is a need for alternate fuels to meet the day to day mounting energy demands. As fossil fuel started depleting, a quest for alternate forms of fuel was initiated and biofuel is one of its promising outcomes. First-generation biofuels are made from edible sources like vegetable oils, starch, and sugars. Second-generation biofuels (SGB) are derived from lignocellulosic crops and the third-generation involves algae for biofuel production. Technical challenges in the production of SGB are hampering its commercialization. Advanced molecular technologies like metagenomics can help in the discovery of novel lignocellulosic biomass-degrading enzymes for commercialization and industrial production of SGB. This review discusses the metagenomic outcomes to enlighten the importance of unexplored habitats for novel cellulolytic gene mining. It also emphasizes the potential of different metagenomic approaches to explore the uncultivable cellulose-degrading microbiome as well as cellulolytic enzymes associated with them. This review also includes effective pre-treatment technology and consolidated bioprocessing for efficient biofuel production.

Biochemical characterization of a novel glucose-tolerant GH3 β-glucosidase (Bgl1973) from Leifsonia sp. ZF2019

Beta-glucosidase (Bgl) is an enzyme with considerable food, beverage, and biofuel processing potential. However, as many Bgls are inhibited by their reaction end product glucose, their industrial applications are greatly limited. In this study, a novel Bgl gene (Bgl1973) was cloned from Leifsonia sp. ZF2019 and heterologously expressed in E. coli. Sequence analysis and structure modeling revealed that Bgl1973 was 748 aa, giving it a molecular weight of 78 kDa, and it showed high similarity with the glycoside hydrolase 3 (GH3) family Bgls with which its active site residues were conserved. By using pNPGlc (p-nitrophenyl-β-D-glucopyranoside) as substrate, the optimum temperature and pH of Bgl1973 were shown to be 50 °C and 7.0, respectively. Bgl1973 was insensitive to most metal ions (12.5 mM), 1% urea, and even 0.1% Tween-80. This enzyme maintained 60% of its original activity in the presence of 20% NaCl, demonstrating its excellent salt tolerance. Furthermore, it still had 83% residual activity in 1 M of glucose, displaying its outstanding glucose tolerance. The K_m, V_max, and k_cat of Bgl1973 were 0.22 mM, 44.44 μmol/min mg, and 57.78 s^-1, respectively. Bgl1973 had a high specific activity for pNPGlc (19.10 ± 0.59 U/mg) and salicin (20.43 ± 0.92 U/mg). Furthermore, molecular docking indicated that the glucose binding location and the narrow and deep active channel geometry might contribute to the glucose tolerance of Bgl1973. Our results lay a foundation for the studying of this glucose-tolerant β-glucosidase and its applications in many industrial settings. KEY POINTS: • A novel β-glucosidase from GH3 was obtained from Leifsonia sp. ZF2019. • Bgl1973 demonstrated excellent glucose tolerance. • The glucose tolerance of Bgl1973 was explained using molecular docking analysis.

Fungal cellulases: protein engineering and post-translational modifications

Enzymatic degradation of lignocelluloses into fermentable sugars to produce biofuels and other biomaterials is critical for environmentally sustainable development and energy resource supply. However, there are problems in enzymatic cellulose hydrolysis, such as the complex cellulase composition, low degradation efficiency, high production cost, and post-translational modifications (PTMs), all of which are closely related to specific characteristics of cellulases that remain unclear. These problems hinder the practical application of cellulases. Due to the rapid development of computer technology in recent years, computer-aided protein engineering is being widely used, which also brings new opportunities for the development of cellulases. Especially in recent years, a large number of studies have reported on the application of computer-aided protein engineering in the development of cellulases; however, these articles have not been systematically reviewed. This article focused on the aspect of protein engineering and PTMs of fungal cellulases. In this manuscript, the latest literatures and the distribution of potential sites of cellulases for engineering have been systematically summarized, which provide reference for further improvement of cellulase properties. KEY POINTS: •Rational design based on virtual mutagenesis can improve cellulase properties. •Modifying protein side chains and glycans helps obtain superior cellulases. •N-terminal glutamine-pyroglutamate conversion stabilizes fungal cellulases.

VTR: A Web Tool for Identifying Analogous Contacts on Protein Structures and Their Complexes

Evolutionarily related proteins can present similar structures but very dissimilar sequences. Hence, understanding the role of the inter-residues contacts for the protein structure has been the target of many studies. Contacts comprise non-covalent interactions, which are essential to stabilize macromolecular structures such as proteins. Here we show VTR, a new method for the detection of analogous contacts in protein pairs. The VTR web tool performs structural alignment between proteins and detects interactions that occur in similar regions. To evaluate our tool, we proposed three case studies: we 1) compared vertebrate myoglobin and truncated invertebrate hemoglobin; 2) analyzed interactions between the spike protein RBD of SARS-CoV-2 and the cell receptor ACE2; and 3) compared a glucose-tolerant and a non-tolerant β-glucosidase enzyme used for biofuel production. The case studies demonstrate the potential of VTR for the understanding of functional similarities between distantly sequence-related proteins, as well as the exploration of important drug targets and rational design of enzymes for industrial applications. We envision VTR as a promising tool for understanding differences and similarities between homologous proteins with similar 3D structures but different sequences. VTR is available at http://bioinfo.dcc.ufmg.br/vtr.

Cloning, expression, biochemical characterization, and molecular docking studies of a novel glucose tolerant β-glucosidase from Saccharomonospora sp. NB11

Most of the presently known β-glucosidases are sensitive to end-product inhibition by glucose, restricting their potential use in many industrial applications. Identification of novel glucose tolerant β-glucosidase can prove a pivotal solution to eliminate end-product inhibition and enhance the overall lignocellulosic saccharification process. In this study, a novel gene encoding β-glucosidase BglNB11 of 1405bp was identified in the genome of Saccharomonospora sp. NB11 and was successfully cloned and heterologously expressed in E. coli BL21 (DE3).The presence of conserved amino acids; NEPW and TENG indicated that BglNB11 belonged to GH1 β-glucosidases. The recombinant enzyme was purified using a Ni-NTA column, with the molecular mass of 51 kDa, using SDS-PAGE analysis. BglNB11 showed optimum activity at 40 °C and pH 7 and did not require any tested co-factors for activation. The kinetic values, K_m, V_max, k_cat, and k_cat/K_m of purified enzyme were 0.4037 mM, 5735.8 μmol/min/mg, 5042.16 s^-1 and 12487.71 s^-1 mM^-1, respectively. The enzyme was not inhibited by glucose to a concentration of 4 M but was slightly stimulated in the presence of glucose. Molecular docking of BglNB11 with glucose suggested that the relative binding position of glucose in the active site channel might be responsible for modulating end product tolerance and stimulation. β-glucosidase from BglNB11 is an excellent enzyme with high catalytic efficiency and enhanced glucose tolerance compared to many known glucose tolerant β-glucosidases. These unique properties of BglNB11 make it a prime candidate to be utilized in many biotechnological applications.

A novel β-glucosidase from a hot-spring metagenome shows elevated thermal stability and tolerance to glucose and ethanol

β-glucosidase causes hydrolysis of β-1,4-glycosidic bond in glycosides and oligosaccharides. It is an industrially important enzyme owing to its potential in biomass processing applications. In this study, computational screening of an extreme temperature aquatic habitat metagenomic resource was done, leading to the identification of a novel gene, bgl_M, encoding a β-glucosidase. The comparative protein sequence and homology structure analyses designated it as a GH1 family β-glucosidase. The bgl_M gene was expressed in a heterologous host, Escherichia coli. The purified protein, Bgl_M, was biochemically characterized for β-glucosidase activity. Bgl_M exhibited noteworthy hydrolytic potential towards cellobiose and lactose. Bgl_M, showed substantial catalytic activity in the pH range of 5.0-7.0 and at the temperature 40 °C-70 °C. The enzyme was found quite stable at 50 °C with a loss of hardly 20% after 40 h of heat exposure. Furthermore, any drastically negative effect was not observed on the enzyme's activity in the presence of metal ions, non-ionic surfactants, metal chelating, and denaturing agents. A significantly high glucose tolerance, retaining 80% relative activity at 1 M, and 40% at 5 M glucose, and ethanol tolerance, exhibiting 80% relative activity in 10% ethanol, enrolled Bgl_M as a promising enzyme for cellulose saccharification. Furthermore, its ability to catalyze the hydrolysis of daidzin and polydatin ascertained it as an admirably suited biocatalyst for enhancement of nutritional values in soya and wine industries.

Proteus: An algorithm for proposing stabilizing mutation pairs based on interactions observed in known protein 3D structures

Background

Protein engineering has many applications for industry, such as the development of new drugs, vaccines, treatment therapies, food, and biofuel production. A common way to engineer a protein is to perform mutations in functionally essential residues to optimize their function. However, the discovery of beneficial mutations for proteins is a complex task, with a time-consuming and high cost for experimental validation. Hence, computational approaches have been used to propose new insights for experiments narrowing the search space and reducing the costs.

Results

In this study, we developed Proteus (an acronym for Protein Engineering Supporter), a new algorithm for proposing mutation pairs in a target 3D structure. These suggestions are based on contacts observed in other known structures from Protein Data Bank (PDB). Proteus’ basic assumption is that if a non-interacting pair of amino acid residues in the target structure is exchanged to an interacting pair, this could enhance protein stability. This trade is only allowed if the main-chain conformation of the residues involved in the contact is conserved. Furthermore, no steric impediment is expected between the proposed mutations and the surrounding protein atoms. To evaluate Proteus, we performed two case studies with proteins of industrial interests. In the first case study, we evaluated if the mutations suggested by Proteus for four protein structures enhance the number of inter-residue contacts. Our results suggest that most mutations proposed by Proteus increase the number of interactions into the protein. In the second case study, we used Proteus to suggest mutations for a lysozyme protein. Then, we compared Proteus’ outcomes to mutations with available experimental evidence reported in the ProTherm database. Four mutations, in which our results agree with the experimental data, were found. This could be initial evidence that changes in the side-chain of some residues do not cause disturbances that harm protein structure stability.

Conclusion

We believe that Proteus could be used combined with other methods to give new insights into the rational development of engineered proteins. Proteus user-friendly web-based tool is available at <http://proteus.dcc.ufmg.br>.

Glutantβase: a database for improving the rational design of glucose-tolerant β-glucosidases

Β-glucosidases are key enzymes used in second-generation biofuel production. They act in the last step of the lignocellulose saccharification, converting cellobiose in glucose. However, most of the β-glucosidases are inhibited by high glucose concentrations, which turns it a limiting step for industrial production. Thus, β-glucosidases have been targeted by several studies aiming to understand the mechanism of glucose tolerance, pH and thermal resistance for constructing more efficient enzymes. In this paper, we present a database of β-glucosidase structures, called Glutantβase. Our database includes 3842 GH1 β-glucosidase sequences collected from UniProt. We modeled the sequences by comparison and predicted important features in the 3D-structure of each enzyme. Glutantβase provides information about catalytic and conserved amino acids, residues of the coevolution network, protein secondary structure, and residues located in the channel that guides to the active site. We also analyzed the impact of beneficial mutations reported in the literature, predicted in analogous positions, for similar enzymes. We suggested these mutations based on six previously described mutants that showed high catalytic activity, glucose tolerance, or thermostability (A404V, E96K, H184F, H228T, L441F, and V174C). Then, we used molecular docking to verify the impact of the suggested mutations in the affinity of protein and ligands (substrate and product). Our results suggest that only mutations based on the H228T mutant can reduce the affinity for glucose (product) and increase affinity for cellobiose (substrate), which indicates an increment in the resistance to product inhibition and agrees with computational and experimental results previously reported in the literature. More resistant β-glucosidases are essential to saccharification in industrial applications. However, thermostable and glucose-tolerant β-glucosidases are rare, and their glucose tolerance mechanisms appear to be related to multiple and complex factors. We gather here, a set of information, and made predictions aiming to provide a tool for supporting the rational design of more efficient β-glucosidases. We hope that Glutantβase can help improve second-generation biofuel production. Glutantβase is available at http://bioinfo.dcc.ufmg.br/glutantbase .

In-Silico Characterization of Glycosyl Hydrolase Family 1 β-Glucosidase from Trichoderma asperellum UPM1

β-glucosidases (Bgl) are widely utilized for releasing non-reducing terminal glucosyl residues. Nevertheless, feedback inhibition by glucose end product has limited its application. A noticeable exception has been found for β-glucosidases of the glycoside hydrolase (GH) family 1, which exhibit tolerance and even stimulation by glucose. In this study, using local isolate Trichoderma asperellum UPM1, the gene encoding β-glucosidase from GH family 1, hereafter designated as TaBgl2, was isolated and characterized via in-silico analyses. A comparison of enzyme activity was subsequently made by heterologous expression in Escherichia coli BL21(DE3). The presence of N-terminal signature, cis-peptide bonds, conserved active site motifs, non-proline cis peptide bonds, substrate binding, and a lone conserved stabilizing tryptophan (W) residue confirms the identity of Trichoderma sp. GH family 1 β-glucosidase isolated. Glucose tolerance was suggested by the presence of 14 of 22 known consensus residues, along with corresponding residues L167 and P172, crucial in the retention of the active site's narrow cavity. Retention of 40% of relative hydrolytic activity on ρ-nitrophenyl-β-D-glucopyranoside (ρNPG) in a concentration of 0.2 M glucose was comparable to that of GH family 1 β-glucosidase (Cel1A) from Trichoderma reesei. This research thus underlines the potential in the prediction of enzymatic function, and of industrial importance, glucose tolerance of family 1 β-glucosidases following relevant in-silico analyses.

Engineering of β-Glucosidase Bgl15 with Simultaneously Enhanced Glucose Tolerance and Thermostability To Improve Its Performance in High-Solid Cellulose Hydrolysis

In this study, a Petri-dish-based double-layer high-throughput screening method was established to improve glucose tolerance of β-glucosidase Bgl15. Two beneficial mutations were identified, and the joint mutant 2R1 improved the half-maximal inhibitory concentration of glucose from 0.04 to 2.1 M. The crystal structure of 2R1 was subsequently determined at 2.7 Å. Structure analysis revealed that enhancement of glucose tolerance may be due to improved transglycosylation activity made possible by a hydrophobic binding site for glucose as an acceptor and more stringent control of a putative water channel. To further ameliorate the application potential of the enzyme, it was engineered to increase the half-life at 50 °C from 0.8 h (Bgl15) to 180 h (mutant 5R1). Furthermore, supplementation of 5R1 to the cellulase cocktail significantly improved glucose production from pretreated sugar cane bagasse by 38%. Consequently, this study provided an efficient approach to enhance glucose tolerance and generated a promising catalyst for cellulose saccharification.

Conformational flexibility correlates with glucose tolerance for point mutations in β-glucosidases - a computational study

β-glucosidases (EC 3.2.1.21) have been described as essential to second-generation biofuel production. They act in the last step of the lignocellulosic saccharification, cleaving the β - 1,4 glycosidic bonds in cellobiose to produce two molecules of glucose. However, β-glucosidases have been described as strongly inhibited by glucose, causing an increment of cellobiose concentration. Also, cellobiose is an inhibitor of other enzymes used in this process, such as exoglucanases and endoglucanases. Hence, the engineering of thermostable and glucose-tolerant β-glucosidases has been targeted by many studies. In this study, we performed high sampling accelerated molecular dynamics for a wild glucose-tolerant GH1 β-glucosidase (Bgl1A), a wild non-tolerant (Bgl1B), and a set of glucose-tolerant Bgl1B's mutants: V302F, N301Q/V302F, F172I, V227M, G246S, T299S, and H228T. Our results suggest that point mutations promissory to induce glucose tolerance trend to enhance the mobility of the flexible loops around the active site. Mutations affected B and C loops regions, and an αβ-hairpin motif between them. Conformational clusters and free energy landscape profiles suggest that the mobility acquired by mutants allows a higher closure of the substrate channel. This closure is compatible with a higher impedance for glucose entrance and stimulus of its withdrawal. Based on mutants' structural analyses, we inferred that both the direct stereochemical effect on the glucose path and the changes in the mobility affect glucose tolerance. We hope these results be useful for the rational design of glucose-tolerant and industrially promising enzymes.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Gives New Insights into the Glucose Tolerance and Inhibition Mechanisms on β-Glucosidases

β-Glucosidases are enzymes with high importance for many industrial processes, catalyzing the last and limiting step of the conversion of lignocellulosic material into fermentable sugars for biofuel production. However, β-glucosidases are inhibited by high concentrations of the product (glucose), which limits the biofuel production on an industrial scale. For this reason, the structural mechanisms of tolerance to product inhibition have been the target of several studies. In this study, we performed in silico experiments, such as molecular dynamics (MD) simulations, free energy landscape (FEL) estimate, Poisson-Boltzmann surface area (PBSA), and grid inhomogeneous solvation theory (GIST) seeking a better understanding of the glucose tolerance and inhibition mechanisms of a representative GH1 β-glucosidase and a GH3 one. Our results suggest that the hydrophobic residues Y180, W350, and F349, as well the polar one D238 act in a mechanism for glucose releasing, herein called "slingshot mechanism", dependent also on an allosteric channel (AC). In addition, water activity modulation and the protein loop motions suggest that GH1 β-Glucosidases present an active site more adapted to glucose withdrawal than GH3, in consonance with the GH1s lower product inhibition. The results presented here provide directions on the understanding of the molecular mechanisms governing inhibition and tolerance to the product in β-glucosidases and can be useful for the rational design of optimized enzymes for industrial interests.

A Computational Method to Propose Mutations in Enzymes Based on Structural Signature Variation (SSV)

With the use of genetic engineering, modified and sometimes more efficient enzymes can be created for different purposes, including industrial applications. However, building modified enzymes depends on several in vitro experiments, which may result in the process being expensive and time-consuming. Therefore, computational approaches could reduce costs and accelerate the discovery of new technological products. In this study, we present a method, called structural signature variation (SSV), to propose mutations for improving enzymes’ activity. SSV uses the structural signature variation between target enzymes and template enzymes (obtained from the literature) to determine if randomly suggested mutations may provide some benefit for an enzyme, such as improvement of catalytic activity, half-life, and thermostability, or resistance to inhibition. To evaluate SSV, we carried out a case study that suggested mutations in β-glucosidases: Essential enzymes used in biofuel production that suffer inhibition by their product. We collected 27 mutations described in the literature, and manually classified them as beneficial or not. SSV was able to classify the mutations with values of 0.89 and 0.92 for precision and specificity, respectively. Then, we used SSV to propose mutations for Bgl1B, a low-performance β-glucosidase. We detected 15 mutations that could be beneficial. Three of these mutations (H228C, H228T, and H228V) have been related in the literature to the mechanism of glucose tolerance and stimulation in GH1 β-glucosidase. Hence, SSV was capable of detecting promising mutations, already validated by in vitro experiments, that improved the inhibition resistance of a β-glucosidase and, consequently, its catalytic activity. SSV might be useful for the engineering of enzymes used in biofuel production or other industrial applications.

Enhancing the Thermostability of Highly Active and Glucose-Tolerant β-Glucosidase Ks5A7 by Directed Evolution for Good Performance of Three Properties

A high-performance β-glucosidase for efficient cellulose hydrolysis needs to excel in thermostability, catalytic efficiency, and resistance to glucose inhibition. However, it is challenging to achieve superb properties in all three aspects in a single enzyme. In this study, a hyperactive and glucose-tolerant β-glucosidase Ks5A7 was employed as the starting point. Four rounds of random mutagenesis were then performed, giving rise to a thermostable mutant 4R1 with five amino acid substitutions. The half-life of 4R1 at 50 °C is 8640-fold that of Ks5A7 (144 h vs 1 min). Meanwhile, 4R1 had a higher specific activity (374.26 vs 243.18 units·mg^-1) than the wild type with a similar glucose tolerance. When supplemented to Celluclast 1.5L, the mutant significantly enhanced the hydrolysis of pretreated sugar cane bagasse, improving the released glucose concentration by 44%. With excellent performance in thermostability, activity, and glucose tolerance, 4R1 will serve as an exceptional catalyst for industrial applications.