Restrained molecular dynamics of solvated duplex DNA using the particle mesh Ewald method

Restrained and unrestrained aqueous solution molecular dynamics simulations applying the particle mesh Ewald (PME) method to DNA duplex structures previously determined via in vacuo restrained molecular dynamics with NMR-derived restraints are reported. Without experimental restraints, the DNA decamer, d(CATTTGCATC)⋅d(GATGCAAATG) and trisdecamer, d(AGCTTGCCTTGAG)⋅d(CTCAAGGCAAGCT), structures are stable on the nanosecond time scale and adopt conformations in the B-DNA family. These free DNA simulations exhibit behavior characteristic of PME simulations previously performed on DNA sequences, including a low helical twist, frequent sugar pucker transitions, BI- BII(ε−ζ) transitions and coupled crankshaft (α−γ) motion. Refinement protocols similar to the original in vacuo restrained molecular dynamics (RMD) refinements but in aqueous solution using the Cornell et al. force field [Cornell et al. (1995) J. Am. Chem. Soc., 117, 5179–5197] and a particle mesh Ewald treatment produce structures which fit the restraints very well and are very similar to the original in vacuo NMR structure, except for a significant difference in the average helical twist. Figures of merit for the average structure found in the RMD PME decamer simulations in solution are equivalent to the original in vacuo NMR structure while the figures of merit for the free MD simulations are significantly higher. The free MD simulations with the PME method, however, lead to some sequence-dependent structural features in common with the NMR structures, unlike free MD calculations with earlier force fields and protocols. There is some suggestion that the improved handling of electrostatics by PME improves long-range structural aspects which are not well defined by the short-range nature of NMR restraints.

[1]  G L Kenyon,et al.  Solution structure of a conserved DNA sequence from the HIV-1 genome: restrained molecular dynamics simulation with distance and torsion angle restraints derived from two-dimensional NMR spectra. , 1994, Biochemistry.

[2]  V. Zhurkin,et al.  Systematic study of nuclear Overhauser effects vis-à-vis local helical parameters, sugar puckers, and glycosidic torsions in B DNA: insensitivity of NOE to local transitions in B DNA oligonucleotides due to internal structural compensations. , 1992, Biochemistry.

[3]  P. Kollman,et al.  COMBINED LOCALLY ENHANCED SAMPLING AND PARTICLE MESH EWALD AS A STRATEGY TO LOCATE THE EXPERIMENTAL STRUCTURE OF A NONHELICAL NUCLEIC ACID , 1998 .

[4]  R. Shafer,et al.  Solution structure of the octamer motif in immunoglobulin genes via restrained molecular dynamics calculations. , 1994, Biochemistry.

[5]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[6]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[7]  D. Beveridge,et al.  A 5-nanosecond molecular dynamics trajectory for B-DNA: analysis of structure, motions, and solvation. , 1997, Biophysical journal.

[8]  T. Darden,et al.  The effect of long‐range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods , 1993 .

[9]  A. Klug,et al.  Helical periodicity of DNA determined by enzyme digestion , 1980, Nature.

[10]  R. Sarma,et al.  Structure, Motion, Interaction and Expression of Biological Macromolecules. , 1998, Journal of biomolecular structure & dynamics.

[11]  E. O'Brien,et al.  The position of tropomyosin in muscle thin filaments , 1980, Nature.

[12]  Ashok Kumar,et al.  Probability assessment of conformational ensembles: sugar repuckering in a DNA duplex in solution. , 1995, Biophysical journal.

[13]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[14]  R. Shafer,et al.  The octamer motif in immunoglobulin genes: extraction of structural constraints from two-dimensional NMR studies. , 1992, Biochemistry.

[15]  P. Kollman,et al.  Molecular Dynamics Simulations on Solvated Biomolecular Systems: The Particle Mesh Ewald Method Leads to Stable Trajectories of DNA, RNA, and Proteins , 1995 .

[16]  T. James,et al.  How to generate accurate solution structures of double-helical nucleic acid fragments using nuclear magnetic resonance and restrained molecular dynamics. , 1995, Methods in enzymology.

[17]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[18]  R. Elber,et al.  Modeling side chains in peptides and proteins: Application of the locally enhanced sampling and the simulated annealing methods to find minimum energy conformations , 1991 .

[19]  P A Kollman,et al.  Observation of the A-DNA to B-DNA transition during unrestrained molecular dynamics in aqueous solution. , 1996, Journal of molecular biology.