Rational Enzyme Design without Structural Knowledge: A Sequence-Based Approach for Efficient Generation of Transglycosylases

Unlocking the potential of enzyme engineering via rational computational design strategies

Enzymes play a pivotal role in various industries by enabling efficient, eco-friendly, and sustainable chemical processes. However, the low turnover rates and poor substrate selectivity of enzymes limit their large-scale applications. Rational computational enzyme design, facilitated by computational algorithms, offers a more targeted and less labor-intensive approach. There has been notable advancement in employing rational computational protein engineering strategies to overcome these issues, it has not been comprehensively reviewed so far. This article reviews recent developments in rational computational enzyme design, categorizing them into three types: structure-based, sequence-based, and data-driven machine learning computational design. Case studies are presented to demonstrate successful enhancements in catalytic activity, stability, and substrate selectivity. Lastly, the article provides a thorough analysis of these approaches, highlights existing challenges and potential solutions, and offers insights into future development directions.

Synthesis of fucosyllactose using α-L-fucosidases GH29 from infant gut microbial metagenome

Fucosyl-oligosaccharides (FUS) provide many health benefits to breastfed infants, but they are almost completely absent from bovine milk, which is the basis of infant formula. Therefore, there is a growing interest in the development of enzymatic transfucosylation strategies for the production of FUS. In this work, the α-L-fucosidases Fuc2358 and Fuc5372, previously isolated from the intestinal bacterial metagenome of breastfed infants, were used to synthesize fucosyllactose (FL) by transfucosylation reactions using p-nitrophenyl-α-L-fucopyranoside (pNP-Fuc) as donor and lactose as acceptor. Fuc2358 efficiently synthesized the major fucosylated human milk oligosaccharide (HMO) 2'-fucosyllactose (2'FL) with a 35% yield. Fuc2358 also produced the non-HMO FL isomer 3'-fucosyllactose (3'FL) and traces of non-reducing 1-fucosyllactose (1FL). Fuc5372 showed a lower transfucosylation activity compared to Fuc2358, producing several FL isomers, including 2'FL, 3'FL, and 1FL, with a higher proportion of 3'FL. Site-directed mutagenesis using rational design was performed to increase FUS yields in both α-L-fucosidases, based on structural models and sequence identity analysis. Mutants Fuc2358-F184H, Fuc2358-K286R, and Fuc5372-R230K showed a significantly higher ratio between 2'FL yields and hydrolyzed pNP-Fuc than their respective wild-type enzymes after 4 h of transfucosylation. The results with the Fuc2358-F184W and Fuc5372-W151F mutants showed that the residues F184 of Fuc2358 and W151 of Fuc5372 could have an effect on transfucosylation regioselectivity. Interestingly, phenylalanine increases the selectivity for α-1,2 linkages and tryptophan for α-1,3 linkages. These results give insight into the functionality of the active site amino acids in the transfucosylation activity of the GH29 α-L-fucosidases Fuc2358 and Fuc5372. KEY POINTS: Two α-L-fucosidases from infant gut bacterial microbiomes can fucosylate glycans Transfucosylation efficacy improved by tailored point-mutations in the active site F184 of Fuc2358 and W151 of Fuc5372 seem to steer transglycosylation regioselectivity.

Improved Enzymatic Production of the Fucosylated Human Milk Oligosaccharide LNFP II with GH29B α-1,3/4-l-Fucosidases

Five GH29B α-1,3/4-l-fucosidases (EC 3.2.1.111) were investigated for their ability to catalyze the formation of the human milk oligosaccharide lacto-N-fucopentaose II (LNFP II) from lacto-N-tetraose (LNT) and 3-fucosyllactose (3FL) via transglycosylation. We studied the effect of pH on transfucosylation and hydrolysis and explored the impact of specific mutations using molecular dynamics simulations. LNFP II yields of 91 and 65% were obtained for the wild-type SpGH29C and CpAfc2 enzymes, respectively, being the highest LNFP II transglycosylation yields reported to date. BbAfcB and BiAfcB are highly hydrolytic enzymes. The results indicate that the effects of pH and buffer systems are enzyme-dependent yet relevant to consider when designing transglycosylation reactions. Replacing Thr284 in BiAfcB with Val resulted in increased transglycosylation yields, while the opposite replacement of Val258 in SpGH29C and Val289 CpAfc2 with Thr decreased the transfucosylation, confirming a role of Thr and Val in controlling the flexibility of the acid/base loop in the enzymes, which in turn affects transglycosylation. The substitution of an Ala residue with His almost abolished secondary hydrolysis in CpAfc2 and BbAfcB. The results are directly applicable in the enhancement of transglycosylation and may have significant implications for manufacturing of LNFP II as a new infant formula ingredient.

