Loss of loop adenines alters human telomere d[AG3(TTAG3)3] quadruplex folding

Abasic (AP) lesions are the most frequent type of damages occurring in cellular DNA. Here we describe the conformational effects of AP sites substituted for 2′-deoxyadenosine in the first (ap7), second (ap13) or third (ap19) loop of the quadruplex formed in K⁺ by the human telomere DNA 5′-d[AG₃(TTAG₃)₃]. CD spectra and electrophoresis reveal that the presence of AP sites does not hinder the formation of intramolecular quadruplexes. NMR spectra show that the structural heterogeneity is substantially reduced in ap7 and ap19 as compared to that in the wild-type. These two (ap7 and ap19) sequences are shown to adopt the hybrid-1 and hybrid-2 quadruplex topology, respectively, with AP site located in a propeller-like loop. All three studied sequences transform easily into parallel quadruplex in dehydrating ethanol solution. Thus, the AP site in any loop region facilitates the formation of the propeller loop. Substitution of all adenines by AP sites stabilizes the parallel quadruplex even in the absence of ethanol. Whereas guanines are the major determinants of quadruplex stability, the presence or absence of loop adenines substantially influences quadruplex folding. The naturally occurring adenine-lacking sites in the human telomere DNA can change the quadruplex topology in vivo with potentially vital biological consequences.     

Kinetics of tetramolecular quadruplexes

The melting of tetramolecular DNA or RNA quadruplexes is kinetically irreversible. However, rather than being a hindrance, this kinetic inertia allows us to study association and dissociation processes independently. From a kinetic point of view, the association reaction is fourth order in monomer and the dissociation first order in quadruplex. The association rate constant k_on, expressed in M⁻³·s⁻¹ decreases with increasing temperature, reflecting a negative activation energy (E_on) for the sequences presented here. Association is favored by an increase in monocation concentration. The first-order dissociation process is temperature dependent, with a very positive activation energy E_off, but nearly ionic strength independent. General rules may be drawn up for various DNA and RNA sequence motifs, involving 3–6 consecutive guanines and 0–5 protruding bases. RNA quadruplexes are more stable than their DNA counterparts as a result of both faster association and slower dissociation. In most cases, no dissociation is found for G-tracts of 5 guanines or more in sodium, 4 guanines or more in potassium. The data collected here allow us to predict the amount of time required for 50% (or 90%) quadruplex formation as a function of strand sequence and concentration, temperature and ionic strength.     

How long is too long? Effects of loop size on G-quadruplex stability

We compared here 80 different sequences containing four tracts of three guanines with loops of variable length (between 1 and 15 bases for unmodified sequences, up to 30 for fluorescently labeled oligonucleotides). All sequences were capable of forming stable quadruplexes, with T_m above physiological temperature in most cases. Unsurprisingly, the melting temperature was systematically lower in sodium than in potassium but the difference between both ionic conditions varied between 1 and >39°C (average difference: 18.3°C). Depending on the sequence context, and especially for G4 sequences involving two very short loops, the third one may be very long without compromising the stability of the quadruplex. A strong inverse correlation between total loop length and T_m was found in K⁺: each added base leads to a 2°C drop in T_m or ∼0.3 kcal/mol loss in ΔG°. The trend was less clear in Na⁺, with a longer than expected optimal loop length (up to 5 nt). This study will therefore extend the sequence repertoire of quadruplex-prone sequences, arguing for a modification of the widely used consensus (maximal loop size of 7 bases).     

Stability of intramolecular quadruplexes: sequence effects in the central loop

