Strategies for Modulating the pH-Dependent Activity of a Family 11 Glycoside Hydrolase

Positive activation entropy of Bacillus circulans xylanase catalyzed ONPX2 hydrolysis: A mechanistic and engineering study

Transition state (TS) stabilization by enzymes greatly accelerates catalytic reactions. For some enzymes, the TS complex has entropy higher than enzyme substrate (ES) complex. But the origin of favorable entropy remains unclear. In this work, we studied the mechanism of Bacillus Circulans xylanase (BCX) 11 catalyzed o-nitrophenyl β-xylobioside (ONPX2) glycoside hydrolysis. The catalytic reaction exhibits a positive activation entropy, and an increase in ionic strength leads to a decrease in entropy without affecting the activation free energy, indicating that the entropy is predominantly influenced by electrostatic forces. Moreover, NMR measurements of electrostatic attractions within the active site demonstrate a positive entropy, aligning with molecular dynamics (MD) simulations showing that electrostatic interactions contribute to the entropic stabilization of the TS complex. These findings suggest that the positive entropy primarily originates from alterations in electrostatic interactions due to the formation of the oxocarbenium ion at C1 in the TS. Differences of electrostatic interactions in ES and TS modify hydrogen bonding of surrounding residues in the active site which causes their side chain dynamics and thus conformational entropy changes. Residues critical for the positive activation entropy are identified. A new BCX mutant with an increased activation entropy and catalytic activity is found.

Metal ions-anchored bacterial outer membrane vesicles for enhanced ferroptosis induction and immune stimulation in targeted antitumor therapy

The activation of ferroptosis presents a versatile strategy for enhancing the antitumor immune responses in cancer therapy. However, developing ferroptosis inducers that combine high biocompatibility and therapeutic efficiency remains challenging. In this study, we propose a novel approach using biological nanoparticles derived from outer membrane vesicles (OMVs) of Escherichia coli for tumor treatment, aiming to activate ferroptosis and stimulate the immune responses. Specifically, we functionalize the OMVs by anchoring them with ferrous ions via electrostatic interactions and loading them with the STING agonist-4, followed by tumor-targeting DSPE-PEG-FA decoration, henceforth referred to as OMV/SaFeFA. The anchoring of ferrous ions endows the OMVs with peroxidase-like activity, capable of inducing cellular lipid peroxidation by catalyzing H2O2 to •OH. Furthermore, OMV/SaFeFA exhibits pH-responsive release of ferrous ions and the agonist, along with tumor-targeting capabilities, enabling tumor-specific therapy while minimizing side effects. Notably, the concurrent activation of the STING pathway and ferroptosis elicits robust antitumor responses in colon tumor-bearing mouse models, leading to exceptional therapeutic efficacy and prolonged survival. Importantly, no acute toxicity was observed in mice receiving OMV/SaFeFA treatments, underscoring its potential for future tumor therapy and clinical translation.

Production of highly water-soluble genistein α-diglucoside using an engineered O-α-glycoligase with enhanced transglycosylation activity and altered substrate specificity

Genistein is one of isoflavones, showing various biological functions for human health. MalA-D416A, termed O-α-glycoligase, is an acid/base catalytic residue-deficient mutant of a α-glucosidase from Sulfolobus solfataricus, synthesizing genistein 7-O-α-glucoside using α-glucosyl fluoride as the donor substrate. Through mutagenesis toward MalA-D416A, an O-α-glycoligase variant with two mutations (D416R and Q450S) was identified as a biocatalyst with a 58.8-fold enhanced catalytic efficiency for genistein compared to the parent enzyme. The use of a 2:1 ratio of α-glucosyl fluoride and genistein at pH 9 facilitated the synthesis of genistein 7,4'-O-α-diglucoside by MalA-D416R/Q450S. The α-diglucoside exhibited 2,459-fold improved water solubility compared to genistein itself as well as facile deglycosylation by the intestinal α-glucosidase from rat, suggesting the potential of the α-diglucoside for improved bioavailability in human intestine. Through molecular docking analyses the modulation of the active site conformation by these mutations was expected for proper binding of both genistein and the monoglucoside.

Design of oscillatory dynamics in numerical simulations of compartment-based enzyme systems

