Published online 30 January 2006
Article |
Similarities and differences in interaction of K+ and Na+ with condensed ordered DNA. A molecular dynamics computer simulation study
School of Biological Sciences, Nanyang Technological University 60 Nanyang Drive, Singapore 637551
*To whom correspondence should be addressed. Tel: +65 6316 2856; Fax: +65 6791 3856; Email: LarsNor{at}ntu.edu.sg
Received November 12, 2005. Revised December 21, 2005. Accepted December 21, 2005.
Four 20 ns molecular dynamics simulations have been performed with two counterions, K+ or Na+, at two water contents, 15 or 20 H2O per nucleotide. A hexagonal simulation cell comprised of three identical DNA decamers [d(5'-ATGCAGTCAG) x d(5'-TGACTGCATC)] with periodic boundary condition along the DNA helix was used. The simulation setup mimics the DNA state in oriented DNA fibers or in crystals of DNA oligomers. Variation of counterion nature and water content do not alter averaged DNA structure. K+ and Na+ binding to DNA are different. K+ binds to the electronegative sites of DNA bases in the major and the minor grooves, while Na+ interacts preferentially with the phosphate groups. Increase of water causes a shift of both K+ and Na+ from the first hydration shell of O1P/O2P and of the DNA bases in the minor groove with lesser influence for the cation binding to the bases in the major groove. Mobility of both water and cations in the KDNA systems is faster than in the NaDNA systems: Na+ organizes and immobilizes water structure around itself and near DNA while for K+ water is less organized and more dynamic.