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.