Differences in Gluco and Galacto Substrate-Binding Interactions in a Dual 6Pß-Glucosidase/6Pß-Galactosidase Glycoside Hydrolase 1 Enzyme from Bacillus licheniformis.

Bacterial glycoside hydrolase 1 (GH1) enzymes with 6-phospho-ß-galactosidase and 6-phospho-ß-glucosidase activities have the important task of releasing phosphorylated and nonphosphorylated monosaccharides into the cytoplasm. Curiously, dual 6-phospho-ß-galactosidase/6-phospho-ß-glucosidase (dual-phospho) enzymes have broad specificity and are able to hydrolyze galacto- and gluco-derived substrates. This study investigates the structure and substrate specificity of a GH family 1 enzyme from Bacillus licheniformis, hereafter known as BlBglC. The enzyme structure has been solved, and sequence analysis, molecular dynamics simulations, and binding free energy calculations offered evidence of dual-phospho activity. Both test ligands p-nitrophenyl-ß-d-galactoside-6-phosphate (PNP6Pgal) and p-nitrophenyl-ß-d-glucoside-6-phosphate (PNP6Pglc) demonstrated strong binding to BlBglC although the pose and interactions of the PNP6Pglc triplicates were slightly more consistent. Interestingly, known specificity-inducing residues, Gln23 and Trp433, bind strongly to the ligand O3 hydroxyl group in the PNP6Pgal-BlBglC complex and to the ligand O4 hydroxyl group in the PNP6Pglc-BlBglC complex. Additionally, the BlBglC-His124 residue is a major contributor of hydrogen bonds to the PNP6Pgal O3 hydroxyl group but does not form any hydrogen bonds with PNP6Pglc. On the other hand, BlBglC residues Tyr173, Tyr301, Gln302, and Thr321 form hydrogen bonds with PNP6Pglc but not PNP6Pgal. These findings provide important details of the broad specificity of dual-phospho activity GH1 enzymes.

Combining Free Energy Simulations and NMR Chemical-Shift Perturbation To Identify Transient Cation-π Contacts in Proteins.

Flexible protein regions containing cationic and aromatic side-chains exposed to solvent may form transient cation-π interactions with structural and functional roles. To evaluate their stability and identify important intramolecular cation-π contacts, a combination of free energy profiles estimated from umbrella sampling with molecular dynamics simulations and chemical shift perturbations (CSP) obtained from nuclear magnetic resonance (NMR) experiments is applied here to the complete catalytic domain of human phosphatase Cdc25B. This protein is a good model system for transient cation-π interactions as it contains only one Trp residue (W550) in the disordered C-terminal segment and a total of 17 Arg residues, many exposed to solvent. Eight putative Arg-Trp pairs were simulated here. Only R482 and R544 show bound profiles corresponding to important transient cation-π interactions, while the others have dissociative or almost flat profiles. These results are corroborated by CSP analysis of three Cdc25B point mutants (W550A, R482A, and R544A) disrupting cation-π contacts. The proposed validation of statistically representative molecular simulations by NMR spectroscopy could be applied to identify transient contacts of proteins in general but carefully, as NMR chemical shifts are sensitive to changes in both molecular contacts and conformational distributions.

X-ray Structure, Bioinformatics Analysis, and Substrate Specificity of a 6-Phospho-ß-glucosidase Glycoside Hydrolase 1 Enzyme from Bacillus licheniformis.

