On the Truncation of Long-Range Electrostatic Interactions in DNA

Abstract

Long-range interactions are known to play an important role in highly polar biomolecules like DNA. In molecular dynamics simulations of nucleic acids and proteins, an accurate treatment of the long-range interactions are crucial for achieving stable nanosecond trajectories. In this report, we evaluate the structural and dynamic effects on a highly charged oligonucleotide in aqueous solution from different long-range truncation methods. Two group-based truncation methods, one with a switching function and one with a force-switching function were found to fail to give accurate stable trajectories close to the crystal structure. For these group-based truncation methods, large root mean square (rms) deviations from the initial structure were obtained and severe distortions of the oligonucleotide were observed. Another group-based truncation scheme, which used an abrupt truncation at 8. 0 A or at 12.0 A was also investigated. For the short cutoff distance, the conformations deviated far away from the initial structure and were significantly distorted. However, for the longer cutoff, where all necessary electrostatic interactions were included, the trajectory was quite stable. For the particle mesh Ewald (PME) truncation method, a stable DNA simulation with a heavy atom rms deviation of 1.5 A was obtained. The atom-based truncation methods also resulted in stable trajectories, according to the rms deviation from the initial B-DNA structure, of between 1.5 and 1.7 A for the heavy atoms. In these stable simulations, the heavy atom rms deviations were approximately 0.6-1.0 A lower for the bases than for the backbone. An increase of the cutoff radius from 8 to 12 A decreased the rms deviation by approximately 0.2 A for the atom-based truncation method with a force-shifting function, but increased the computational time by a factor of 2. Increasing the cutoff from 12 to 18 A for the atom-based truncation method with a force-shifting function requires 2-3 times more computational time, but did not significantly change the rms deviation. Similar rms deviations from the initial structure were found for the atom-based method with a force-shifting function and for the PME method. The computational cost was longer for the PME method with a cutoff of 12. 0 A for the direct space nonbonded calculations than for the atom-based truncation method with a force-shifting function and a cutoff of 12.0 A. If a nonperiodic boundary, e.g., a spherical boundary, was used, a considerable speedup could be achieved. From the rms fluctuations, the terminal nucleotides and especially the cytidines were found to be more flexible than the nonterminal nucleotides. The B-DNA form of the oligonucleotide was maintained throughout the simulations and is judged to depend on the parameters of the energy function and not on the truncation method used to handle the long-range electrostatic interactions. To perform accurate and stable simulations of highly charged biological macromolecules, we recommend that the atom-based force-shift method or the PME method...

Ultrafast Dynamics in DNA: “Fraying” at the End of the Helix

“…Previous NMR measurements have detected an increased amplitude of motion at the ends of oligonucleotides in the picosecond time range, but could not isolate a time constant for this motion. [4][5][6] X-ray, 3 NMR 4-6 and computer [7][8][9] experiments agree that fraying is largely confined to the last two base-pairs, with most of the motion on the terminal base pair. Two general types of motion can be envisioned, each of which could be detected in a TRSS experiment.…”

Orchestration of cooperative events in DNA synthesis and repair mechanism unraveled by transition path sampling of DNA polymerase β's closing

“…The truncation of long-range electrostatic interactions makes the computation of our 40,238-atom system tractable; studies on DNA (e.g., ref. 55) conclude that this approach is satisfactory.…”

Empirical force fields for biological macromolecules: Overview and issues

“…[222][223][224] These successes were initially attributed to the use of the Ewald method, 4 typically in the particle-mesh Ewald (PME) formalism, 129 for the treatment of long-range electrostatic interactions; however, it was subsequently shown that use of adequate lengths for truncation of electrostatic interactions combined with the appropriate smoothing functions 123 led to stable MD simulations of DNA. [225][226][227] In addition, improvements in second-generation force fields made significant contributions to the ability to perform stable simulations of oligonucleotides. Second-generation forces fields for nucleic acids included the Cornell et al AMBER (PARM94) 17 and CHARMM all-atom 163 force fields.…”