Computational modeling and small-angle X-ray scattering based structure analysis and identifying ligand cleavage mechanism by processive endocellulase of family 9 glycoside hydrolase (HtGH9) from Hungateiclostridium thermocellum ATCC 27405

Molecular dynamics-based structural insights of the first putative endoglucanase, PsGH5A of glycoside hydrolase family 5 from Pseudopedobacter saltans

CONTEXT

The putative endoglucanase, PsGH5A, from Pseudopedobacter saltans of family GH5 contains a catalytic module, PsGH5 (β/α)₈-TIM barrel), at N-terminal followed by a family 6 carbohydrate-binding module (CBM6, β-sandwich). Superposition of PsGH5A with PDB homologs revealed Glu220 and Glu318 as evolutionarily conserved and catalytic residues performing the hydrolysis through retaining-type mechanism, a canonical property of GH5 family. PsGH5A showed higher affinity for longer cellooligosaccharides, as long as cellodecaose with binding free energy (∆G) of - 13.72 kcal/mol upon the molecular docking, thereby indicating the endo-mode of hydrolysis. The radius of gyration, Rg (2.7 nm), and solvent accessible surface area, SASA (229.6 nm²), of PsGH5A-Cellotetraose complex were determined by MD simulation which was lower than that of PsGH5A (Rg, 2.8 nm, SASA, 267 nm²) demonstrating the compactness and affinity of PsGH5A with the cellulosic ligands. Cellulose compatibility of PsGH5A was further confirmed by MMPBSA and per-residue decomposition analysis, where notable ∆G of - 54.38 kcal/mol for PsGH5A-Cellotetraose complex was observed. Thus, PsGH5A could be potentially an efficient endoglucanase as it accommodated larger cellooligosaccharides at its active-site. PsGH5A is the first putative endoglucanase studied here from P. saltans which could be genome-mined for lignocellulosic biomass saccharification in the renewable energy sector.

METHODS

The 3-D structure of PsGH5A generated by AlphaFold2, RaptorX, SwissModel, Phyre2 and Robetta tool; YASARA was used for energy minimization of built models. UCLA SAVES-v6 was used for quality assessment of models. Molecular Docking was performed using SWISS-DOCK server and Chimera software. Molecular Dynamics simulations and MMPBSA analysis of PsGH5A and PsGH5A-Cellotetraose complex were performed on GROMACS 2019.6.

Contribution of calcium ligands in substrate binding and product release in the Acetovibrio thermocellus glycoside hydrolase family 9 cellulase CelR

Enzymatic deconstruction of lignocellulosic biomass is crucial to establishment of the renewable biofuel and bioproduct economy. Better understanding of these enzymes, including their catalytic and binding domains, and other features offer potential avenues for improvement. Glycoside hydrolase family 9 (GH9) enzymes are attractive targets because they have members that exhibit exo- and endo-cellulolytic activity, processivity of reaction, and thermostability. This study examines a GH9 from Acetovibrio thermocellus ATCC 27405, AtCelR containing a catalytic domain and a carbohydrate binding module (CBM3c). Crystal structures of the enzyme without substrate, bound to cellohexaose (substrate) or cellobiose (product), show the positioning of ligands to calcium and adjacent residues in the catalytic domain that may contribute to substrate binding and facilitate product release. We also investigated the properties of the enzyme engineered to contain an additional carbohydrate binding module (CBM3a). Relative to the catalytic domain alone, CBM3a gave improved binding for Avicel (a crystalline form of cellulose), and catalytic efficiency (kcat/KM) was improved 40× with both CBM3c and CBM3a present. However, because of the molecular weight added by CBM3a, the specific activity of the engineered enzyme was not increased relative to the native construct consisting of only the catalytic and CBM3c domains. This work provides new insight into a potential role of the conserved calcium in the catalytic domain and identifies contributions and limitations of domain engineering for AtCelR and perhaps other GH9 enzymes.

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.

Computational and SAXS-based structure insights of pectin acetyl esterase (CtPae12B) of family 12 carbohydrate esterase from Clostridium thermocellum ATCC 27405

Pectin is a complex form of polysaccharide and is composed of several structural components that require the concerted action of several pectinases for its complete degradation. In this study, in silico and solution structure of a pectin acetyl esterase (CtPae12B) of family 12 carbohydrate esterase (CE12) from Clostridium thermocellum was determined. The CtPae12B modelled structure, showed a new α/β hydrolase fold, similar to the fold found in the crystal structures of its nearest homologues from CE12 family, which differed from α/β hydrolase fold found in glycoside hydrolases. In the active site of CtPae12B, two loops (loop1 and loop6) play an important role in the formation of a catalytic triad Ser15-Asp187-His190, where Ser15 acts as a nucleophile. The structural stability of CtPae12B and its catalytic site was detected by performing molecular dynamic (MD) simulation which showed stable and compact conformation of the structure. Molecular docking method was employed to analyse the conformations of various suitable ligands docked at the active site of CtPae12B. The stability and structural specificity of the catalytic residues with the ligand, 4-nitrophenyl acetate (4-NPA) was confirmed by MD simulation of CtPae12B-4NPA docked complex. Moreover, it was found that the nucleophile Ser15, forms hydrophobic interaction with 4-NPA in the active site to complete covalent catalysis. Small angle X-ray scattering analysis of CtPae12B at 3 mg/mL displayed elongated, compact and monodispersed nature in solution. The ab initio derived dummy model showed that CtPae12B exists as a homotrimer at 3 mg/mL which was also confirmed by dynamic light scattering.Communicated by Ramaswamy H. Sarma.