QM/MM Studies on the β-Galactosidase Catalytic Mechanism: Hydrolysis and Transglycosylation Reactions

QM/MM study of the catalytic reaction of aphid myrosinase

Brevicoryne brassicae, an aphid species, exclusively consumes plants from the Brassicaceae family and employs a sophisticated defense mechanism involving a myrosinase enzyme that breaks down glucosinolates obtained from its host plants. In this work, we employed combined quantum mechanical and molecular mechanical (QM/MM) calculations and molecular dynamics (MD) simulations to study the catalytic reaction of aphid myrosinase. A proper QM region to study the myrosinase reaction should contain the whole substrate, models of Gln-19, His-122, Asp-124, Asn-166, Glu-167, Lys-173, Tyr-180, Val-228, Tyr-309, Tyr-346, Ile-347, Glu-374, Glu-423, Trp-424, and a water molecule. The calculations show that Asp-124 and Glu-423 must be charged, His-122 must be protonated on NE2, and Glu-167 must be protonated on OE2. Our model reproduces the anomeric retaining characteristic of myrosinase and indicates that the deglycosylation reaction is the rate-determining step of the reaction. Based on the calculations, we propose a reaction mechanism for aphid myrosinase-mediated hydrolysis of glucosinolates with an overall barrier of 15.2 kcal/mol. According to the results, removing a proton from Arg-312 or altering it to valine or methionine increases glycosylation barriers but decreases the deglycosylation barrier.

Invited review: Properties of β-galactosidases derived from Lactobacillaceae species and their capacity for galacto-oligosaccharide production

β-galactosidase (enzymatic class 3.2.1.23) is one of the dairy industry's most important and widely used enzymes. The enzyme is part of a large family known to catalyze hydrolysis and transglycosylation reactions. Its hydrolytic activity is commonly used to decrease lactose content in dairy products, while its transglycosylase activity has recently been used to synthesize galacto-oligosaccharides (GOS). During the past couple of years, researchers have focused on studying β-galactosidase isolated and purified from lactic acid bacteria. This review will focus on β-galactosidase purified and characterized from what used to be the Lactobacillus genera. Furthermore, particular emphasis is given to its kinetics, biochemical characteristics, GOS production, market, and utilization by Lactobacilllaceae species.

Mangosteen Metabolites as Promising Alpha-Amylase Inhibitor Candidates: In Silico and In Vitro Evaluations

Diabetes is a chronic metabolic disorder characterized by raised glucose levels in the blood, resulting in grave damage over time to various body organs, including the nerves, heart, kidneys, eyes, and blood vessels. One of its therapeutic treatment approaches involves the inhibition of enzymes accountable for carbohydrate digestion and absorption. The present work is aimed at evaluating the potential of some reported metabolites from Garcinia mangostana (mangosteen, Guttiferae) as alpha-amylase inhibitors. Forty compounds were assessed for their capacity to inhibit alpha-amylase using in silico studies as well as in vitro assays. Molecular docking was carried out to analyze their binding capacities in the 3D structure of alpha-amylase (PDB ID: 4GQR). Among the tested compounds, 6-O-β-D-glucopyranosyl-2,4,6,3',4',6'-hexahydroxybenzophenone (8), aromadendrin-8-C-glucoside (5), epicatechin (6), rhodanthenone (4), and garcixanthone D (40) had a high XP G.score and a Glide G.score of -12.425, -11.855, -11.135, and -11.048 Kcal/mol, respectively. Compound 8 possessed the XP and Glide docking score of -12.425 Kcal/mol compared to the reference compounds myricetin and acarbose which had an XP and Glide docking score of -12.319 and 11.201 Kcal/mol, respectively. It interacted through hydrogen bond formations between its hydroxyl groups and the residues His 101, Asp 197, Glu 233, Asp 300, and His 305, in addition to water bridges and hydrophobic interactions. Molecular mechanics-generalized born surface area (MM-GBSA) was used to calculate the binding free energy and molecular dynamic studies that indicated the stability of the alpha-amylase-compound 8 complex during the 100 ns simulation in comparison with myricetin- and acarbose-alpha-amylase complexes. Additionally, the in vitro alpha-amylase inhibition assay findings validated the in silico study's findings. This could further validate the potential of G. mangostana as a candidate for diabetes management.

Role of Enzyme and Active Site Conformational Dynamics in the Catalysis by α-Amylase Explored with QM/MM Molecular Dynamics

