Ian R. Greig

Computational mutagenesis reveals the role of active-site tyrosine in stabilising a boat conformation for the substrate: QM/MM molecular dynamics studies of wild-type and mutant xylanases

Molecular dynamics simulations have been performed for non-covalent complexes of phenyl beta-xylobioside with the retaining endo-beta-1,4-xylanase from B. circulans (BCX) and its Tyr69Phe mutant using a hybrid QM/MM methodology. A trajectory initiated for the wild-type enzyme-substrate complex with the proximal xylose ring bound at the -1 subsite (adjacent to the scissile glycosidic bond) in the (4)C(1) chair conformation shows spontaneous transformation to the (2,5)B boat conformation, and potential of mean force calculations indicate that the boat is approximately 30 kJ mol(-1) lower in free energy than the chair. Analogous simulations for the mutant lacking one oxygen atom confirm the key role of Tyr69 in stabilizing the boat in preference to the (4)C(1) chair conformation, with a relative free energy difference of about 20 kJ mol(-1), by donating a hydrogen bond to the endocyclic oxygen of the proximal xylose ring. QM/MM MD simulations for phenyl beta-xyloside in water, with and without a propionate/propionic acid pair to mimic the catalytic glutamate/glutamic acid pair of the enzyme, show the (4)C(1) chair to be stable, although a hydrogen bond between the OH group at C2 of xylose and the propionate moiety seems to provide some stabilization for the (2,5)B conformation.

Probing Synergy between Two Catalytic Strategies in the Glycoside Hydrolase O-GlcNAcase Using Multiple Linear Free Energy Relationships

Human O-GlcNAcase plays an important role in regulating the post-translational modification of serine and threonine residues with beta-O-linked N-acetylglucosamine monosaccharide unit (O-GlcNAc). The mechanism of O-GlcNAcase involves nucleophilic participation of the 2-acetamido group of the substrate to displace a glycosidically linked leaving group. The tolerance of this enzyme for variation in substrate structure has enabled us to characterize O-GlcNAcase transition states using several series of substrates to generate multiple simultaneous free-energy relationships. Patterns revealing changes in mechanism, transition state, and rate-determining step upon concomitant variation of both nucleophilic strength and leaving group abilities are observed. The observed changes in mechanism reflect the roles played by the enzymic general acid and the catalytic nucleophile. Significantly, these results illustrate how the enzyme synergistically harnesses both modes of catalysis; a feature that eludes many small molecule models of catalysis. These studies also suggest the kinetic significance of an oxocarbenium ion intermediate in the O-GlcNAcase-catalyzed hydrolysis of glucosaminides, probing the limits of what may be learned using nonatomistic investigations of enzymic transition-state structure and offering general insights into how the superfamily of retaining glycoside hydrolases act as efficient catalysts.

Mechanism of glycoside hydrolysis: A comparative QM/MM molecular dynamics analysis for wild type and Y69F mutant retaining xylanases

Computational simulations have been performed using hybrid quantum-mechanical/molecular-mechanical potentials to investigate the catalytic mechanism of the retaining endo-beta-1, 4-xylanase (BCX) from B. circulans. Two-dimensional potential-of-mean-force calculations based upon molecular dynamics with the AM1/OPLS method for wild-type BCX with a p-nitrophenyl xylobioside substrate in water clearly indicates a stepwise mechanism for glycosylation: the rate-determining step is nucleophilic substitution by Glu78 to form the covalently bonded enzyme-substrate intermediate without protonation of the leaving group by Glu172. The geometrical configuration of the transition state for the enzymic reaction is essentially the same as found for a gas-phase model involving only the substrate and a propionate/propionic acid pair to represent the catalytic glutamate/glutamic acid groups. In addition to stabilizing the (2,5)B boat conformation of the proximal xylose in the non-covalent reactant complex of the substrate with BCX, Tyr69 lowers the free-energy barrier for glycosylation by 42 kJ mol(-1) relative to that calculated for the Y69F mutant, which lacks the oxygen atom O(Y). B3LYP/6-31+G* energy corrections reduce the absolute height of the barrier to reaction. In the oxacarbenium ion-like transition state O(Y) approaches closer to the endocyclic oxygen O(ring) of the sugar ring but donates its hydrogen bond not to O(ring) but rather to the nucleophilic oxygen of Glu78. Comparison of the average atomic charge distributions for the wild-type and mutant indicates that charge separation along the bond between the anomeric carbon and O(ring) is matched in the former by a complementary separation of charge along the O(Y)-H(Y) bond, corresponding to a pair of roughly antiparallel bond dipoles, which is not present in the latter.

