Force Field Dependence of NMR-Based, Restrained Molecular Dynamics DNA Structure Calculations Including an Analysis of the Influence of Residual Dipolar Coupling Restraints

Abstract

Restrained molecular dynamics is widely used to calculate DNA structures from NMR data. Here, results of an in silico experiment show that the force field can be significant compared to the NMR restraints in driving the final structures to converge. Specifically, we observed that i) the influence of the force field leads to artificially tight convergence within final families of structures and ii) the precision and character of resulting structures depend on the choice of force field used in the calculations. A canonical B-DNA model was used as a target structure. Distances, dihedral angles, and simulated residual dipolar couplings were measured in the target structure and used as restraints. X-PLOR and Discover, which use force fields developed for CHARMM and AMBER programs, respectively, were tested and found to produce different final structures despite the use of identical distance and dihedral restraints. Incorporation of residual dipolar coupling restraints in X-PLOR improves convergence with the target structure and between families of structures indicating that the force field dependence can potentially be overcome if residual dipolar coupling restraints are employed.