3020–3030 Nucleic Acids Research, 2016, Vol. 44, No. 7 Published online 14 March 2016 doi: 10.1093/nar/gkw160 Molecular dynamics simulations reveal the balance of forces governing the formation of a guanine tetrad––a common structural unit of G-quadruplex DNA Mateusz Kogut, Cyprian Kleist and Jacek Czub*
Department of Physical Chemistry, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
Received January 20, 2016; Revised February 26, 2016; Accepted March 1, 2016
ABSTRACT mation of G4 motifs at the single-stranded telomeric 3 - overhang leads to telomerase inhibition (10) and competes G-quadruplexes (G4) are nucleic acid conformations with the binding of the ssDNA overhang by POT1, a com- of guanine-rich sequences, in which guanines are ponent of shelterin protein complex protecting chromoso- arranged in the square-planar G-tetrads, stacked on mal ends (11–13). Accordingly, G4-stabilizing agents, such one another. G4 motifs form in vivo and are impli- as telomestatin, have been shown to impair telomere home- cated in regulation of such processes as gene ex- ostasis and thereby to induce growth arrest or apoptosis in pression and chromosome maintenance. The struc- cultured cells (14–16). As a result, G4 DNA is now consid- ture and stability of various G4 topologies were de- ered a promising target for anticancer therapy (17–19). G4- termined experimentally; however, the driving forces forming sequences are also found to be over-represented in for their formation are not fully understood at the other regulatory regions of the genome, including promot- molecular level. Here, we used all-atom molecular ers (20,21), introns (22), mitotic and meiotic double-strand break sites (6), immunoglobulin switch regions (9) and 5 - dynamics to probe the microscopic origin of the G4 UTRs (23). To fully understand the role of G4 in cell biol- motif stability. By computing the free energy profiles ogy and to allow the efficient use of G4-interactive ligands, governing the dissociation of the 3 -terminal G-tetrad it is necessary to understand structural diversity, conforma- in the telomeric parallel-stranded G4, we examined tional dynamics and molecular determinants of the stability the thermodynamic and kinetic stability of a single of these non-canonical DNA structural motifs. G-tetrad, as a common structural unit of G4 DNA. A common structural feature of G4 DNA is the pres- Our results indicate that the energetics of guanine ence of at least two stacked guanine tetrads (G-tetrads) with association alone does not explain the overall sta- four co-planar guanine bases, interconnected through cyclic bility of the G-tetrad and that interactions involving Hoogsteen-type hydrogen bonding and additionally stabi- sugar–phosphate backbone, in particular, the con- lized by coordination of monovalent cations to guanine O6 strained minimization of the phosphate–phosphate carbonyl groups (Figure 1A) (24). In intramolecular G4, the stacked G-tetrads are formed by (typically) contiguous gua- repulsion energy, are crucial in providing the ob- nine runs separated by intervening loops of variable length served enthalpic stabilization. This enthalpic gain and sequence. Based on the relative G-strand orientation, is largely compensated by the unfavorable entropy G4 motifs are classified into parallel, antiparallel and hy- change due to guanine association and optimization brid conformational classes, in which a variety of different of the backbone topology. folding topologies are possible (24). In particular, human telomeric G4 has been shown by x-ray crystallography and nuclear magnetic resonance spectroscopy to adopt six dis- INTRODUCTION tinct conformations containing three (or in one case, two) Certain guanine-rich DNA (and RNA) sequences have long G-tetrads and the three-nucleotide TTA loops, and differing been known to fold into intra- or inter-molecular four- in strand orientation, glycosidic conformation of guanines, stranded secondary structures called G-quadruplexes (G4) loop arrangement, etc. (25–30). (1,2). Recently, G4 have been also detected in ciliate and Calorimetric and thermal denaturation measurements mammalian cells (3–5), where they are believed to be in- combined with spectroscopic data show, on the one hand, volved in regulation of DNA replication and transcription, that biologically relevant G4 motifs are characterized by a genetic recombination, maintaining chromosome stability high thermodynamic stability (24,31,32), similar in magni- and other cellular processes (6–9). In particular, the for- tude to the stability of the competing Watson–Crick pairing
*To whom correspondence should be addressed. Tel: +48 583472092; Fax: +48 583472694; Email: [email protected]