Anharmonic vibrational self-consistent field calculations as an approach to improving force fields for monosaccharides

S. K. Gregurick, Jennifer H.Y. Liu, David A. Brant, R. B. Gerber*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

34 Scopus citations

Abstract

An Amber-type of force field, based on experimental vibrational frequencies which is suitable for monosaccharides, is presented. In the present force field, the atomic partial charges and some torsional parameters are derived from a fit of calculated vibrational energy levels to known experimental spectra for α-D-glucose. The vibrational spectra were calculated using the vibrational self-consistent field (VSCF) method, which includes contributions from both anharmonic and mode-coupled terms. We find that with a reparametrization of the force field the agreement between the experimental and calculated vibrational spectra is ±3.3 cm-1 for α-glucose and ±5.1 cm-1 for β-glucose. Using the VSCF method, we are also able to lend support to the idea that the COH bending motion is strongly coupled to the methylene and methine modes, as well as the other internal modes. We then test our spectroscopically derived force field by calculating the anomeric effect for α → β glucose. Molecular dynamics simulations are performed separately for both anomers (α and β) in order to evaluate configurational entropy, and hence free energy. We find that out of six simulations, half correctly predict the anomeric free energy ΔGα→β = -0.3 kcal/mol, while two simulations yield a ΔGα→β = +0.2 kcal/mol, and in one simulation ΔGα→β ∼ 0.0 kcal/mol. We also calculate the atomic pair distribution function, g(r), and show that in most simulations, the β-conformer is slightly more adept at structuring the surrounding water molecules, resulting in better hydrogen bonding for this anomer. However, we believe that our force field, which is static, is unable to represent adequately the dynamic interactions between the pyranose sugar and the surrounding water molecules. This resulted in a large fluctuation about the average calculated anomeric free energy of ΔGα→β ∼ -0.1 kcal/mol. Thus, while the current all atom force field is well suited for spectroscopic studies of monosaccharides, it is not yet well suited for dynamical studies.

Original languageEnglish
Pages (from-to)3476-3488
Number of pages13
JournalJournal of Physical Chemistry B
Volume103
Issue number17
DOIs
StatePublished - 29 Apr 1999

Fingerprint

Dive into the research topics of 'Anharmonic vibrational self-consistent field calculations as an approach to improving force fields for monosaccharides'. Together they form a unique fingerprint.

Cite this