Ionic network at the C-terminus of the beta-glycosidase from the hyperthermophilic archaeon Sulfolobus solfataricus: Functional role in the quaternary structure thermal stabilization

Residue Effect-Guided Design: Engineering of S. Solfataricus β-Glycosidase to Enhance Its Thermostability and Bioproduction of Ginsenoside Compound K

β-Glycosidase from Sulfolobus solfataricus (SS-BGL) is a highly effective biocatalyst for the synthesis of compound K (CK) from glycosylated protopanaxadiol ginsenosides. In order to improve the thermal stability of SS-BGL, molecular dynamics simulations were used to determine the residue-level binding energetics of ginsenoside Rd in the SS-BGL-Rd docked complex and to identify the top ten critical contributors. Target sites for mutations were determined using dynamic cross-correlation mapping of residues via the Ohm server to identify networks of distal residues that interact with the key binding residues. Target mutations were determined rationally based on site characteristics. Single mutants and then recombination of top hits led to the two most promising variants SS-BGL-Q96E/N97D/N302D and SS-BGL-Q96E/N97D/N128D/N302D with 2.5-fold and 3.3-fold increased half-lives at 95 °C, respectively. The enzyme activities relative to those of wild-type for ginsenoside conversion were 161 and 116%, respectively..

Transcript Regulation of the Recoded Archaeal α-l-Fucosidase In Vivo

Genetic decoding is flexible, due to programmed deviation of the ribosomes from standard translational rules, globally termed "recoding". In Archaea, recoding has been unequivocally determined only for termination codon readthrough events that regulate the incorporation of the unusual amino acids selenocysteine and pyrrolysine, and for -1 programmed frameshifting that allow the expression of a fully functional α-l-fucosidase in the crenarchaeon Saccharolobus solfataricus, in which several functional interrupted genes have been identified. Increasing evidence suggests that the flexibility of the genetic code decoding could provide an evolutionary advantage in extreme conditions, therefore, the identification and study of interrupted genes in extremophilic Archaea could be important from an astrobiological point of view, providing new information on the origin and evolution of the genetic code and on the limits of life on Earth. In order to shed some light on the mechanism of programmed -1 frameshifting in Archaea, here we report, for the first time, on the analysis of the transcription of this recoded archaeal α-l-fucosidase and of its full-length mutant in different growth conditions in vivo. We found that only the wild type mRNA significantly increased in S. solfataricus after cold shock and in cells grown in minimal medium containing hydrolyzed xyloglucan as carbon source. Our results indicated that the increased level of fucA mRNA cannot be explained by transcript up-regulation alone. A different mechanism related to translation efficiency is discussed.

Biomass-degrading glycoside hydrolases of archaeal origin

During the last decades, the impact of hyperthermophiles and their enzymes has been intensively investigated for implementation in various high-temperature biotechnological processes. Biocatalysts of hyperthermophiles have proven to show extremely high thermo-activities and thermo-stabilities and are identified as suitable candidates for numerous industrial processes with harsh conditions, including the process of an efficient plant biomass pretreatment and conversion. Already-characterized archaea-originated glycoside hydrolases (GHs) have shown highly impressive features and numerous enzyme characterizations indicated that these biocatalysts show maximum activities at a higher temperature range compared to bacterial ones. However, compared to bacterial biomass-degrading enzymes, the number of characterized archaeal ones remains low. To discover new promising archaeal GH candidates, it is necessary to study in detail the microbiology and enzymology of extremely high-temperature habitats, ranging from terrestrial to marine hydrothermal systems. State-of-the art technologies such as sequencing of genomes and metagenomes and automated binning of genomes out of metagenomes, combined with classical microbiological culture-dependent approaches, have been successfully performed to detect novel promising biomass-degrading hyperthermozymes. In this review, we will focus on the detection, characterization and similarities of archaeal GHs and their unique characteristics. The potential of hyperthermozymes and their impact on high-temperature industrial applications have not yet been exhausted.

Discovery of hyperstable carbohydrate-active enzymes through metagenomics of extreme environments

