Toward the design of efficient transglycosidases: the case of the GH1 of Thermus thermophilus

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.

The Role of a Loop in the Non-catalytic Domain B on the Hydrolysis/Transglycosylation Specificity of the 4-α-Glucanotransferase from Thermotoga maritima

The mechanism by which glycoside hydrolases control the reaction specificity through hydrolysis or transglycosylation is a key element embedded in their chemical structures. The determinants of reaction specificity seem to be complex. We looked for structural differences in domain B between the 4-α-glucanotransferase from Thermotoga maritima (TmGTase) and the α-amylase from Thermotoga petrophila (TpAmylase) and found a longer loop in the former that extends towards the active site carrying a W residue at its tip. Based on these differences we constructed the variants W131G and the partial deletion of the loop at residues 120-124/128-131, which showed a 11.6 and 11.4-fold increased hydrolysis/transglycosylation (H/T) ratio relative to WT protein, respectively. These variants had a reduction in the maximum velocity of the transglycosylation reaction, while their affinity for maltose as the acceptor was not substantially affected. Molecular dynamics simulations allow us to rationalize the increase in H/T ratio in terms of the flexibility near the active site and the conformations of the catalytic acid residues and their associated pKas.

Modulating Glycoside Hydrolase Activity between Hydrolysis and Transfer Reactions Using an Evolutionary Approach

The proteins within the CAZy glycoside hydrolase family GH13 catalyze the hydrolysis of polysaccharides such as glycogen and starch. Many of these enzymes also perform transglycosylation in various degrees, ranging from secondary to predominant reactions. Identifying structural determinants associated with GH13 family reaction specificity is key to modifying and designing enzymes with increased specificity towards individual reactions for further applications in industrial, chemical, or biomedical fields. This work proposes a computational approach for decoding the determinant structural composition defining the reaction specificity. This method is based on the conservation of coevolving residues in spatial contacts associated with reaction specificity. To evaluate the algorithm, mutants of α-amylase (TmAmyA) and glucanotransferase (TmGTase) from Thermotoga maritima were constructed to modify the reaction specificity. The K98P/D99A/H222Q variant from TmAmyA doubled the transglycosydation/hydrolysis (T/H) ratio while the M279N variant from TmGTase increased the hydrolysis/transglycosidation ratio five-fold. Molecular dynamic simulations of the variants indicated changes in flexibility that can account for the modified T/H ratio. An essential contribution of the presented computational approach is its capacity to identify residues outside of the active center that affect the reaction specificity.

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.

Tuning the Transglycosylation Reaction of a GH11 Xylanase by a Delicate Enhancement of its Thumb Flexibility

Glycoside hydrolases (GHs) are attractive tools for multiple biotechnological applications. In conjunction with their hydrolytic function, GHs can perform transglycosylation under specific conditions. In nature, oligosaccharide synthesis is performed by glycosyltransferases (GTs); however, the industrial use of GTs is limited by their instability in solution. A key difference between GTs and GHs is the flexibility of their binding site architecture. We have used the xylanase from Bacillus circulans (BCX) to study the interplay between active-site flexibility and transglycosylation. Residues of the BCX "thumb" were substituted to increase the flexibility of the enzyme binding site. Replacement of the highly conserved residue P116 with glycine shifted the balance of the BCX enzymatic reaction toward transglycosylation. The effects of this point mutation on the structure and dynamics of BCX were investigated by NMR spectroscopy. The P116G mutation induces subtle changes in the configuration of the thumb and enhances the millisecond dynamics of the active site. Based on our findings, we propose the remodelling of the GH enzymes glycon site flexibility as a strategy to improve the transglycosylation efficiency of these biotechnologically important catalysts.