Structure-guided engineering of the substrate specificity of a fungal β-glucuronidase toward triterpenoid saponins

Efficient glycyrrhetinic acid biomanufacturing through protein engineering and dual-GUS combination strategy with novel β-glucuronidase from Aspergillus calidoustus CLH-22

Glycyrrhetinicacid (GA) is a high-value pentacyclic triterpenoid with broad applications. However, the industrial production of GA is hindered by low yield and the accumulation of the intermediate product GlycyrrhetinicAcid3-O-Mono-β-D-Glucuronide (GAMG). This study first identified a novel β-glucuronidase (AcGUS) from Aspergillus calidoustus CLH-22 through transcriptomic analysis, demonstrating a substrate preference for GAMG. Subsequently, mutant AcGUS3G461C/Q462H/I575K with significantly improved activity (kcat/Km of 11.02-fold) was obtained via computer-aided engineering. Furthermore, the dual-GUS combination strategy was employed for the first timeto construct engineered Pichia pastoris for GA production, offering multiple advantages of enhanced conversion efficiency and reduced fermentation viscosity. Finally, under systematically optimized conditions and employing Glycyrrhizin (GL) as the substrate, the final concentration of GA was 48.73 g/L with a conversion of 97.26 % in a 1000-L fermenter, representing the optimal biocatalytic performance reported to date. This study provides new ideas and insights for industrial GA production.

Astraoleanosides E-P, oleanane-type triterpenoid saponins from the aerial parts of Astragalus membranaceus Bunge and their β-glucuronidase inhibitory activity

Historically, Astragalus membranaceus Bunge has been used as a beneficial medicinal plant, particularly in the Asian traditional medical systems, for the treatment of various human diseases such as stomach ulcers, diarrhea, and respiratory issues associated with phlegm. In this study, a phytochemical characterization of the aerial parts of A. membranaceusled to the isolation of 29 oleanane-type triterpenoid saponins, including 11 new compounds named astraoleanosides E-P (6-9, 13, 14, 18-22), as well as 18 known ones. The structures of these compounds were elucidated using nuclear magnetic resonance (NMR) spectroscopy and high-resolution mass spectrometry. Among them, astraoleanoside H (9) and cloversaponin III (15) demonstrated the most potent β-glucuronidase inhibitory activities, with IC50 values of 21.20 ± 0.75 and 9.05 ± 0.47 µM, respectively, compared to the positive control d-saccharic acid 1,4-lactone (IC50 = 20.62 ± 1.61 µM). Enzyme kinetics studies were then conducted to investigate the type of inhibition exhibited by these active compounds. In addition, the binding mechanism, key interactions, binding stability, and dynamic behavior of protein-ligand complexes were investigated through in silico approaches, such as molecular docking and molecular dynamics simulations. These findings highlight the promising potential of triterpenoid saponins from A. membranaceus as lead compounds for β-glucuronidase inhibitors, offering new possibilities for the development of therapeutic agents targeting various diseases where β-glucuronidase plays a crucial role.

Site-specific crosslinking and assembly of tetrameric β-glucuronidase improve glycyrrhizin hydrolysis

In this study, eight nonconserved residues with exposed surfaces and flexible conformations of the homotetrameric PGUS (β-glucuronidase from Aspergillus oryzae Li-3) were identified. Single-point mutation into cysteine enabled the thiol-maleimide reaction and site-specific protein assembly using a two-arm polyethylene glycol (PEG)-maleimide crosslinker (Mal2 ). The Mal2 (1k) (with 1 kDa PEG spacer)-crosslinked PGUS assemblies showed low crosslinking efficiency and unimproved thermostability except for G194C-Mal2 (1k). To improve the crosslinking efficiency, a lengthened crosslinker Mal2 (2k) (with 2 kDa PEG spacer) was used to produce PGUS assembly and a highly improved thermostability was achieved with a half-life of 47.2-169.2 min at 70°C, which is 1.04-3.74 times that of wild type PGUS. It is found that the thermostability of PGUS assembly was closely associated with the formation of inter-tetramer assembly and intratetramer crosslinking, rather than the PEGylation of the enzyme. Therefore, the four-arm PEG-maleimide crosslinker Mal4 (2k) (with 2 kDa PEG spacer) was employed to simultaneously increase the inter-tetramer assembly and intratetramer crosslinking, and the resulting PGUS assemblies showed further improved thermostabilities compared with Mal2 (2k)-crosslinked assemblies. Finally, the application of PGUS assemblies with significantly improved thermostability to the bioconversion of GL proved that the PGUS assembly is a strong catalyst for glycyrrhizin (GL) hydrolysis in industrial applications.