The enzymes from hyperthermophilic microorganisms populating volcanic sites represent interesting cases of protein adaptation and biotransformations under conditions where conventional enzymes quickly denature. The difficulties in cultivating extremophiles severely limit access to this class of biocatalysts. To circumvent this problem, we embarked on the exploration of the biodiversity of the solfatara Pisciarelli, Agnano (Naples, Italy), to discover hyperthermophilic carbohydrate-active enzymes (CAZymes) and to characterize the entire set of such enzymes in this environment (CAZome). Here, we report the results of the metagenomic analysis of two mud/water pools that greatly differ in both temperature and pH (T = 85 °C and pH 5.5; T = 92 °C and pH 1.5, for Pool1 and Pool2, respectively). DNA deep sequencing and following in silico analysis led to 14 934 and 17 652 complete ORFs in Pool1 and Pool2, respectively. They exclusively belonged to archaeal cells and viruses with great genera variance within the phylum Crenarchaeota, which reflected the difference in temperature and pH of the two Pools. Surprisingly, 30% and 62% of all of the reads obtained from Pool1 and 2, respectively, had no match in nucleotide databanks. Genes associated with carbohydrate metabolism were 15% and 16% of the total in the two Pools, with 278 and 308 putative CAZymes in Pool1 and 2, corresponding to ~ 2.0% of all ORFs. Biochemical characterization of two CAZymes of a previously unknown archaeon revealed a novel subfamily GH5_19 β-mannanase/β-1,3-glucanase whose hemicellulose specificity correlates with the vegetation surrounding the sampling site, and a novel NAD⁺ -dependent GH109 with a previously unreported β-N-acetylglucosaminide/β-glucoside specificity. DATABASES: The sequencing reads are available in the NCBI Sequence Read Archive (SRA) database under the accession numbers SRR7545549 (Pool1) and SRR7545550 (Pool2). The sequences of GH5_Pool2 and GH109_Pool2 are available in GenBank database under the accession numbers MK869723 and MK86972, respectively. The environmental data relative to Pool1 and Pool2 (NCBI BioProject PRJNA481947) are available in the Biosamples database under the accession numbers SAMN09692669 (Pool1) and SAMN09692670 (Pool2).

Effects of Random Mutagenesis and In Vivo Selection on the Specificity and Stability of a Thermozyme

Factors that give enzymes stability, activity, and substrate recognition result from the combination of few weak molecular interactions, which can be difficult to study through rational protein engineering approaches. We used irrational random mutagenesis and in vivo selection to test if a β-glycosidase from the thermoacidophile Saccharolobus solfataricus (Ssβ-gly) could complement an Escherichia coli strain unable to grow on lactose. The triple mutant of Ssβ-gly (S26L, P171L, and A235V) was more active than the wild type at 85 °C, inactivated at this temperature almost 300-fold quicker, and showed a 2-fold higher kcat on galactosides. The three mutations, which were far from the active site, were analyzed to test their effect at the structural level. Improved activity on galactosides was induced by the mutations. The S26L and P171L mutations destabilized the enzyme through the removal of a hydrogen bond and increased flexibility of the peptide backbone, respectively. However, the flexibility added by S26L mutation improved the activity at T > 60 °C. This study shows that random mutagenesis and biological selection allowed the identification of residues that are critical in determining thermal activity, stability, and substrate recognition.

Identification and characterization of a core fucosidase from the bacterium Elizabethkingia meningoseptica

All reported α-l-fucosidases catalyze the removal of nonreducing terminal l-fucoses from oligosaccharides or their conjugates, while having no capacity to hydrolyze core fucoses in glycoproteins directly. Here, we identified an α-fucosidase from the bacterium Elizabethkingia meningoseptica with catalytic activity against core α-1,3-fucosylated substrates, and we named it core fucosidase I (cFase I). Using site-specific mutational analysis, we found that three acidic residues (Asp-242, Glu-302, and Glu-315) in the predicted active pocket are critical for cFase I activity, with Asp-242 and Glu-315 acting as a pair of classic nucleophile and acid/base residues and Glu-302 acting in an as yet undefined role. These findings suggest a catalytic mechanism for cFase I that is different from known α-fucosidase catalytic models. In summary, cFase I exhibits glycosidase activity that removes core α-1,3-fucoses from substrates, suggesting cFase I as a new tool for glycobiology, especially for studies of proteins with core fucosylation.

Engineering the stability and the activity of a glycoside hydrolase