Hundreds of thousands of putative quadruplex sequences have been found in the human genome. It is important to understand the rules that govern the stability of these intramolecular structures. In this report, we analysed sequence effects in a 3-base-long central loop, keeping the rest of the quadruplex unchanged. A first series of 36 different sequences were compared; they correspond to the general formula GGGTTTGGGHNHGGGTTTGGG. One clear rule emerged from the comparison of all sequence motifs: the presence of an adenine at the first position of the loop was significantly detrimental to stability. In contrast, adenines have no detrimental effect when present at the second or third position of the loop. Cytosines may either have a stabilizing or destabilizing effect depending on their position. In general, the correlation between the T_m or ΔG° in sodium and potassium was weak. To determine if these sequence effects could be generalized to different quadruplexes, specific loops were tested in different sequence contexts. Analysis of 26 extra sequences confirmed the general destabilizing effect of adenine as the first base of the loop(s). Finally, analysis of some of the sequences by microcalorimetry (DSC) confirmed the differences found between the sequence motifs.     

Influence of loop size on the stability of intramolecular DNA quadruplexes

We have determined the stability of intramolecular DNA quadruplexes in which the four G₃-tracts are connected by non-nucleosidic linkers containing propanediol, octanediol or hexaethylene glycol, replacing the TTA loops in the human telomeric repeat sequence. We find that these sequences all fold to form intramolecular complexes, which are stabilized by lithium < sodium < potassium. Quadruplex stability increases in the order propanediol < hexaethylene glycol < octanediol. The shallower shape of the melting profile with propanediol linkers and its lower dependency on potassium concentration suggests that this complex contains fewer stacks of G-quartets. The sequence with octanediol linkers displays a biphasic melting profile, suggesting that it can adopt more than one stable structure. All these complexes display melting temperatures above 310 K in the presence of 10 mM lithium, without added potassium, in contrast to the telomeric repeat sequence. These complexes also fold much faster than the telomeric repeat and there is little or no hysteresis between their melting and annealing profiles. In contrast, the human telomeric repeat sequence and a complex containing two hexaethylene glycol groups in each loop, are less stable and fold more slowly. The melting and annealing profiles for the latter sequence show significant differences, even when heated at 0.2°C min^–1. CD spectra for the oligonucleotides containing non-nucleosidic linkers show positive maxima at 264 nm, with negative minima ~244 nm, which are characteristic of parallel quadruplex structures. These results show that the structure and stability of intramolecular quadruplexes is profoundly influenced by the length and composition of the loops.     

Improved thrombin binding aptamer by incorporation of a single unlocked nucleic acid monomer

A 15-mer DNA aptamer (named TBA) adopts a G-quadruplex structure that strongly inhibits fibrin-clot formation by binding to thrombin. We have performed thermodynamic analysis, binding affinity and biological activity studies of TBA variants modified by unlocked nucleic acid (UNA) monomers. UNA-U placed in position U3, U7 or U12 increases the thermodynamic stability of TBA by 0.15–0.50 kcal/mol. In contrast, modification of any position within the two G-quartet structural elements is unfavorable for quadruplex formation. The intramolecular folding of the quadruplexes is confirmed by T_m versus ln c analysis. Moreover, circular dichroism and thermal difference spectra of the modified TBAs displaying high thermodynamic stability show bands that are characteristic for antiparallel quadruplex formation. Surface plasmon resonance studies of the binding of the UNA-modified TBAs to thrombin show that a UNA monomer is allowed in many positions of the aptamer without significantly changing the thrombin-binding properties. The biological effect of a selection of the modified aptamers was tested by a thrombin time assay and showed that most of the UNA-modified TBAs possess anticoagulant properties, and that the construct with a UNA-U monomer in position 7 is a highly potent inhibitor of fibrin-clot formation.     

A sequence-independent study of the influence of short loop lengths on the stability and topology of intramolecular DNA G-quadruplexes

