Toward the design of new DNA G-quadruplex ligands through rational analysis of polymorphism and binding data

Human telomeres play a key role in protecting chromosomal ends from fusion events; they are composed of d(TTAGGG) repeats, ranging in size from 3 to 15 kb. They form G-quadruplex DNA structures, stabilized by G-quartets in the presence of cations, and are involved in several biological processes. In particular, a telomere maintenance mechanism is provided by a specialized enzyme called telomerase, a reverse transcriptase able to add multiple copies of the 5'-GGTTAG-3' motif to the end of the G-strand of the telomere and which is over-expressed in the majority of cancer cells. The central cation has a crucial role in maintaining the stability of the structure. Based on its nature, it can be associated with different topological telomeric quadruplexes, which depend also on the orientation of the DNA strands and the syn/anti conformation of the guanines. Such a polymorphism, confirmed by the different structures deposited in the Protein Data Bank (PDB), prompted us to apply a computational protocol in order to investigate the conformational properties of a set of known G-quadruplex ligands and their molecular recognition against six different experimental models of the human telomeric sequence d[AG3(T2AG3)3]. The average AutoDock correlation between theoretical and experimental data yielded an r2 value equal to 0.882 among all the studied models. Such a result was always improved with respect to those of the single folds, with the exception of the parallel structure (r2 equal to 0.886), thus suggesting a key role of this G4 conformation in the stacking interaction network. Among the studied binders, a trisubstituted acridine and a dibenzophenanthroline derivative were well recognized by the parallel and the mixed G-quadruplex structures, allowing the identification of specific key contacts with DNA and the further design of more potent or target specific G-quadruplex ligands.