Elucidating the Nature of the Streptomyces plicatus β-Hexosaminidase-Bound Intermediate Using ab initio Molecular Dynamics Simulations

By using all-atom ab initio molecular dynamics simulations, the solution pK(a) of the oxazolinium ion intermediate formed during the Streptomyces plicatus beta-hexosaminidase (SpHex)-catalyzed hydrolysis of beta-D-N-acetylglucosaminides is estimated as pK(a) = 7.7. The structure and protonation state of the enzyme-bound intermediate have also been investigated, using hybrid QM/MM methods. The protonation state and conformational properties of the enzyme bound intermediate are found to be sensitive to the protonation state of a number of ionisable residues (other than the aspartate-glutamate catalytic dyad) suggesting that the microscopic electrostatic environment of SpHex not only perturbs the relative magnitudes of the pK(a) values of the Asp side chain carboxylate and oxazolinium ion but also that SpHex binds its intermediate in a distorted conformation with respect to its ground-state conformation in solution.

Structural, mechanistic, and computational analysis of the effects of anomeric fluorines on anomeric fluoride departure in 5-fluoroxylosyl fluorides

The effects of fluorine substitution at the C-5 center of pyranosyl fluorides on the reactivity at the C-1 anomeric center was probed by studying a series of 5-fluoroxylosyl fluoride derivatives. X-ray structures of their per-O-acetates detailed the effects on the ground-state structures. First-order rate constants for spontaneous hydrolysis, in conjunction with computational studies, revealed that changes in the stereochemistry of the 5-fluorine had minimal effects on the solvolysis rate constants and that the observed rate reductions were broadly similar to those caused by additional fluorine substitution at C-1 but significantly less than those due to substitution at C-2. Differences in the trapping behavior of 5- versus 2-fluoro-substituted glycosyl fluorides with α- and β-glycosidases arise more from differences in steric effects and hydrogen-bonding interactions than from intrinsic reactivity differences.

Tuning Mechanism‐Based Inactivators of Neuraminidases: Mechanistic and Structural Insights

3-Fluorosialosyl fluorides are inhibitors of sialidases that function by the formation of a long-lived covalent active-site adduct and have potential as therapeutics if made specific for the pathogen sialidase. Surprisingly, human Neu2 and the Trypanosoma cruzi trans-sialidase are inactivated more rapidly by the reagent with an equatorial fluorine at C3 than by its axial epimer, with reactivation following the same pattern. To explore a possible stereoelectronic basis for this, rate constants for spontaneous hydrolysis of the full series of four 3-fluorosialosyl fluorides were measured, and ground-state energies for each computed. The alpha (equatorial) anomeric fluorides hydrolyze more rapidly than their beta anomers, consistent with their higher ground-state energies. However ground-state energies do not explain the relative spontaneous reactivities of the 3-fluoro-epimers. The three-dimensional structures of the two 3-fluoro-sialosyl enzyme intermediates of human Neu2 were solved, revealing key stabilizing interactions between Arg21 and the equatorial, but not the axial, fluorine. Because of changes in geometry these interactions will increase at the transition state, likely explaining the difference in reaction rates.

Solvent Effects on Isotope Effects: Methyl Cation as a Model System

The isotopic sensitivity (CH3(+) vs CD3(+)) of the equilibrium between the methyl cation in vacuum and in solution has been investigated. Two alternative options for describing the shape of the solute cavity within the widely used polarized continuum model for implicit solvation were compared; the UFF and UA0 methods give equilibrium isotope effects (EIEs) that vary as a function of the dielectric constant in opposite directions. The same isotope effect was also obtained as the average over 40 structures from a hybrid quantum mechanical/molecular mechanical molecular dynamics simulation for the methyl cation explicitly solvated by many water molecules; the inverse value of the EIE agrees with UFF but not UA0. The opposing trends may be satisfactorily explained in terms of the different degrees of exposure of the atomic charges to the dielectric continuum in cavities of different shapes.