G-Rich sequences found within biologically important regions of the genome have been shown to form intramolecular G-quadruplexes with varied loop lengths and sequences. Many of these quadruplexes will be distinguishable from each other on the basis of their thermodynamic stabilities and folded conformations. It has been proposed that loop lengths can strongly influence the topology and stability of intramolecular G-quadruplexes. Previous studies have been limited to the analysis of quadruplex sequences with particular loop sequences, making it difficult to make generalizations. Here, we describe an original study that aimed to elucidate the effect of loop length on the biophysical properties of G-quadruplexes in a sequence-independent context. We employed UV melting and circular dichroism spectroscopy to examine and compare the properties of 21 DNA quadruplex libraries, each comprising partially randomized loop sequences with lengths ranging from one to three nucleotides. Our work supports a number of general predictions that can be made solely on the basis of loop lengths. In particular, the results emphasize the strong influence of single-nucleotide loops on quadruplex properties. This study provides a predictive framework that may help identify or classify biologically relevant G-quadruplex-forming sequences.     

Long promoter sequences form higher-order G-quadruplexes: an integrative structural biology study of c-Myc,k-Ras and c-Kit promoter sequences

We report on higher-order G-quadruplex structures adopted by long promoter sequences obtained by an iterative integrated structural biology approach. Our approach uses quantitative biophysical tools (analytical ultracentrifugation, small-angle X-ray scattering, and circular dichroism spectroscopy) combined with modeling and molecular dynamics simulations, to derive self-consistent structural models. The formal resolution of our approach is 18 angstroms, but in some cases structural features of only a few nucleotides can be discerned. We report here five structures of long (34–70 nt) wild-type sequences selected from three cancer-related promoters: c-Myc, c-Kit and k-Ras. Each sequence studied has a unique structure. Three sequences form structures with two contiguous, stacked, G-quadruplex units. One longer sequence from c-Myc forms a structure with three contiguous stacked quadruplexes. A longer c-Kit sequence forms a quadruplex-hairpin structure. Each structure exhibits interfacial regions between stacked quadruplexes or novel loop geometries that are possible druggable targets. We also report methodological advances in our integrated structural biology approach, which now includes quantitative CD for counting stacked G-tetrads, DNaseI cleavage for hairpin detection and SAXS model refinement. Our results suggest that higher-order quadruplex assemblies may be a common feature within the genome, rather than simple single quadruplex structures.     

Synthesis,biophysical characterization and anti-HIV activity of d(TG3AG) Quadruplexes bearing hydrophobic tails at the 5′-end

Novel conjugated G-quadruplex-forming d(TG₃AG) oligonucleotides, linked to hydrophobic groups through phosphodiester bonds at 5′-end, have been synthesized as potential anti-HIV aptamers, via a fully automated, online phosphoramidite-based solid-phase strategy. Conjugated quadruplexes showed pronounced anti-HIV activity with some preference for HIV-1, with inhibitory activity invariably in the low micromolar range. The CD and DSC monitored thermal denaturation studies on the resulting quadruplexes, indicated the insertion of lipophilic residue at the 5′-end, conferring always improved stability to the quadruplex complex (20 < ΔTm < 40 °C). The data suggest no direct functional relationship between the thermal stability and anti-HIV activity of the folded conjugated G-quartets. It would appear that the nature of the residue at 5′ end of the d(TG₃AG) quadruplexes plays an important role in the thermodynamic stabilization but a minor influence on the anti-HIV activity. Moreover, a detailed CD and DSC analyses indicate a monophasic behaviour for sequences I and V, while for ODNs (II–IV) clearly show that these quadruplex structures deviate from simple two-state melting, supporting the hypothesis that intermediate states along the dissociation pathway may exist.     

Effects of abasic sites on structural,thermodynamic and kinetic properties of quadruplex structures

