Pyranose Dehydrogenase Ligand Promiscuity: A Generalized Approach to Simulate Monosaccharide Solvation, Binding, and Product Formation

The flavoenzyme pyranose dehydrogenase (PDH) from the litter decomposing fungus Agaricus meleagris oxidizes many different carbohydrates occurring during lignin degradation. This promiscuous substrate specificity makes PDH a promising catalyst for bioelectrochemical applications. A generalized approach to simulate all 32 possible aldohexopyranoses in the course of one or a few molecular dynamics (MD) simulations is reported. Free energy calculations according to the one-step perturbation (OSP) method revealed the solvation free energies (ΔGsolv) of all 32 aldohexopyranoses in water, which have not yet been reported in the literature. The free energy difference between β- and α-anomers (ΔGβ-α) of all d-stereoisomers in water were compared to experimental values with a good agreement. Moreover, the free-energy differences (ΔG) of the 32 stereoisomers bound to PDH in two different poses were calculated from MD simulations. The relative binding free energies (ΔΔGbind) were calculated and, where available, compared to experimental values, approximated from K m values. The agreement was very good for one of the poses, in which the sugars are positioned in the active site for oxidation at C1 or C2. Distance analysis between hydrogens of the monosaccharide and the reactive N5-atom of the flavin adenine dinucleotide (FAD) revealed that oxidation is possible at HC1 or HC2 for pose A, and at HC3 or HC4 for pose B. Experimentally detected oxidation products could be rationalized for the majority of monosaccharides by combining ΔΔGbind and a reweighted distance analysis. Furthermore, several oxidation products were predicted for sugars that have not yet been tested experimentally, directing further analyses. This study rationalizes the relationship between binding free energies and substrate promiscuity in PDH, providing novel insights for its applicability in bioelectrochemistry. The results suggest that a similar approach could be applied to study promiscuity of other enzymes.

[1]  Wilfred F van Gunsteren,et al.  Free energies of binding of polychlorinated biphenyls to the estrogen receptor from a single simulation , 2003, Proteins.

[2]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[3]  Chris Oostenbrink,et al.  Free‐energy differences between states with different conformational ensembles , 2013, J. Comput. Chem..

[4]  Riccardo Baron,et al.  Computational Drug Discovery and Design , 2012, Methods in Molecular Biology.

[5]  Tim N. Heinz,et al.  Comparison of four methods to compute the dielectric permittivity of liquids from molecular dynamics simulations , 2001 .

[6]  J. Volc,et al.  Pyranose dehydrogenases: biochemical features and perspectives of technological applications , 2009, Applied Microbiology and Biotechnology.

[7]  Dan S. Tawfik,et al.  Enzyme promiscuity: evolutionary and mechanistic aspects. , 2006, Current opinion in chemical biology.

[8]  Harel Weinstein,et al.  Thermodynamic basis for promiscuity and selectivity in protein-protein interactions: PDZ domains, a case study. , 2006, Journal of the American Chemical Society.

[9]  W. V. van Gunsteren,et al.  Prediction of folding equilibria of differently substituted peptides using one-step perturbation. , 2010, Journal of the American Chemical Society.

[10]  C. Oostenbrink.  Free energy calculations from one-step perturbations. , 2012, Methods in molecular biology.

[11]  D. Haltrich,et al.  Regioselective control of β-d-glucose oxidation by pyranose 2-oxidase is intimately coupled to conformational degeneracy. , 2011, Journal of molecular biology.

[12]  Markus Christen,et al.  Architecture, implementation and parallelisation of the GROMOS software for biomolecular simulation , 2012, Comput. Phys. Commun..

[13]  L. Wackett.  Evolution of Enzymes for the Metabolism of New Chemical Inputs into the Environment* , 2004, Journal of Biological Chemistry.

[14]  Rengaswami Chandrasekaran,et al.  Conformation of Carbohydrates , 1998 .

[15]  Alexei Kurakin,et al.  Self‐organization versus Watchmaker: ambiguity of molecular recognition and design charts of cellular circuitry , 2007, Journal of molecular recognition : JMR.

[16]  U. Bornscheuer,et al.  Catalytic promiscuity in biocatalysis: using old enzymes to form new bonds and follow new pathways. , 2004, Angewandte Chemie.

[17]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[18]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[19]  R. Jensen.  Enzyme recruitment in evolution of new function. , 1976, Annual review of microbiology.

[20]  Dan S. Tawfik,et al.  Conformational diversity and protein evolution--a 60-year-old hypothesis revisited. , 2003, Trends in biochemical sciences.

[21]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[22]  R. Mariuzza.  Multiple paths to multispecificity. , 2006, Immunity.

[23]  S. Copley.  Enzymes with extra talents: moonlighting functions and catalytic promiscuity. , 2003, Current opinion in chemical biology.

[24]  B. Henrissat,et al.  Recent structural insights into the expanding world of carbohydrate-active enzymes. , 2005, Current opinion in structural biology.

[25]  D. Haltrich,et al.  Engineering of Pyranose Dehydrogenase for Increased Oxygen Reactivity , 2014, PloS one.