In bacteria, mono- and disaccharides are phosphorylated during the uptake processes through the vastly spread transport system phosphoenolpyruvate-dependent phosphotransferase. As an initial step in the phosphorylated disaccharide metabolism pathway, 6-phospho-ß-glucosidases and 6-phospho-ß-galactosidases play a crucial role by releasing phosphorylated and nonphosphorylated monosaccharides. However, structural determinants for the specificity of these enzymes still need to be clarified. Here, an X-ray structure of a glycoside hydrolase family 1 enzyme from Bacillus licheniformis, hereafter known as BlBglH, was determined at 2.2 Å resolution, and its substrate specificity was investigated. The sequence of BlBglH was compared to the sequences of 58 other GH1 enzymes using sequence alignments, sequence identity calculations, phylogenetic analysis, and motif discovery. Through these various analyses, BlBglH was found to have sequence features characteristic of the 6-phospho-ß-glucosidase activity enzymes. Motif and structural observations highlighted the importance of loop L8 in 6-phospho-ß-glucosidase activity enzymes. To further affirm enzyme specificity, molecular docking and molecular dynamics simulations were performed using the crystallographic structure of BlBglH. Docking was carried out with a 6-phospho-ß-glucosidase enzyme activity positive and negative control ligand, followed by 400 ns of MD simulations. The positive and negative control ligands were PNP6Pglc and PNP6Pgal, respectively. PNP6Pglc maintained favorable interactions within the active site until the end of the MD simulation, while PNP6Pgal exhibited instability. The favorable binding of substrate stabilized the loops that surround the active site. Binding free energy calculations showed that the PNP6Pglc complex had a substantially lower binding energy compared to the PNP6Pgal complex. Altogether, the findings of this study suggest that BlBglH possesses 6-phospho-ß-glucosidase enzymatic activity and revealed sequence and structural differences between bacterial GH1 enzymes of various activities.

Conformational flexibility of the complete catalytic domain of Cdc25B phosphatases.

Cdc25B phosphatases are involved in cell cycle checkpoints and have become a possible target for developing new anticancer drugs. A more rational design of Cdc25B ligands would benefit from detailed knowledge of its tertiary structure. The conformational flexibility of the C-terminal region of the Cdc25B catalytic domain has been debated recently and suggested to play an important structural role. Here, a combination of experimental NMR measurements and molecular dynamics simulations for the complete catalytic domain of the Cdc25B phosphatase is presented. The stability of the C-terminal α-helix is confirmed, but the last 20 residues in the complete catalytic domain are very flexible, partially occlude the active site and may establish transient contacts with the protein core. This flexibility in the C-terminal tail may modulate the molecular recognition of natural substrates and competitive inhibitors by Cdc25B. Proteins 2016; 84:1567-1575. © 2016 Wiley Periodicals, Inc.

Mechanism of imidazole inhibition of a GH1 ß-glucosidase.

Imidazole is largely employed in recombinant protein purification, including GH1 ß-glucosidases, but its effect on the enzyme activity is rarely taken into consideration. Computational docking suggested that imidazole interacts with residues forming the active site of the GH1 ß-glucosidase from Spodoptera frugiperda (Sfßgly). We confirmed this interaction by showing that imidazole reduces the activity of Sfßgly, which does not result from enzyme covalent modification or promotion of transglycosylation reactions. Instead, this inhibition occurs through a partial competitive mechanism. Imidazole binds to the Sfßgly active site, reducing the substrate affinity by about threefold, whereas the rate constant of product formation remains unchanged. The binding of imidazole within the active site was further confirmed by enzyme kinetic experiments in which imidazole and cellobiose competed to inhibit the hydrolysis of p-nitrophenyl ß-glucoside. Finally, imidazole interaction in the active site was also demonstrated by showing that it hinders access of carbodiimide to the Sfßgly catalytic residues, protecting them from chemical inactivation. In conclusion, imidazole binds in the Sfßgly active site, generating a partial competitive inhibition. Considering that GH1 ß-glucosidases share conserved active sites, this inhibition phenomenon is probably widespread among these enzymes, and this should be taken into account when considering the characterization of their recombinant forms.

Single mutations outside the active site affect the substrate specificity in a ß-glycosidase.