Qualitative metabolomics-based characterization of a phenolic UDP-xylosyltransferase with a broad substrate spectrum from Lentinus brumalis

Wood-decaying fungi are the major decomposers of plant litter. Heavy sequencing efforts on genomes of wood-decaying fungi have recently been made due to the interest in their lignocellulolytic enzymes; however, most parts of their proteomes remain uncharted. We hypothesized that wood-decaying fungi would possess promiscuous enzymes for detoxifying antifungal phytochemicals remaining in the dead plant bodies, which can be useful biocatalysts. We designed a computational mass spectrometry-based untargeted metabolomics pipeline for the phenotyping of biotransformation and applied it to 264 fungal cultures supplemented with antifungal plant phenolics. The analysis identified the occurrence of diverse reactivities by the tested fungal species. Among those, we focused on O-xylosylation of multiple phenolics by one of the species tested, Lentinus brumalis. By integrating the metabolic phenotyping results with publicly available genome sequences and transcriptome analysis, a UDP-glycosyltransferase designated UGT66A1 was identified and validated as an enzyme catalyzing O-xylosylation with broad substrate specificity. We anticipate that our analytical workflow will accelerate the further characterization of fungal enzymes as promising biocatalysts.

Tetramerization of GH2 β-Glucuronidases is Essential for Catalyzing the Hydrolysis of the Large Substrate Glycyrrhizin

In this study, structural analysis was employed to identify three hotspot residues that contribute most to the tetramer formation of a glycoside hydrolase family 2 (GH2) β-glucuronidase (GUS) from Aspergillus oryzae Li-3. Single-point mutation at these sites completely disrupted the tetramer structure and abolished the glycyrrhizin (GL)-hydrolyzing activity. Then, the W522A dimer was refactored into a tetramer by disulfide bonding, and partial GL activity was restored. Further saturated mutation showed a strong correlation between the GL activity of the mutants and their tetramer ratios. Molecular simulations were employed to illustrate the critical role of the tetramer interface in maintaining a functional active-site structure. The three highly conserved tetramer-forming residues were finally applied to two other GH2 GUSs for tetramer dissociation and demonstrated the significance of the homotetramerization for GL-hydrolyzing activity of GH2 GUSs. This study lays foundation for engineering GL-hydrolyzing GUSs at the quaternary structure level for function regulations.

O-glycosyltransferases from Homo sapiens contributes to the biosynthesis of and Glycyrrhizin in Saccharomyces cerevisiae

Glycyrrhizin (GL) and Glycyrrhetic Acid 3-O-mono-β-D-glucuronide (GAMG) are the typical triterpenoid glycosides found in the root of licorice, a popular medicinal plant that exhibits diverse physiological effects and pharmacological manifestations. However, only few reports are available on the glycosylation enzymes involved in the biosynthesis of these valuable compounds with low conversion yield so far. In mammals, glycosyltransferases are involved in the phase II metabolism and may provide new solutions for us to engineer microbial strains to produce high valued compounds due to the substrate promiscuity of these glycosyltransferases. In this study, we mined the genomic databases of mammals and evaluated 22 candidate genes of O-glycosyltransferases by analyzing their catalytic potential for O-glycosylation of the native substrate, glycyrrhetinic acid (GA) for its glycodiversification. Out of 22 selected glycosyltransferases, only UGT1A1 exhibited high catalytic performance for biosynthesis of the key licorice compounds GL and GAMG. Molecular docking results proposed that the enzymatic activity of UGT1A1 was likely owing to the stable hydrogen bonding interactions and favorite conformations between the amino acid residues around substrate channels (P82~R85) and substrates. Furthermore, the complete biosynthesis pathway of GL was reconstructed in Saccharomyces cerevisiae for the first time, resulting in the production of 5.98 ± 0.47 mg/L and 2.31 ± 0.21 mg/L of GL and GAMG, respectively.

Improving the activity and thermostability of GH2 β-glucuronidases via domain reassembly