Glycosidases, the enzymes responsible in nature for the catabolism of carbohydrates, are well-studied catalysts widely used in industrial biotransformations and oligosaccharide synthesis, which are also attractive targets for drug development. Glycosidases from hyperthermophilic organisms (thriving at temperatures > 85 °C) are also interesting models to understand the molecular basis of protein stability and to produce robust tools for industrial applications. Here, we review the results obtained in the last two decades by our group on a β-glycosidase from the hyperthermophilic Archaeon Sulfolobus solfataricus. Our findings will be presented in the general context of the stability of proteins from hyperthermophiles and of the chemo-enzymatic synthesis of oligosaccharides.

Comparative study and mutational analysis of distinctive structural elements of hyperthermophilic enzymes

Comparison of the three-dimensional structure of hyperthermophilic and mesophilic beta-glycosidases shows differences in secondary structure composition. The enzymes from hyperthermophilic archaea have a significantly larger number of beta-strands arranged in supernumerary beta-sheets compared to mesophilic enzymes from bacteria and other organisms. Amino acid replacements designed to alter the structure of the supernumerary beta-strands were introduced by site directed mutagenesis into the sequence encoding the beta-glycosidase from Sulfolobus solfataricus. Most of the replacements caused almost complete loss of activity but some yielded enzyme variants whose activities were affected specifically at higher temperatures. Far-UV CD spectra recorded as a function of temperature for both wild type beta-glycosidase and mutant V349G, one of the mutants with reduced activity at higher temperatures, were similar, showing that the protein structure of the mutant was stable at the highest temperatures assayed. The properties of mutant V349G show a difference between thermostability (stability of the protein structure at high temperatures) and thermophilicity (optimal activity at high temperatures).

A comparative infrared spectroscopic study of glycoside hydrolases from extremophilic archaea revealed different molecular mechanisms of adaptation to high temperatures

The identification of the determinants of protein thermal stabilization is often pursued by comparing enzymes from hyperthermophiles with their mesophilic counterparts while direct structural comparisons among proteins and enzymes from hyperthermophiles are rather uncommon. Here, oligomeric beta-glycosidases from the hyperthermophilic archaea Sulfolobus solfataricus (Ss beta-gly), Thermosphaera aggregans (Ta beta-gly), and Pyrococcus furiosus (Pf beta-gly), have been compared. Studies of FTIR spectroscopy and kinetics of thermal inactivation showed that the three enzymes had similar secondary structure composition, but Ss beta-gly and Ta beta-gly (temperatures of melting 98.1 and 98.4 degrees C, respectively) were less stable than Pf beta-gly, which maintained its secondary structure even at 99.5 degrees C. The thermal denaturation of Pf beta-gly, followed in the presence of SDS, suggested that this enzyme is stabilized by hydrophobic interactions. A detailed inspection of the 3D-structures of these enzymes supported the experimental results: Ss beta-gly and Ta beta-gly are stabilized by a combination of ion-pairs networks and intrasubunit S-S bridges while the increased stability of Pf beta-gly resides in a more compact protein core. The different strategies of protein stabilization give experimental support to recent theories on thermophilic adaptation and suggest that different stabilization strategies could have been adopted among archaea.

Structural basis of the destabilization produced by an amino-terminal tag in the β-glycosidase from the hyperthermophilic archeon Sulfolobus solfataricus

We have previously shown that the major ion-pairs network of the tetrameric beta-glycosidase from the hyperthermophilic archeon Sulfolobus solfataricus involves more than 16 ion-pairs and hydrogen bonds between several residues from the four subunits and protects the protein from thermal unfolding by sewing the carboxy-termini of the enzyme. We show here that the amino-terminal of the enzyme also plays a relevant role in the thermostabilization of the protein. In fact, the addition of four extra amino acids at the amino-terminal of the beta-glycosidase, though not affecting the catalytic machinery of the enzyme and its thermophilicity, produced a faster enzyme inactivation in the temperature range 85-95 degrees C and decreased the Tm of the protein of 6 degrees C, measured by infrared spectroscopy. In addition, detailed two-dimensional IR correlation analysis revealed that the quaternary structure of the tagged enzyme is destabilized at 85 degrees C whilst that of the wild type enzyme is stable up to 98 degrees C. Molecular models allowed the rationalization of the experimental data indicating that the longer amino-terminal tail may destabilize the beta-glycosidase by enhancing the molecular fraying of the polypeptide and loosening the dimeric interfaces. The data support the hypothesis that fraying of the polypeptide chain termini is a relevant event in protein unfolding.

