How sugars pucker: electronic structure calculations map the kinetic landscape of five biologically paramount...

Glycoside hydrolases (GHs) distort carbohydrate ring geometry along particular "catalytic itineraries" during the cleavage of glycosidic bonds, illustrating the relationship between substrate conformation and reactivity. Previous theoretical studies of thermodynamics of isolated monosaccharides offer insights into the catalytic itineraries of particular sugars. However, kinetic accessibility of carbohydrate puckering conformations and the role of exocyclic groups have not yet been thoroughly addressed. Here we present the first complete library of low-energy local minima and puckering interconversion transition states for five biologically relevant pyranose sugars: β-xylose, β-mannose, α-glucose, β-glucose, and β-N-acetylglucosamine. These were obtained by a thorough theoretical investigation each of the 38 IUPAC designated puckering geometries and all possible conformations of the exocyclic groups. These calculations demonstrate that exocyclic groups must be explicitly considered when examining these interconversion pathways. Furthermore, these data enable evaluation of previous hypotheses of why enzymes perturb ring geometries from the low-energy equatorial chair ((4)C1) conformation. They show that the relative thermodynamics alone do not universally correlate with GH catalytic itineraries. For some sugars, particular puckers offer both catalytically favorable electronic structure properties, such as anomeric carbon partial charge, and low kinetic barriers to achieve a given puckering conformation. However, different factors correlate with catalytic itineraries for other sugars; for β-N-acetylglucosamine, the key N-acetyl arm confounds the puckering landscape and appears to be the crucial factor. Overall, this study reveals a more comprehensive understanding of why particular puckering geometries are favored in carbohydrate catalysis concomitant with the complexity of glycobiology.