Conformational Free Energy Surface of α-N-Acetylneuraminic Acid: An Interplay Between Hydrogen Bonding and Solvation

The conformational free energy surface of α-N-acetylneuraminic acid (Neu5Ac, sialic acid) in the space of ring-puckering coordinates was calculated using the metadynamics method. Free energy surfaces in vacuum and with an explicit solvent were calculated in GLYCAM 06 force field. In vacuum three structures are almost equivalently populated, namely, the ²C₅ chair and the B_3,6/²S₆ and ^OS₃ boat/skew-boat conformations. The B_3,6/²S₆ structure is stabilized by an ionic hydrogen bond between the amide N−H bond and the carboxylic group. However, this structure is unfavorable in a water environment in which the experimentally observed ²C₅ chair conformation is predicted to be more stable than the other structures. These results indicate that environment significantly influences conformation of Neu5Ac and that Neu5Ac-processing enzymes might modify a conformation of their substrates solely by a changing polarity of the environment. The structure of Neu5Ac bound in influenza neuraminidase (⁴S₂/B_2,5) belongs to conformations preferred in a water environment. The free energy penalty of this conformational change was calculated (relative to ²C₅) as 10.2 ± 2.0 and 17.3 ± 2.0 kJ/mol for ^4,OB/^OS₃ and ⁴S₂, respectively. This result indicates that mimicking of the enzyme-bound conformation is likely to be a viable strategy for the design of neuraminidase inhibitors.