Glycosidic-Bond Hydrolysis Mechanism Catalyzed by Cellulase Cel7A from Trichoderma reesei: A Comprehensive Theoretical Study by Performing MD, QM, and QM/MM Calculations

Combination of engineering the substrate and Ca2+ binding domains of heparinase I to improve the catalytic activity

Heparinase I (EC 4.2.2.7), is an enzyme that cleaves heparin, showing great potential for eco-friendly production of low molecular weight heparin (LMWH). However, owing to its poor catalytic activity and thermal stability, the industrial application of heparinase I has been severely hindered. To improve the catalytic activity, we proposed to engineer both the substrate and Ca2+ binding domains of heparinase I. Several heparinases I from different organisms were selected for multiple sequence alignment and molecular docking to screen the key residues in the binding domain. Nine single-point mutations were selected to enhance the catalytic activity of heparinase I. Among them, T250D was the most highly active one, whereas mutations around Ca2+ binding domain yielded two active mutants. Mutant D152S/R244K/T250D with significantly increased catalytic activity was obtained by combined mutation. The catalytic efficiency of the mutant was 118,875.8 min-1·µM-1, which was improved 5.26 times. Molecular modeling revealed that the improved activity and stability of the mutants were probably attributed to the formation of new hydrogen bonds. The highly active mutant had great potential applications in industry and the strategy could be used to improve the performance of other enzymes.

CelS-Catalyzed Processive Cellulose Degradation and Cellobiose Extraction for the Production of Bioethanol

Bacterial cellulase enzymes are potent candidates for the efficient production of bioethanol, a promising alternative to fossil fuels, from cellulosic biomass. These enzymes catalyze the breakdown of cellulose in plant biomass into simple sugars and then to bioethanol. In the absence of the enzyme, the cellulosic biomass is recalcitrant to decomposition due to fermentation-resistant lignin and pectin coatings on the cellulose surface, which make them inaccessible for hydrolysis. Cellobiohydrolase CelS is a microbial enzyme that binds to cellulose fiber and efficiently cleaves it into a simple sugar (cellobiose) by a repeated processive chopping mechanism. The two contributing factors to the catalytic reaction rate and the yield of cellobiose are the efficient product expulsion from the product binding site of CelS and the movement of the substrate or cellulose chain into the active site. Despite progress in understanding product expulsion in other cellulases, much remains to be understood about the molecular mechanism of processive action of these enzymes. Here, nonequilibrium molecular dynamics simulations using suitable reaction coordinates are carried out to investigate the energetics and mechanism of the substrate dynamics and product expulsion in CelS. The calculated free energy barrier for the product expulsion is three times lower than that for the processive action indicating that product removal is relatively easier and faster than the sliding of the substrate to the catalytic active site. The water traffic near the active site in response to the product expulsion and the processive action is also explored.

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 of Sugar Ring Contraction and Closure Catalyzed by UDP-d-apiose/UDP-d-xylose Synthase (UAXS)

Uridine diphosphate (UDP)-apiose/UDP-xylose synthase (UAXS) is a member of the short-chain dehydrogenase/reductase superfamily (SDR), which catalyzes the ring contraction and closure of UDP-d-glucuronic acid (UDP-GlcA), affording UDP-apiose and UDP-xylose. UAXS is a special enzyme that integrates ring-opening, decarboxylation, rearrangement, and ring closure/contraction in a single active site. Recently, the ternary complex structure of UAXS was crystallized from Arabidopsis thaliana. In this work, to gain insights into the detailed formation mechanism of UDP-apiose and UDP-xylose, an enzyme-substrate reactant model has been constructed and quantum mechanical/molecular mechanical (QM/MM) calculations have been performed. Our calculation results reveal that the reaction starts from the C4-OH oxidation, which is accompanied by the conformational transformation of the sugar ring from chair type to boat type. The sugar ring-opening is prior to decarboxylation, and the deprotonation of the C2-OH group is the prerequisite for sugar ring-opening. Moreover, the keto-enol tautomerization of the decarboxylated intermediate is a necessary step for ring closure/contraction. Based on our calculation results, more UDP-apiose product was expected, which is in line with the experimental observation. Three titratable residues, Tyr185, Cys100, and Cys140, steer the reaction by proton transfer from or to UDP-GlcA, and Arg182, Glu141, and D337 constitute a proton conduit for sugar C2-OH deprotonation. Although Thr139 and Tyr105 are not directly involved in the enzymatic reaction, they are responsible for promoting the catalysis by forming hydrogen-bonding interactions with GlcA. Our calculations may provide useful information for understanding the catalysis of the SDR family.

Enhancing the catalytic activity of a GH5 processive endoglucanase from Bacillus subtilis BS-5 by site-directed mutagenesis

