Microsecond-scale sucrose conformational dynamics in aqueous solution via molecular dynamics methods

Molecular dynamics methods have proven their applicability for the reproduction and prediction of molecular conformations during the past decades. However, most of works considered dilute solutions and relatively short trajectories that limit insights into conformational dynamics. In this study, we investigate the conformational dynamics of sucrose in aqueous solution using microsecond-scale molecular dynamics simulations. For the most of the calculations, we use the OPLS-AA/1.14*CM1A-LBCC force field, but we also utilize OPLS-AA/1.14*CM1A and GLYCAM06 for the comparison. We focused on the glycosidic linkage conformers and their lifetimes, glucopyranose and fructofuranose ring puckering. Our findings indicate that the 1C4 glucopyranose ring conformation can stabilize the sucrose conformer, appeared only in the GLYCAM06 and OPLS-AA/1.14*CM1A force fields. All the results are strengthened by comparison with the available experimental data on nuclear magnetic resonance J-coupling constants and ultrasonic spectra.