If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
To whom correspondence may be addressed: Dept. of Chemical and Materials Engineering, University of Kentucky, 177 F. Paul Anderson Tower, Lexington, KY 40506-0046. Tel.: 703-292-2895;
Department of Chemical and Materials Engineering, University of Kentucky, Lexington, Kentucky 40506-0046Department of Molecular Sciences, Swedish University of Agricultural Sciences, Box 7015, 750 07 Uppsala, Sweden
β-Glucosidases enhance enzymatic biomass conversion by relieving cellobiose inhibition of endoglucanases and cellobiohydrolases. However, the susceptibility of these enzymes to inhibition and transglycosylation at high glucose or cellobiose concentrations severely limits their activity and, consequently, the overall efficiency of enzyme mixtures. We determined the impact of these two processes on the hydrolytic activity of the industrially relevant family 3 β-glucosidases from Hypocrea jecorina, HjCel3A and HjCel3B, and investigated the underlying molecular mechanisms through kinetic studies, binding free energy calculations, and molecular dynamics (MD) simulations. HjCel3B had a 7-fold higher specificity for cellobiose than HjCel3A but greater tendency for glucose inhibition. Energy decomposition analysis indicated that cellobiose has relatively weak electrostatic interactions with binding site residues, allowing it to be easily displaced by glucose and free to inhibit other hydrolytic enzymes. HjCel3A is, thus, preferable as an industrial β-glucosidase despite its lower activity caused by transglycosylation. This competing pathway to hydrolysis arises from binding of glucose or cellobiose at the product site after formation of the glycosyl-enzyme intermediate. MD simulations revealed that binding is facilitated by hydrophobic interactions with Trp-37, Phe-260, and Tyr-443. Targeting these aromatic residues for mutation to reduce substrate affinity at the product site would therefore potentially mitigate transglycosidic activity. Engineering improved variants of HjCel3A and other structurally similar β-glucosidases would have a significant economic effect on enzymatic biomass conversion in terms of yield and production cost as the process can be consequently conducted at higher substrate loadings.
Plant biomass is a rich, sustainable carbon source that presents an attractive alternative to fossil reserves for the generation of liquid fuel products (
). Current industrial practices for converting biomass to small, fermentable sugar precursors often include an enzymatic hydrolysis process that utilizes the secreted hydrolytic enzymes of filamentous fungi, most commonly Hypocrea jecorina (teleomorph of Trichoderma reesei) (
). The fungal endoglucanases, lytic polysaccharide monooxygenases, processive cellobiohydrolases, and β-glucosidases work synergistically to decompose cellulose, the most abundant polysaccharide of plant biomass, to glucose (
). β-Glucosidases play a particularly important role in maintaining efficient enzyme mixture turnover as they are responsible for relieving endoglucanase and cellobiohydrolase product inhibition through hydrolysis of cellobiose to glucose (
free energy perturbation with Hamiltonian replica-exchange MD
Aa
A. aculeatus.
enzymes are among the best cellobiose-hydrolyzing enzymes, and accordingly, they comprise the majority of β-glucosidases in industrial enzyme mixtures (
). For instance, the optimum level of Cel3A from H. jecorina (HjCel3A) in cellulase mixtures enhanced the conversion of various cellulosic substrates by nearly 10% (
Cloning and amplification of the gene encoding an extracellular β-glucosidase from Trichoderma reesei: evidence for improved rates of saccharification of cellulosic substrates.
) exhibited broad substrate specificity, hydrolyzing not only cellobiose but also sophorose (β-1,2), laminaribiose (β-1,3), and gentiobiose (β-1,6). However, these GH3 β-glucosidases had the highest specificity for longer cello-oligomers (cellotriose and cellotetraose).
The solution of the first HjCel3A crystal structure revealed new insights into catalytic activity and substrate affinity (
). Namely, a network of hydrogen bonds, originating from surrounding charged residues, stabilizes a glucose molecule in the −1 site (substrate site). Aromatic residues flank the +1 site (product site) where the reducing end of cellobiose would be situated. Additionally, the ability of HjCel3A to catalyze hydrolysis of longer cello-oligomers points toward the existence of a putative +2 product site.
Despite these very important findings, many of the fundamental molecular mechanisms of HjCel3A activity, along with that of related GH3s, have yet to be defined. Decomposition of polysaccharides by GH3 β-glucosidases is proposed to follow the retaining mechanism (see Fig. 1). In the initial glycosylation step, the glycosidic oxygen is protonated by the acid/base catalyst, Glu, whereas the nucleophile, Asp, attacks the anomeric carbon of the glycone moiety to form a glycosyl-enzyme intermediate (GEI). In the subsequent hydrolysis reaction, the deprotonated Glu activates a water molecule that carries out another nucleophilic attack, yielding a product with the same stereochemistry as the substrate (
Figure 1.Proposed retaining mechanism of GH3 β-glucosidases. The acid/base catalyst (Glu) protonates the glycosidic oxygen, and the nucleophile (Asp) attacks the anomeric carbon to form the glycosyl-enzyme intermediate in the glycosylation step. Two competing pathways subsequently occur: (a) hydrolysis, in which water attacks the intermediate to release the sugar, and (b) transglycosylation, in which a second sugar molecule acts as an acceptor to form a new β-glucoside.
The activity of β-glucosidases is seriously impeded at high cellobiose or glucose concentration because of their susceptibility to product and/or substrate inhibition and transglycosylation (
). Transglycosylation is a competing catalytic pathway to hydrolysis wherein glucose or another cellobiose molecule, instead of water, attacks the GEI to form a di- or trisaccharide (Fig. 1). Accumulation of cellobiose and glucose inevitably occurs under high substrate loadings, a necessary condition to increase the yield and lower the production costs of industrial biomass conversion (
Trp-49 of the family 3 β-glucosidase from Aspergillus niger is important for its transglucosidic activity: creation of novel β-glucosidases with low transglucosidic efficiencies.
Transglucosidic reactions of the Aspergillus niger family 3 β-glucosidase: qualitative and quantitative analyses and evidence that the transglucosidic rate is independent of pH.
Trp-262 is a key residue for the hydrolytic and transglucosidic reactivity of the Aspergillus niger family 3 β-glucosidase: substitution results in enzymes with mainly transglucosidic activity.
Kinetics of substrate transglycosylation by glycoside hydrolase family 3 glucan (1→3)-β-glucosidase from the white-rot fungus Phanerochaete chrysosporium.
) that allow quantitative determination of the effects of these two processes on hydrolytic activity. One of the notable findings is that transglycosylation is the main cause of decay in activity at high cellobiose concentration but not at high glucose concentration (
The critical role of β-glucosidases in the overall biomass conversion efficiency provides the impetus to engineer variants that overcome limitations arising from environmental process conditions. Protein engineering requires knowledge of the molecular factors underlying inhibition and transglycosylation. An analysis of published Michaelis (Km) and inhibition (Ki) constants for cellobiose and glucose, respectively, suggests positive correlation between affinity for cellobiose and strength of glucose inhibition (
). Modulation of the transglycosylation activity of the A. niger β-glucosidase upon mutation of tryptophan residues was linked to affinity of glucose for the −1 and +1 sites (
Trp-49 of the family 3 β-glucosidase from Aspergillus niger is important for its transglucosidic activity: creation of novel β-glucosidases with low transglucosidic efficiencies.
Trp-262 is a key residue for the hydrolytic and transglucosidic reactivity of the Aspergillus niger family 3 β-glucosidase: substitution results in enzymes with mainly transglucosidic activity.
). The results of these studies imply that reducing affinity at the −1 and +1 sites of β-glucosidases will result in lower product inhibition and decrease the incidence of transglycosylation, respectively.
We investigated this hypothesis on HjCel3A and HjCel3B, which exhibit the highest activity against cellobiose among the GH3 β-glucosidases in H. jecorina (
). Kinetic studies and binding free energy calculations were performed to determine which process, inhibition or transglycosylation, has a more significant impact on hydrolytic activity. Protein residues responsible for cellobiose and glucose affinity at the different binding sites were then identified using molecular dynamics (MD) simulations.
Results
Physical properties and model structure of HjCel3B
SDS-PAGE analysis under nonreducing and nondenaturing conditions suggested that HjCel3B is a dimer with a molecular mass of 100 kDa. In contrast, HjCel3A is a monomer with a molecular mass of 60 kDa (
). Both enzymes have the same pH optimum of 5.0 at 37 °C (Fig. S1).
The HjCel3B homology model reveals an overall fold similar to HjCel3A, consisting of 1) a triose-phosphate isomerase (TIM) barrel-like domain, 2) α/β sandwich domain, and 3) fibronectin type III (FnIII) domain (Fig. 2A). However, HjCel3A lacks the insertion region in the FnIII domain. The active site is located at the interface of domains 1 and 2. Glu-516 (domain 2) acts as the acid/base catalyst, and Asp-287 (domain 1) acts as the nucleophile. A glucose substrate is shown bound at the −1 site (Fig. 2B). It is held in place by hydrogen bond interactions with Asp-99, Arg-163, Lys-196, and His-197 (Asp-61, Arg-125, Lys-158, and His-159 in HjCel3A). Glucose also has a hydrophobic interaction with Trp-288 (Trp-237 in HjCel3A) at this site.
Figure 2.A, cartoon representation of the homodimeric HjCel3B homology model. The TIM barrel–like domain is shown in blue, α/β sandwich domain is in purple, FnIII domain is in green, linker regions are in yellow, and loop-insertion region is in red. The acid/base catalyst (Glu-516) and nucleophile (Asp-287) are shown as spheres. B, stick representation of the HjCel3B active site with bound glucose (tan carbon atoms) in the −1 site. Charged residues forming hydrogen bonds with the substrate are shown with gray carbon atoms; interactions are represented as dashed lines. C and D, substrate pockets of HjCel3A and HjCel3B, respectively, defined by gray surface with residues shown in stick representation (green carbon atoms) and glucose in ball-and-stick representation (cyan carbon atoms). In HjCel3A, Trp-37 and Tyr-443 form the +1 site, and Tyr-68 and Phe-260 form the putative +2 site. In HjCel3B, Trp-75, Tyr-518, and Phe-312 form the +1 site, and Gln-208, Gln-289, and Glu-211 form the putative +2 site. Gln-208 and Gln-289 correspond to Glu-170 and Asn-238 in HjCel3A. Glu-211 has no equivalent as it is located in an α-helix (highlighted in yellow) that is absent in HjCel3A. Phe-106 (corresponding to Tyr-68 in HjCel3A) and Trp-366 at the HjCel3B substrate pocket entrance are also shown.
Aromatic residues Trp-75 and Tyr-518 (Trp-37 and Tyr-443 in HjCel3A) are central to the +1 product site. Phe-312 is also included in this site, but the corresponding residue in HjCel3A, Phe-260, is oriented away from the substrate pocket. Based on the predicted shape of the substrate pocket (Fig. 2C), the putative +2 site of HjCel3A is formed by Phe-260 and Tyr-68. Although the two aromatic residues at the entrance of the HjCel3B substrate pocket, Phe-106 and Trp-366, may form the corresponding site, the predicted +2 site instead consists of Gln-208, Glu-211, and Gln-289 (Fig. 2D). Gln-208 and Gln-289 correspond to Glu-170 and Asn-238 in HjCel3A, whereas Glu-211 is located in an α-helix in domain 1 that is missing in HjCel3A. The calculated solvent-excluded volume of the HjCel3B substrate pocket is smaller (431 Å3) than that of HjCel3A (519 Å3).
2-Chloro-4-nitrophenyl-β-glucoside (CNPG), 2-chloro-4-nitrophenyl-β-xyloside (CNPX), and cellobiose kinetics
Fig. 3A shows the plot of initial velocity as a function of the concentration of the chromogenic model substrate CNPG. For HjCel3A, there is a hyperbolic increase in initial velocity at low substrate concentrations followed by a decrease from about 1 mm CNPG. In contrast, no decrease in initial velocity is observed with HjCel3B within the same concentration range. The same trend was observed with cellobiose as substrate. The two β-glucosidases are also active on CNPX, which is not unexpected considering the prevalence of xylosidase activity in structurally similar GHs (
Km and kcat values for the hydrolysis of CNPG by HjCel3A (Table 1) were determined at low substrate concentrations where inhibition and/or transglycosylation is minimal and the kinetics obey the Michaelis–Menten equation. As the kinetics of CNPX hydrolysis did not follow the Michaelis–Menten equation, the parameters were not calculated. HjCel3A has a higher specificity for CNPG than cellobiose. Interestingly, HjCel3B not only has a higher specificity for cellobiose but also hydrolyzes this substrate more efficiently than HjCel3A. This is contrary to the comparative study of HjCel3A and HjCel3B expressed in A. oryzae, showing that the latter has a slightly lower kcat/Km (21.6 against 37.3 mm−1 s−1) (
The inhibition constants (Ki) determined from the Lineweaver–Burk plot show that, for both HjCel3A and HjCel3B, glucose inhibition increases as the pH increases (Table 2). HjCel3B is more inhibited than HjCel3A at acidic pH, although the Ki values of the two β-glucosidases are similar at higher pH.
Table 2Inhibition constants (Ki; mm) for glucose inhibition of HjCel3A and HjCel3B at different pH values
Identification of transglycosylation products with cellobiose as substrate
A variety of transglycosylation products were obtained for both HjCel3A and HjCel3B (Fig. 4). These include disaccharides (gentiobiose, laminaribiose, and sophorose), formed with the glucose product as acceptor, and cellotriose, formed with the substrate as acceptor. Additional transglycosylation products were formed, likely longer glucosides with mixed β-linkages; however, we were unable to identify these products without appropriate comparative standards.
Figure 4.Transglycosylation products obtained with cellobiose as substrate and acceptor (pH 5. 0 and 5.7 and 37 °C).A, glucose; B, cellobiose; C, cellotriose; D, cellotetraose; E, cellopentaose; F, cellohexaose; G, gentiobiose; H, laminaribiose; I, sophorose. The chromatogram was obtained after removal of glucose and concentration of the products. The units of response are given in nanocoulombs (nC).
Binding free energies from free energy perturbation with Hamiltonian replica-exchange MD (FEP/λ-REMD)
The free energies for binding glucose and cellobiose in the active site models (Fig. 5) of HjCel3A and HjCel3B are summarized in Table 3. Glucose has similar affinity for both reactive forms of HjCel3A and HjCel3B (Model I). Comparison of individual free energy contributions when the substrate is in the binding site and in bulk solution shows that electrostatic interactions provide the most favorable contribution (−18 kcal/mol) in each case. Cellobiose (Model II) has a much lower affinity for HjCel3B dominated by dispersion (−18 kcal/mol) in contrast to HjCel3A where the binding energy is largely electrostatic (−28 kcal/mol). The energy for glucose binding in the product site (Model III) is essentially zero in both enzymes. Cellobiose binds more stably in the product site when the O4 atom of the nonreducing end is positioned for attack by Glu (Model IV) with comparable energies for HjCel3A and HjCel3B. The binding of both compounds to the GEI is mainly facilitated by dispersion interactions, which contribute between −8 and −14 kcal/mol to the binding energy.
Figure 5.Enzyme–substrate complexes modeled using MD simulations. The nucleophile is Asp (Asp-236 and Asp-287 for HjCel3A and HjCel3B, respectively), and the acid/base catalyst is Glu (Glu-441 and Glu-516 for HjCel3A and HjCel3B, respectively). Binding sites are labeled −1 and +1 according to standard nomenclature (
). Models I and II represent the reactive form of the enzyme (protonated Glu) with bound glucose and cellobiose, respectively. Models III and IV(V) represent the glycosyl-enzyme intermediate (glucose covalently bound to Asp) with bound glucose and cellobiose, respectively. The ligand inserted into the enzyme in the FEP/λ-REMD calculations is shown in blue.
Hydrogen bond occupancies are listed in Tables S1–S7, and van der Waals interaction energies are listed in Tables S8–S10. The active site residues of HjCel3A and HjCel3B form several hydrogen bonds with the glucose substrate in the −1 site (Table S1). In HjCel3A, the residues involved are Asp-61, Lys-158, His-159, and the nucleophile Asp-236 (Fig. 6A). There is also a hydrogen bond with Arg-67, which was not observed in the crystal structure (Protein Data Bank (PDB) code 3ZYZ (
)). Glucose interacts with the corresponding residues in HjCel3B with the exception of the arginine residue; however, an additional interaction with Arg-163 is formed (Fig. 6B). Among the hydrophobic residues, only the tryptophan adjacent to the nucleophile has significant van der Waals interaction with glucose (−5.0 kcal/mol with Trp-237 in HjCel3A and −3.5 kcal/mol with Trp-288 in HjCel3B). There are also relatively weak interactions with Leu-110 and Met-201 in HjCel3A (Table S8). Glucose is oriented almost perpendicular to Trp-288 in HjCel3B, and its O6 atom also briefly formed a hydrogen bond with the indole side chain (Table S1).
Figure 6.Snapshots from MD simulations of β-glucosidases with glucose in the −1 site (Model I) or cellobiose across the −1 and +1 sites (Model II). Glu-441 (Glu-516) is protonated in the models. Hydrogen atoms are omitted for clarity. Substrate pocket residues that formed hydrogen bond and van der Waals interactions with the ligand during the simulation are also shown, where dashed lines illustrate the interactions. A, HjCel3A–glucose, 90 ns. B, HjCel3B–glucose, 60 ns. C, HjCel3A–cellobiose, 60 ns. D, HjCel3B–cellobiose, 60 ns.
The nonreducing end of cellobiose forms hydrogen bonds with many of the same residues as glucose in HjCel3A (HjCel3B) (Table S2). The reducing end, however, interacts with the side chains of Arg-67 and Asp-236, and the backbone oxygen atom of Trp-237 in HjCel3A (Fig. 6C). The glycosidic oxygen also formed a hydrogen bond with Glu-441, albeit a short-lived one (Table S2). In HjCel3B, the reducing end only forms hydrogen bonds with Arg-105 and Arg-207 (Fig. 6D). The van der Waals interaction of cellobiose with Trp-288 in HjCel3B (−6.9 kcal/mol) is slightly stronger compared with that of cellobiose with Trp-237 in HjCel3A (−4.9 kcal/mol) (Table S8). This substrate additionally interacts with Tyr-443 (Tyr-518) in HjCel3A (HjCel3B) with an energy of about −3 kcal/mol. The van der Waals interactions with Trp-37 (Trp-75) and Phe-260 (Phe-312) are negligible.
In HjCel3A– and HjCel3B–GEI, hydrogen bonding with the acid/base catalyst (Glu-441 and Glu-516) positions glucose in the +1 site slightly above the covalently bound glucose instead of next to it where the reducing end of cellobiose would be (Fig. 7). In addition to Glu-441, hydrogen bonds are also formed with Arg-67, Asp-370, and, briefly, Arg-169 and Tyr-204 in HjCel3A (Fig. 7A and Table S3). The van der Waals interactions with +1 site aromatic residues are negligible except with Tyr-443 (−3 kcal/mol; Table S9). Glucose has fewer hydrogen bond interactions in HjCel3B–GEI; the only residues forming hydrogen bonds with greater than 10% occupancy are Glu-516, Arg-207, and Tyr-255 (Fig. 7B and Table S3). Unlike in HjCel3A, glucose is located almost centrally between Trp-288 and Tyr-518, resulting in comparable van der Waals interaction energies (−3 kcal/mol; Table S9).
Figure 7.Snapshots from MD simulations of glycosyl-enzyme intermediate with glucose in the product site (Model III). Glu-441 (Glu-516) is unprotonated in the models. Hydrogen atoms are omitted for clarity. Substrate pocket residues that formed hydrogen bond and van der Waals interactions with the ligand during the simulation are also shown, where dashed lines illustrate the interactions. A, HjCel3A, 65 ns. B, HjCel3B, 60 ns.
Only the results of the stable binding orientation of cellobiose (Model IV) in the substrate pocket of the GEI will be discussed here in detail. Results for Model V simulations are included in Tables S6, S7, and S10 and Fig. S2. Cellobiose is initially above the covalently bound glucose (similar to glucose in Model III) but becomes very mobile during the simulations (Fig. S3). The hydrogen bond between cellobiose and Glu-441 of HjCel3A (Fig. 8A) lasts for only ∼55 ns as Glu-441 eventually forms hydrogen bonds with Arg-169, Tyr-443, and Ser-384 instead of with cellobiose. Consequently, the substrate moves farther from the covalently bound glucose and forms interactions with Glu-170 and Asn-238 (Fig. 8B and Table S4). The strongest van der Waals interactions are with Tyr-443 and Phe-260 (3–4 kcal/mol) (Table S9). Interaction with Trp-37 at the +1 site is rather short-lived, whereas that with Tyr-68 at the putative +2 is negligible. In HjCel3B–GEI, cellobiose initially forms hydrogen bonds with the side chains of Glu-516, Asp-287, Arg-207, and Glu-211 and the backbone oxygen atom of Trp-288 (Fig. 8C and Table S5). Hydrogen bonds with Gln-208 and Gln-289 in the putative +2 site have occupancies less than 10%. The substrate flips, and the reducing end eventually forms hydrogen bonds with Glu-516, Asp-443, and Ser-457 after ∼130 ns (Fig. 8D). Only Tyr-518 forms significant van der Waals interaction with cellobiose (−3 kcal/mol; Table S9).
Figure 8.Snapshots from MD simulations of glycosyl-enzyme intermediate with cellobiose in the product sites (Model IV). Glu-441 (Glu-516) is unprotonated in the models. Hydrogen atoms are omitted for clarity. Substrate pocket residues that formed hydrogen bond and van der Waals interactions with the ligand during the simulation are also shown, where dashed lines illustrate the interactions. A, HjCel3A, 35 ns. The nonreducing end of cellobiose is hydrogen-bonded to Glu-441. B, HjCel3A, 70 ns. The hydrogen bond with Glu-441 is broken, and cellobiose moves farther to interact with Glu-170, Asn-238, and Phe-260. C, HjCel3B, 60 ns. The nonreducing end of cellobiose is hydrogen-bonded to Glu-516. D, HjCel3B, 140 ns. Cellobiose flips during the simulation, and the reducing end now forms hydrogen bonds with Glu-516 as well as Asp-443 and Ser-457.
Optimizing the hydrolytic activity of β-glucosidases is critical to enhancing enzymatic biomass conversion efficiency. At high substrate loadings, product inhibition and transglycosylation tend to be the main limitations to efficient catalytic turnover and can be linked to the affinity of cellobiose and glucose for the product site(s). Thus, kinetic data of industrially important GH3 β-glucosidases from H. jecorina, HjCel3A and HjCel3B, were interpreted with the aid of binding free energy calculations and MD simulations. Decomposition of the binding free energy provides insight on the nature of the interaction governing substrate affinity. The residues involved can then be identified from the MD simulations and targeted for mutation to reduce inhibition and transglycosylation.
HjCel3A is susceptible to substrate inhibition and/or transglycosylation as evidenced by the decrease in initial velocity at both high CNPG and cellobiose concentrations; in contrast, the activity of HjCel3B did not decay within the same concentration range. Substrate inhibition arises from binding to a nonproductive position (i.e. +1 and +2 sites) instead of across the −1 and +1 sites. Affinity at the +1 site can be estimated from the binding energy difference between HjCel3A Models I and II (∼13 kcal/mol). The existence of a +2 site is implied by the reported activity of HjCel3A toward cellotriose and cellotetraose (
). MD simulations of Models IV and V suggested that this site is the pocket formed by Glu-170 and Asn-238 (instead of the predicted pocket formed by Tyr-68 and Phe-260), and its affinity is relatively low (∼−3 kcal/mol). Binding of cellobiose across the −1 and +1 sites, with a free energy of −21 kcal/mol, would, therefore, be more favorable, rendering substrate inhibition less likely.
Moreover, a previous study on various GH3 β-glucosidases demonstrated through quantitative determination of products that reduced glucose production at high substrate concentrations can be attributed to transglycosylation (
Transglucosidic reactions of the Aspergillus niger family 3 β-glucosidase: qualitative and quantitative analyses and evidence that the transglucosidic rate is independent of pH.
). Cellobiose, either from solution or the original cello-oligomer substrate, binds to the product site following formation of the GEI to yield longer glucosides of various β-linkages. In the HjCel3A Model IV simulation, the O4 atom of the nonreducing end of cellobiose formed a hydrogen bond with Glu-441, positioning it for formation of the experimentally observed cellotriose product. In comparison, the main trisaccharide product obtained for the A. niger β-glucosidase was 6-O-glucosyl-cellobiose (
Transglucosidic reactions of the Aspergillus niger family 3 β-glucosidase: qualitative and quantitative analyses and evidence that the transglucosidic rate is independent of pH.
Free energy calculations for HjCel3A Model IV showed that binding is favorable mainly due to dispersion effects. The nonreducing end of cellobiose interacted with Tyr-443, whereas the reducing end interacted with Phe-260. Although cellobiose was still hydrogen-bonded to Glu-441 in the early stages of simulation, there was also hydrophobic interaction between the reducing end and Trp-37. A study on the A. niger β-glucosidase suggested that this residue might also be critical to transglycosylation activity. Mutation of the corresponding residue, Trp-49, to Gly, Ala, Asp, or Asn resulted in lower efficiency of the transglycosylation pathway due to higher Km. However, the authors suggested that Trp-49 only has an indirect conformational effect on cellobiose binding on the basis of the similarity of the Km values regardless of the nature of the substituted residue (
Trp-49 of the family 3 β-glucosidase from Aspergillus niger is important for its transglucosidic activity: creation of novel β-glucosidases with low transglucosidic efficiencies.
Mutation of Trp-262 (adjacent to the nucleophile Asp-261) of this A. niger β-glucosidase led to preferential transglycosidic activity, which was attributed to a lower hydrolysis rate constant and decreased affinity at the −1 site (whereas product sites were unaffected) (
Trp-262 is a key residue for the hydrolytic and transglucosidic reactivity of the Aspergillus niger family 3 β-glucosidase: substitution results in enzymes with mainly transglucosidic activity.
). The latter was inferred by comparing Ki values for inhibitors of increasing length (glucose, cellobiose, gentiobiose, and cellotriose). This hypothesis is consistent with MD simulations of HjCel3A Models I and II showing that the corresponding residue, Trp-237, forms a strong van der Waals interaction with either glucose or the nonreducing end of cellobiose at the −1 site. In the case of cellobiose, the O2 atom of the reducing end even formed a hydrogen bond with the backbone oxygen atom of Trp-237 toward the end of the simulation (from ∼175 ns).
Accumulation of glucose from hydrolysis can also cause inhibition as glucose competes with the substrate for binding to the −1 or +1 site. The pH dependence of Ki indicated that glucose preferentially binds at the −1 site where it forms hydrogen bonds with charged residues. Comparison of Km(CNPG) and Ki(glucose) values obtained at pH 5.7 using CNPG as substrate suggested that HjCel3A is more glucose-tolerant than HjCel3B. However, it is well known that Ki values depend on the substrate used (
). Although product inhibition studies with cellobiose were not performed, the relative affinities for cellobiose and glucose can be evaluated from the binding free energy calculations. Glucose bound at the −1 site of HjCel3A (as observed in the crystal structure) had a much more positive binding energy than cellobiose. The two substrates had comparable binding energies in HjCel3B; thus, glucose can compete with cellobiose for the binding site. The more positive binding energy of cellobiose in HjCel3B compared with that in HjCel3A may be attributed to the weaker affinity at the +1 site where electrostatic effects do not make a significant contribution unlike the case in HjCel3A. Thus, MD simulations provide further support that HjCel3B is more product-inhibited than HjCel3A.
The significant formation of disaccharide products suggested that transglycosylation, to some extent, relieves product inhibition in HjCel3A. Model III represents the case where glucose is retained after the glycosylation step. The binding free energy was essentially zero, even despite the strong affinity at the HjCel3A +1 site; this reduced affinity may be due to the glucose positioning for nucleophilic attack, which was no longer exactly in the +1 product site but slightly above the covalently bound glucose. This indicates that glucose is easily released following glycosidic bond cleavage, making it more likely that the acceptor for the transglycosylation reaction is another molecule that reentered from solution. Because protein–substrate interactions at or near the +1 site are primarily hydrophobic, glucose can freely adopt different orientations, which may have more favorable binding free energies than Models III–V. This would account for the formation of various β-linked disaccharides. During the simulation of the GEI–glucose complex of HjCel3A, Glu-441 eventually formed hydrogen bonds with the O2 and O3 atoms (instead of O4 as in the initial simulation coordinates), which would be consistent with formation of sophorose and laminaribiose, respectively. However, the major disaccharide product observed from previous studies of other GH3 β-glucosidases was gentiobiose. This was attributed to the higher reactivity of O6 compared with other hydroxyl groups (
Transglucosidic reactions of the Aspergillus niger family 3 β-glucosidase: qualitative and quantitative analyses and evidence that the transglucosidic rate is independent of pH.
Kinetics of substrate transglycosylation by glycoside hydrolase family 3 glucan (1→3)-β-glucosidase from the white-rot fungus Phanerochaete chrysosporium.
Although HjCel3A is more prone to transglycosylation than HjCel3B, the calculated binding free energies indicated that the GEIs of both β-glucosidases have similar affinities at the product site(s). Water availability within the first and second solvation shells of the GEI is also comparable (Table S11). The more accessible and larger substrate pocket of HjCel3A may account for the significant occurrence of transglycosylation as glucose or cellobiose from solution can easily enter and bind at the product site(s). Another explanation, suggested by the study on the A. niger β-glucosidase, is that transglycosylation activity depends on not only the product–site binding affinity but also the hydrolysis rate (
Trp-49 of the family 3 β-glucosidase from Aspergillus niger is important for its transglucosidic activity: creation of novel β-glucosidases with low transglucosidic efficiencies.
Transglucosidic reactions of the Aspergillus niger family 3 β-glucosidase: qualitative and quantitative analyses and evidence that the transglucosidic rate is independent of pH.
Trp-262 is a key residue for the hydrolytic and transglucosidic reactivity of the Aspergillus niger family 3 β-glucosidase: substitution results in enzymes with mainly transglucosidic activity.
Transglucosidic reactions of the Aspergillus niger family 3 β-glucosidase: qualitative and quantitative analyses and evidence that the transglucosidic rate is independent of pH.
). A constant-pH MD study of HjCel3A–GEI (with or without glucose in the +1 site) indicated that Glu-441 has a very low pKa (∼2) due to hydrogen bond networks with Arg-125 and Arg-169. The low basicity of Glu-441 would reduce its ability to deprotonate water, which has a higher deprotonation enthalpy than sugars (
). In contrast, HjCel3B homology modeling showed negatively charged residues (Glu-211, Glu-454, and Glu-525) situated within 10 Å of Glu-516 (Fig. S4). MD simulations of HjCel3B–GEI without cellobiose or glucose in the product site(s), moreover, showed that Glu-516 only forms relatively short-lived hydrogen bonds with neighboring residues (Tables S12 and S13 and Fig. S5). These conditions favor an elevated pKa (
), which would account for the predominance of hydrolysis over transglycosylation in HjCel3B.
Based solely on the specificity constant kcat/Km, HjCel3B is a better β-glucosidase than HjCel3A. However, a survey of various family 1 and 3 β-glucosidases showed that those with high kcat/Km (>105m−1 s−1) tend to have almost equal affinity for cellobiose and glucose (
) as demonstrated to be the case for HjCel3B by free energy calculations. This may prove to be detrimental to the efficiency of cellulase mixtures as cellobiose will eventually accumulate and inhibit endoglucanases and cellobiohydrolases. From this perspective, the more glucose-tolerant HjCel3A remains industrially preferable despite the lower kcat/Km. Transglycosylation is a major limiting factor for HjCel3A, but this can be mitigated by point mutations of residues responsible for cellobiose and glucose affinity at the product site, specifically Trp-37, Phe-260, and Tyr-443. This is another approach in addition to elevating the pKa of Glu-441 as we proposed in the recent constant-pH MD study (
). Finally, HjCel3A’s sequence similarity with xylosidic GH3s suggests that xylosidase activity can be conferred with relatively small modifications to the active site; this is advantageous in potentially reducing the number of enzyme components required in cellulase mixtures. Protein engineering of HjCel3A variants with reduced transglycosidic activity and broad substrate specificity, therefore, presents a preferential path forward for developing highly efficient industrial biomass conversion mixtures.
Experimental procedures
Reagents
Cellobiose was purchased from CU Chemie Uetikon (Lahr, Germany), glucose was from Janssen Pharmaceuticals (New Brunswick, NJ), laminaribiose was from Megazyme (Wicklow, Ireland), gentiobiose and sophorose were from Sigma, and cello-oligosaccharides were from Merck. CNPG and CNPX were synthesized as described previously (
Epoxyalkyl glycosides of d-xylose and xylo-oligosaccharides are active-site markers of xylanases from glycoside hydrolase family 11, not from family 10.
). The purified β-glucosidases HjCel3A and HjCel3B were a gift from DuPont Industrial Biosciences.
Enzyme purification
HjCel3A was further purified by affinity chromatography using a p-aminobenzyl-thio-β-glucopyranoside (coupled to Sepharose) affinity column as described previously (
). The column was equilibrated and washed with 100 mm acetate buffer (pH 5.0) containing 200 mm NaCl. The bound fraction was eluted with 100 mm glucose in 100 mm acetate buffer (pH 5.0). HjCel3B did not bind to the column and was used as received. Protein concentrations were determined according to the Bradford (
Kinetic studies of HjCel3B with CNPG and CNPX as substrates were carried out in 100 mm phosphate buffer (pH 5.7) in a microtiter plate. The substrate (190 μl) was added to 10 μl of enzyme at an appropriate concentration, and the release of 2-chloro-4-nitrophenol (CNP) was followed continuously for 10 min at 37 °C by spectrophotometry at 405 nm on a EL808 microplate reader (Bio-Tek Instruments, Inc.). The initial velocity ([CNP] (μm/min)) was calculated using a standard curve of CNP in the range of 0–200 μm. The kinetic parameters were calculated by fitting the data to the Michaelis–Menten equation using the software KaleidaGraphTM 3.0. The procedure for HjCel3A kinetics has been described previously by Karkehabadi et al. (
). For HjCel3A with CNPG as substrate, the kinetic parameters were calculated at low substrate concentrations (up to 800 μm) because of the inhibition and subsequent decrease of activity observed at higher substrate concentrations. In the kinetic experiments of the β-glucosidases, the substrate was incubated with the enzyme at 37 °C at the optimum pH. At timed intervals (2 min), an aliquot was withdrawn and heated to 100 °C for 5 min to stop the reaction. The released glucose was determined using the glucose oxidase–peroxidase assay.
Identification of transglycosylation products
Experiments were performed using 100 and 50 mm cellobiose for HjCel3A and HjCel3B, respectively. The reaction was incubated at 37 °C and optimum pH for 19 h for HjCel3A and 24 h for HjCel3B. The reaction mixture was loaded to a Carbograph column (1000 mg/15 ml; Alltech), which was subsequently washed with water to remove excess glucose. The transglycosylation products with cellobiose were eluted with 30% (v/v) acetonitrile (
). The eluent was removed on a Rotavapor, and the residue was redissolved in water before analysis by high-performance anion-exchange chromatography with pulsed amperometric detection (ICS-3000, Dionex, Sunnyvale, CA).
Enzyme–substrate complex models
The crystal structures of HjCel3A with and without a glucose molecule bound in the −1 substrate site (PDB codes 3ZYZ and 3ZZ1, respectively) were reported previously by our group (
). The structure of HjCel3B (UniProt ID Q7Z9M5) was modeled based on the crystal structure of A. aculeatus β-glucosidase 1 (AaBGL1; PDB code 4IIB; 58.6% sequence identity) (
). Residues 1–28 were not included in the final model owing to the unavailability of a structural template. Side chain conformations of residues in the active site region were manually refined using Coot (
Structural and functional studies of the glycoside hydrolase family 3 β-glucosidase Cel3A from the moderately thermophilic fungus Rasamsonia emersonii.
Acta Crystallogr. D Struct. Biol.2016; 72 (27377383): 860-870
). Residues 641–655 of HjCel3B, constituting an inserted loop region when compared with the latter two enzymes, had an implausible conformation in the generated homology model. This 15-residue peptide was manually adjusted using Coot. First, an α-helix at the N-terminal end of the loop was extended such that the length agreed with both R. emersonii Cel3A and AaBGL1. Second, the 10-residue insertion was adjusted such that a relaxed loop conformation was adopted, stretching from the end of the helix to residue 656; each side chain in the manually adjusted loop region was initially placed in the most commonly occurring rotamer. The CASTp server was used to identify the substrate pocket of HjCel3B as well as HjCel3A and calculate the solvent-excluded volume (
Five simulation models each were constructed for HjCel3A and HjCel3B (Fig. 5). HjCel3B was modeled as a monomer to reduce computational resource requirements. This is not expected to affect the results as the active site is far from the homodimeric interface (∼40 Å) (Fig. 2A). In Models I and II, the enzyme is in the reactive state during the glycosylation step wherein the acid/base catalyst Glu is protonated (Fig. 1). Glucose is bound in the −1 site (Model I), and cellobiose is bound across the −1 and +1 sites (Model II). Models III–V were used to study the binding of glucose or cellobiose as an acceptor for the transglycosylation reaction. The enzyme is in the GEI state with a glucose molecule covalently bound to the nucleophile Asp (Fig. 1). Model III represents the case in which either glucose is retained in the +1 product site after cleavage of the cellobiose glycosidic bond or glucose has reentered the active site from solution. Models IV and V represent the scenario in which either cellobiose, under high substrate conditions, enters the enzyme product site or the enzyme has acted upon a cellotriose molecule; two different orientations of cellobiose are considered, both with the nonreducing end positioned in the +1 site.
MD simulations
System preparation and equilibration are described in detail in the supporting information. Protein and substrates were modeled using the CHARMM36 force field (
Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone ϕ, ψ and side-chain χ1 and χ2 dihedral angles.
J. Chem. Theory Comput.2012; 8 (23341755): 3257-3273
) in the NVT ensemble at 300 K with a 2-fs time step. The simulation time varied between 110 and 250 ns as, in some cases, the substrate diffused from the relatively shallow binding site pockets. The Langevin thermostat was used for temperature control (
) was used to treat long-range electrostatic interactions. All other simulation parameters are the same as those used in equilibration. Hydrogen bond occupancies, with a distance cutoff of 3.0 Å and angle cutoff of 135° (
). Pairwise interaction energies were calculated using NAMD.
FEP/λ-REMD calculations
Free energy calculations were performed by dividing the binding process according to a thermodynamic cycle wherein 1) the bound ligand (shown in blue in Fig. 5) is decoupled from the enzyme and 2) the solvated ligand is decoupled from bulk solution. The binding free energy, ΔGb0, for the solvated enzyme–substrate complex is the difference between the free energies associated with 1) and 2). Insertion of the ligand into the binding pocket was done in three stages using the thermodynamic coupling parameters λrep, λdisp, and λelec, giving the repulsive (ΔGrep), dispersive (ΔGdisp), and electrostatic (ΔGelec) contributions, respectively. A fourth parameter, λrstr, controls the translational and orientational restraints on the ligand and gives ΔGrstr. Insertion of the ligand into the bulk was calculated similarly but without the restraint (
) calculations were taken from the 25-ns snapshot of each MD simulation of the enzyme–substrate complexes. At least 20 sequential, 0.1-ns simulations were performed using 128 replicas (72 repulsive, 24 dispersive, and 32 electrostatic) with an exchange frequency of 1/100 steps. Further details of the calculation are as described in a previous publication (
). Repulsive, dispersive, and electrostatic free energy contributions were obtained from the last 1 ns of data using the multistate Bennett acceptance ratio (
), and the standard deviation was subsequently calculated.
Author contributions
I. G., P. N., K. P., M. G., H. H., and M. S. data curation; I. G., P. N., K. P., M. G., and C. M. P. formal analysis; I. G., P. N., K. P., M. G., and H. H. investigation; I. G. visualization; I. G., H. H., M. S., and C. M. P. methodology; I. G. and C. M. P. writing-original draft; I. G., M. G., H. H., and C. M. P. writing-review and editing; M. S. and C. M. P. conceptualization; M. S. and C. M. P. resources; M. S. and C. M. P. supervision; M. S. and C. M. P. funding acquisition; C. M. P. project administration.
Acknowledgments
We thank Dr. Wim Nerinckx of Ghent University for assistance in synthesizing CNPG and CNPX and constructing the p-aminobenzyl-thio-β-glucopyranoside affinity column. Computing resources were provided by the University of Kentucky (DLX and Kingsley clusters) and the National Science Foundation (NSF) Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by NSF Grant ACI-1548562 (Gordon cluster under allocation MCB090159).
Cloning and amplification of the gene encoding an extracellular β-glucosidase from Trichoderma reesei: evidence for improved rates of saccharification of cellulosic substrates.
Trp-49 of the family 3 β-glucosidase from Aspergillus niger is important for its transglucosidic activity: creation of novel β-glucosidases with low transglucosidic efficiencies.
Transglucosidic reactions of the Aspergillus niger family 3 β-glucosidase: qualitative and quantitative analyses and evidence that the transglucosidic rate is independent of pH.
Trp-262 is a key residue for the hydrolytic and transglucosidic reactivity of the Aspergillus niger family 3 β-glucosidase: substitution results in enzymes with mainly transglucosidic activity.
Kinetics of substrate transglycosylation by glycoside hydrolase family 3 glucan (1→3)-β-glucosidase from the white-rot fungus Phanerochaete chrysosporium.
Epoxyalkyl glycosides of d-xylose and xylo-oligosaccharides are active-site markers of xylanases from glycoside hydrolase family 11, not from family 10.
Structural and functional studies of the glycoside hydrolase family 3 β-glucosidase Cel3A from the moderately thermophilic fungus Rasamsonia emersonii.
Acta Crystallogr. D Struct. Biol.2016; 72 (27377383): 860-870
Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone ϕ, ψ and side-chain χ1 and χ2 dihedral angles.
J. Chem. Theory Comput.2012; 8 (23341755): 3257-3273
This work was supported in part by the Swedish University of Agricultural Sciences August T. Larsson Guest Researcher Programme (to C. M. P.), National Science Foundation Experimental Program to Stimulate Competitive Research Grant 1355438 (to I. G.), and Swedish Governmental Agency for Innovation Systems (VINNOVA) Grant 2015-04835. The authors declare that they have no conflicts of interest with the contents of this article.