Abasic sites represent the most frequent lesion in DNA. Since several events generating abasic sites concern guanines, this damage is particularly important in quadruplex forming G-rich sequences, many of which are believed to be involved in several biological roles. However, the effects of abasic sites in sequences forming quadruplexes have been poorly studied. Here, we investigated the effects of abasic site mimics on structural, thermodynamic and kinetic properties of parallel quadruplexes. Investigation concerned five oligodeoxynucleotides based on the sequence d(TGGGGGT), in which all guanines have been replaced, one at a time, by an abasic site mimic (dS). All sequences preserve their ability to form quadruplexes; however, both spectroscopic and kinetic experiments point to sequence-dependent different effects on the structural flexibility and stability. Sequences d(TSGGGGT) and d(TGGGGST) form quite stable quadruplexes; however, for the other sequences, the introduction of the dS in proximity of the 3′-end decreases the stability more considerably than the 5′-end. Noteworthy, sequence d(TGSGGGT) forms a quadruplex where dS does not hamper the stacking between the G-tetrads adjacent to it. These results strongly argue for the central role of apurinic/apyrimidinic site damages and they encourage the production of further studies to better delineate the consequences of their presence in the biological relevant regions of the genome.     

The formation pathway of tetramolecular G-quadruplexes

Oligonucleotides containing guanosine stretches associate into tetrameric structures stabilized by monovalent ions. In order to describe the sequence of reactions leading to association of four identical strands, we measured by NMR the formation and dissociation rates of (TGnT)₄ quadruplexes (n = 3–6), their dissociation constants and the reaction orders for quadruplex formation. The quadruplex formation rates increase with the salt concentration but weakly depend on the nature (K⁺, Na⁺ or Li⁺) of the counter ions. The activation energies for quadruplex formation are negative. The quadruplex lifetimes strongly increase with the G-tract length and are much more longer in K⁺ solution than in Na⁺ or Li⁺ solutions. The reaction order for quadruplex formation is 3 in 0.125 M KCl and 4 in LiCl solutions. The kinetics measurements suggest that quadruplex formation proceeds step by step via sequential strand association into duplex and triplex intermediate species. Triplex formation is rate limiting in 0.125 M KCl solution. In LiCl, each step of the association process depends on the strand concentration. Parallel reactions to formation of the fully matched canonical quadruplex may result in kinetically trapped mismatched quadruplexes making the canonical quadruplex practically inaccessible in particular at low temperature in KCl solution.     

Quadruplex-forming properties of FRAXA (CGG) repeats interrupted by (AGG) triplets

The (CGG) repeats associated with X-chromosome fragility are generally believed to form quadruplexes. This notion has persisted although it had been shown that only very short (CGG)_n sequences form quadruplexes and that this quadruplex formation occurs in conditions far from physiological. We have now studied, using CD and absorption spectroscopies, quadruplex formation of (CGG)_n (n = 4, 7, 8, or 16) and their analogs interrupted by (AGG) triplets under various solvent conditions. In healthy individuals, (AGG) triplets are interspersed throughout the (CGG) repeat regions and appear to hinder (CGG)_n motif expansion. Here we show that (CGG) repeats do not form quadruplexes under physiological conditions in aqueous solution but, interestingly, quadruplexes are readily formed in water–ethanol solutions. The presence of (AGG) triplets markedly stabilized quadruplex formation. Quadruplexes may thus hinder rather than support (CGG)_n motif expansion.     

Distance-dependent duplex DNA destabilization proximal to G-quadruplex/i-motif sequences

G-quadruplexes and i-motifs are complementary examples of non-canonical nucleic acid substructure conformations. G-quadruplex thermodynamic stability has been extensively studied for a variety of base sequences, but the degree of duplex destabilization that adjacent quadruplex structure formation can cause has yet to be fully addressed. Stable in vivo formation of these alternative nucleic acid structures is likely to be highly dependent on whether sufficient spacing exists between neighbouring duplex- and quadruplex-/i-motif-forming regions to accommodate quadruplexes or i-motifs without disrupting duplex stability. Prediction of putative G-quadruplex-forming regions is likely to be assisted by further understanding of what distance (number of base pairs) is required for duplexes to remain stable as quadruplexes or i-motifs form. Using oligonucleotide constructs derived from precedented G-quadruplexes and i-motif-forming bcl-2 P1 promoter region, initial biophysical stability studies indicate that the formation of G-quadruplex and i-motif conformations do destabilize proximal duplex regions. The undermining effect that quadruplex formation can have on duplex stability is mitigated with increased distance from the duplex region: a spacing of five base pairs or more is sufficient to maintain duplex stability proximal to predicted quadruplex/i-motif-forming regions.     