A library of random mutants of the ß-glycosidase Sfßgly was screened for mutations that affect its specificity for the substrate glycone (ß-d-fucoside versus ß-d-glucoside). Among mutations selected (T35A, R189G, Y345C, P348L, S358F, S378G, N400D, S424F, F460L, and R474H), eight occurred in the C-terminal half of Sfßgly and only two were at the active site (R189G and N400D). Tryptophan fluorescence spectra and thermal inactivation showed that the selected mutants and wild-type Sfßgly are similarly folded. Enzyme kinetics confirmed that these mutations resulted in broadening or narrowing of the preference for the substrate glycone. Structural modeling and interaction maps revealed contact pathways that connect the sites of the selected mutations through up to three interactions to the active site residues E399, W444, and E187, which are involved in substrate binding and catalysis. Interestingly, independently selected mutations (Y345C, P348L, and R189G; S424F and N400D) were placed on the same contact pathway. Moreover, (k(cat)/K(m) fucoside)/(k(cat)/K(m) glucoside) ratios showed that mutations at intermediate residues of the same contact pathway often had similar effects on substrate specificity. Finally mutations in the same contact pathway caused similar structural disturbance as evidenced by acrylamide quenching of the Sfßgly fluorescence. Based on these data, it is proposed that the effects of the selected mutations were propagated into the active site through groups of interacting residues (contact pathways) changing the Sfßgly substrate specificity.

A simple method to determine changes in the affinity between HisF and HisH in the Imidazole Glycerol Phosphate Synthase heterodimer.

The bi-enzyme HisF-HisH heterodimer is part of the pathway that produces histidine and purines in bacteria and lower eukaryotes, but it is absent in mammals. This heterodimer has been largely studied probing the basis of the allosteric effects and the structural stability in proteins. It is also a potential target for antibacterial drugs. In this work, we developed a simple method to evaluate changes in the affinity between HisF and HisH in the heterodimer of the bacteria Thermotoga maritima. HisH contains a single tryptophan residue, which is exposed in the free protein, but buried in the heterodimer interface. Hence, the intrinsic fluorescence maximum of this residue changes to shorter wavelengths upon dimerization. Thus, we used the fluorescence intensity at this shorter wavelength to monitor heterodimer accumulation when HisH was combined with sub-stoichiometric HisF. Under conditions where the HisF-HisH heterodimer is in equilibrium with the free states of these enzymes, when [HisH] > [HisF], we deduced a linear function connecting [HisF-HisH] to [HisF], in which the slope depends on the heterodimer dissociation constant (Kd). Based on this equation, taking fluorescence intensities as proxies of the heterodimer and HisF concentrations, we experimentally determined the Kd at four different temperatures. These Kd values were compared to those evaluated using ITC. Both methods revealed an increase in the HisF and HisH binding affinity as the temperature increases. In spite of differences in their absolute values, the Kd determined using these methods presented an evident linear correlation. To demonstrate the effectiveness of the fluorescence method we determined the effect on the Kd caused by 12 single mutations in HisF. Coherently, this test singled out the only mutation in the binding interface. In brief, the method described here effectively probes qualitative effects on the Kd, can be carried out using common laboratory equipment and is scalable.

Homodimerization of a glycoside hydrolase family GH1 ß-glucosidase suggests distinct activity of enzyme different states.

In this work, we investigated how activity and oligomeric state are related in a purified GH1 ß-glucosidase from Spodoptera frugiperda (Sfßgly). Gel filtration chromatography coupled to a multiple angle light scattering detector allowed separation of the homodimer and monomer states and determination of the dimer dissociation constant (KD ), which was in the micromolar range. Enzyme kinetic parameters showed that the dimer is on average 2.5-fold more active. Later, we evaluated the kinetics of homodimerization, scanning the changes in the Sfßgly intrinsic fluorescence over time when the dimer dissociates into the monomer after a large dilution. We described how the rate constant of monomerization (koff ) is affected by temperature, revealing the enthalpic and entropic contributions to the process. We also evaluated how the rate constant (kobs ) by which equilibrium is reached after dimer dilution behaves when varying the initial Sfßgly concentration. These data indicated that Sfßgly dimerizes through the conformational selection mechanism, in which the monomer undergoes a conformational exchange and then binds to a similar monomer, forming a more active homodimer. Finally, we noted that conformational selection reports and experiments usually rely on a ligand whose concentration is in excess, but for homodimerization, this approach does not hold. Hence, since our approach overcomes this limitation, this study not only is a new contribution to the comprehension of GH1 ß-glucosidases, but it can also help to elucidate protein interaction pathways.