Enzymatic reactions that yield non-neutral products are known to involve feedback due to the bell-shaped pH-rate curve of the enzyme. Compartmentalizing the reaction has been shown to lead to transport-driven oscillations in theory; however, there have been few reproducible experimental examples. Our objective was to determine how the conditions could be optimized to achieve pH oscillations. We employed numerical simulations to investigate the hydrolysis of ethyl acetate in a confined esterase enzyme system, examining the influence of key factors on its behavior. Specific parameter ranges that lead to bistability and self-sustained pH oscillations and the importance of fast base transport for oscillations in this acid-producing system are highlighted. Suggestions are made to expand the parameter space for the occurrence of oscillations, including modifying the maximum of the enzyme pH-rate curve and increasing the negative feedback rate. This research not only sheds light on the programmable nature of enzyme-driven pH regulation but also furthers knowledge on the optimal design of such feedback systems for experimentalists.

Strategies for tailoring pH performances of glycoside hydrolases

Glycoside hydrolases (GHs) exhibit high activity and stability under harsh conditions, such as high temperatures and extreme pHs, given their wide use in industrial biotechnology. However, strategies for improving the acidophilic and alkalophilic adaptations of GHs are poorly summarized due to the complexity of the mechanisms of these adaptations. This review not only highlights the adaptation mechanisms of acidophilic and alkalophilic GHs under extreme pH conditions, but also summarizes the recent advances in engineering the pH performances of GHs with a focus on four strategies of protein engineering, enzyme immobilization, chemical modification, and medium engineering (additives). The examples described here summarize the methods used in modulating the pH performances of GHs and indicate that methods integrated in different protein engineering techniques or methods are efficient to generate industrial biocatalysts with the desired pH performance and other adapted enzyme properties.

Redirecting substrate regioselectivity using engineered ΔN123-GBD-CD2 branching sucrases for the production of pentasaccharide repeating units of S. flexneri 3a, 4a and 4b haptens

The (chemo-)enzymatic synthesis of oligosaccharides has been hampered by the lack of appropriate enzymatic tools with requisite regio- and stereo-specificities. Engineering of carbohydrate-active enzymes, in particular targeting the enzyme active site, has notably led to catalysts with altered regioselectivity of the glycosylation reaction thereby enabling to extend the repertoire of enzymes for carbohydrate synthesis. Using a collection of 22 mutants of ΔN₁₂₃-GBD-CD2 branching sucrase, an enzyme from the Glycoside Hydrolase family 70, containing between one and three mutations in the active site, and a lightly protected chemically synthesized tetrasaccharide as an acceptor substrate, we showed that altered glycosylation product specificities could be achieved compared to the parental enzyme. Six mutants were selected for further characterization as they produce higher amounts of two favored pentasaccharides compared to the parental enzyme and/or new products. The produced pentasaccharides were shown to be of high interest as they are precursors of representative haptens of Shigella flexneri serotypes 3a, 4a and 4b. Furthermore, their synthesis was shown to be controlled by the mutations introduced in the active site, driving the glucosylation toward one extremity or the other of the tetrasaccharide acceptor. To identify the molecular determinants involved in the change of ΔN₁₂₃-GBD-CD2 regioselectivity, extensive molecular dynamics simulations were carried out in combination with in-depth analyses of amino acid residue networks. Our findings help to understand the inter-relationships between the enzyme structure, conformational flexibility and activity. They also provide new insight to further engineer this class of enzymes for the synthesis of carbohydrate components of bacterial haptens.

Studying the Role of a Single Mutation of a Family 11 Glycoside Hydrolase Using High-Resolution X-ray Crystallography

XynII is a family 11 glycoside hydrolase that uses the retaining mechanism for catalysis. In the active site, E177 works as the acid/base and E86 works as the nucleophile. Mutating an uncharged residue (N44) to an acidic residue (D) near E177 decreases the enzyme's optimal pH by ~ 1.0 unit. D44 was previously suggested to be a second proton carrier for catalysis. To test this hypothesis, we abolished the activity of E177 by mutating it to be Q, and mutated N44 to be D or E. These double mutants have dramatically decreased activities. Our high-resolution crystallographic structures and the microscopic pK_a calculations show that D44 has similar position and pK_a value during catalysis, indicating that D44 changes electrostatics around E177, which makes it prone to rotate as the acid/base in acidic conditions, thus decreases the pH optimum. Our results could be helpful to design enzymes with different pH optimum.

