Abstract
Molecular dynamics simulations were carried out on two conformations of the dinucleoside monophosphate guanylyl-3',5'-uridine (GpU) in aqueous solution with one sodium counterion. One stacked conformation and one with the C3'-O3'-P-O5' backbone torsion angle twisted 180 degrees to create an unstacked conformation. We observed a relatively stable behavior of the stacked conformation, which remained stacked throughout the simulation, whereas the unstacked conformation showed major changes in the backbone torsion and glycosidic angles. During the simulation the unstacked conformation transformed into a more stacked form and then back again to an unstacked one. The calculated correlation times for rotational diffusion from the molecular dynamics simulations are in agreement with fluorescence anisotropy and nuclear magnetic resonance data. As expected, the correlation times for rotational diffusion of the unstacked conformation were observed to be longer than for the stacked conformation. The 2'OH group may contribute in stabilizing the stacked conformation, where the O2'-H...O4' hydrogen bond occurred in 82.7% of the simulation.