Molecular Dynamics Simulations of the d(TAT) Triple Helix

Molecular dynamics (MD) simulations have been used to study the dynamical and time-averaged characteristics of the DNA triple helix d(T)10âd(A)10âd(T)10. The structures sampled during the trajectory resemble closely the B-type model for the DNA triplex proposed on the basis of NMR data, although there are some subtle differences. Alternative P- and A-type conformations for the triplex, suggested from X-ray experiments, are not predicted to contribute significantly to the structure of the DNA triplex in solution. Comparison with the best available experimental data supports the correctnes of the MD-generated structures. The analysis of the collected data gives a detailed picture of the characteristics of triple-helix DNA. A new and interesting pattern of hydration, specific for triplex DNA, is an important observation. The results suggest that molecular dynamics can be useful for the study of novel nucleic acid structures.

Ramachandran-type Plots for Glycosidic Linkages: Examples from Molecular Dynamics Simulations using the Glycam06 Force Field

The goals of this article are to (1) provide further validation of the Glycam06 force field, specifically for its use in implicit solvent molecular dynamic (MD) simulations, and (2) to present the extension of G.N. Ramachandran's idea of plotting amino acid phi and psi angles to the glycosidic phi, psi, and omega angles formed between carbohydrates. As in traditional Ramachandran plots, these carbohydrate Ramachandran-type (carb-Rama) plots reveal the coupling between the glycosidic angles by displaying the allowed and disallowed conformational space. Considering two-bond glycosidic linkages, there are 18 possible conformational regions that can be defined by (α, ϕ, ψ) and (β, ϕ, ψ), whereas for three-bond linkages, there are 54 possible regions that can be defined by (α, ϕ, ψ, ω) and (β, ϕ, ψ, ω). Illustrating these ideas are molecular dynamic simulations on an implicitly hydrated oligosaccharide (700 ns) and its eight constituent disaccharides (50 ns/disaccharide). For each linkage, we compare and contrast the oligosaccharide and respective disaccharide carb-Rama plots, validate the simulations and the Glycam06 force field through comparison to experimental data, and discuss the general trends observed in the plots.