Mutational Effects on the Catalytic Mechanism of Cellobiohydrolase I from Trichoderma reesei

A Case Study of the Glycoside Hydrolase Enzyme Mechanism Using an Automated QM-Cluster Model Building Toolkit

Glycoside hydrolase enzymes are important for hydrolyzing the β-1,4 glycosidic bond in polysaccharides for deconstruction of carbohydrates. The two-step retaining reaction mechanism of Glycoside Hydrolase Family 7 (GH7) was explored with different sized QM-cluster models built by the Residue Interaction Network ResidUe Selector (RINRUS) software using both the wild-type protein and its E217Q mutant. The first step is the glycosylation, in which the acidic residue 217 donates a proton to the glycosidic oxygen leading to bond cleavage. In the subsequent deglycosylation step, one water molecule migrates into the active site and attacks the anomeric carbon. Residue interaction-based QM-cluster models lead to reliable structural and energetic results for proposed glycoside hydrolase mechanisms. The free energies of activation for glycosylation in the largest QM-cluster models were predicted to be 19.5 and 31.4 kcal mol⁻¹ for the wild-type protein and its E217Q mutant, which agree with experimental trends that mutation of the acidic residue Glu217 to Gln will slow down the reaction; and are higher in free energy than the deglycosylation transition states (13.8 and 25.5 kcal mol⁻¹ for the wild-type protein and its mutant, respectively). For the mutated protein, glycosylation led to a low-energy product. This thermodynamic sink may correspond to the intermediate state which was isolated in the X-ray crystal structure. Hence, the glycosylation is validated to be the rate-limiting step in both the wild-type and mutated enzyme.

Mechanism and biomass association of glucuronoyl esterase: an α/β hydrolase with potential in biomass conversion

Glucuronoyl esterases (GEs) are α/β serine hydrolases and a relatively new addition in the toolbox to reduce the recalcitrance of lignocellulose, the biggest obstacle in cost-effective utilization of this important renewable resource. While biochemical and structural characterization of GEs have progressed greatly recently, there have yet been no mechanistic studies shedding light onto the rate-limiting steps relevant for biomass conversion. The bacterial GE OtCE15A possesses a classical yet distinctive catalytic machinery, with easily identifiable catalytic Ser/His completed by two acidic residues (Glu and Asp) rather than one as in the classical triad, and an Arg side chain participating in the oxyanion hole. By QM/MM calculations, we identified deacylation as the decisive step in catalysis, and quantified the role of Asp, Glu and Arg, showing the latter to be particularly important. The results agree well with experimental and structural data. We further calculated the free-energy barrier of post-catalysis dissociation from a complex natural substrate, suggesting that in industrial settings non-catalytic processes may constitute the rate-limiting step, and pointing to future directions for enzyme engineering in biomass utilization.

Zong and coworkers combine computational and experimental methods to decipher in detail the mechanism of action of glucuronoyl esterases, enzymes with significant biotechnological potential for decoupling lignin from polysaccharides in biomass.

A structured understanding of cellobiohydrolase I binding to poplar lignin fractions after dilute acid pretreatment

BACKGROUND

Cellulase adsorption to lignin is considered a cost barrier for bioethanol production; however, its detailed association mechanism is still not fully understood. In this study, two natural poplar variants with high and low sugar release performance were selected as the low and high recalcitrant raw materials (named L and H, respectively). Three different lignin fractions were extracted using ethanol, followed by p-dioxane and then cellulase treatment from the dilute acid pretreated poplar solids (fraction 1, 2, and 3, respectively).

RESULTS

Each lignin fraction had different physicochemical properties. Ethanol-extracted lignin had the lowest weight average molecular weight, while the molecular weights for the other two lignin fractions were similar. ³¹P NMR analysis revealed that lignin fraction with higher molecular weight contained more aliphatic hydroxyl groups and less phenolic hydroxyl groups. Semi-quantitative analysis by 2D HSQC NMR indicated that the lignin fractions isolated from the natural variants had different contents of syringyl (S), guaiacyl (G) and interunit linkages. Lignin extracted by ethanol contained the largest amount of S units, the smallest amounts of G and p-hydroxybenzoate (PB) subunits, while the contents of these lignin subunits in the other two lignin fractions were similar. The lignin fraction obtained after cellulase treatment was primarily comprised of β-O-4 linkages with small amounts of β-5 and β-β linkages. The binding strength of these three lignin fractions obtained by Langmuir equations were in the order of L₁ > L₃ > L₂ for the low recalcitrance poplar and H₁ > H₂ > H₃ for the high recalcitrance poplar.

CONCLUSIONS

Overall, adsorption ability of lignin was correlated with the sugar release of poplar. Structural features of lignin were associated with its binding to CBH. For natural poplar variants, lignin fractions with lower molecular weight and polydispersity index (PDI) exhibited more CBH adsorption ability. Lignins with more phenolic hydroxyl groups had higher CBH binding strength. It was also found that lignin fractions with more condensed aromatics adsorbed more CBH likely attributed to stronger hydrophobic interactions.

Solvent effect on hydrogen bonded Tyr⋯Asp⋯Arg triads: Enzymatic catalyzed model system