Two-dimensional IR correlation spectroscopy of mutants of the beta-glycosidase from the hyperthermophilic archaeon Sulfolobus solfataricus identifies the mechanism of quaternary structure stabilization and unravels the sequence of thermal unfolding events

Beta-glycosidase from the hyperthermophilic archaeon Sulfolobus solfataricus is a homotetramer with a higher number of ion pairs compared with mesophilic glycoside hydrolases. The ion pairs are arranged in large networks located mainly at the tetrameric interface of the molecule. In the present study, the structure and thermal stability of the wild-type beta-glycosidase and of three mutants in residues R488 and H489 involved in the C-terminal ionic network were examined by FTIR (Fourier-transform IR) spectroscopy. The FTIR data revealed small differences in the secondary structure of the proteins and showed a lower thermostability of the mutant proteins with respect to the wild-type. Generalized 2D-IR (two-dimensional IR correlation spectroscopy) at different temperatures showed different sequences of thermal unfolding events in the mutants with respect to the wild-type, indicating that punctual mutations affect the unfolding and aggregation process of the protein. A detailed 2D-IR analysis of synchronous maps of the proteins allowed us to identify the temperatures at which the ionic network that stabilizes the quaternary structure of the native and mutant enzymes at the C-terminal breaks down. This evidence gives support to the current theories on the mechanism of ion-pair stabilization in proteins from hyperthermophilic organisms.

Large improvement in the thermal stability of a tetrameric malate dehydrogenase by single point mutations at the dimer-dimer interface

The stability of tetrameric malate dehydrogenase from the green phototrophic bacterium Chloroflexus aurantiacus (CaMDH) is at least in part determined by electrostatic interactions at the dimer-dimer interface. Since previous studies had indicated that the thermal stability of CaMDH becomes lower with increasing pH, attempts were made to increase the stability by removal of (excess) negative charge at the dimer-dimer interface. Mutation of Glu165 to Gln or Lys yielded a dramatic increase in thermal stability at pH 7.5 (+23.6 -- + 23.9 degrees C increase in apparent t(m)) and a more moderate increase at pH 4.4 (+4.6 -- + 5.4 degrees C). The drastically increased stability at neutral pH was achieved without forfeiture of catalytic performance at low temperatures. The crystal structures of the two mutants showed only minor structural changes close to the mutated residues, and indicated that the observed stability effects are solely due to subtle changes in the complex network of electrostatic interactions in the dimer-dimer interface. Both mutations reduced the concentration dependency of thermal stability, suggesting that the oligomeric structure had been reinforced. Interestingly, the two mutations had similar effects on stability, despite the charge difference between the introduced side-chains. Together with the loss of concentration dependency, this may indicate that both E165Q and E165K stabilize CaMDH to such an extent that disruption of the inter-dimer electrostatic network around residue 165 no longer limits kinetic thermal stability.

X-ray structure of a membrane-bound β-glycosidase from the hyperthermophilic archaeon Pyrococcus horikoshii

The beta-glycosidase of the hyperthermophilic Archaeon Pyrococcus horikoshii is a membrane-bound enzyme with the preferred substrate of alkyl-beta-glycosides. In this study, the unusual structural features that confer the extreme thermostability and substrate preferences of this enzyme were investigated by X-ray crystallography and docking simulation. The enzyme was crystallized in the presence of a neutral surfactant, and the crystal structure was solved by the molecular replacement method and refined at 2.5 A. The main-chain fold of the enzyme belongs to the (betaalpha)8 barrel structure common to the Family 1 glycosyl hydrolases. The active site is located at the center of the C-termini of the barrel beta-strands. The deep pocket of the active site accepts one sugar unit, and a hydrophobic channel extending radially from there binds the nonsugar moiety of the substrate. The docking simulation for oligosaccharides and alkylglucosides indicated that alkylglucosides with a long aliphatic chain are easily accommodated in the hydrophobic channel. This sparingly soluble enzyme has a cluster of hydrophobic residues on its surface, situated at the distal end of the active site channel and surrounded by a large patch of positively charged residues. We propose that this hydrophobic region can be inserted into the membrane while the surrounding positively charged residues make favorable contacts with phosphate groups on the inner surface of the membrane. The enzyme could thus adhere to the membrane in the proximity of its glycolipid substrate.