Processive endoglucanases possess both endo- and exoglucanase activity, making them attractive discovery and engineering targets. Here, a processive endoglucanase EG5C-1 from Bacillus subtilis was employed as the starting point for enzyme engineering. Referring to the complex structure information of EG5C-1 and cellohexaose, the amino acid residues in the active site architecture were identified and subjected to alanine scanning mutagenesis. The residues were chosen for a saturation mutagenesis since their variants showed similar activities to EG5C-1. Variants D70Q and S235W showed increased activity towards the substrates CMC and Avicel, an increase was further enhanced in D70Q/S235W double mutant, which displayed a 2.1- and 1.7-fold improvement in the hydrolytic activity towards CMC and Avicel, respectively. In addition, kinetic measurements showed that double mutant had higher substrate affinity (K_m) and a significantly higher catalytic efficiency (k_cat/K_m). The binding isotherms of wild-type EG5C-1 and double mutant D70Q/S235W suggested that the binding capability of EG5C-1 for the insoluble substrate was weaker than that of D70Q/S235W. Molecular dynamics simulations suggested that the collaborative substitutions of D70Q and S235W altered the hydrogen bonding network within the active site architecture and introduced new hydrogen bonds between the enzyme and cellohexaose, thus enhancing both substrate affinity and catalytic efficiency.

Binding Process and Free Energy Characteristics of Cellulose Chain into the Catalytic Domain of Cellobiohydrolase TrCel7A

It was observed in experiments that the catalytic domain (CD) of Trichoderma reesei Cel7A (TrCel7A) hydrolyzes crystalline cellulose in a processive manner, but the underlying binding mechanism is still unknown. Here, through replica-exchange molecular dynamics simulations, we find that the loading and sucking-in process of the cellulose chain into CD is entropy-driven and enthalpy-unfavorable, which firmly relate to the desolvation of the binding channel of CD. During the loading process, hydrophobic interactions play a dominant role because several aromatic residues have been identified to guide the cellulose chain processing. At the active site, a transition from enthalpy- to entropy-driven is detected for the driving force. Such a finding reveals the indispensability of the catalytic reaction of the glycosidic bond to provide the energy to drive the movements of the cellulose chain. Our study reveals the interaction pictures between the cellulose chain and TrCel7A at the atomic level, which helps better understand the catalytic mechanism of TrCel7A.

Substrate binding versus escape dynamics in a pH-affected fungal beta-glucosidase revealed by molecular dynamics simulations

The cellulolytic ability of fungal species is important to both natural and engineered biocycling of plant matter. One essential step is the conversion of cellobiose into glucose catalyzed by beta-glucosidases. Mutagenesis studies have implicated altering the substrate binding pocket to influence the pH-activity profile of this enzyme. However, structural understanding of the pH-affected substrate binding environment is lacking. Here we conducted molecular dynamics simulations of fully hydrated TrBgl2, a beta-glucosidase of Trichoderma reesei, equilibrated at its optimal pH (pH 6) and two unfavorable pHs (pH 5 and pH 7.5). We identified structural arrangement of specific residues that facilitated substrate escape from the catalytic site at pH 5 but locked the bound substrate in an unfavorable orientation at pH 7.5. For comparative analysis, we also performed simulations of a mutated TrBgl2 with previously demonstrated improved catalysis as a function of pH. We captured the responsible conformational changes in the engineered substrate binding pocket.

Effects of Tyr555 and Trp678 on the processivity of cellobiohydrolase A from Ruminiclostridium thermocellum: A simulation study

Cellobiohydrolase A from Ruminiclostridium thermocellum (Cbh9A) is a processive exoglucanase from family 9 and is an important cellobiohydrolase that hydrolyzes cello-oligosaccharide into cellobiose. Residues Tyr555 and Trp678 considerably affect catalytic activity, but their mechanisms are still unknown. To investigate how the Tyr555 and Trp678 affect the processivity of Cbh9A, conventional molecular dynamics, steered molecular dynamics, and free energy calculation were performed to simulate the processive process of wild type (WT)-Cbh9A, Y555S mutant, and W678G mutant. Analysis of simulation results suggests that the binding free energies between the substrate and WT-Cbh9A are lower than those of Y555S and W678G mutants. The pull forces and energy barrier in Y555S and W678G mutants also reduced significantly during the steered molecular dynamics (SMD) simulation compared with that of the WT-Cbh9A. And the potential mean force calculations showed that the pulling energy barrier of Y555S and W678G mutants is much lower than that of WT-Cbh9A.

Molecular dynamics investigations of cello-oligosaccharide recognition by Cel9G-CBM3c from Clostridium cellulovorans

The processive mechanism of cellulases against cellulose represents one of the key mechanisms in the conversion of biomass. A reliable model of substrate binding in a multidomain cellulase is a prerequisite for fully understanding this mechanism. In this study, the specificity of the recognition of the polysaccharide by the multidomain endoglucanase Cel9G from Clostridium cellulovorans was investigated by molecular dynamics simulations. Aromatic ring-containing residues were found to be critical for stabilizing the substrate. The calculated subtotal contributions of polar residues close to the active site, e.g., D58, E244, R315 and D420, also have some critical functions in substrate binding. Unlike other members of the carbohydrate-binding module family, CBM3c alone is shown not to bind cellulose very well, which is also consistent with experimental conclusions.

