The alkali metal ion affinity of guanine quadruplexes has been studied using dispersion-corrected density functional theory (DFT-D). We have done computational investigations 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. In both cases, an alkali metal cation is needed to assist self-assembly. Our quantum chemical computations on these supramolecular systems are able to reproduce the experimental order of affinity of the guanine quadruplexes for the cations Li(+), Na(+), K(+), Rb(+), and Cs(+). The strongest binding is computed between the potassium cation and the quadruplex as it occurs in nature. The desolvation and the size of alkali metal cations are thought to be responsible for the order of affinity. Until now, the relative importance of these two factors has remained unclear and debated. By assessing the quantum chemical 'size' of the cation, determining the amount of deformation of the quadruplex needed to accommodate the cation and through the energy decomposition analysis (EDA) of the interaction energy between the cation and the guanines, we reveal that the desolvation and size of the alkali metal cation are both almost equally responsible for the order of affinity.

We show that the cooperative reinforcement between hydrogen bonds in guanine quartets is not caused by resonance-assisted hydrogen bonding (RAHB). This follows from extensive computational analyses of guanine quartets (G(4)) and xanthine quartets (X(4)) based on dispersion-corrected density functional theory (DFT-D). Our investigations cover the situation of quartets in the gas phase, in aqueous solution as well as in telomere-like stacks. A new mechanism for cooperativity between hydrogen bonds in guanine quartets emerges from our quantitative Kohn-Sham molecular orbital (MO) and corresponding energy decomposition analyses (EDA). Our analyses reveal that the intriguing cooperativity originates from the charge separation that goes with donor-acceptor orbital interactions in the σ-electron system, and not from the strengthening caused by resonance in the π-electron system. The cooperativity mechanism proposed here is argued to apply, beyond the present model systems, also to other hydrogen bonds that show cooperativity effects.

Our computational studies suggest that 3-substituted xanthines are good candidates for tetrad and quadruplex structures. 3-Methylxanthine (3MX) has been synthesized from 7-benzylxanthine, and the existence of tetrameric and octameric aggregates of 3MX with NH 4 + , Na + and K + ions in the gas phase (MS) and in DMSO-d 6 solution (NMR) has been observed. The ''internal'' H-bonds (N1HÁ Á ÁO6) are stronger than the ''external'' ones (N7HÁ Á ÁO2) in these clusters (NMR). Fig. 1 Tetrads composed of 7H-3-methylxanthine (3MX) and different cations [cat + = none, Na + , K + , NH 4 + ; purine atom and spectroscopic (underlined) numbering shown in upper left ring].