Pyranose Dehydrogenase Ligand Promiscuity: A Generalized Approach to Simulate Monosaccharide Solvation

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 (DGsolv) of all 32 aldohexopyranoses in water, which have not yet been reported in the literature. The free energy difference between b- and a-anomers (DGb-a) of all D-stereoisomers in water were compared to experimental values with a good agreement. Moreover, the free-energy differences (DG) of the 32 stereoisomers bound to PDH in two different poses were calculated from MD simulations. The relative binding free energies (DDGbind) were calculated and, where available, compared to experimental values, approximated from Km 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 DDGbind and a reweighted distance analysis. Furthermore, several oxidation products were predicted for sugars that hav