NMR structure refinement and dynamics of the K+-[d(G3T4G3)]2 quadruplex via particle mesh Ewald molecular dynamics simulations

The solution structure and dynamical properties of the potassium-stabilized, hairpin dimer quadruplex formed by the oligonucleotide d(G(3)T(4)G(3)) have been elucidated by a combination of high-resolution NMR and molecular dynamics simulations. Refinement calculations were carried out both in vacuo, without internally coordinated K+ cations, and in explicit water, with internally coordinated K+ cations, In the latter case, the electrostatic interactions were calculated using the particle mesh Ewald (PME) method. The NMR restraints indicate that the K+ quadruplex has a folding arrangement similar to that formed by the same oligonucleotide in the presence of sodium, but with significant local differences. Unlike the Nat quadruplex, the thymine loops found in K+ exhibit considerable flexibility, and appear to interconvert between two preferred conformations. Furthermore, the NMR evidence points toward K+-stabilized guanine quartets of slightly larger diameter relative to the Na+-stabilized structure. The characteristics of the quartet stem are greatly affected by the modeling technique employed: caged cations alter the size and symmetry of the quartets, and explicit water molecules form hydration spines within the grooves. These results provide insight into those factors that determine the overall stability of hairpin dimer quadruplexes and the effects of different cations in modulating the relative stability of the dimeric hairpin and linear, four-stranded, quadruplex forms.