We assessed enzyme:substrate conformational dynamics and the rate-limiting glycosylation step of a human pancreatic α-amylase:maltopentose complex. Microsecond molecular dynamics simulations suggested that the distance of the catalytic Asp197 nucleophile to the anomeric carbon of the buried glucoside is responsible for most of the enzyme active site fluctuations and that both Asp197 and Asp300 interact the most with the buried glucoside unit. The buried glucoside binds either in a ⁴C₁ chair or ²S_O skew conformations, both of which can change to TS-like conformations characteristic of retaining glucosidases. Starting from four distinct enzyme:substrate complexes, umbrella sampling quantum mechanics/molecular mechanics simulations (converged within less than 1 kcal·mol^-1 within a total simulation time of 1.6 ns) indicated that the reaction occurrs with a Gibbs barrier of 13.9 kcal·mol ^-1, in one asynchronous concerted step encompassing an acid-base reaction with Glu233 followed by a loose S_N2-like nucleophilic substitution by the Asp197. The transition state is characterized by a ²H₃ half-chair conformation of the buried glucoside that quickly changes to the E₃ envelope conformation preceding the attack of the anomeric carbon by the Asp197 nucleophile. Thermodynamic analysis of the reaction supported that a water molecule tightly hydrogen bonded to the glycosidic oxygen of the substrate at the reactant state (∼1.6 Å) forms a short hydrogen bond with Glu233 at the transition state (∼1.7 Å) and lowers the Gibbs barrier in over 5 kcal·mol^-1. The resulting Asp197-glycosyl was mostly found in the ⁴C₁ conformation, although the more endergonic B_3,O conformation was also observed. Altogether, the combination of short distances for the acid-base reaction with the Glu233 and for the nucleophilic attack by the Asp197 nucleophile and the availability of water within hydrogen bonding distance of the glycosidic oxygen provides a reliable criteria to identify reactive conformations of α-amylase complexes.

A novel salt-tolerant GH42 β-galactosidase with transglycosylation activity from deep-sea metagenome

β-Galactosidase is a widely adopted enzyme in the food and pharmaceutical industries. Metagenome techniques have the advantage of discovering novel functional genes, particularly potential genes from uncultivated microbes. In this study, a novel GH42 β-galactosidase isolated from a deep-sea metagenome was overexpressed in Escherichia coli BL21 (DE3) and purified by affinity chromatography. The optimal temperatures and pH of the enzyme for o-nitrophenyl-β-D-galactopyranoside (oNPG) and lactose were 40 ℃, 6.5 and 50 ℃, 7, respectively. The enzyme was stable at temperatures between 4 and 30 ℃ and within the pH range of 6-9. Moreover, it was highly tolerant to salt and inhibited by Zn²⁺ and Cu²⁺. The kinetic values of K_m and k_cat of the enzyme against oNPG were 1.1 mM and 57.8 s^-1, respectively. Furthermore, it showed hydrolysis and transglycosylation activity to lactose and the extra monosaccharides could improve the productivity of oligosaccharides. Overall, this recombinant β-galactosidase is a potential biocatalyst for the hydrolysis of milk lactose and synthesis of functional oligosaccharides.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Structure based virtual screening of natural product molecules as glycosidase inhibitors

Objective of the present investigation comprised of the application of in silico methods to discover novel natural product (NP) based potential inhibitors for carbohydrate mediated diseases. Structure based drug design studies (molecular docking and structure based pharmacophore analysis)  were carried out on a series of natural product compounds to identify significant bioactive molecules to inhibit α-mannosidase (I and II) and β-galactosidase enzymes. Furthermore, protein ligand interaction fingerprint analysis, molecular dynamics simulations and molecular access system (MACCS) fingerprint analysis were performed to understand the binding behaviors of the studied molecules. The results derived from these analyses showed that the identified compounds exhibit significant binding interactions with the active site residues. The compounds, NP-51, NP-81 and NP-165 have shown significant docking score against the studied enzymes (α-mannosidases-I, α-mannosidases-II and β-galactosidases). The fingerprint studies showed that the presence of rings (aromatic or aliphatic) with sulfur atoms, nitrogen atoms, methyl groups, etc. have favorable effects on the α-mannosidase II inhibitory activity. However, the presence of halogen atoms substituted in the molecules have reduced inhibitory ability against α-mannosidase II. The compound, NP-165 has significant activity against both enzymes (α-mannosidases and β-galactosidases). These studies accomplished that the compounds identified through in silico methodologies can be used to develop semisynthetic derivatives of the glycosidase inhibitors and can be screened for the treatment of different carbohydrate mediated diseases.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-021-00115-9.

Galacto-oligosaccharides as infant prebiotics: production, application, bioactive activities and future perspectives

Galacto-oligosaccharides (GOS) are non-digestible oligosaccharides characterized by a mix of structures that vary in their degree of polymerization (DP) and glycosidic linkage between the galactose moieties or between galactose and glucose. They have enjoyed extensive scientific scrutiny, and their health-promoting effects are supported by a large number of scientific and clinical studies. A variety of GOS-associated health-promoting effects have been reported, such as growth promotion of beneficial bacteria, in particular bifidobacteria and lactobacilli, inhibition of pathogen adhesion and improvement of gut barrier function. GOS have attracted significant interest from food industries for their versatility as a bioactive ingredient and in particular as a functional component of infant formulations. These oligosaccharides are produced in a kinetically-controlled reaction involving lactose transgalactosylation, being catalyzed by particular β-galactosidases of bacterial or fungal origin. Despite the well-established technology applied for GOS production, this process may still meet with technological challenges when employed at an industrial scale. The current review will cover relevant scientific literature on the beneficial physiological properties of GOS as a prebiotic for the infant gut microbiota, details of GOS structures, the associated reaction mechanism of β-galactosidase, and its (large-scale) production.