Altering the Y137-K164-K166 triad of mandelate racemase and its effect on the observed pKa of the Brønsted base catalysts

Mandelate racemase (MR) catalyzes the interconversion of the enantiomers of mandelate using a two-base mechanism with Lys 166 acting as the Brønsted base to abstract the α-proton from (S)-mandelate. The resulting intermediate is subsequently re-protonated by the conjugate acid of His 297 to yield (R)-mandelate. The roles of these amino acids are reversed when (R)-mandelate is the substrate. The side chains of Tyr 137, Lys 164, and Lys 166 form a H-bonding network and the proximity of the two ε-NH₃⁺ groups is believed to lower the pK_a of Lys 166. We used site-directed mutagenesis, kinetics, and pH-rate studies to explore the roles of Lys 164 (K164 C/M) and Tyr 137 (Y137  L/F/S/T) in catalysis. The efficiency (k_cat/K_m) was reduced ∼3.5 × 10⁵-fold for K164C MR, relative to wild-type MR, indicating a major role for this residue in catalysis. The efficiency of Y137F MR, however, was reduced only 25-30-fold. pH-Rate profiles (log k_cat vs. pH) revealed that substitution of Tyr 137 by Phe increased the kinetic pK_a of Lys 166 from 5.88 ± 0.02 to 7.3 ± 0.2. Hence, Tyr 137 plays an important role in facilitating the reduction of the pK_a of the Brønsted base Lys 166 by ∼1.4 units. Interestingly, the Phe substitution also increased the kinetic pK_a of His 297 from 5.97 ± 0.04 to 7.1 ± 0.1. Thus, the Tyr 137-Lys 164-Lys 166 H-bonding network plays a broader role in modulating the pK_a of catalytic residues by influencing the electrostatic character of the entire active site, not only by decreasing the observed pK_a value of Lys 166, but also by decreasing the pK_a of His 297 by 1.1 units.

Characterization of the two conformations adopted by the T3SS inner-membrane protein PrgK

The pathogenic bacterium Salmonella enterica serovar Typhimurium utilizes two type III secretion systems (T3SS) to inject effector proteins into target cells upon infection. The T3SS secretion apparatus (the injectisome) is a large macromolecular assembly composed of over twenty proteins, many in highly oligomeric states. A sub-structure of the injectisome, termed the basal body, spans both membranes and the periplasmic space of the bacterium. It is primarily composed of three integral membranes proteins, InvG, PrgH, and PrgK, that form ring structures through which components are secreted. In particular, PrgK possesses a periplasmic region consisting of two globular domains joined by a linker polypeptide. We showed previously that in isolation, this region adopts two distinct conformations, of with only one is observed in the assembled basal body complex. Here, using NMR spectroscopy, we further characterize these two conformations. In particular, we demonstrate that the interaction of the linker region with the first globular domain, as found in the intact basal body, is dependent upon the cis conformation of the Leu77-Pro78 peptide. Furthermore, this interaction is pH-dependent due to coupling with hydrogen bond formation between Tyr75 and His42 in its neutral N^δ1 H tautomeric form. This pH-dependent interaction may play a role in the regulation of the secretion apparatus disassembly in the context of bacterial infection.

Rationalising pKa shifts in Bacillus circulans xylanase with computational studies

Bacillus circulans xylanase (BcX), a family 11 glycoside hydrolase, catalyses the hydrolysis of xylose polymers with a net retention of stereochemistry. Glu172 in BcX is believed to act as a general acid by protonating the aglycone during glycosylation, and then as a general base to facilitate the deglycosylation step. The key to the dual role of this general acid/base lies in its protonation states, which depend on its intrinsic pK_a value and the specific environment which it resides within. To fully understand the detailed molecular features in BcX to establish the dual role of Glu172, we present a combined study based on both atomistic simulations and empirical models to calculate pK_a shifts for the general acid/base Glu172 in BcX at different functional states. Its pK_a values and those of nearby residues, obtained based on QM/MM free energy calculations, MCCE and PROPKA, show a good agreement with available experimental data. Additionally, our study provides additional insights into the effects of structural and electrostatic perturbations caused by mutations and chemical modifications, suggesting that the local solvation environment and mutagenesis of the residues adjacent to Glu172 establish its dual role during hydrolysis. The strengths and limitations of various methods for calculating pK_as and pK_a shifts have also been discussed.

