A nanosecond molecular dynamics study of antiparallel d(G)₇ quadruplex structures: effect of the coordinated cations

Abstract

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G₇) quadruplex structures with different coordinated ions, namely Na⁺ and K⁺ ion, water and Na⁺ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of ˜1.5 Å from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na⁺ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K⁺ ions are quite immobile. MD studies at 400K indicate that K⁺ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K⁺ ion coordinated structure is more favourable than that for the Na⁺ ion coordinated antiparallel quadruplex structure.