Infant-Associated Bifidobacterial β-Galactosidases and Their Ability to Synthesize Galacto-Oligosaccharides

Galacto-oligosaccharides (GOS) represent non-digestible glycans that are commercially produced by transgalactosylation of lactose, and that are widely used as functional food ingredients in prebiotic formulations, in particular in infant nutrition. GOS consumption has been reported to enhance growth of specific bacteria in the gut, in particular bifidobacteria, thereby supporting a balanced gut microbiota. In a previous study, we assessed the hydrolytic activity and substrate specificity of seventeen predicted β-galactosidases encoded by various species and strains of infant-associated bifidobacteria. In the current study, we further characterized seven out of these seventeen bifidobacterial β-galactosidases in terms of their kinetics, enzyme stability and oligomeric state. Accordingly, we established whether these β-galactosidases are capable of synthesizing GOS via enzymatic transgalactosylation employing lactose as the feed substrate. Our findings show that the seven selected enzymes all possess such transgalactosylation activity, though they appear to differ in their efficiency by which they perform this reaction. From chromatography analysis, it seems that these enzymes generate two distinct GOS mixtures: GOS with a relatively short or long degree of polymerization profile. These findings may be the stepping stone for further studies aimed at synthesizing new GOS variants with novel and/or enhanced prebiotic activities and potential for industrial applications.

Transglycosylation abilities of β-d-galactosidases from GH family 2

The ability to predict the transglycosylation activity of glycosidases by in silico analysis was investigated. The transglycosylation abilities of 7 different β-d-galactosidases from GH family 2 were tested experimentally using 7 different acceptors and p-nitrophenyl-β-d-galactopyranoside as a donor of galactosyl moiety. Similar transglycosylation abilities were confirmed for all enzymes originating from bacteria belonging to Enterobacteriaceae, which were able to use all tested acceptor molecules. Higher acceptor selectivity was observed for all others used bacterial strains. Structure models of all enzymes were constructed using homology modeling. Ligand-docking method was used for enzymes-transglycosylation products models construction and evaluation. Results obtained by in silico analysis were compared with results arisen out of experimental testing. The experiments confirmed that significant differences in transglycosylation abilities are caused by small differences in active sites composition of analyzed enzymes. According to obtained result, it is possible to conclude that homology modeling may serve as a quick starting point for detection or exclusion of enzymes with defined transglycosylation abilities, which can be used for subsequent synthesis of e.g., pharmaceutically interesting glycosides.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02715-w.

QM/MM Study of the Catalytic Reaction of Myrosinase; Importance of Assigning Proper Protonation States of Active-Site Residues

Myrosinase from Sinapis alba hydrolyzes glycosidic bonds of β-d-S-glucosides. The enzyme shows an enhanced activity in the presence of l-ascorbic acid. In this work, we employed combined quantum mechanical and molecular mechanical (QM/MM) calculations and molecular dynamics simulations to study the catalytic reaction of wild-type myrosinase and its E464A, Q187A, and Q187E mutants. Test calculations show that a proper QM region to study the myrosinase reaction must contain the whole substrate, models of Gln-187, Glu-409, Gln-39, His-141, Asn-186, Tyr-330, Glu-464, Arg-259, and a water molecule. Furthermore, to make the deglycosylation step possible, Arg-259 must be charged, Glu-464 must be protonated on OE2, and His-141 must be protonated on the NE2 atom. The results indicate that assigning proper protonation states of the residues is more important than the size of the model QM system. Our model reproduces the anomeric retaining characteristic of myrosinase and also reproduces the experimental fact that ascorbate increases the rate of the reaction. A water molecule in the active site, positioned by Gln-187, helps the aglycon moiety of the substrate to stabilize the buildup of negative charge during the glycosylation reaction and this in turn makes the moiety a better leaving group. The water molecule also lowers the glycosylation barrier by ∼9 kcal/mol. The results indicate that the Q187E and E464A mutants but not the Q187A mutant can perform the glycosylation step. However, the energy profiles for the deglycosylation step of the mutants are not similar to that of the wild-type enzyme. The Glu-464 residue lowers the barriers of the glycosylation and deglycosylation steps. The ascorbate ion can act as a general base in the reaction of the wild-type enzyme only if the Glu-464 and His-141 residues are properly protonated.

Cysteine Nucleophiles in Glycosidase Catalysis: Application of a Covalent β-l-Arabinofuranosidase Inhibitor