Optimum temperature may be a misleading parameter in enzyme characterization and application.

The optimum temperature is commonly determined in enzyme characterization. A search in the PubMed database for "optimum temperature" and "enzymes" yielded more than 1,700 manuscripts reporting this parameter over the last five years. Here, we show that the optimum temperature is not a constant. The catalytic activity of the mesophylic ß-glucosidase Sfßgly was determined at different temperatures using different assay times and enzyme concentrations. We observed that the optimum temperature for Sfßgly changed by 5°C simply by modifying the assay length, and it was inversely correlated with enzyme concentration. These observations rely on the fact that close to the melting temperature, thermal denaturation continuously decreases the active enzyme concentration as the assay progresses. Thus, as the denaturation rate increases with increasing temperature, the bell-shaped curves observed in "activity versus temperature plots" occur only if the enzyme is denatured at and above the optimum temperature, which was confirmed using the thermostable ß-glucosidase bglTm. Thus, the optimum temperature hardly reflects an intrinsic enzyme property and is actually a mere consequence of the assay condition. Thus, adoption of the "optimum temperature" determined under bench conditions for large-scale uses, which differ in assay length and enzyme concentration, may result in lower yields and financial losses.

The interplay of processivity, substrate inhibition and a secondary substrate binding site of an insect exo-beta-1,3-glucanase.

Abracris flavolineata midgut contains a processive exo-beta-glucanase (ALAM) with lytic activity against Saccharomyces cerevisiae, which was purified (yield, 18%; enrichment, 37 fold; specific activity, 1.89 U/mg). ALAM hydrolyses fungal cells or callose from the diet. ALAM (45 kDa; pI 5.5; pH optimum 6) major products with 0.6 mM laminarin as substrate are beta-glucose (61%) and laminaribiose (39%). Kinetic data obtained with laminaridextrins and methylumbelliferyl glucoside suggest that ALAM has an active site with at least six subsites. The best fitting of kinetic data to theoretical curves is obtained using a model where one laminarin molecule binds first to a high-affinity accessory site, causing active site exposure, followed by the transference of the substrate to the active site. The two-binding-site model is supported by results from chemical modifications of amino acid residues and by ALAM action in MUbetaGlu plus laminarin. Low laminarin concentrations increase the modification of His, Tyr and Asp or Glu residues and MUbetaGlu hydrolysis, whereas high concentrations abolish modification and inhibit MUbetaGlu hydrolysis. Our data indicate that processivity results from consecutive transferences of substrate between accessory and active site and that substrate inhibition arises when both sites are occupied by substrate molecules abolishing processivity.

The role in the substrate specificity and catalysis of residues forming the substrate aglycone-binding site of a beta-glycosidase.

The relative contributions to the specificity and catalysis of aglycone, of residues E190, E194, K201 and M453 that form the aglycone-binding site of a beta-glycosidase from Spodoptera frugiperda (EC 3.2.1.21), were investigated through site-directed mutagenesis and enzyme kinetic experiments. The results showed that E190 favors the binding of the initial portion of alkyl-type aglycones (up to the sixth methylene group) and also the first glucose unit of oligosaccharidic aglycones, whereas a balance between interactions with E194 and K201 determines the preference for glucose units versus alkyl moieties. E194 favors the binding of alkyl moieties, whereas K201 is more relevant for the binding of glucose units, in spite of its favorable interaction with alkyl moieties. The three residues E190, E194 and K201 reduce the affinity for phenyl moieties. In addition, M453 favors the binding of the second glucose unit of oligosaccharidic aglycones and also of the initial portion of alkyl-type aglycones. None of the residues investigated interacted with the terminal portion of alkyl-type aglycones. It was also demonstrated that E190, E194, K201 and M453 similarly contribute to stabilize ES(double dagger). Their interactions with aglycone are individually weaker than those formed by residues interacting with glycone, but their joint catalytic effects are similar. Finally, these interactions with aglycone do not influence glycone binding.