Glycoside hydrolase family 2 (GH2) enzymes are generally composed of three domains: TIM-barrel domain (TIM), immunoglobulin-like β-sandwich domain (ISD), and sugar-binding domain (SBD). The combination of these three domains yields multiple structural combinations with different properties. Theoretically, the drawbacks of a given GH2 fold may be circumvented by efficiently reassembling the three domains. However, very few successful cases have been reported. In this study, we used six GH2 β-glucuronidases (GUSs) from bacteria, fungi, or humans as model enzymes and constructed a series of mutants by reassembling the domains from different GUSs. The mutants PGUS-At, GUS-PAA, and GUS-PAP, with reassembled domains from fungal GUSs, showed improved expression levels, activity, and thermostability, respectively. Specifically, compared to the parental enzyme, the mutant PGUS-At displayed 3.8 times higher expression, the mutant GUS-PAA displayed 1.0 time higher catalytic efficiency (k_cat /K_m ), and the mutant GUS-PAP displayed 7.5 times higher thermostability at 65°C. Furthermore, two-hybrid mutants, GUS-AEA and GUS-PEP, were constructed with the ISD from a bacterial GUS and SBD and TIM domain from fungal GUSs. GUS-AEA and GUS-PEP showed 30.4% and 23.0% higher thermostability than GUS-PAP, respectively. Finally, molecular dynamics simulations were conducted to uncover the molecular reasons for the increased thermostability of the mutant.

Variants of glycosyl hydrolase family 2 β-glucuronidases have increased activity on recalcitrant substrates

Glucuronidated drug metabolites can be quantified from urine samples by first hydrolyzing conjugates with β-glucuronidase (β-GUS) and then separating free drug molecules by liquid chromatography and mass spectrometry detection (LC-MS). To improve the activity and specificity of various β-GUS, we designed enzyme chimeras and generated site-saturation variants based on structural analyses, then screened them for improved activity on drug metabolites important to clinical and forensic drug-testing laboratories. Often, an increase of activity on one substrate of interest was countered by loss of activity against another, and there was no strong correlation of activity on standard β-glucuronidase substrates to activity on recalcitrant drug glucuronides. However, we discovered a chimera of two enzymes from different species of Aspergillus that displays a 27 % increase in activity on morphine-3-glucuronide than the parent proteins. Furthermore, mutations in the M-loop, which is a loop near the active site, resulted in numerous variants with dramatically increased rates of hydrolysis on drug glucuronides. Specifically, the M-loop variant Q451D/A452E of a β-GUS from Brachyspira pilosicoli has a 50-fold and 25-fold increase in activity on the recalcitrant substrates codeine-6-glucuronide and dihydrocodeine-6-glucuronide, respectively, compared to the parent enzyme.

Bioengineering oligomerization and monomerization of enzymes: learning from natural evolution to matching the demands for industrial applications

It is generally accepted that oligomeric enzymes evolve from their monomeric ancestors, and the evolution process generates superior structural benefits for functional advantages. Furthermore, adjusting the transition between different oligomeric states is an important mechanism for natural enzymes to regulate their catalytic functions for adapting environmental fluctuations in nature, which inspires researchers to mimic such a strategy to develop artificially oligomerized enzymes through protein engineering for improved performance under specific conditions. On the other hand, transforming oligomeric enzymes into their monomers is needed in fundamental research for deciphering catalytic mechanisms as well as exploring their catalytic capacities for better industrial applications. In this article, strategies for developing artificially oligomerized and monomerized enzymes are reviewed and highlighted by their applications. Furthermore, advances in the computational prediction of oligomeric structures are introduced, which would accelerate the systematic design of oligomeric and monomeric enzymes. Finally, the current challenges and future directions in this field are discussed.

Construction of a CaHPO4-PGUS1 hybrid nanoflower through protein-inorganic self-assembly, and its application in glycyrrhetinic acid 3-O-mono-β-d-glucuronide preparation

Glycyrrhetinic acid 3-O-mono-β-d-glucuronide (GAMG), an important pharmaceutical intermediate and functional sweetener, has broad applications in the food and medical industries. A green and cost-effective method for its preparation is highly desired. Using site-directed mutagenesis, we previously obtained a variant of β-glucuronidase from Aspergillus oryzae Li-3 (PGUS1), which can specifically transform glycyrrhizin (GL) into GAMG. In this study, a facile method was established to prepare a CaHPO4-PGUS1 hybrid nanoflower for enzyme immobilization, based on protein-inorganic hybrid self-assembly. Under optimal conditions, 1.2 mg of a CaHPO4-PGUS1 hybrid nanoflower precipitate with 71.2% immobilization efficiency, 35.60 mg·g-1 loading capacity, and 118% relative activity was obtained. Confocal laser scanning microscope and scanning electron microscope results showed that the enzyme was encapsulated in the CaHPO4-PGUS1 hybrid nanoflower. Moreover, the thermostability of the CaHPO4-PGUS1 hybrid nanoflower at 55°C was improved, and its half-life increased by 1.3 folds. Additionally, the CaHPO4-PGUS1 hybrid nanoflower was used for the preparation of GAMG through GL hydrolysis, with the conversion rate of 92% in 8 h, and after eight consecutive runs, it had 60% of its original activity.