The recent discovery of zinc-dependent retaining glycoside hydrolases (GHs), with active sites built around a Zn(Cys)₃ (Glu) coordination complex, has presented unresolved mechanistic questions. In particular, the proposed mechanism, depending on a Zn-coordinated cysteine nucleophile and passing through a thioglycosyl enzyme intermediate, remains controversial. This is primarily due to the expected stability of the intermediate C-S bond. To facilitate the study of this atypical mechanism, we report the synthesis of a cyclophellitol-derived β-l-arabinofuranosidase inhibitor, hypothesised to react with the catalytic nucleophile to form a non-hydrolysable adduct analogous to the mechanistic covalent intermediate. This β-l-arabinofuranosidase inhibitor reacts exclusively with the proposed cysteine thiol catalytic nucleophiles of representatives of GH families 127 and 146. X-ray crystal structures determined for the resulting adducts enable MD and QM/MM simulations, which provide insight into the mechanism of thioglycosyl enzyme intermediate breakdown. Leveraging the unique chemistry of cyclophellitol derivatives, the structures and simulations presented here support the assignment of a zinc-coordinated cysteine as the catalytic nucleophile and illuminate the finely tuned energetics of this remarkable metalloenzyme clan.

Technical integrative approaches to cheese whey valorization towards sustainable environment

Whey, a byproduct of cheese production, is often treated as an industrial dairy waste. A large volume of this product is disposed of annually due to inadequate bioconversion approaches. With its high pollutant load, disposal without pretreatment has raised a lot of environmental concerns alerting the need to seek optimal methods for adequately extracting and utilizing its organic content. In recent years, several techniques for whey valorization have emerged which may serve as interventionary measures against its environmental effects after disposal. In this review, we discuss five major approaches, by which whey can be converted into eco-friendly products, to significantly cut whey wastage. The approaches to whey valorization are therefore examined under the following perspectives: whey as a raw material for the production of bioethanol and prebiotic oligosaccharides via β-galactosidase and microbe catalyzed reactions, for the production of refined lactose as an excipient for pharmaceutical purposes, and the clinical significance of whey hydrolysates and their antifungal activity in food processing.

Mapping the Transglycosylation Relevant Sites of Cold-Adapted β-d-Galactosidase from Arthrobacter sp. 32cB

β-Galactosidase from Arthrobacter sp. 32cB (ArthβDG) is a cold-adapted enzyme able to catalyze hydrolysis of β-d-galactosides and transglycosylation reaction, where galactosyl moiety is being transferred onto an acceptor larger than a water molecule. Mutants of ArthβDG: D207A and E517Q were designed to determine the significance of specific residues and to enable formation of complexes with lactulose and sucrose and to shed light onto the structural basis of the transglycosylation reaction. The catalytic assays proved loss of function mutation E517 into glutamine and a significant drop of activity for mutation of D207 into alanine. Solving crystal structures of two new mutants, and new complex structures of previously presented mutant E441Q enables description of introduced changes within active site of enzyme and determining the importance of mutated residues for active site size and character. Furthermore, usage of mutants with diminished and abolished enzymatic activity enabled solving six complex structures with galactose, lactulose or sucrose bounds. As a result, not only the galactose binding sites were mapped on the enzyme's surface but also the mode of lactulose, product of transglycosylation reaction, and binding within the enzyme's active site were determined and the glucopyranose binding site in the distal of active site was discovered. The latter two especially show structural details of transglycosylation, providing valuable information that may be used for engineering of ArthβDG or other analogous galactosidases belonging to GH2 family.

Recent advances in β-galactosidase and fructosyltransferase immobilization technology

The highly demanding conditions of industrial processes may lower the stability and affect the activity of enzymes used as biocatalysts. Enzyme immobilization emerged as an approach to promote stabilization and easy removal of enzymes for their reusability. The aim of this review is to go through the principal immobilization strategies addressed to achieve optimal industrial processes with special care on those reported for two types of enzymes: β-galactosidases and fructosyltransferases. The main methods used to immobilize these two enzymes are adsorption, entrapment, covalent coupling and cross-linking or aggregation (no support is used), all of them having pros and cons. Regarding the support, it should be cost-effective, assure the reusability and an easy recovery of the enzyme, increasing its stability and durability. The discussion provided showed that the type of enzyme, its origin, its purity, together with the type of immobilization method and the support will affect the performance during the enzymatic synthesis. Enzymes' immobilization involves interdisciplinary knowledge including enzymology, nanotechnology, molecular dynamics, cellular physiology and process design. The increasing availability of facilities has opened a variety of possibilities to define strategies to optimize the activity and re-usability of β-galactosidases and fructosyltransferases, but there is still great place for innovative developments.

Active Site Architecture and Reaction Mechanism Determination of Cold Adapted β-d-galactosidase from Arthrobacter sp. 32cB

ArthβDG is a dimeric, cold-adapted β-d-galactosidase that exhibits high hydrolytic and transglycosylation activity. A series of crystal structures of its wild form, as well as its ArthβDG_E441Q mutein complexes with ligands were obtained in order to describe the mode of its action. The ArthβDG_E441Q mutein is an inactive form of the enzyme designed to enable observation of enzyme interaction with its substrate. The resulting three-dimensional structures of complexes: ArthβDG_E441Q/LACs and ArthβDG/IPTG (ligand bound in shallow mode) and structures of complexes ArthβDG_E441Q/LACd, ArthβDG/ONPG (ligands bound in deep mode), and galactose ArthβDG/GAL and their analysis enabled structural characterization of the hydrolysis reaction mechanism. Furthermore, comparative analysis with mesophilic analogs revealed the most striking differences in catalysis mechanisms. The key role in substrate transfer from shallow to deep binding mode involves rotation of the F581 side chain. It is worth noting that the 10-aa loop restricting access to the active site in mesophilic GH2 βDGs, in ArthβDG is moved outward. This facilitates access of substrate to active site. Such a permanent exposure of the entrance to the active site may be a key factor for improved turnover rate of the cold adapted enzyme and thus a structural feature related to its cold adaptation.