[26]  Dan S. Tawfik,et al.  Enzyme promiscuity: a mechanistic and evolutionary perspective. , 2010, Annual review of biochemistry.

[27]  Gábor Nagy,et al.  Entropic and Enthalpic Contributions to Stereospecific Ligand Binding from Enhanced Sampling Methods , 2014, J. Chem. Inf. Model..

[28]  D. Haltrich,et al.  Characterization of pyranose dehydrogenase from Agaricus meleagris and its application in the C-2 specific conversion of D-galactose. , 2008, Journal of biotechnology.

[29]  D. Haltrich,et al.  Properties of pyranose dehydrogenase purified from the litter‐degrading fungus Agaricus xanthoderma , 2007, The FEBS journal.

[30]  Jozef Hritz,et al.  Malleability and versatility of cytochrome P450 active sites studied by molecular simulations. , 2012, Current drug metabolism.

[31]  Alan E. Mark,et al.  Estimating the Relative Free Energy of Different Molecular States with Respect to a Single Reference State , 1996 .

[32]  D. Wood,et al.  Pyranose 2-dehydrogenase, a novel sugar oxidoreductase from the basidiomycete fungus Agaricus bisporus , 1997, Archives of Microbiology.

[33]  C. Oostenbrink,et al.  On the use of one-step perturbation to investigate the dependence of NOE-derived atom–atom distance bound violations of peptides upon a variation of force-field parameters , 2014, European Biophysics Journal.

[34]  Andrea Amadei,et al.  Molecular dynamics simulations with constrained roto-translational motions: Theoretical basis and statistical mechanical consistency , 2000 .

[35]  R. Zwanzig.  High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[36]  Ian A. Wilson,et al.  Molecular basis of crossreactivity and the limits of antibody–antigen complementarity , 1993, Nature.

[37]  K. Hult,et al.  Enzyme promiscuity: mechanism and applications. , 2007, Trends in biotechnology.

[38]  Jozef Hritz,et al.  Efficient free energy calculations for compounds with multiple stable conformations separated by high energy barriers. , 2009, The journal of physical chemistry. B.

[39]  Angelo D. Favia,et al.  Protein promiscuity and its implications for biotechnology , 2009, Nature Biotechnology.

[40]  D. Haltrich,et al.  New biotransformations of some reducing sugars to the corresponding (di)dehydro(glycosyl) aldoses or aldonic acids using fungal pyranose dehydrogenase , 2006 .

[41]  R. Raines,et al.  Identifying latent enzyme activities: substrate ambiguity within modern bacterial sugar kinases. , 2004, Biochemistry.

[42]  D. Haltrich,et al.  The 1.6 Å Crystal Structure of Pyranose Dehydrogenase from Agaricus meleagris Rationalizes Substrate Specificity and Reveals a Flavin Intermediate , 2013, PloS one.

[43]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[44]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[45]  Chris Oostenbrink,et al.  Molecular dynamics simulations and free energy calculations on the enzyme 4‐hydroxyphenylpyruvate dioxygenase , 2011, J. Comput. Chem..

[46]  Juhan Kim,et al.  Why metabolic enzymes are essential or nonessential for growth of Escherichia coli K12 on glucose. , 2007, Biochemistry.

[47]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[48]  D. Haltrich,et al.  Screening of basidiomycete fungi for the quinone-dependent sugar C-2/C-3 oxidoreductase, pyranose dehydrogenase, and properties of the enzyme from Macrolepiota rhacodes , 2001, Archives of Microbiology.

[49]  Dietmar Haltrich,et al.  Molecular dynamics simulations give insight into d-glucose dioxidation at C2 and C3 by Agaricus meleagris pyranose dehydrogenase , 2013, Journal of Computer-Aided Molecular Design.

[50]  Chris Oostenbrink,et al.  Free Energy Calculations Give Insight into the Stereoselective Hydroxylation of α-Ionones by Engineered Cytochrome P450 BM3 Mutants , 2012, J. Chem. Inf. Model..

[51]  Ken A Dill,et al.  Use of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations. , 2007, Journal of chemical theory and computation.

[52]  D. Cavener,et al.  GMC oxidoreductases. A newly defined family of homologous proteins with diverse catalytic activities. , 1992, Journal of molecular biology.

[53]  H. Gabius.  The sugar code : fundamentals of glycosciences , 2009 .

[54]  Philippe H. Hünenberger,et al.  Using the local elevation method to construct optimized umbrella sampling potentials: Calculation of the relative free energies and interconversion barriers of glucopyranose ring conformers in water , 2010, J. Comput. Chem..

[55]  W. Whelan.  Carbohydrate Chemistry and Biochemistry , 1963, Nature.

[56]  D. Haltrich,et al.  Structural Basis for Substrate Binding and Regioselective Oxidation of Monosaccharides at C3 by Pyranose 2-Oxidase* , 2006, Journal of Biological Chemistry.

[57]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[58]  Roberto D. Lins,et al.  A new GROMOS force field for hexopyranose‐based carbohydrates , 2005, J. Comput. Chem..