The hydrogen bond plays a vital role in structural arrangement, intermediate state stabilization, materials function, and biological activity of certain enzymatic reactions. The solvent and electronic effects on hydrogen bonds are illustrated employing the polarizable contimuum model at B3LYP/6-311++G(d,p) level. Geometry optimizations reflect the significant solvent and electronic effect. The proton departs spontaneously upon oxidation from the hydroxyl group of tyrosyl in hydrogen bonded Tyr⋯Asp⋯Arg triads in both gas phase and solvents. The electron transfer isomers are observed for anionic triads, no matter what the solvent is. The difference of distance between two hydrogen bonds is enlarged in solvent as compared to that in gas phase. The electronic effect on IR spectra is distinctive. The tyrosyl fragment tends to be oxidized and the arginine moiety is easier to bind an excess electron. The variations of chemical shift and spin-spin coupling constant are more significant upon electron transfer than upon solvent dielectric constant. The augmentation of solvent dielectric constant stabilizes the system, enhances the difference of isomers, and increases the vertical ionization potential and vertical electron affinity values.

Improving the activity of Trichoderma reesei cel7B through stabilizing the transition state

Trichoderma reesei (Tr.) cellulases, which convert cellulose to reducing sugars, are a promising catalyst used in the lignocellulosic biofuel production. Improving Tr. cellulases activity, though very difficult, is highly desired due to the recalcitrance of lignocellulose. Meanwhile, it is preferable to enhance the cellulase's promiscuity so that substrates other than cellulose can also be hydrolyzed. In this work, an attempt is made to improve the catalytic activity of a major endogluanase Tr. Cel7B against xylan which crosslinks with cellulose in lignocellulose. By using quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations, the transition state of the xylo-oligosaccharide hydrolysis is identified. Then, mutations are introduced and their effect on the transition state stabilization is ranked based on the free energy calculations. Seven top ranked mutants are evaluated experimentally. Three mutants A208Q, A222D, and G230R show a higher activity than the wild-type Tr. Cel7B in the hydrolysis of xylan (by up to 47%) as well as filter paper (by up to 50%). The combination of the single mutants can further improve the enzyme activity. Our work demonstrates that the free energy method is effective in engineering the Tr. Cel7B activity against xylan and cellulose, and thus may also be useful for improving the activity of other Tr. cellulases. Biotechnol. Bioeng. 2016;113: 1171-1177. © 2015 Wiley Periodicals, Inc.

Free Energy Diagram for the Heterogeneous Enzymatic Hydrolysis of Glycosidic Bonds in Cellulose

Kinetic and thermodynamic data have been analyzed according to transition state theory and a simplified reaction scheme for the enzymatic hydrolysis of insoluble cellulose. For the cellobiohydrolase Cel7A from Hypocrea jecorina (Trichoderma reesei), we were able to measure or collect relevant values for all stable and activated complexes defined by the reaction scheme and hence propose a free energy diagram for the full heterogeneous process. For other Cel7A enzymes, including variants with and without carbohydrate binding module (CBM), we obtained activation parameters for the association and dissociation of the enzyme-substrate complex. The results showed that the kinetics of enzyme-substrate association (i.e. formation of the Michaelis complex) was almost entirely entropy-controlled and that the activation entropy corresponded approximately to the loss of translational and rotational degrees of freedom of the dissolved enzyme. This implied that the transition state occurred early in the path where the enzyme has lost these degrees of freedom but not yet established extensive contact interactions in the binding tunnel. For dissociation, a similar analysis suggested that the transition state was late in the path where most enzyme-substrate contacts were broken. Activation enthalpies revealed that the rate of dissociation was far more temperature-sensitive than the rates of both association and the inner catalytic cycle. Comparisons of one- and two-domain variants showed that the CBM had no influence on the transition state for association but increased the free energy barrier for dissociation. Hence, the CBM appeared to promote the stability of the complex by delaying dissociation rather than accelerating association.

Characterization of the Dielectric Constant in the Trichoderma reesei Cel7B Active Site

An attempt is made to evaluate the dielectric constant of the Trichoderma reesei Cel7B active site. Through kinetic measurements, the pKa value of the catalytic acid E201 is determined. Mutations (away from E201) with net charge changes are introduced to perturb the E201 pKa. It is shown that the mutation with a +1 charge change (including G225R, G230R, and A335R) decreases the pKa of E201, whereas the mutation with a -1 charge change (including Q149E, A222D, G225D, and G230D) increases the pKa. This effect is consistent with the electrostatic interaction between the changed charge and the E201 side chain. The fitting of the experimental data yields an apparent dielectric constant of 25-80. Molecular dynamics simulations with explicit water molecules indicate that the high solvent accessibility of the active site contributes largely to the high dielectric constant. ONIOM calculations show that high dielectric constant benefits the catalysis through decreasing the energy of the transition state relative to that of the enzyme substrate complex.

The impact of active site protonation on substrate ring conformation in Melanocarpus albomyces cellobiohydrolase Cel7B

Understanding how the protonation state of active site residues impacts the enzyme's structure and substrate conformation is important for improving the efficiency and economic viability of the degradation of cellulosic materials as feedstock for liquid fuel and value-added chemicals.

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.