Multifunctional properties of Lactobacillus plantarum strains WiKim83 and WiKim87 as a starter culture for fermented food

This study aimed to evaluate the safety (hemolysis and enzyme activity), probiotic properties (gastrointestinal tract tolerance, adhesion, hydrophobicity, and auto-aggregation), and functional characteristics (antimicrobial, antioxidant, and β-galactosidase activities) of lactic acid bacteria (LAB), isolated from kimchi, in order to select a multifunctional LAB strain for starter culture in fermented food. The five isolated strains included Lactobacillus plantarum WiKim83, L. plantarum WiKim84, Pediococcus pentosaceus WiKim85, P. pentosaceus WiKim86, and L. plantarum WiKim87, as identified by 16S rRNA gene sequence analysis; they were confirmed to be nonhemolytic and not able to produce β-glucuronidase, a carcinogenic enzyme. Probiotic properties of the five LAB strains were evaluated relative to those of commercial Lactobacillus rhamnosus GG, and results revealed probiotic potential of three strains (L. plantarum WiKim83, L. plantarum WiKim84, and L. plantarum WiKim87) to be superior. L. plantarum WiKim84 showed high antimicrobial activity against pathogens, and L. plantarum WiKim83 exhibited the highest antioxidant and β-galactosidase activities. Based on the probiotic and functional properties, the main characteristics of each strain were highlighted and two of them, L. plantarum WiKim83 and L. plantarum WiKim87, were selected as the most potent by principal component analysis. These strains showed antimicrobial, β-galactosidase, and antioxidant activities, which recommend their suitability as starter culture in various fermented foods.

Intrinsic dynamic behavior of enzyme:substrate complexes govern the catalytic action of β-galactosidases across clan GH-A

The conformational itineraries taken by carbohydrate residues in the catalytic subsite of retaining glycoside hydrolases (GHs), harness the link between substrate conformation and reactivity. GHs' active sites may be described as a combination of subsites dedicated to the binding of individual sugar residues and to catalysis. The three-dimensional structure of GH:carbohydrate complexes has demonstrated that carbohydrate ring conformation changes in an ordered manner during catalysis. Here we demonstrate in silico that a link exists between subsite binding dynamics and substrate specificity for β-galactosidases from clan GH-A families GH1, GH2, GH35, GH42 and GH59. Different oligosaccharides were docked in the active site of reference β-galactosidase structures using Vina-Carb. Subsequent molecular dynamics (MD) simulations revealed that these enzymes favor a high degree of flexibility and ring distortion of the substrate the lytic subsite -1. Although the β-galactosidase families examined are structurally and mechanistically related, distinct patterns of ring distortion were unveiled for the different families. For β-galactosidases, three different family-dependent reaction itineraries (1S3 → 4H3‡ → 4C1, 1,4B → 4H3/ 4E‡ → 4C1, and 1S5 → 4E/ 4H5‡ → 4C1) were identified, all compatible with the antiperiplanar lone pair hypothesis (ALPH) for the hydrolysis of β-glycosides. This comparative study reveals the fuzzy character of the changes in carbohydrate ring geometry prior to carbohydrate hydrolysis.

Comparing Hydrolysis and Transglycosylation Reactions Catalyzed by Thermus thermophilus β-Glycosidase. A Combined MD and QM/MM Study

The synthesis of oligosaccharides and other carbohydrate derivatives is of relevance for the advancement of glycosciences both at the fundamental and applied level. For many years, glycosyl hydrolases (GHs) have been explored to catalyze the synthesis of glycosidic bonds. In particular, retaining GHs can catalyze a transglycosylation (T) reaction that competes with hydrolysis (H). This has been done either employing controlled conditions in wild type GHs or by engineering new mutants. The goal, which is to increase the T/H ratio, has been achieved with moderate success in several cases despite the fact that the molecular basis for T/H modulation are unclear. Here we have used QM(DFT)/MM calculations to compare the glycosylation, hydrolysis and transglycosylation steps catalyzed by wild type Thermus thermophilus β-glycosidase (family GH1), a retaining glycosyl hydrolase for which a transglycosylation yield of 36% has been determined experimentally. The three transition states have a strong oxocarbenium character and ring conformations between ⁴H₃ and ⁴E. The atomic charges at the transition states for hydrolysis and transglycosylation are very similar, except for the more negative charge of the oxygen atom of water when compared to that of the acceptor Glc. The glycosylation transition state has a stronger S_N2 character than the deglycosylation ones and the proton transfer is less advanced. At the QM(PBE0/TZVP)/MM level, the TS for transglycosylation has shorter O4_GLC-C1_FUC (forming bond) distance and longer OE2_GLU338-C1_FUC (breaking) distance than the hydrolysis one, although the H_ACC proton is closer to the Glu164 base in the hydrolysis TS. The QM(SCC-DFTB)/MM free energy maxima show the inverted situation, although the hydrolysis TS presents significant structural fluctuations. The 3-OH_GLC group of the acceptor Glc (transglycosylation) and WAT432 (neighbor water in hydrolysis) are identified to stabilize the oxocarbenium transition states through interaction with O5_FUC and O4_FUC. The analysis of interaction suggests that perturbing the Glu392-Fuc interaction could increase the T/H ratio, either by direct mutation of this residue or indirectly as reported experimentally in the Asn390I and Phe401S cases. The molecular understanding of similarities and differences between hydrolysis and transglycosylation steps may be of help in the design of new biocatalysts for glycan synthesis.

