Stability and migration of metal ions in G4-wires by molecular dynamics simulations

We present a molecular dynamics investigation of guanine quadruple helices based on classical force fields. We analyze the dependence of the helical conformation on various compositional factors, such as the length of the G4-wire, as well as the incorporation into the helix channel of alkali ions of different species and in different amounts. In compliance with previous indications, our results suggest that monovalent alkali cations assist the stability of the quadruplex arrangement against disruption on the few nanoseconds time scale in the order of increasing van der Waals radius. Whereas very short G4-wire fragments immediately unfold in the absence of coordinating metal ions or in the presence of tiny ions (e.g., Li+) in agreement with the experimental evidence that empty short guanine quadruplexes are not formed in any synthetic conditions, our simulations show that longer empty helices do not discompose. This finding supports the possibility of producing long G4-wires with different guanine-cation stoichiometries than those routinely known. The classical trajectories allow us to identify different stationary axial sites for the different metal species, which are confirmed by complementary quantum calculations.