Tuning the pH profile of β-glucuronidase by rational site-directed mutagenesis for efficient transformation of glycyrrhizin

In this study, we aimed to shift the optimal pH of acidic β-glucuronidase from Aspergillus oryzae Li-3 (PGUS) to the neutral region by site-directed mutagenesis, thus allowing high efficient biotransformation of glycyrrhizin (GL) into glycyrrhetinic acid (GA) under higher pH where the solubility of GL could be greatly enhanced. Based on PGUS structure analysis, five critical aspartic acid and glutamic acid residues were replaced with arginine on the surface to generate a variant 5Rs with optimal pH shifting from 4.5 to 6.5. The catalytic efficiency (k_cat /K_m) value of 5Rs at pH 6.5 was 10.7-fold higher than that of PGUS wild-type at pH 6.5, even 1.4-fold higher than that of wild-type at pH 4.5. Molecular dynamics simulation was performed to explore the molecular mechanism for the shifted pH profile and enhanced pH stability of 5Rs.

Enhanced production, overexpression and characterization of a hyperthermophilic multimodular GH family 2 β‑glucuronidase (TpGUS) cloned from Thermotoga petrophila RKU-1T in a mesophilic host

A multimodular hyperthermophilic β‑glucuronidase (TpGUS) from Thermotoga petrophila RKU-1^T, belongs to glycoside hydrolase family 2 (GH2), was cloned and overexpressed in Escherichia coli BL21 CodonPlus (DE3)-RIPL. Expression and production of extracellular TpGUS was enhanced through various specific cultivation and induction strategies. Extracellular TpGUS activity was improved by 3.44 and 7 fold in 4 × ZB medium induced with 0.5 mM IPTG and 100 mM lactose, respectively. The enzyme was purified to homogeneity with a single band of 65.6 kDa on SDS-PAGE, using two subsequent steps of anion exchange and hydrophobic interaction chromatography after heat precipitation (70 °C, 1 h). Optimal activity of TpGUS was observed at 95 °C and pH 6.0; and it displayed prodigious thermal stability over a temperature range of 50-85 °C for 12 h at pH 6.0-7.5. K_m, V_max, V_maxK_m^-1, k_cat, and k_catK_m^-1 were calculated to be 0.7 mM, 227 mmol mg^-1 min^-1, 324.3 min^-1, 164,492.7 s^-1 and 234,989.6 mM^-1 s^-1, respectively using pNPGU as a substrate. Recombinant TpGUS exhibited favorable properties which make this a promising candidate for various biotechnological and pharmacological applications.

β-Glucuronidases of opportunistic bacteria are the major contributors to xenobiotic-induced toxicity in the gut

Gut bacterial β-D-glucuronidases (GUSs) catalyze the removal of glucuronic acid from liver-produced β-D-glucuronides. These reactions can have deleterious consequences when they reverse xenobiotic metabolism. The human gut contains hundreds of GUSs of variable sequences and structures. To understand how any particular bacterial GUS(s) contributes to global GUS activity and affects human health, the individual substrate preference(s) must be known. Herein, we report that representative GUSs vary in their ability to produce various xenobiotics from their respective glucuronides. To attempt to explain the distinct substrate preference, we solved the structure of a bacterial GUS complexed with coumarin-3-β-D-glucuronide. Comparisons of this structure with other GUS structures identified differences in loop 3 (or the α2-helix loop) and loop 5 at the aglycone-binding site, where differences in their conformations, hydrophobicities and flexibilities appear to underlie the distinct substrate preference(s) of the GUSs. Additional sequence, structural and functional analysis indicated that several groups of functionally related gut bacterial GUSs exist. Our results pinpoint opportunistic gut bacterial GUSs as those that cause xenobiotic-induced toxicity. We propose a structure-activity relationship that should allow both the prediction of the functional roles of GUSs and the design of selective inhibitors.

Computation-Aided Rational Deletion of C-Terminal Region Improved the Stability, Activity, and Expression Level of GH2 β-Glucuronidase

In this study, computation-aided design on the basis of structural analysis was employed to rationally identify a highly dynamic C-terminal region that regulates the stability, expression level, and activity of a GH2 fungal glucuronidase from Aspergillus oryzae Li-3 (PGUS). Then, four mutants with a precisely truncated C-terminal region in different lengths were constructed; among them, mutant D591-604 with a 3.8-fold increase in half-life at 65 °C and a 6.8 kJ/mol increase in Gibbs free energy showed obviously improved kinetic and thermodynamic stability in comparison to PGUS. Mutants D590-604 and D591-604 both showed approximately 2.4-fold increases in the catalytic efficiency k_cat/ K_m and 1.8-fold increases in the expression level. Additionally, the expression level of PGUS was doubled through a C-terminal region swap with bacterial GUS from E. coli (EGUS). Finally, the robust PGUS mutants D590-604 and D591-604 were applied in the preparation of glycyrrhetinic acid with 4.0- and 4.4-fold increases in concentration through glycyrrhizin hydrolysis by a fed-batch process.