The Driving Force for the Acylation of β-Lactam Antibiotics by L,D-Transpeptidase 2: Quantum Mechanics/Molecular Mechanics (QM/MM) Study

β-lactam antibiotics, which are used to treat infectious diseases, are currently the most widely used class of antibiotics. This study focused on the chemical reactivity of five- and six-membered ring systems attached to the β-lactam ring. The ring strain energy (RSE), force constant (FC) of amide (C-N), acylation transition states and second-order perturbation stabilization energies of 13 basic structural units of β-lactam derivatives were computed using the M06-2X and G3/B3LYP multistep method. In the ring strain calculations, an isodesmic reaction scheme was used to obtain the total energies. RSE is relatively greater in the five-(1a-2c) compared to the six-membered ring systems except for 4b, which gives a RSE that is comparable to five-membered ring lactams. These variations were also observed in the calculated inter-atomic amide bond distances (C-N), which is why the six-membered ring lactams C-N bond are more rigid than those with five-membered ring lactams. The calculated ΔG^# values from the acylation reaction of the lactams (involving the S-H group of the cysteine active residue from L,D transpeptidase 2) revealed a faster rate of C-N cleavage in the five-membered ring lactams especially in the 1-2 derivatives (17.58 kcal mol^-1 ). This observation is also reflected in the calculated amide bond force constant (1.26 mDyn/A) indicating a weaker bond strength, suggesting that electronic factors (electron delocalization) play more of a role on reactivity of the β-lactam ring, than ring strain.

Computational prediction of the bioactivity potential of proteomes based on expert knowledge

Advances in the field of genome sequencing have enabled a comprehensive analysis and annotation of the dynamics of the protein inventory of cells. This has been proven particularly rewarding for microbial cells, for which the majority of proteins are already accessible to analysis through automatic metagenome annotation. The large-scale in silico screening of proteomes and metaproteomes is key to uncover bioactivities of translational, clinical and biotechnological interest, and to help assign functions to certain proteins, such as those predicted as hypothetical. This work introduces a new method for the prediction of the bioactivity potential of proteomes/metaproteomes, supporting the discovery of functionally relevant proteins based on prior knowledge. This methodology complements functional annotation enrichment methods by allowing the assignment of functions to proteins annotated as hypothetical/putative/uncharacterised, as well as and enabling the detection of specific bioactivities and the recovery of proteins from defined taxa. This work shows how the new method can be applied to screen proteome and metaproteome sets to obtain predictions of clinical or biotechnological interest based on reference datasets. Notably, with this methodology, the large information files obtained after DNA sequencing or protein identification experiments can be associated for translational purposes that, in cases such as antibiotic-resistance pathogens or foodborne diseases, may represent changes in how these important and global health burdens are approached in the clinical practice. Finally, the Sequence-based Expert-driven pRoteome bioactivity Prediction EnvironmENT, a public Web service implemented in Scala functional programming style, is introduced as means to ensure broad access to the method as well as to discuss main implementation issues, such as modularity, extensibility and interoperability.

Rational design of the beta-galactosidase from Aspergillus oryzae to improve galactooligosaccharide production

The β-galactosidase from Aspergillus oryzae produces galactooligosaccharides with beneficial prebiotic effects through the transglycosylation of lactose. However, its low galactooligosaccharide yield greatly limits its use in industrial applications. Rational design and site-directed mutagenesis were used to produce three mutants of A. oryzae β-galactosidase: N140C, W806F, and N140C/W806F. The galactooligosaccharide yields of N140C (50.7%), W806F (49.3%), and N140C/W806F (59.8%) represent substantial improvements over that of the wild-type (35.7%). The galactooligosaccharide yield of N140C/W806F is the highest reported to date. Our rational design approach suggests novel strategies for further study of the β-galactosidase reaction mechanism and has produced mutants may be more useful in industrial applications than wild-type A. oryzae β-galactosidase.

α-d-Gal-cyclophellitol cyclosulfamidate is a Michaelis complex analog that stabilizes therapeutic lysosomal α-galactosidase A in Fabry disease