Strategy in manipulating transglycosylation activity of glycosyl hydrolase for oligosaccharide production

BACKGROUND

The increasing market demand for oligosaccharides has intensified the need for efficient biocatalysts. Glycosyl hydrolases (GHs) are still gaining popularity as biocatalyst for oligosaccharides synthesis owing to its simple reaction and high selectivity.

PURPOSE

Over the years, research has advanced mainly directing to one goal; to reduce hydrolysis activity of GHs for increased transglycosylation activity in achieving high production of oligosaccharides.

DESIGN AND METHODS

This review concisely presents the strategies to increase transglycosylation activity of GHs for oligosaccharides synthesis, focusing on controlling the reaction equilibrium, and protein engineering. Various modifications of the subsites of GHs have been demonstrated to significantly modulate the hydrolysis and transglycosylation activity of the enzymes. The clear insight of the roles of each amino acid in these sites provides a platform for designing an enzyme that could synthesize a specific oligosaccharide product.

CONCLUSIONS

The key strategies presented here are important for future improvement of GHs as a biocatalyst for oligosaccharide synthesis.

Ligand Binding Enhances Millisecond Conformational Exchange in Xylanase B2 from Streptomyces lividans

Xylanases catalyze the hydrolysis of xylan, an abundant carbon and energy source with important commercial ramifications. Despite tremendous efforts devoted to the catalytic improvement of xylanases, success remains limited because of our relatively poor understanding of their molecular properties. Previous reports suggested the potential role of atomic-scale residue dynamics in modulating the catalytic activity of GH11 xylanases; however, dynamics in these studies was probed on time scales orders of magnitude faster than the catalytic time frame. Here, we used nuclear magnetic resonance titration and relaxation dispersion experiments ((15)N-CPMG) in combination with X-ray crystallography and computational simulations to probe conformational motions occurring on the catalytically relevant millisecond time frame in xylanase B2 (XlnB2) and its catalytically impaired mutant E87A from Streptomyces lividans 66. Our results show distinct dynamical properties for the apo and ligand-bound states of the enzymes. The apo form of XlnB2 experiences conformational exchange for residues in the fingers and palm regions of the catalytic cleft, while the catalytically impaired E87A variant displays millisecond dynamics only in the fingers, demonstrating the long-range effect of the mutation on flexibility. Ligand binding induces enhanced conformational exchange of residues interacting with the ligand in the fingers and thumb loop regions, emphasizing the potential role of residue motions in the fingers and thumb loop regions for recognition, positioning, processivity, and/or stabilization of ligands in XlnB2. To the best of our knowledge, this work represents the first experimental characterization of millisecond dynamics in a GH11 xylanase family member. These results offer new insights into the potential role of conformational exchange in GH11 enzymes, providing essential dynamic information to help improve protein engineering and design applications.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

The complete conformational free energy landscape of β-xylose reveals a two-fold catalytic itinerary for β-xylanases

Unraveling the conformational catalytic itinerary of glycoside hydrolases (GHs) is a growing topic of interest in glycobiology, with major impact in the design of GH inhibitors. β-xylanases are responsible for the hydrolysis of glycosidic bonds in β-xylans, a group of hemicelluloses of high biotechnological interest that are found in plant cell walls. The precise conformations followed by the substrate during catalysis in β-xylanases have not been unambiguously resolved, with three different pathways being proposed from structural analyses. In this work, we compute the conformational free energy landscape (FEL) of β-xylose to predict the most likely catalytic itineraries followed by β-xylanases. The calculations are performed by means of ab initio metadynamics, using the Cremer-Pople puckering coordinates as collective variables. The computed FEL supports only two of the previously proposed itineraries, 2SO → [2,5B]ǂ → 5S1 and 1S3 → [4H3]ǂ → 4C1, which clearly appear in low energy regions of the FEL. Consistently, 2SO and 1S3 are conformations preactivated for catalysis in terms of free energy/anomeric charge and bond distances. The results however exclude the OE → [OS2]ǂ → B2,5 itinerary that has been recently proposed for a family 11 xylanase. Classical and ab initio QM/MM molecular dynamics simulations reveal that, in this case, the observed OE conformation has been enforced by enzyme mutation. These results add a word of caution on using modified enzymes to inform on catalytic conformational itineraries of glycoside hydrolases.

