Repetitive guanine-rich nucleic acid sequences play a crucial role in maintaining genome stability and the cell life cycle and represent potential targets for regulatory drugs. Recently, it has been demonstrated that guanine-based ligands with a porphyrin core can be used as markers of G-quadruplex assemblies in cell tissues. Herein, model systems of guanine-based ligands are explored by DFT methods. The energies of formation of modified guanine tetrads and those of modified tetrads stacked on the top of natural guanine tetrads have been calculated. The interaction energy has been decomposed into contributions from hydrogen bonding, stacking, and ion coordination and a twist-rise potential energy scan has been performed to find the individual local minima. Energy decomposition analysis reveals the impact of various substituents (F, Cl, Br, I, Me, NMe2 ) on individual energy terms. In addition, cooperative reinforcement in forming the modified and stacked tetrads, as well as the frontier orbitals participating in the hydrogen-bonding framework involving the HOMO-LUMO gap between the occupied σHOMO on the proton-accepting C=O and =N- groups and unoccupied σLUMO on the N-H groups, has been studied. The investigated systems are demonstrated to have a potential in ligand development, mainly due to stacking enhancement compared with natural guanine, which is used as a reference.

CEITEC Central European Institute of Technology Masaryk University Kamenice 5A4 62500 Brno Czech Republic

Department of Theoretical Chemistry Amsterdam Center for Multiscale Modeling Vrije Universiteit Amsterdam De Boelelaan 1083 1081 HV Amsterdam The Netherlands

National Center for Biomolecular Research Faculty of Science Masaryk University Kamenice 5A4 62500 Brno Czech Republic

References provided by Crossref.org