α-d-Gal-cyclophellitol cyclosulfamidate is a new class of neutral, conformationally-constrained competitive glycosidase inhibitor that stabilizes α-gal A and prevents its degradation both in vitro and in cellulo by mimicry of the Michaelis complex conformation.

Heterologous Expression of a Thermostable β-1,3-Galactosidase and Its Potential in Synthesis of Galactooligosaccharides

A thermostable β-1,3-galactosidase from Marinomonas sp. BSi20414 was successfully heterologously expressed in Escherichia coli BL21 (DE3), with optimum over-expression conditions as follows: the recombinant cells were induced by adding 0.1 mM of IPTG to the medium when the OD600 of the culture reached between 0.6 and 0.9, followed by 22 h incubation at 20 °C. The recombinant enzyme β-1,3-galactosidase (rMaBGA) was further purified to electrophoretic purity by immobilized metal affinity chromatography and size exclusion chromatography. The specific activity of the purified enzyme was 126.4 U mg-1 at 37 °C using ONPG (o-nitrophenyl-β-galactoside) as a substrate. The optimum temperature and pH of rMaBGA were determined as 60 °C and 6.0, respectively, resembling with its wild-type counterpart, wild type (wt)MaBGA. However, rMaBGA and wtMaBGA displayed different thermal stability and steady-state kinetics, although they share identical primary structures. It is postulated that the stability of the enzyme was altered by heterologous expression with the absence of post-translational modifications such as glycosylation, as well as the steady-state kinetics. To evaluate the potential of the enzyme in synthesis of galactooligosaccharides (GOS), the purified recombinant enzyme was employed to catalyze the transgalactosylation reaction at the lab scale. One of the transgalactosylation products was resolved as 3'-galactosyl-lactose, which had been proven to be a better bifidogenic effector than GOS with β-1,4 linkage and β-1,6 linkages. The results indicated that the recombinant enzyme would be a promising alternative for biosynthesis of GOS mainly with β-1,3 linkage.

The catalytic role of water in the binding site of l,d-transpeptidase 2 within acylation mechanism: A QM/MM (ONIOM) modelling

Mycobacterium tuberculosis is the causative agent of Tuberculosis. Formation of 3 → 3 crosslinks in the peptidoglycan layer of M. tuberculosis is catalyzed by l,d-transpeptidases. These enzymes can confer resistance against classical β-lactams that inhibit enzymes that generate 4 → 3 peptidoglycan crosslinks. The focus of this study is to investigate the catalytic role of water molecules in the acylation mechanism of the β-lactam ring within two models; 4- and 6-membered ring systems using two-layered our Own N-layer integrated Molecular Mechanics ONIOM (B3LYP/6-311++G(2d,2p): AMBER) model. The obtained thermochemical parameters revealed that the 6-membered ring model best describes the inhibition mechanism of acylation which indicates the role of water in the preference of 6-membered ring reaction pathway. This finding is in accordance with experimental data for the rate-limiting step of cysteine protease with the same class of inhibitor and binding affinity for both inhibitors. As expected, the ΔG^# results also reveal that the 6-membered ring reaction pathway is the most favourable. The electrostatic potential (ESP) and the natural bond orbital analysis (NBO) showed stronger interactions in 6-membered ring transition state (TS-6) mechanism involving water in the active site of the enzyme. This study could be helpful in the development of novel antibiotics against l,d-transpeptidase.

Oxalyltransferase, a plant cell-wall acyltransferase activity, transfers oxalate groups from ascorbate metabolites to carbohydrates

In the plant apoplast, ascorbate is oxidised, via dehydroascorbic acid, to O-oxalyl esters [oxalyl-l-threonate (OxT) and cyclic oxalyl-l-threonate (cOxT)]. We tested whether OxT and cOxT can donate the oxalyl group in transacylation reactions to form oxalyl-polysaccharides, potentially modifying the cell wall. [oxalyl-14 C]OxT was incubated with living spinach (Spinacia oleracea) and Arabidopsis cell-suspension cultures in the presence or absence of proposed acceptor substrates (carbohydrates). In addition, [14 C]OxT and [14 C]cOxT were incubated in vitro with cell-wall enzyme preparations plus proposed acceptor substrates. Radioactive products were monitored electrophoretically. Oxalyltransferase activity was detected. Living cells incorporated oxalate groups from OxT into cell-wall polymers via ester bonds. When sugars were added, [14 C]oxalyl-sugars were formed, in competition with OxT hydrolysis. Preferred acceptor substrates were carbohydrates possessing primary alcohols e.g. glucose. A model transacylation product, [14 C]oxalyl-glucose, was relatively stable in vivo (half-life >24 h), whereas [14 C]OxT underwent rapid turnover (half-life ~6 h). Ionically wall-bound enzymes catalysed similar transacylation reactions in vitro with OxT or cOxT as oxalyl donor substrates and any of a range of sugars or hemicelluloses as acceptor substrates. Glucosamine was O-oxalylated, not N-oxalylated. We conclude that plants possess apoplastic acyltransferase (oxalyltransferase) activity that transfers oxalyl groups from ascorbate catabolites to carbohydrates, forming relatively long-lived O-oxalyl-carbohydrates. The findings increase the range of known metabolites whose accumulation in vivo indicates vitamin C catabolism. Possible signalling roles of the resulting oxalyl-sugars can now be investigated, as can the potential ability of polysaccharide oxalylation to modify the wall's physical properties.