Characterization of a GH20 β-N-Acetylhexosaminidase from Flavobacterium algicola Suitable to Synthesize Lacto-N-triose II

β-N-Acetylhexosaminidases have attracted much attention in the enzymatic synthesis of lacto-N-triose II (LNT2) as a backbone precursor of human milk oligosaccharides (HMOs). In this study, a novel glycoside hydrolase (GH) 20 family β-N-acetylhexosaminidase, FlaNag2353, from Flavobacterium algicola was biochemically characterized and applied to synthesize LNT2. FlaNag2353 displayed optimal activity to p-nitrophenyl N-acetyl-β-d-glucosaminide (pNP-GlcNAc) at 40 °C and pH 8.0. In addition to its excellent hydrolysis activity toward pNP-GlcNAc and chitooligosaccharides, FlaNag2353 showed trans-glycosylation activity. Under conditions of pH 9.0 and 55 °C for 2 h and utilizing 200 mM lactose and 10 mM pNP-GlcNAc, FlaNag2353 synthesized LNT2 with a conversion ratio of 4.15% calculated from pNP-GlcNAc. Moreover, when applied to LNT2 synthesis with 10 mM pNP-GlcNAc and 9.7% (w/v) industrial waste whey powder, FlaNag2353 achieved a conversion ratio of 2.39%. This study has significant implications for broadening the applications of GH20 β-N-acetylhexosaminidases and promoting the high-value utilization of whey powder.

Strategies for modulating transglycosylation activity, substrate specificity, and product polymerization degree of engineered transglycosylases

Glycosides are widely used in many fields due to their favorable biological activity. The traditional plant extractions and chemical methods for glycosides production are limited by environmentally unfriendly, laborious protecting group strategies and low yields. Alternatively, enzymatic glycosylation has drawn special attention due to its mild reaction conditions, high catalytic efficiency, and specific stereo-/regioselectivity. Glycosyltransferases (GTs) and retaining glycoside hydrolases (rGHs) are two major enzymes for the formation of glycosidic linkages. Therein GTs generally use nucleotide phosphate activated donors. In contrast, GHs can use broader simple and affordable glycosyl donors, showing great potential in industrial applications. However, most rGHs mainly show hydrolysis activity and only a few rGHs, namely non-Leloir transglycosylases (TGs), innately present strong transglycosylation activities. To address this problem, various strategies have recently been developed to successfully tailor rGHs to alleviate their hydrolysis activity and obtain the engineered TGs. This review summarizes the current modification strategies in TGs engineering, with a special focus on transglycosylation activity enhancement, substrate specificity modulation, and product polymerization degree distribution, which provides a reference for exploiting the transglycosylation potentials of rGHs.

Lignocellulosic biomass-based glycoconjugates for diverse biotechnological applications

Glycoconjugates are the ubiquitous components of mammalian cells, mainly synthesized by covalent bonds of carbohydrates to other biomolecules such as proteins and lipids, with a wide range of potential applications in novel vaccines, therapeutic peptides and antibodies (Ab). Considering the emerging developments in glycoscience, renewable production of glycoconjugates is of importance and lignocellulosic biomass (LCB) is a potential source of carbohydrates to produce synthetic glycoconjugates in a sustainable pathway. In this review, recent advances in glycobiology aiming on glycoconjugates production is presented together with the recent and cutting-edge advances in the therapeutic properties and application of glycoconjugates, including therapeutic glycoproteins, glycosaminoglycans (GAGs), and nutraceuticals, emphasizing the integral role of glycosylation in their function and efficacy. Special emphasis is given towards the potential exploration of carbon neutral feedstocks, in which LCB has an emerging role. Techniques for extraction and recovery of mono- and oligosaccharides from LCB are critically discussed and influence of the heterogeneous nature of the feedstocks and different methods for recovery of these sugars in the development of the customized glycoconjugates is explored. Although reports on the use of LCB for the production of glycoconjugates are scarce, this review sets clear that the potential of LCB as a source for the production of valuable glycoconjugates cannot be underestimated and encourages that future research should focus on refining the existing methodologies and exploring new approaches to fully realize the potential of LCB in glycoconjugate production.