Search for independent (ß/α)4 subdomains in a (ß/α)8 barrel ß-glucosidase.

Proteins that fold as (ß/α)8 barrels are thought to have evolved from half-barrels that underwent duplication and fusion events. The evidence is particularly clear for small barrels, which have almost identical halves. Additionally, computational calculations of the thermodynamic stability of these structures in the presence of denaturants have revealed that (ß/α)8 barrels contain two subunits or domains corresponding to half-barrels. Hence, within (ß/α)8 barrels, half-barrels are self-contained units. Here, we tested this hypothesis using ß-glucosidase from the bacterium Thermotoga maritima (bglTm), which has a (ß/α)8 barrel structure. Mutations were introduced to disrupt the noncovalent contacts between its halves and reveal the presence of two domains within bglTm, thus resulting in the creation of mutants T1 (containing W12A and I217A mutations) and T2 (containing W12A, H195A, I217A and F404A mutations). Mutants T1 and T2 were properly folded, as indicated by their fluorescence spectra and enzyme kinetic parameters. T1 and wild-type bglTm were equally stable, as shown by the results of thermal inactivation, differential scanning fluorimetry and guanidine hydrochloride denaturation experiments. However, T2 showed a first-order inactivation at 80°C, a single melting temperature of 82°C and only one transition concentration (c50) in 2.4 M guanidine hydrochloride. Additionally, T1 and T2 exhibited a cooperative denaturation process that followed a two-state model (m-values equal to 1.4 and 1.6 kcal/mol/M, respectively), similar to that of wild-type bglTm (1.2 kcal/mol/M). Hence, T1 and T2 each denatured as a single unit, although they contained different degrees of disruption between their halves. In conclusion, bglTm halves are equivalent in terms of their thermal and chemical stability; thus, their separate contributions to (ß/α)8 barrel unfolding cannot be disentangled.

Mutations close to a hub residue affect the distant active site of a GH1 ß-glucosidase.

The tertiary structure of proteins has been represented as a network, in which residues are nodes and their contacts are edges. Protein structure networks contain residues, called hubs or central, which are essential to form short connection pathways between any pair of nodes. Hence hub residues may effectively spread structural perturbations through the protein. To test whether modifications nearby to hub residues could affect the enzyme active site, mutations were introduced in the ß-glycosidase Sfßgly (PDB-ID: 5CG0) directed to residues that form an α-helix (260-265) and a ß-strand (335-337) close to one of its main hub residues, F251, which is approximately 14 Å from the Sfßgly active site. Replacement of residues A263 and A264, which side-chains project from the α-helix towards F251, decreased the rate of substrate hydrolysis. Mutation A263F was shown to weaken noncovalent interactions involved in transition state stabilization within the Sfßgly active site. Mutations placed on the opposite side of the same α-helix did not show these effects. Consistently, replacement of V336, which side-chain protrudes from a ß-strand face towards F251, inactivated Sfßgly. Next to V336, mutation S337F also caused a decrease in noncovalent interactions involved in transition state stabilization. Therefore, we suggest that mutations A263F, A264F, V336F and S337F may directly perturb the position of the hub F251, which could propagate these perturbations into the Sfßgly active site through short connection pathways along the protein network.

Structural and biochemical characterization of a GH3 ß-glucosidase from the probiotic bacteria Bifidobacterium adolescentis.