Data characterizing the energetics of enzyme-catalyzed hydrolysis and transglycosylation reactions by DFT cluster model calculations

The data presented in this paper are related to the research article entitled “QM/MM modeling of the hydrolysis and transfructosylation reactions of fructosyltransferase from Aspergillus japonicas, an enzyme that produces prebiotic fructooligosaccharide” (Jitonnom et al., 2018) [1]. This paper presents the procedure and data for characterizing the whole relative energy profiles of hydrolysis and transglycosylation reactions whose elementary steps differ in chemical composition. The data also reflects the choices of the QM cluster model, the functional/basis set method and the equations in determining the reaction energetics.

The role of catalytic residue pKa on the hydrolysis/transglycosylation partition in family 3 β-glucosidases

The catalytic acid/base residue (E441) pK_a, as modulated by its electrostatic environment, has potential impact on the hydrolysis/transglycosylation partition in β-glucosidases.

Establishing the catalytic mechanism of human pancreatic α-amylase with QM/MM methods

In this work, we studied the catalytic mechanism of human pancreatic α-amylase (HPA). Our goal was to determine the catalytic mechanism of HPA with atomic detail using computational methods. We demonstrated that the HPA catalytic mechanism consists of two steps, the first of which (glycosylation step) involves breaking the glycosidic bond to culminate in the formation of a covalent intermediate. The second (deglycosylation step) consists of the addition of a water molecule to release the enzyme/substrate covalent intermediate, completing the hydrolysis of the sugar. The active site was very open to the solvent. Our mechanism basically differs from the previously proposed mechanism by having two water molecules instead of only one near the active site that participate in the mechanism. We also demonstrate the relevant role of the three catalytic amino acids, two aspartate residues and a glutamate (D197, E233, and D300), during catalysis. It was also shown that the rate limiting step was glycosylation, and its activation energy was in agreement with experimental values obtained for HPA. The experimental activation energy was 14.4 kcal mol(-1), and the activation energy obtained computationally was 15.1 kcal mol(-1).

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Enzymatic Synthesis of Galactosylated Serine/Threonine Derivatives by β-Galactosidase from Escherichia coli

The transgalactosylations of serine/threonine derivatives were investigated using β-galactosidase from Escherichia coli as biocatalyst. Using ortho-nitrophenyl-β-d-galactoside as donor, the highest bioconversion yield of transgalactosylated N-carboxy benzyl l-serine benzyl ester (23.2%) was achieved in heptane:buffer medium (70:30), whereas with the lactose, the highest bioconversion yield (3.94%) was obtained in the buffer reaction system. The structures of most abundant galactosylated serine products were characterized by MS/MS. The molecular docking simulation revealed that the binding of serine/threonine derivatives to the enzyme’s active site was stronger (−4.6~−7.9 kcal/mol) than that of the natural acceptor, glucose, and mainly occurred through interactions with aromatic residues. For N-tert-butoxycarbonyl serine methyl ester (6.8%) and N-carboxybenzyl serine benzyl ester (3.4%), their binding affinities and the distances between their hydroxyl side chain and the 1′-OH group of galactose moiety were in good accordance with the quantified bioconversion yields. Despite its lower predicted bioconversion yield, the high experimental bioconversion yield obtained with N-carboxybenzyl serine methyl ester (23.2%) demonstrated the importance of the thermodynamically-driven nature of the transgalactosylation reaction.

Glycosylation of phenolic compounds by the site-mutated β-galactosidase from Lactobacillus bulgaricus L3

β-Galactosidases can transfer the galactosyl from lactose or galactoside donors to various acceptors and thus are especially useful for the synthesis of important glycosides. However, these enzymes have limitations in the glycosylation of phenolic compounds that have many physiological functions. In this work, the β-galactosidase from Lactobacillus bulgaricus L3 was subjected to site-saturation mutagenesis at the W980 residue. The recombinant pET-21b plasmid carrying the enzyme gene was used as the template for mutation. The mutant plasmids were transformed into Escherichia coli cells for screening. One recombinant mutant, W980F, exhibited increased yield of glycoside when using hydroquinone as the screening acceptor. The enzyme was purified and the effects of the mutation on enzyme properties were determined in detail. It showed improved transglycosylation activity on novel phenolic acceptors besides hydroquinone. The yields of the glycosides produced from phenol, hydroquinone, and catechol were increased by 7.6% to 53.1%. Moreover, it generated 32.3% glycosides from the pyrogallol that could not be glycosylated by the wild-type enzyme. Chemical structures of these glycoside products were further determined by MS and NMR analysis. Thus, a series of novel phenolic galactosides were achieved by β-galactosidase for the first time. This was a breakthrough in the enzymatic galactosylation of the challenging phenolic compounds of great values.