Conformational Switch of the 250s Loop Enables the Efficient Transglycosylation in GH Family 77

Amylomaltases (AMs) play important roles in glycogen and maltose metabolism. However, the molecular mechanism is elusive. Here, we investigated the conformational dynamics of the 250s loop and catalytic mechanism of Thermus aquaticus TaAM using path-metadynamics and QM/MM MD simulations. The results demonstrate that the transition of the 250s loop from an open to closed conformation promotes polysaccharide sliding, leading to the ideal positioning of the acid/base. Furthermore, the conformational dynamics can also modulate the selectivity of hydrolysis and transglycosylation. The closed conformation of the 250s loop enables the tight packing of the active site for transglycosylation, reducing the energy penalty and efficiently preventing the penetration of water into the active site. Conversely, the partially closed conformation for hydrolysis results in a loosely packed active site, destabilizing the transition state. These computational findings guide mutation experiments and enable the identification of mutants with an improved disproportionation/hydrolysis ratio. The present mechanism is in line with experimental data, highlighting the critical role of conformational dynamics in regulating the catalytic reactivity of GHs.

The Engineered Hexosaminidase TtOGA-D120N Is an Efficient O-/N-/S-Glycoligase That Also Catalyzes Formation and Release of Oxazoline Donors for Cascade Syntheses with Glycosynthases or Transglycosylases

Engineering glycoside hydrolases is a major route to obtaining catalysts forming glycosidic bonds. Glycosynthases, thioglycoligases, and transglycosylases represent the main strategies, each having advantages and drawbacks. Here, we show that an engineered enzyme from the GH84 family, the acid-base mutant TtOGA-D120N, is an efficient O-, N-, and S-glycoligase, able to use Ssp3, Osp3, Nsp2, and Nsp nucleophiles. Moreover, TtOGA-D120N catalyzes the formation and release of N-acetyl-d-glucosamine 1,2-oxazoline, the intermediate of hexosaminidases displaying substrate-assisted catalysis. This release of an activated intermediate allows cascade synthesis by combination with transglycosylases or glycosynthases, here exemplified by synthesis of the human milk oligosaccharide lacto-N-triose II.

Second-Shell Residues Contribute to Catalysis by Predominately Preorganizing the Apo State in PafA

Residues beyond the first coordination shell are often observed to make considerable cumulative contributions in enzymes. Due to typically indirect perturbations of multiple physicochemical properties of the active site, however, their individual and specific roles in enzyme catalysis and disease-causing mutations remain difficult to predict and understand at the molecular level. Here we analyze the contributions of several second-shell residues in phosphate-irrepressible alkaline phosphatase of flavobacterium (PafA), a representative system as one of the most efficient enzymes. By adopting a multifaceted approach that integrates quantum-mechanical/molecular-mechanical free energy computations, molecular-mechanical molecular dynamics simulations, and density functional theory cluster model calculations, we probe the rate-limiting phosphoryl transfer step and structural properties of all relevant enzyme states. In combination with available experimental data, our computational results show that mutations of the studied second-shell residues impact catalytic efficiency mainly by perturbation of the apo state and therefore substrate binding, while they do not affect the ground state or alter the nature of phosphoryl transfer transition state significantly. Several second-shell mutations also modulate the active site hydration level, which in turn influences the energetics of phosphoryl transfer. These mechanistic insights also help inform strategies that may improve the efficiency of enzyme design and engineering by going beyond the current focus on the first coordination shell.

Exploration of the Transglycosylation Activity of Barley Limit Dextrinase for Production of Novel Glycoconjugates

