TY - JOUR
T1 - Anharmonic vibrational self-consistent field calculations as an approach to improving force fields for monosaccharides
AU - Gregurick, S. K.
AU - Liu, Jennifer H.Y.
AU - Brant, David A.
AU - Gerber, R. B.
PY - 1999/4/29
Y1 - 1999/4/29
N2 - 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.
AB - 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.
UR - http://www.scopus.com/inward/record.url?scp=0001109647&partnerID=8YFLogxK
U2 - 10.1021/jp9826221
DO - 10.1021/jp9826221
M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???
AN - SCOPUS:0001109647
SN - 1520-6106
VL - 103
SP - 3476
EP - 3488
JO - Journal of Physical Chemistry B
JF - Journal of Physical Chemistry B
IS - 17
ER -