Microwave-assisted ionic liquid-mediated rapid catalytic conversion of non-edible lignocellulosic Sunn hemp fibres to biofuels

Sunn hemp fibre - a cellulose-rich crystalline non-food energy crop, containing 75.6% cellulose, 10.05% hemicellulose, 10.32% lignin, with high crystallinity (80.17%) and degree of polymerization (650) - is identified as a new non-food substrate for lignocellulosic biofuel production. Microwave irradiation is employed to rapidly rupture the cellulose's glycosidic bonds and enhance glucose yield to 78.7% at 160 °C in only 46 min. The reactants - long-chain cellulose, ionic liquid, transition metal catalyst, and water - form a polar supramolecular complex that rotates under the microwave's alternating polarity and rapidly dissipates the electromagnetic energy through molecular collisions, thus accelerating glycosidic bond breakage. In 46 min, 1 kg of Sunn hemp fibres containing 756 g of cellulose produces 595 g of glucose at 160 °C, and 203 g of hydroxymethyl furfural (furanic biofuel precursor) at 180 °C. Yeast mediated glucose fermentation produces 75.6% bioethanol yield at 30 °C, and the ionic liquid is recycled for cost-effectiveness.

Enhancing the catalytic activity of a novel GH5 cellulase GtCel5 from Gloeophyllum trabeum CBS 900.73 by site-directed mutagenesis on loop 6

BACKGROUND

Cellulases of glycosyl hydrolase (GH) family 5 share a (β/α)₈ TIM-barrel fold structure with eight βα loops surrounding the catalytic pocket. These loops exposed on the surface play a vital role in protein functions, primarily due to the interactions of some key amino acids with solvent and ligand molecules. It has been reported that motions of these loops facilitate substrate access and product release, and loops 6 and 7 located at the substrate entrance of the binding pocket promote proton transfer reaction at the catalytic site motions. However, the role of these flexible loops in catalysis of GH5 cellulase remains to be explored.

RESULTS

In the present study, an acidic, mesophilic GH5 cellulase (with optimal activity at pH 4.0 and 70 °C), GtCel5, was identified in Gloeophyllum trabeum CBS 900.73. The specific activities of GtCel5 toward CMC-Na, barley β-glucan, and lichenan were 1117 ± 43, 6257 ± 26 and 5318 ± 54 U/mg, respectively. Multiple sequence alignment indicates that one amino acid residue at position 233 on the loop 6 shows semi-conservativeness and might contribute to the great catalytic performance. Saturation mutagenesis at position 233 was then conducted to reveal the vital roles of this position in enzyme properties. In comparison to the wild type, variants N233A and N233G showed decreased optimal temperature (- 10 °C) but increased activities (27 and 70%) and catalytic efficiencies (k_cat/K_m; 45 and 52%), respectively. The similar roles of position 233 in catalytic performance were also verified in the other two GH5 homologs, TeEgl5A and PoCel5, by reverse mutation. Further molecular dynamics simulations suggested that the substitution of asparagine with alanine or glycine may introduce more hydrogen bonds, increase the flexibility of loop 6, enhance the interactions between enzyme and substrate, and thus improve the substrate affinity and catalytic efficiency.

CONCLUSION

This study proposed a novel cellulase with potentials for industrial application. A specific position was identified to play key roles in cellulase-substrate interactions and enzyme catalysis. It is of great importance for understanding the binding mechanism of GH5 cellulases, and provides an effective strategy to improve the catalytic performance of cellulases.

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.

Hemolytic activity of venom from crown-of-thorns starfish Acanthaster planci spines

Background

The crown-of-thorns starfish Acanthaster planci is a venomous species from Taiwan whose venom provokes strong hemolytic activity. To understand the hemolytic properties of A. planci venom, samples were collected from A. planci spines in the Penghu Islands, dialyzed with distilled water, and lyophilized into A. planci spine venom (ASV) powder.

Results

Both crude venom and ASV cause 50% hemolysis at a concentration of 20 μg/mL. The highest hemolytic activity of ASV was measured at pH 7.0-7.4; ASV-dependent hemolysis was sharply reduced when the pH was lower than 3 or greater than 8. There was almost no hemolytic activity when the Cu²⁺ concentration was increased to 10 mM. Furthermore, incubation at 100°C for 30 to 60 minutes sharply decreased the hemolytic activity of ASV. After treatment with the protease α-chymotrypsin, the glycoside hydrolase cellulase, and the membrane component cholesterin, the hemolytic activity of ASV was significantly inhibited.

Conclusions

The results of this study provide fundamental information about A. planci spine venom. The hemolytic activity was affected by pH, temperature, metal ions, EDTA, cholesterin, proteases, and glycoside hydrolases. ASV hemolysis was inhibited by Cu²⁺, cholesterin, α-chymotrypsin, and cellulose, factors that might prevent the hemolytic activity of venom and provide the medical treatment for sting.

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.