G-quadruplexes (GQs) have become valid targets for anticancer studies in recent decades due to their multifaceted biological function. Herewith, we aim to quantify interactions of potential heterocyclic ligands (Ls) with model GQs. For seven 4-aminoquinazolines and three 2-heteroaryl perimidines, seven of this ten-membered group so far unknown, we use routine quantum chemical modeling. As shown in the literature, a preferred mode of interaction of heterocycles with cellular structures is stacking to exposable faces of G-quadruplexes. To exploit the energy of this interaction as a molecular descriptor and achieve the necessary chemical precision, we use state of the art large-scale density functional theory (DFT) calculations of stacked heterocycles to a GQ. Actually, the GQ has been simplified for the computation by stripping it off all pentose phosphate residues into a naked model of stacked guanine quartets. The described model thus becomes computable. The obtained heterocyclic ligand GQ.L stacking energies, that is, their GQ affinities, are the necessary ligand descriptors. Using the ligand biological inhibitory activities (IC₅₀) on a human malignant melanoma A375 cell line, we obtain a good linear relationship between computed ligand stacking affinities to GQ, and experimental log (IC₅₀) values. Based on the latter relationship, we discuss a putative mechanism of anticancer activity of heterocyclic ligands via stacking interactions with GQs and thereby controlling cell regulatory activity. This mechanism may tentatively be applied to other condensed five- and six-membered small heterocycles as well.

Quadruplexes (GQs), peculiar DNA/RNA motifs concentrated in specific genomic regions, play a vital role in biological processes including telomere stability and, hence, represent promising targets for anticancer therapy. GQs are formed by folding guanine-rich sequences into square planar G-tetrads which stack onto one another. Metal cations, most often potassium, further stabilize the architecture by coordinating the lone electron pairs of the O atoms. The presence of additional nucleic acid bases, however, has been recently observed experimentally and contributes substantially to the structural heterogeneity of quadruplexes. Therefore, it is of paramount significance to understand the factors governing the underlying complex processes in these structures. The current study employs DFT calculations to model the interactions between metal cations (K⁺, Na⁺, Sr²⁺) and diverse tetrads composed of a guanine layer in combination with a guanine (G)-, adenine (A)-, cytosine (C)-, thymine (T)-, or uracil (U)-based tetrad layer. Moreover, the addition of 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline to the modeled quadruplexes as a possible mechanism of its well-exerted antitumor effect is assessed. The calculations imply that the metal cation competition and ligand complexation are influenced by the balance between electronic and implicit/explicit solvation effects, the composition of the tetrad layers, as well as by the solvent exposure to the surrounding environment expressed in terms of different dielectric constant values. The provided results significantly enhance our understanding of quadruplex diversity, ligand recognition, and the underlying mechanisms of stabilization at an atomic level.

Small bioactive molecules show significant propensity to form noncovalent addition complexes with guanine quadruplexes, G4. The stabilization energies of these complexes have been computed precisely at the sufficiently high 6-31G** basis set level of density functional quantum chemical theory, DFT. A decisive factor in present model computations is the adopted size of G4 models, whether these consist simply of stacked quanine quartets, or also involve (deoxy)ribose-phosphate fragments of proper nucleic acids. The challenge is in the preservation of physico-chemical accuracy of DFT computations with increasing sizes of models, involving upwards of 120 atoms for the simplest two-layer G4, plus at least 60 pentose-phosphate linker atoms per each pair of guanine quartets. Bioactive ligand sizes add to the requirements for further rigorous analyses of the roles of G4 complexes in biological processes, which thus remain necessarily open-ended.