Multipolar Description of Atom-Atom Electrostatic Interaction Energies in Single/Double-Stranded DNAs

Area of Science:

• Computational chemistry
• Molecular dynamics
• Biophysics

Background:

• Current DNA force fields often use simplified electrostatic models.
• Atomic charge and low-rank multipole moments are typically used, neglecting electron distribution complexities.
• Accurate electrostatics are crucial for understanding DNA behavior.

Purpose of the Study:

• To develop a more accurate molecular force field for DNA simulations.
• To incorporate high-rank atomic multipole moments and account for electron anisotropy.
• To assess the transferability of the developed models for both single-stranded and double-stranded DNA.

Main Methods:

• Applied a multipole expansion method combining polarizability and anisotropy.
• Used deoxynucleotides from pilot DNA systems, capped to mimic their environment.
• Integrated atomic multipole moments instead of fixed-point charges for electrostatic energy calculations.

Main Results:

• The new method significantly improves the accuracy of electrostatic interactions in DNA simulations.
• Models developed for pilot systems demonstrated applicability and transferability to test systems.
• The approach successfully modeled behaviors of both single-stranded and double-stranded DNA.

Conclusions:

• The developed multipole expansion method offers enhanced accuracy for DNA simulations.
• The transferability of these models is a significant advancement for future DNA force field development.
• This work paves the way for more precise computational studies of DNA properties.