Electrostatic interactions across the dimer-dimer interface contribute to the pH-dependent stability of a tetrameric malate dehydrogenase

Malate dehydrogenase (MDH) from the moderately thermophilic bacterium Chloroflexus aurantiacus (CaMDH) is a tetrameric enzyme, while MDHs from mesophilic bacteria usually are dimers. Using site-directed mutagenesis, we show here that a network of electrostatic interactions across the extra dimer-dimer interface in CaMDH is important for thermal stability and oligomeric integrity. Stability effects of single point mutations (E25Q, E25K, D56N, D56K) varied from -1.2 degrees C to -26.8 degrees C, and depended strongly on pH. Gel-filtration experiments indicated that the 26.8 degrees C loss in stability observed for the D56K mutant at low pH was accompanied by a shift towards a lower oligomerization state.

Identification of the catalytic nucleophile of the family 29 alpha-L-fucosidase from Sulfolobus solfataricus via chemical rescue of an inactive mutant

We have recently reported that a functional alpha-L-fucosidase could be expressed by a single insertional mutation in the region of overlap between the ORFs SSO11867 and SSO3060 of the hyperthermophilic Archaeon Sulfolobus solfataricus [Cobucci-Ponzano et al. J. Biol. Chem. (2003) 278, 14622-14631]. This enzyme, belonging to glycoside hydrolase family 29 (GH29), showed micromolar specificity for p-nitrophenyl-alpha-L-fucoside (pNp-Fuc) and promoted transfucosylation reactions by following a reaction mechanism in which the products retained the anomeric configuration of the substrate. The active site residues in GH29 enzymes are still unknown. We describe here the identification of the catalytic nucleophile of the reaction in the alpha-L-fucosidase from S. solfataricus by reactivation with sodium azide of the mutant Asp242Gly that shows a 10(3)-fold activity reduction on pNp-Fuc. The detailed stereochemical analysis of the fucosyl-azide produced by the mutant reactivated on pNp-Fuc revealed its inverted (beta-fucosyl azide) configuration compared with the substrate. This allows for the first time the unambiguous assignment of Asp242, and its homologous residues, as the nucleophilic catalytic residues of GH29 alpha-L-fucosidases. This is the first time that this approach is used for alpha-L-glycosidases, widening the applicability of this method.

New understandings of thermostable and peizostable enzymes

Recent large-scale studies illustrate the importance of electrostatic interactions near the surface of proteins as a major factor in enhancing thermal stability. Mutagenesis studies have also demonstrated the importance of optimized charge interactions on the surface of the protein, which can significantly augment enzyme thermal stability. Directed evolution studies show that increased stability may be obtained by different routes, which may not mimic those used by nature. Despite observations that some of the most thermotolerant organisms grow under conditions of high pressure, little effort has been made to understand the correlation between pressure and temperature stability. One recent study demonstrates that the active-site volume may be important in increasing pressure stability.

Identification of an archaeal alpha-L-fucosidase encoded by an interrupted gene. Production of a functional enzyme by mutations mimicking programmed -1 frameshifting

The analysis of the complete genome of the thermoacidophilic Archaeon Sulfolobus solfataricus revealed two open reading frames (ORF), named SSO11867 and SSO3060, interrupted by a -1 frameshift and encoding for the N- and the C-terminal fragments, respectively, of an alpha-l-fucosidase. We report here that these ORFs are actively transcribed in vivo, and we confirm the presence of the -1 frameshift between them at the cDNA level, explaining why we could not find alpha-fucosidase activity in S. solfataricus extracts. Detailed analysis of the region of overlap between the two ORFs revealed the presence of the consensus sequence for a programmed -1 frameshifting. Two specific mutations, mimicking this regulative frameshifting event, allow the expression, in Escherichia coli, of a fully active thermophilic and thermostable alpha-l-fucosidase (EC ) with micromolar substrate specificity and showing transfucosylating activity. The analysis of the fucosylated products of this enzyme allows, for the first time, assigning a retaining reaction mechanism to family 29 of glycosyl hydrolases. The presence of an alpha-fucosidase putatively regulated by programmed -1 frameshifting is intriguing both with respect to the regulation of gene expression and, in post-genomic era, for the definition of gene function in Archaea.