Predictive modelling of topology and loop variations in dimeric DNA quadruplex structures

We have used a combination of simulated annealing (SA), molecular dynamics (MD) and locally enhanced sampling (LES) methods in order to predict the favourable topologies and loop conformations of dimeric DNA quadruplexes with T2 or T3 loops. This follows on from our previous MD simulation studies on the influence of loop lengths on the topology of intramolecular quadruplex structures [P. Hazel et al. (2004) J. Am. Chem. Soc., 126, 16405–16415], which provided results consistent with biophysical data. The recent crystal structures of d(G4T3G4)2 and d(G4BrUT2G4) (P. Hazel et al. (2006) J. Am. Chem. Soc., in press) and the NMR-determined topology of d(TG4T2G4T)2 [A.T. Phan et al. (2004) J. Mol. Biol., 338, 93–102] have been used in the present study for comparison with simulation results. These together with MM-PBSA free-energy calculations indicate that lateral T3 loops are favoured over diagonal loops, in accordance with the experimental structures; however, distinct loop conformations have been predicted to be favoured compared to those found experimentally. Several lateral and diagonal loop conformations have been found to be similar in energy. The simulations suggest an explanation for the distinct patterns of observed dimer topology for sequences with T3 and T2 loops, which depend on the loop lengths, rather than only on G-quartet stability.     

Binding properties of human telomeric quadruplex multimers: a new route for drug design

Human telomeric G-quadruplex structures are known to be promising targets for an anticancer therapy. In the past decade, several research groups have been focused on the design of new ligands trying to optimize the interactions between these small molecules and the G-quadruplex motif. In most of these studies, the target structures were the single quadruplex units formed by short human DNA telomeric sequences (typically 21-26 nt). However, the 3′-terminal single-stranded human telomeric DNA is actually 100-200 bases long and can form higher-order structures by clustering several consecutive quadruplex units (multimers). Despite the increasing number of structural information on longer DNA telomeric sequences, very few data are available on the binding properties of these sequences compared with the shorter DNA telomeric sequences.In this paper we use a combination of spectroscopic (CD, UV and fluorescence) and calorimetric techniques (ITC) to compare the binding properties of the (TTAGGG)₈TT structure formed by two adjacent quadruplex units with the binding properties of the (AG₃TT)₄ single quadruplex structure. The three side-chained triazatruxene derivative azatrux and TMPyP4 cationic porphyrin were used as quadruplex ligands. We found that, depending on the drug, the number of binding sites per quadruplex unit available in the multimer structure was smaller or greater than the one expected on the basis of the results obtained from individual quadruplex binding studies. This work suggests that the quadruplex units along a multimer structure do not behave as completely independent. The presence of adjacent quadruplexes results in a diverse binding ability not predictable from single quadruplex binding studies. The existence of quadruplex-quadruplex interfaces in the full length telomeric overhang may provide an advantageous factor in drug design to enhance both affinity and selectivity for DNA telomeric quadruplexes.     

Using principal component analysis to find correlations between loop-related and thermodynamic variables for G-quadruplex-forming sequences

The application of Principal Component Analysis (PCA) is proposed here as a simple means of revealing correlations between thermodynamic variables corresponding to folding equilibria of intramolecular G-quadruplexes and Watson–Crick duplexes, and the length of loops in the corresponding guanine-rich DNA sequences. To this end, two previously studied data sets were analyzed (Arora and Maiti, J. Phys. Chem. B. 2009 and Kumar and Maiti, Nucleic Acids. Res. 2008). All of the sequences considered shared the common structure 5’- GGG - loop1 - GGG - loop2 - GGG - loop3 - GGG -3’. PCA of these data sets supported a series of correlations between the variables studied. First, the association of loop length with thermodynamic stability and quadruplex structure was corroborated. Secondly, it is proposed that the addition of ethylene glycol produces a stronger stabilization on those sequences showing long loop1 and/or loop3. Thirdly, it is proposed that a low content of adenine in loop1 and/or loop3 will produce an increase in the stability of G-quadruplex and its related Watson–Crick duplex.     