Remarkable Reactivity Differences between Glucosides with Identical Leaving Groups

Two isomeric aryl 2-deoxy-2-fluoro-β-glucosides react with a β-glucosidase at rates differing by 106-fold, despite the fact that they release the same aromatic aglycone. In contrast, the equivalent glucoside substrates react with essentially identical rate constants. Insight into the source of these surprising rate differences was obtained through a comprehensive study of the nonenzymatic (spontaneous) hydrolysis of these same substrates, wherein an approximate 105-fold difference in rates was measured, clarifying that the differences were inherent rather than being due to specific interactions with the enzyme. The possibility that an alternate nucleophilic aryl substitution mechanism was responsible for the rapid reaction of the faster substrate was excluded through 18O-labeling studies. Further exploration of the origins of these rate differences involved analysis of X-ray crystal structures as well as quantum chemical calculations, which surprisingly revealed that ground state destabilization and transition state stabilizing effects contribute almost equally to the observed reactivity differences. These studies highlight the dangers of using simple reference equilibria such as pKa values as measures of leaving group ability.

Carbohydrate Chemistry and Biochemistry, by Michael L. Sinnott

Carbohydrate Chemistry and Biochemistry represents a significant, single-author, synthesis of knowledge weighing in at 748 pages and containing only seven chapters. The author, Professor Michael Sinnott, has contributed many classic research papers in the field of the carbohydrate reactivity pertaining to the details of glycosyl hydrolase mechanism (1), the kinetic characterization of enzyme activity in evolving systems (2) and the preservation of knowledge in its physical forms (3). This text extends his prior wide-ranging reviews of the mechanisms of carbohydrate-active enzymes (4, 5), adopting a more tutorial approach. It represents, in the main, a detailed consideration of the literature till the end of 2006. Chapters 1 and 2 contain a conventional introduction to the structure and conformation of monosaccharides, with mutarotation being dealt with in detail. Chapter 3 focuses on the physical organic chemistry of nucleophilic substitution at the anomeric centre. It presages Chapter 5’s extensive treatment of enzyme-catalyzed glycosyl transfer. In between, Chapter 4 deals with primary structures and conformations of polysaccharides and includes a discussion of physical methods, including X-ray crystallography. Chapter 6 is a necessary collection of material, covering both synthetic and biologically-relevant heterolytic chemistry not involving the anomeric centre; in some topics (e.g. hydride transfer from NAD(P)H) discussed in this chapter, carbohydrates play only a very minor role. Chapter 7 deals with reactions involving homolytic chemistry. As noted in the author’s preface, this book was originally intended as a contribution to an RSC-sponsored series of volumes on physical organic chemistry but the untimely death of the series editor, Professor Andrew Williams, resulted in its publication as a stand-alone volume. This provenance has produced a book which extensively draws on the approaches of physical organic chemistry, applying them to the study of carbohydrate-active enzymecatalyzed reactions. The enzyme-catalyzed chemistry of the anomeric centre represents the author’s particular area of expertise, and Chapter 5 contains extensive, personal, valuable commentary on numerous studies in this area. The author provides a detailed, tutorial overview of the area, and so the systems chosen, though extensive, are not necessarily comprehensive. The historical overview of the enzymology of glycosyl hydrolases provided by this chapter is particulary interesting, highlighting several significant and perhaps under-appreciated papers. Sinnott writes with a fluid style throughout, and where he disagrees with the conclusions drawn by authors of the primary research papers, he makes clear why he does so. This encourages the reader to question the conclusions that can reasonably be drawn from experimental data and which experiments may be used to choose amongst competing mechanisms. The general approach adopted towards all enzyme-catalyzed glycosyl group transfers that occur with retention of chemistry should be readily approachable to those involved in the structural elucidation of stable enzyme-intermediate complexes. It assumes that if an