A Novel β-Glucuronidase from Talaromyces pinophilus Li-93 Precisely Hydrolyzes Glycyrrhizin into Glycyrrhetinic Acid 3-O-Mono-β-d-Glucuronide

Glycyrrhetinic acid 3-O-mono-β-d-glucuronide (GAMG), which possesses a higher sweetness and stronger pharmacological activity than those of glycyrrhizin (GL), can be obtained by removal of the distal glucuronic acid (GlcA) from GL. In this study, we isolated a β-glucuronidase (TpGUS79A) from the filamentous fungus Talaromyces pinophilus Li-93 that can specifically and precisely convert GL to GAMG without the formation of the by-product glycyrrhetinic acid (GA) from the further hydrolysis of GAMG. First, TpGUS79A was purified and identified through matrix-assisted laser desorption ionization-tandem time of flight mass spectrometry (MALDI-TOF-TOF MS) and deglycosylation, indicating that TpGUS79A is a highly N-glycosylated monomeric protein with a molecular mass of around 85 kDa, including around 25 kDa of glycan moiety. The gene for TpGUS79A was then cloned and verified by heterologous expression in Pichia pastoris TpGUS79A belonged to glycoside hydrolase family 79 (GH79) but shared low amino acid sequence identity (<35%) with the available GH79 GUS enzymes. TpGUS79A had strict specificity toward the glycan moiety but poor specificity toward the aglycone moiety. Interestingly, TpGUS79A recognized and hydrolyzed the distal glucuronic bond of GL but could not cleave the glucuronic bond in GAMG. TpGUS79A showed a much higher catalytic efficiency on GL (k_cat/K_m of 11.14 mM^-1 s^-1) than on the artificial substrate pNP β-glucopyranosiduronic acid (k_cat/K_m of 0.01 mM^-1 s^-1), which is different from the case for most GUSs. Homology modeling, substrate docking, and sequence alignment were employed to identify the key residues for substrate recognition. Finally, a fed-batch fermentation in a 150-liter fermentor was established to prepare GAMG through GL hydrolysis by T. pinophilus Li-93. Therefore, TpGUS79A is potentially a powerful biocatalyst for environmentally friendly and cost-effective production of GAMG.IMPORTANCE Compared to chemical methods, the biotransformation of glycyrrhizin (GL) into glycyrrhetinic acid 3-O-mono-β-d-glucuronide (GAMG), which has a higher sweetness and stronger pharmacological activity than those of GL, via catalysis by β-glucuronidase is an environmentally friendly approach due to the mild reaction conditions and the high yield of GAMG. However, currently available GUSs show low substrate specificity toward GL and further hydrolyze GAMG to glycyrrhetinic acid (GA) as a by-product, increasing the difficulty of subsequent separation and purification. In the present study, we succeeded in isolating a novel β-glucuronidase (named TpGUS79A) from Talaromyces pinophilus Li-93 that specifically hydrolyzes GL to GAMG without the formation of GA. TpGUS79A also shows higher activity on GL than those of the previously characterized GUSs. Moreover, the gene for TpGUS79A was cloned and its function verified by heterologous expression in P. pastoris Therefore, TpGUS79A can serve as a powerful biocatalyst for the cost-effective production of GAMG through GL transformation.

Design of Glyco-Linkers at Multiple Structural Levels to Modulate Protein Stability

N-glycosylation has critical roles in regulating protein stability, but the molecular basis is poorly understood. In this study, we integrated experimental and computational techniques to investigate the mechanism by which full-length N-glycans modulate protein stability from quaternary structure perspective. We found the two inherent N-glycans of β-glucuronidase expressed in Pichia pastoris function as "glyco-linkers" that hold spatially proximal motifs together to compact the local protein structure. We further designed and placed glyco-linkers in the unusual form of glyco-bridge and glyco-hairpin at the interfaces between domains and monomers with higher structural level, respectively, which conferred dramatically higher kinetic stability and thermodynamic stability than the inherent N-glycans. Our study not only provides unique insight into the interactions between glycans and proteins from a quaternary structure perspective but also facilitates the rational design of N-glycans as general tools that can enhance protein stability.