Effect of loop length variation on quadruplex-Watson Crick duplex competition

The effect of loop length on quadruplex stability has been studied when the G-rich strand is present along with its complementary C-rich strand, thereby resulting in competition between quadruplex and duplex structures. Using model sequences with loop lengths varying from T to T5, we carried out extensive FRET to discover the influence of loop length on the quadruplex-Watson Crick duplex competition. The binding data show an increase in the binding affinity of quadruplexes towards their complementary strands upon increasing the loop length. Our kinetic data reveal that unfolding of the quadruplex in presence of a complementary strand involves a contribution from a predominant slow and a small population of fast opening conformer. The contribution from the fast opening conformer increases upon increasing the loop length leading to faster duplex formation. FCS data show an increase in the interconversion between the quadruplex conformers in presence of the complementary strand, which shifts the equilibrium towards the fast opening conformer with an increase in loop length. The relative free-energy difference (Delta DeltaG(o)) between the duplex and quadruplex indicates that an increase in loop length favors duplex formation and out competes the quadruplex.     

Re-evaluation of G-quadruplex propensity with G4Hunter

Critical evidence for the biological relevance of G-quadruplexes (G4) has recently been obtained in seminal studies performed in a variety of organisms. Four-stranded G-quadruplex DNA structures are promising drug targets as these non-canonical structures appear to be involved in a number of key biological processes. Given the growing interest for G4, accurate tools to predict G-quadruplex propensity of a given DNA or RNA sequence are needed. Several algorithms such as Quadparser predict quadruplex forming propensity. However, a number of studies have established that sequences that are not detected by these tools do form G4 structures (false negatives) and that other sequences predicted to form G4 structures do not (false positives). Here we report development and testing of a radically different algorithm, G4Hunter that takes into account G-richness and G-skewness of a given sequence and gives a quadruplex propensity score as output. To validate this model, we tested it on a large dataset of 392 published sequences and experimentally evaluated quadruplex forming potential of 209 sequences using a combination of biophysical methods to assess quadruplex formation in vitro. We experimentally validated the G4Hunter algorithm on a short complete genome, that of the human mitochondria (16.6 kb), because of its relatively high GC content and GC skewness as well as the biological relevance of these quadruplexes near instability hotspots. We then applied the algorithm to genomes of a number of species, including humans, allowing us to conclude that the number of sequences capable of forming stable quadruplexes (at least in vitro) in the human genome is significantly higher, by a factor of 2–10, than previously thought.     

NMR solution structures of LNA (locked nucleic acid) modified quadruplexes

We have determined the NMR solution structures of the quadruplexes formed by d(TGLGLT) and d(TL₄T), where L denotes LNA (locked nucleic acid) modified G-residues. Both structures are tetrameric, parallel and right-handed and the native global fold of the corresponding DNA quadruplex is retained upon introduction of the LNA nucleotides. However, local structural alterations are observed owing to the locked LNA sugars. In particular, a distinct change in the sugar–phosphate backbone is observed at the G2pL3 and L2pL3 base steps and sequence dependent changes in the twist between tetrads are also seen. Both the LNA modified quadruplexes have raised thermostability as compared to the DNA quadruplex. The quadruplex-forming capability of d(TGLGLT) is of particular interest as it expands the design flexibility for stable parallel LNA quadruplexes and shows that LNA nucleotides can be mixed with DNA or other modified nucleic acids. As such, LNA-based quadruplexes can be decorated by a variety of chemical modifications. Such LNA quadruplex scaffolds might find applications in the developing field of nanobiotechnology.