A few α-glucan debranching enzymes (DBEs) of the large glycoside hydrolase family 13 (GH13), also known as the α-amylase family, have been shown to catalyze transglycosylation as well as hydrolysis. However, little is known about their acceptor and donor preferences. Here, a DBE from barley, limit dextrinase (HvLD), is used as a case study. Its transglycosylation activity is studied using two approaches; (i) natural substrates as donors and different p-nitrophenyl (pNP) sugars as well as different small glycosides as acceptors, and (ii) α-maltosyl and α-maltotriosyl fluorides as donors with linear maltooligosaccharides, cyclodextrins, and GH inhibitors as acceptors. HvLD showed a clear preference for pNP maltoside both as acceptor/donor and acceptor with the natural substrate pullulan or a pullulan fragment as donor. Maltose was the best acceptor with α-maltosyl fluoride as donor. The findings highlight the importance of the subsite +2 of HvLD for activity and selectivity when maltooligosaccharides function as acceptors. However, remarkably, HvLD is not very selective when it comes to aglycone moiety; different aromatic ring-containing molecules besides pNP could function as acceptors. The transglycosylation activity of HvLD can provide glycoconjugate compounds with novel glycosylation patterns from natural donors such as pullulan, although the reaction would benefit from optimization.

Plant glycosyltransferases for expanding bioactive glycoside diversity

Glycosylation is a successful strategy to alter the pharmacological properties of small molecules, and it has emerged as a unique approach to expand the chemical space of natural products that can be explored in drug discovery. Traditionally, most glycosylation events have been carried out chemically, often requiring many protection and deprotection steps to achieve a target molecule. Enzymatic glycosylation by glycosyltransferases could provide an alternative strategy for producing new glycosides. In particular, the glycosyltransferase family has greatly expanded in plants, representing a rich enzymatic resource to mine and expand the diversity of glycosides with novel bioactive properties. This article highlights previous and prospective uses for plant glycosyltransferases in generating bioactive glycosides and altering their pharmacological properties.

Data-driven enzyme engineering to identify function-enhancing enzymes

Identifying function-enhancing enzyme variants is a 'holy grail' challenge in protein science because it will allow researchers to expand the biocatalytic toolbox for late-stage functionalization of drug-like molecules, environmental degradation of plastics and other pollutants, and medical treatment of food allergies. Data-driven strategies, including statistical modeling, machine learning, and deep learning, have largely advanced the understanding of the sequence-structure-function relationships for enzymes. They have also enhanced the capability of predicting and designing new enzymes and enzyme variants for catalyzing the transformation of new-to-nature reactions. Here, we reviewed the recent progresses of data-driven models that were applied in identifying efficiency-enhancing mutants for catalytic reactions. We also discussed existing challenges and obstacles faced by the community. Although the review is by no means comprehensive, we hope that the discussion can inform the readers about the state-of-the-art in data-driven enzyme engineering, inspiring more joint experimental-computational efforts to develop and apply data-driven modeling to innovate biocatalysts for synthetic and pharmaceutical applications.

Building mutational bridges between carbohydrate-active enzymes

The commercial value of specialty carbohydrates and glycosylated compounds has sparked considerable interest in the synthetic potential of carbohydrate-active enzymes (CAZymes). Protein engineering methods have proven to be highly successful in expanding the range of glycosylation reactions that these enzymes can perform efficiently and cost-effectively. The past few years have witnessed meaningful progress in this area, largely due to a sharper focus on the understanding of structure-function relationships and mechanistic intricacies. Here, we summarize recent studies that demonstrate how protein engineers have become much better at traversing the fitness landscape of CAZymes through mutational bridges that connect the different activity types.

A Single Hydrogen Bond Controls the Selectivity of Transglycosylation vs Hydrolysis in Family 13 Glycoside Hydrolases

Converting glycoside hydrolases (GHs) from hydrolytic to synthetic enzymes via transglycosylation is a long-standing goal for the biosynthesis of complex carbohydrates. However, the molecular determinants for the selectivity of transglycosylation (T) vs hydrolysis (H) are still not fully unraveled. Herein, we show experimentally that mutation of one active site residue can switch the enzyme activity between hydrolysis and transglycosylation in two highly homologous GHs. Further QM/MM simulations reveal that the mutation modulates the T vs H reaction barriers via the presence/absence of a single H-bond with the nucleophile Asp. Such a H-bond controls the product selectivity via a dual effect: on one hand, it facilitates the breaking of the glycosyl-enzyme intermediate. On the other, it displaces the sugar acceptor, resulting in a reduced affinity and significant steric repulsion for transglycosylation. These findings expand our understanding of the molecular mechanisms that modulate the T/H balance in GHs.

