Ion motions in molecular dynamics simulations on DNA

Area of Science:

• Biophysics
• Computational Chemistry
• Molecular Biology

Background:

• Counterions are essential for DNA structure and function.
• Molecular dynamics (MD) simulations provide detailed insights into ion dynamics at the molecular level.
• Convergence of MD simulations for mobile counterions on DNA is challenging.

Purpose of the Study:

• To determine the necessary sampling for accurate MD simulations of DNA-protein interactions.
• To assess the convergence of DNA-Na+ interactions and DNA structural parameters.
• To investigate sequence-dependent ion occupancies and their relation to DNA minor groove width.

Main Methods:

• Performed MD simulations of a d(CGCGAATTCGCG) DNA duplex with 22 Na+ counterions in aqueous solution.
• Calculated first shell ion occupancies and DNA-Na+ radial distribution functions over time.
• Compared ion dynamics convergence with relaxation times of DNA structural parameters (shift, slide, rise, tilt, roll, twist).

Main Results:

• Achieved convergence for DNA-Na+ interactions and DNA structural parameters.
• Observed sequence-dependent fractional occupancies of Na+ ions in DNA major and minor grooves.
• Investigated potential correlations between ion proximity and DNA minor groove width.

Conclusions:

• Established convergence criteria for MD simulations of DNA-counterion systems.
• Demonstrated sequence-specific ion binding patterns influencing DNA conformation.
• Highlighted the importance of ion dynamics in understanding DNA structure-function relationships.