pH-dependent random coil (1)H, (13)C, and (15)N chemical shifts of the ionizable amino acids: a guide for protein pK a measurements

The pK a values and charge states of ionizable residues in polypeptides and proteins are frequently determined via NMR-monitored pH titrations. To aid the interpretation of the resulting titration data, we have measured the pH-dependent chemical shifts of nearly all the (1)H, (13)C, and (15)N nuclei in the seven common ionizable amino acids (X = Asp, Glu, His, Cys, Tyr, Lys, and Arg) within the context of a blocked tripeptide, acetyl-Gly-X-Gly-amide. Alanine amide and N-acetyl alanine were used as models of the N- and C-termini, respectively. Together, this study provides an essentially complete set of pH-dependent intra-residue and nearest-neighbor reference chemical shifts to help guide protein pK a measurements. These data should also facilitate pH-dependent corrections in algorithms used to predict the chemical shifts of random coil polypeptides. In parallel, deuterium isotope shifts for the side chain (15)N nuclei of His, Lys, and Arg in their positively-charged and neutral states were also measured. Along with previously published results for Asp, Glu, Cys, and Tyr, these deuterium isotope shifts can provide complementary experimental evidence for defining the ionization states of protein residues.

X-ray crystallographic studies of family 11 xylanase Michaelis and product complexes: implications for the catalytic mechanism

Xylanases catalyze the hydrolysis of plant hemicellulose xylan into oligosaccharides by cleaving the main-chain glycosidic linkages connecting xylose subunits. To study ligand binding and to understand how the pH constrains the activity of the enzyme, variants of the Trichoderma reesei xylanase were designed to either abolish its activity (E177Q) or to change its pH optimum (N44H). An E177Q-xylohexaose complex structure was obtained at 1.15 Å resolution which represents a pseudo-Michaelis complex and confirmed the conformational movement of the thumb region owing to ligand binding. Co-crystallization of N44H with xylohexaose resulted in a hydrolyzed xylotriose bound in the active site. Co-crystallization of the wild-type enzyme with xylopentaose trapped an aglycone xylotriose and a transglycosylated glycone product. Replacing amino acids near Glu177 decreased the xylanase activity but increased the relative activity at alkaline pH. The substrate distortion in the E177Q-xylohexaose structure expands the possible conformational itinerary of this xylose ring during the enzyme-catalyzed xylan-hydrolysis reaction.

A computational method for the systematic screening of reaction barriers in enzymes: searching for Bacillus circulans xylanase mutants with greater activity towards a synthetic substrate

We present a semi-empirical (PM6-based) computational method for systematically estimating the effect of all possible single mutants, within a certain radius of the active site, on the barrier height of an enzymatic reaction. The intent of this method is not a quantitative prediction of the barrier heights, but rather to identify promising mutants for further computational or experimental study. The method is applied to identify promising single and double mutants of Bacillus circulans xylanase (BCX) with increased hydrolytic activity for the artificial substrate ortho-nitrophenyl β-xylobioside (ONPX2). The estimated reaction barrier for wild-type (WT) BCX is 18.5 kcal/mol, which is in good agreement with the experimental activation free energy value of 17.0 kcal/mol extracted from the observed k cat using transition state theory (Joshi et al., 2001). The PM6 reaction profiles for eight single point mutations are recomputed using FMO-MP2/PCM/6-31G(d) single points. PM6 predicts an increase in barrier height for all eight mutants while FMO predicts an increase for six of the eight mutants. Both methods predict that the largest change in barrier occurs for N35F, where PM6 and FMO predict a 9.0 and 15.8 kcal/mol increase, respectively. We thus conclude that PM6 is sufficiently accurate to identify promising mutants for further study. We prepared a set of all theoretically possible (342) single mutants in which every amino acid of the active site (except for the catalytically active residues E78 and E172) was mutated to every other amino acid. Based on results from the single mutants we construct a set of 111 double mutants consisting of all possible pairs of single mutants with the lowest barrier for a particular position and compute their reaction profile. None of the mutants have, to our knowledge, been prepared experimentally and therefore present experimentally testable predictions.