DFT studies on the β-glycosidase catalytic mechanism: The deglycosylation step

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.

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.

Impaired coordination of nucleophile and increased hydrophobicity in the +1 subsite shift levansucrase activity towards transfructosylation

Bacterial levansucrases produce β(2,6)-linked levan-type polysaccharides using sucrose or sucrose analogs as donor/acceptor substrates. However, the dominant reaction of Bacillus megaterium levansucrase (Bm-LS) is hydrolysis. Single domain levansucrases from Gram-positive bacteria display a wide substrate-binding pocket with open access to water, challenging engineering for transfructosylation-efficient enzymes. We pursued a shift in reaction specificity by either modifying the water distribution in the active site or the coordination of the catalytic acid/base (E352) and the nucleophile (D95), thus affecting the fructosyl-transfer rate and allowing acceptors other than water to occupy the active site. Two serine (173/422) and two water-binding tyrosine (421/439) residues located in the first shell of the catalytic pocket were modified. Library variants of S173, Y421 and S422, which coordinate the position of D95 and E352, show increased transfructosylation (30–200%) and modified product spectra. Substitutions at position 422 have a higher impact on sucrose affinity, while changes at position 173 and 421 have a strong effect on the overall catalytic rate. As most retaining glycoside hydrolases (GHs) Bm-LS catalyzes hydrolysis and transglycosylation via a double displacement reaction involving two-transition states (TS1 and TS2). Hydrogen bonds of D95 with the side chains of S173 and S422 contribute a total of 2.4 kcal mol⁻¹ to TS1 stabilization, while hydrogen bonds between invariant Y421, E352 and the glucosyl C-2 hydroxyl-group of sucrose contribute 2.15 kcal mol⁻¹ stabilization. Changes at Y439 render predominantly hydrolytic variants synthesizing shorter oligosaccharides.

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.

Glycosynthesis in a waterworld: new insight into the molecular basis of transglycosylation in retaining glycoside hydrolases

Carbohydrates are ubiquitous in Nature and play vital roles in many biological systems. Therefore the synthesis of carbohydrate-based compounds is of considerable interest for both research and commercial purposes. However, carbohydrates are challenging, due to the large number of sugar subunits and the multiple ways in which these can be linked together. Therefore, to tackle the challenge of glycosynthesis, chemists are increasingly turning their attention towards enzymes, which are exquisitely adapted to the intricacy of these biomolecules. In Nature, glycosidic linkages are mainly synthesized by Leloir glycosyltransferases, but can result from the action of non-Leloir transglycosylases or phosphorylases. Advantageously for chemists, non-Leloir transglycosylases are glycoside hydrolases, enzymes that are readily available and exhibit a wide range of substrate specificities. Nevertheless, non-Leloir transglycosylases are unusual glycoside hydrolases in as much that they efficiently catalyse the formation of glycosidic bonds, whereas most glycoside hydrolases favour the mechanistically related hydrolysis reaction. Unfortunately, because non-Leloir transglycosylases are almost indistinguishable from their hydrolytic counterparts, it is unclear how these enzymes overcome the ubiquity of water, thus avoiding the hydrolytic reaction. Without this knowledge, it is impossible to rationally design non-Leloir transglycosylases using the vast diversity of glycoside hydrolases as protein templates. In this critical review, a careful analysis of literature data describing non-Leloir transglycosylases and their relationship to glycoside hydrolase counterparts is used to clarify the state of the art knowledge and to establish a new rational basis for the engineering of glycoside hydrolases.

The catalytic mechanism of mouse renin studied with QM/MM calculations

Hypertension is a chronic condition that affects nearly 25% of adults worldwide. As the Renin-Angiotensin-Aldosterone System is implicated in the control of blood pressure and body fluid homeostasis, its combined blockage is an attractive therapeutic strategy currently in use for the treatment of several cardiovascular conditions. We have performed QM/MM calculations to study the mouse renin catalytic mechanism in atomistic detail, using the N-terminal His6-Asn14 segment of angiotensinogen as substrate. The enzymatic reaction (hydrolysis of the peptidic bond between residues in the 10th and 11th positions) occurs through a general acid/base mechanism and, surprisingly, it is characterized by three mechanistic steps: it begins with the creation of a first very stable tetrahedral gem-diol intermediate, followed by protonation of the peptidic bond nitrogen, giving rise to a second intermediate. In a final step the peptidic bond is completely cleaved and both gem-diol hydroxyl protons are transferred to the catalytic dyad (Asp32 and Asp215). The final reaction products are two separate peptides with carboxylic acid and amine extremities. The activation energy for the formation of the gem-diol intermediate was calculated as 23.68 kcal mol(-1), whereas for the other steps the values were 15.51 kcal mol(-1) and 14.40 kcal mol(-1), respectively. The rate limiting states were the reactants and the first transition state. The associated barrier (23.68 kcal mol(-1)) is close to the experimental values for the angiotensinogen substrate (19.6 kcal mol(-1)). We have also tested the influence of the density functional on the activation and reaction energies. All eight density functionals tested (B3LYP, B3LYP-D3, X3LYP, M06, B1B95, BMK, mPWB1K and B2PLYP) gave very similar results.