A Comparison of the Transglycosylation Capacity between the Guar GH27 Aga27A and Bacteroides GH36 BoGal36A α-Galactosidases

The transglycosylation behavior and capacity of two clan GH-D α-galactosidases, BoGal36A from the gut bacterium Bacteroides ovatus and Aga27A from the guar plant, was investigated and compared. The enzymes were screened for the ability to use para-nitrophenyl-α-galactoside (pNP-Gal), raffinose and locust bean gum (LBG) galactomannan as glycosyl donors with the glycosyl acceptors methanol, propanol, allyl alcohol, propargyl alcohol and glycerol using mass spectrometry. Aga27A was, in general, more stable in the presence of the acceptors. HPLC analysis was developed and used as a second screening method for reactions using raffinose or LBG as a donor substrate with methanol, propanol and glycerol as acceptors. Time-resolved reactions were set up with raffinose and methanol as the donor and acceptor, respectively, in order to develop an insight into the basic transglycosylation properties, including the ratio between the rate of transglycosylation (methyl galactoside synthesis) and rate of hydrolysis. BoGal36A had a somewhat higher ratio (0.99 compared to 0.71 for Aga27A) at early time points but was indicated to be more prone to secondary (product) hydrolysis in prolonged incubations. The methyl galactoside yield was higher when using raffinose (48% for BoGal36A and 38% for Aga27A) compared to LBG (27% for BoGal36A and 30% for Aga27A).

Engineered Glycosidases for the Synthesis of Analogs of Human Milk Oligosaccharides

Enzymatic synthesis is an elegant biocompatible approach to complex compounds such as human milk oligosaccharides (HMOs). These compounds are vital for healthy neonatal development with a positive impact on the immune system. Although HMOs may be prepared by glycosyltransferases, this pathway is often complicated by the high price of sugar nucleotides, stringent substrate specificity, and low enzyme stability. Engineered glycosidases (EC 3.2.1) represent a good synthetic alternative, especially if variations in the substrate structure are desired. Site-directed mutagenesis can improve the synthetic process with higher yields and/or increased reaction selectivity. So far, the synthesis of human milk oligosaccharides by glycosidases has mostly been limited to analytical reactions with mass spectrometry detection. The present work reveals the potential of a library of engineered glycosidases in the preparative synthesis of three tetrasaccharides derived from lacto-N-tetraose (Galβ4GlcNAcβ3Galβ4Glc), employing sequential cascade reactions catalyzed by β3-N-acetylhexosaminidase BbhI from Bifidobacterium bifidum, β4-galactosidase BgaD-B from Bacillus circulans, β4-N-acetylgalactosaminidase from Talaromyces flavus, and β3-galactosynthase BgaC from B. circulans. The reaction products were isolated and structurally characterized. This work expands the insight into the multi-step catalysis by glycosidases and shows the path to modified derivatives of complex carbohydrates that cannot be prepared by standard glycosyltransferase methods.

Computer Simulation to Rationalize “Rational” Engineering of Glycoside Hydrolases and Glycosyltransferases

Glycoside hydrolases and glycosyltransferases are the main classes of enzymes that synthesize and degrade carbohydrates, molecules essential to life that are a challenge for classical chemistry. As such, considerable efforts have been made to engineer these enzymes and make them pliable to human needs, ranging from directed evolution to rational design, including mechanism engineering. Such endeavors fall short and are unreported in numerous cases, while even success is a necessary but not sufficient proof that the chemical rationale behind the design is correct. Here we review some of the recent work in CAZyme mechanism engineering, showing that computational simulations are instrumental to rationalize experimental data, providing mechanistic insight into how native and engineered CAZymes catalyze chemical reactions. We illustrate this with two recent studies in which (i) a glycoside hydrolase is converted into a glycoside phosphorylase and (ii) substrate specificity of a glycosyltransferase is engineered toward forming O-, N-, or S-glycosidic bonds.

Improvement of the Transglycosylation Efficiency of a Lacto-N-Biosidase from Bifidobacterium bifidum by Protein Engineering