Bifidobacterium is an important genus of probiotic bacteria colonizing the human gut. These bacteria can uptake oligosaccharides for the fermentative metabolism of hexoses and pentoses, producing lactate, acetate as well as short-chain fatty acids and propionate. These end-products are known to have important effects on human health. ß-glucosidases (EC 3.2.1.21) are pivotal enzymes for the metabolism and homeostasis of Bifidobacterium, since they hydrolyze small and soluble saccharides, typically producing glucose. Here we describe the cloning, expression, biochemical characterization and the first X-ray structure of a GH3 ß-glucosidase from the probiotic bacteria Bifidobacterium adolescentis (BaBgl3). The purified BaBgl3 showed a maximal activity at 45 °C and pH 6.5. Under the optimum conditions, BaBgl3 is highly active on 4-nitrophenyl-ß-d-glucopyranoside (pNPG) and, at a lesser degree, on 4-nitrophenyl-ß-d-xylopyranoside (pNPX, about 32% of the activity observed for pNPG). The 2.4 Šresolution crystal structure of BaBgl3 revealed a three-domain structure composed of a TIM barrel domain, which together with α/ß sandwich domain accommodate the active site and a third C-terminal fibronectin type III (FnIII) domain with unknown function. Modeling of the substrate in the active site indicates that an aspartate interacts with the hydroxyl group of the C6 present in pNPG but absent in pNPX, which explains the substrate preference. Finally, the enzyme is significantly stabilized by glycerol and galactose, resulting in considerable increase in the enzyme activity and its lifetime. The structural and biochemical studies presented here provide a deeper understanding of the molecular mechanisms of complex carbohydrates degradation utilized by probiotic bacteria as well as for the development of new prebiotic oligosaccharides.

Structural insights into ß-glucosidase transglycosylation based on biochemical, structural and computational analysis of two GH1 enzymes from Trichoderma harzianum.

ß-glucosidases are glycoside hydrolases able to cleave small and soluble substrates, thus producing monosaccharides. These enzymes are distributed among families GH1, GH2, GH3, GH5, GH9, GH30 and GH116, with GH1 and GH3 being the most relevant families with characterized enzymes to date. A recent transcriptomic analysis of the fungus Trichoderma harzianum, known for its increased ß-glucosidase activity as compared to Trichoderma reesei, revealed two enzymes from family GH1 with high expression levels. Here we report the cloning, recombinant expression, purification and crystallization of these enzymes, ThBgl1 and ThBgl2. A close inspection of the enzymatic activity of these enzymes surprisingly revealed a marked difference between them despite the sequence similarity (53%). ThBgl1 has an increased tendency to catalyze transglycosylation reaction while ThBgl2 acts more as a hydrolyzing enzyme. Detailed comparison of their crystal structures and the analysis of the molecular dynamics simulations reveal the presence of an asparagine residue N186 in ThBgl2, which is replaced by the phenylalanine F180 in ThBgl1. This single amino acid substitution seems to be sufficient to create a polar environment that culminates with an increased availability of water molecules in ThBgl2 as compared to ThBgl1, thus conferring stronger hydrolyzing character to the former enzyme.

Structural and biochemical data of Trichoderma harzianum GH1 ß-glucosidases.

Here the statistics concerning X-ray data processing and structure refinement are given, together with the substrate preference analysis for ThBgl1 and ThBgl2. Finally, the analysis of the influence of temperature and pH on the activities of both enzymes are shown.

Protein thermal denaturation is modulated by central residues in the protein structure network.

Network structural analysis, known as residue interaction networks or graphs (RIN or RIG, respectively) or protein structural networks or graphs (PSN or PSG, respectively), comprises a useful tool for detecting important residues for protein function, stability, folding and allostery. In RIN, the tertiary structure is represented by a network in which residues (nodes) are connected by interactions (edges). Such structural networks have consistently presented a few central residues that are important for shortening the pathways linking any two residues in a protein structure. To experimentally demonstrate that central residues effectively participate in protein properties, mutations were directed to seven central residues of the ß-glucosidase Sfßgly (ß-D-glucoside glucohydrolase; EC 3.2.1.21). These mutations reduced the thermal stability of the enzyme, as evaluated by changes in transition temperature (Tm ) and the denaturation rate at 45 °C. Moreover, mutations directed to the vicinity of a central residue also caused significant decreases in the Tm of Sfßgly and clearly increased the unfolding rate constant at 45 °C. However, mutations at noncentral residues or at surrounding residues did not affect the thermal stability of Sfßgly. Therefore, the data reported in the present study suggest that the perturbation of the central residues reduced the stability of the native structure of Sfßgly. These results are in agreement with previous findings showing that networks are robust, whereas attacks on central nodes cause network failure. Finally, the present study demonstrates that central residues underlie the functional properties of proteins.

