To gain better understanding of the stabilizing interactions between metal ions and DNA quadruplexes, dispersion-corrected density functional theory (DFT-D) based calculations were performed on double-, triple- and four-layer guanine tetrads interacting with alkali metal cations. All computations were performed in aqueous solution that mimics artificial supramolecular conditions where guanine bases assemble into stacked quartets as well as biological environments in which telomeric quadruplexes are formed. To facilitate the computations on these significant larger systems, optimization of the DFT description was performed first by evaluating the performance of partial reduced basis sets. Analysis of the stabilizing interactions between alkali cations and the DNA bases in double and triple-layer guanine quadruplex DNA reproduced the experimental affinity trend of the order Li+< Rb+ < Na+ < K+. The desolvation and the size of alkali metal cations are thought to be responsible for the order of affinity. Nevertheless, for the alkali metal cation species individually, the magnitude of the bond energy stays equal for binding as first, second or third cation in double, triple and four-layer guanine quadruplexes, respectively. This is the result of an interplay between a decreasingly stabilizing interaction energy and increasingly stabilizing solvation effects, along the consecutive binding events. This diminished interaction energy is the result of destabilizing electrostatic repulsion between the hosted alkali metal cations. This work emphasizes the stabilizing effect of aqueous solvent on large highly charged biomolecules.