RESUMO
G-Quadruplex (GQ) is a secondary structural unit of DNA, formed at the telomere region of the chromosome with a high guanine content. It is reported that the GQs can hinder many biological processes. Thus, research thrives to explore the structural stability of GQs. Studies based on circular dichroism (CD) and nuclear magnetic resonance (NMR) experiments established the vital role of cations such as K+ and Mg2+ in the stability of antiparallel G-quadruplexes (AGQs). However, there is a need to understand how stability in AGQ is attained in the presence of cations. Here, we employed molecular dynamics (MD) simulations, steered MD (SMD) simulations, and QM/MM calculations to understand the biophysical and electronic bases of the stability imparted to AGQs via cation binding. Our results showed that Mg2+ prefers to bind in the plane with the guanine tetrad, whereas K+ binds in between the AGQ tetrads. Thus, three Mg2+ cations or two K+ ions are needed to stabilize an AGQ molecule, where each and every tetrad binds to Mg2+ more robustly with a higher binding affinity. SMD revealed that the traversal of K+ through the AGQ central channel required less force than that of Mg2+, illustrating the presence of more strong interactions between Mg2+ and AGQ tetrads compared to K+. The stabilization in the AGQ tetrads due to cation binding is reassessed by employing ab initio simulations. Mixed QM/MM calculations confirmed that Mg2+ binds strongly to AGQ compared to K+, and it induces higher interactions between the guanine tetrads. However, K+ binding to AGQ induces a higher stabilization energy than Mg2+ binding to AGQ tetrads. Despite the higher binding energy, Mg2+ binding imparts lower stabilization to AGQ due to its unfavorable fermionic quantum energy.