A mutant ß-glucosidase increases the rate of the cellulose enzymatic hydrolysis.

The enzymatic hydrolysis of cellulose and lignocellulosic materials is marked by a rate decrease along the reaction time. Cellobiohydrolase slow dissociation from the substrate and its inhibition by the cellobiose produced are relevant factors associated to the rate decrease. In that sense, addition of ß-glucosidases to the enzyme cocktails employed in cellulose enzymatic hydrolysis not only produces glucose as final product but also reduces the cellobiohydrolase inhibition by cellobiose. The digestive ß-glucosidase GH1 from the fall armyworm Spodoptera frugiperda, hereafter called Sfßgly, containing the mutation L428V showed an increased kcat for cellobiose hydrolysis. In comparison to assays conducted with the wild-type Sfßgly and cellobiohydrolase TrCel7A, the presence of the mutant L428V increased in 5 fold the initial rate of crystalline cellulose hydrolysis and reduced to one quarter the time needed to TrCel7A produce the maximum glucose yield. As our results show that mutant L428V complement the action of TrCel7A, the introduction of the equivalent replacement in ß-glucosidases is a promising strategy to reduce costs in the enzymatic hydrolysis of lignocellulosic materials.

Using the Amino Acid Network to Modulate the Hydrolytic Activity of ß-Glycosidases.

The active site residues in GH1 ß-glycosidases are compartmentalized into 3 functional regions, involved in catalysis or binding of glycone and aglycone motifs from substrate. However, it still remains unclear how residues outside the active site modulate the enzymatic activity. To tackle this question, we solved the crystal structure of the GH1 ß-glycosidase from Spodoptera frugiperda (Sfßgly) to systematically map its residue contact network and correlate effects of mutations within and outside the active site. External mutations neighbouring the functional residues involved in catalysis and glycone-binding are deleterious, whereas mutations neighbouring the aglycone-binding site are less detrimental or even beneficial. The large dataset of new and previously characterized Sfßgly mutants supports that external perturbations are coherently transmitted to active site residues possibly through contacts and specifically disturb functional regions they interact to, reproducing the effects observed for direct mutations of functional residues. This allowed us to suggest that positions related to the aglycone-binding site are preferential targets for introduction of mutations aiming to further improve the hydrolytic activity of ß-glycosidases.

Half-Barrels Derived from a (ß/α)8 Barrel ß-Glycosidase Undergo an Activation Process.

The evolution of (ß/α)8 barrel proteins is currently thought to have involved the fusion of two (ß/α)4 half-barrels, thereby conferring stability on the protein structure. After the formation of a whole (ß/α)8 barrel, this structure could evolve and diverge to form fully active enzymes. Interestingly, we show here that isolated (ß/α)4 half-barrels derived from the N- and C-terminal domains of the ß-glucosidase Sfßgly (Sfßgly-N: residues 1 to 265; Sfßgly-C: residues 266 to 509) undergo an activation process, which renders them catalytically active. The rate constants of the activation process were calculated to be 0.029 and 0.032 h-1 for Sfßgly-N and Sfßgly-C, respectively. Moreover, the Sfßgly-N and Sfßgly-C activation processes were simultaneous with modifications in their initial structure, which reduced the exposure of their tryptophan residues. Importantly, this activation was also coincident with an increase in the sizes of Sfßgly-N and Sfßgly-C particles. These novel observations suggest that the change in catalytic activity associated with the transition from a half to whole (ß/α)8 barrel might also have driven such an evolutionary process.