The lacto-N-biosidase LnbB from Bifidobacterium bifidum JCM 1254 was engineered to improve its negligible transglycosylation efficiency with the purpose of enzymatically synthesizing lacto-N-tetraose (LNT; Gal-β1,3-GlcNAc-β1,3-Gal-β1,4-Glc) in one enzymatic step. LNT is a prebiotic human milk oligosaccharide in itself and constitutes the structural core of a range of more complex human milk oligosaccharides as well. Thirteen different LnbB variants were expressed and screened for transglycosylation activity by monitoring transglycosylation product formation using lacto-N-biose 1,2-oxazoline as donor substrate and lactose as acceptor substrate. LNT was the major reaction product, yet careful reaction analysis revealed the formation of three additional LNT isomers, which we identified to have a β1,2-linkage, a β1,6-linkage, and a 1,1-linkage, respectively, between lacto-N-biose (Gal-β1,3-GlcNAc) and lactose. Considering both maximal transglycosylation yield and regioselectivity as well as minimal product hydrolysis, the best variant was LnbB W394H, closely followed by W465H and Y419N. A high transglycosylation yield was also obtained with W394F, yet the substitution of W394 and W465 of the subsite −1 hydrophobic platform in the enzyme with His dramatically impaired the undesirable product hydrolysis as compared to substitution with Phe; the effect was most pronounced for W465. Using p-nitrophenyl-β-lacto-N-bioside as donor substrate manifested W394 as an important target position. The optimization of the substrate concentrations confirmed that high initial substrate concentration and high acceptor-to-donor ratio both favor transglycosylation.

Modulation of the Activity and Regioselectivity of a Glycosidase: Development of a Convenient Tool for the Synthesis of Specific Disaccharides

The synthesis of disaccharides, particularly those containing hexofuranoside rings, requires a large number of steps by classical chemical means. The use of glycosidases can be an alternative to limit the number of steps, as they catalyze the formation of controlled glycosidic bonds starting from simple and easy to access building blocks; the main drawbacks are the yields, due to the balance between the hydrolysis and transglycosylation of these enzymes, and the enzyme-dependent regioselectivity. To improve the yield of the synthesis of β-d-galactofuranosyl-(1→X)-d-mannopyranosides catalyzed by an arabinofuranosidase, in this study we developed a strategy to mutate, then screen the catalyst, followed by a tailored molecular modeling methodology to rationalize the effects of the identified mutations. Two mutants with a 2.3 to 3.8-fold increase in transglycosylation yield were obtained, and in addition their accumulated regioisomer kinetic profiles were very different from the wild-type enzyme. Those differences were studied in silico by docking and molecular dynamics, and the methodology revealed a good predictive quality in regards with the regioisomer profiles, which is in good agreement with the experimental transglycosylation kinetics. So, by engineering CtAraf51, new biocatalysts were enabled to obtain the attractive central motif from the Leishmania lipophosphoglycan core with a higher yield and regioselectivity.

Synthesis of α-l-Araf and β-d-Galf series furanobiosides using mutants of a GH51 α-l-arabinofuranosidase

The GH-51 α-l-arabinofuranosidase from Thermobacillus xylanilyticus (TxAbf) possesses versatile catalytic properties, displaying not only the ability to hydrolyze glycosidic linkages but also to synthesize furanobiosides in α-l-Araf and β-d-Galf series. Herein, mutants are investigated to evaluate their ability to perform self-condensation, assessing both yield improvements and changes in regioselectivity. Overall yields of oligo-α-l-arabino- and oligo-β-d-galactofuranosides were increased up to 4.8-fold compared to the wild-type enzyme. In depth characterization revealed that the mutants exhibit increased transfer rates and thus a hydrolysis/self-condensation ratio in favor of synthesis. The consequence of the substitution N216W is the creation of an additional binding subsite that provides the basis for an alternative acceptor substrate binding mode. As a result, mutants bearing N216W synthesize not only (1,2)-linked furanobiosides, but also (1,3)- and even (1,5)-linked furanobiosides. Since the self-condensation is under kinetic control, the yield of homo-disaccharides was maximized using higher substrate concentrations. In this way, the mutant R69H-N216W produced oligo-β-d-galactofuranosides in > 70% yield. Overall, this study further demonstrates the potential usefulness of TxAbf mutants for glycosynthesis and shows how these might be used to synthesize biologically-relevant glycoconjugates.