Kinetic Partitioning Modulates Human Telomere DNA G-Quadruplex Structural Polymorphism

Telomeres are specialized chromatin structures found at the end of chromosomes and are crucial to the maintenance of eukaryotic genome stability. Human telomere DNA is comprised of the repeating sequence (T₂AG₃)_n, which is predominantly double-stranded but terminates with a 3’ single-stranded tail. The guanine-rich tail can fold into secondary structures known as a G-quadruplexes (GQs) that may exist as a polymorphic mixture of anti-parallel, parallel, and several hybrid topological isomers. Using single-molecule Förster resonance energy transfer (smFRET), we have reconstructed distributions of telomere DNA GQ conformations generated by an in situ refolding protocol commonly employed in single-molecule studies of GQ structure, or using a slow cooling DNA annealing protocol typically used in the preparation of GQ samples for ensemble biophysical analyses. We find the choice of GQ folding protocol has a marked impact on the observed distributions of DNA conformations under otherwise identical buffer conditions. A detailed analysis of the kinetics of GQ folding over timescales ranging from minutes to hours revealed the distribution of GQ structures generated by in situ refolding gradually equilibrates to resemble the distribution generated by the slow cooling DNA annealing protocol. Interestingly, conditions of low ionic strength, which promote transient GQ unfolding, permit the fraction of folded DNA molecules to partition into a distribution that more closely approximates the thermodynamic folding equilibrium. Our results are consistent with a model in which kinetic partitioning occurs during in situ folding at room temperature in the presence of K⁺ ions, producing a long-lived non-equilibrium distribution of GQ structures in which the parallel conformation predominates on the timescale of minutes. These results suggest that telomere DNA GQ folding kinetics, and not just thermodynamic stability, likely contributes to the physiological ensemble GQ structures.     

The first crystal structures of hybrid and parallel four-tetrad intramolecular G-quadruplexes

G-quadruplexes (GQs) are non-canonical DNA structures composed of stacks of stabilized G-tetrads. GQs play an important role in a variety of biological processes and may form at telomeres and oncogene promoters among other genomic locations. Here, we investigate nine variants of telomeric DNA from Tetrahymena thermophila with the repeat (TTGGGG)_n. Biophysical data indicate that the sequences fold into stable four-tetrad GQs which adopt multiple conformations according to native PAGE. Excitingly, we solved the crystal structure of two variants, TET25 and TET26. The two variants differ by the presence of a 3′-T yet adopt different GQ conformations. TET25 forms a hybrid [3 + 1] GQ and exhibits a rare 5′-top snapback feature. Consequently, TET25 contains four loops: three lateral (TT, TT, and GTT) and one propeller (TT). TET26 folds into a parallel GQ with three TT propeller loops. To the best of our knowledge, TET25 and TET26 are the first reported hybrid and parallel four-tetrad unimolecular GQ structures. The results presented here expand the repertoire of available GQ structures and provide insight into the intricacy and plasticity of the 3D architecture adopted by telomeric repeats from T. thermophila and GQs in general.     

G-quadruplex conformation and dynamics are determined by loop length and sequence

The quadruplex forming G-rich sequences are unevenly distributed throughout the human genome. Their enrichment in oncogenic promoters and telomeres has generated interest in targeting G-quadruplex (GQ) for an anticancer therapy. Here, we present a quantitative analysis on the conformations and dynamics of GQ forming sequences measured by single molecule fluorescence. Additionally, we relate these properties to GQ targeting ligands and G4 resolvase 1 (G4R1) protein binding. Our result shows that both the loop (non-G components) length and sequence contribute to the conformation of the GQ. Real time single molecule traces reveal that the folding dynamics also depend on the loop composition. We demonstrate that GQ-stabilizing small molecules, N-methyl mesoporphyrin IX (NMM), its analog, NMP and the G4R1 protein bind selectively to the parallel GQ conformation. Our findings point to the complexity of GQ folding governed by the loop length and sequence and how the GQ conformation determines the small molecule and protein binding propensity.     

Frustrated folding of guanine quadruplexes in telomeric DNA

Human chromosomes terminate in long, single-stranded, DNA overhangs of the repetitive sequence (TTAGGG)n. Sets of four adjacent TTAGGG repeats can fold into guanine quadruplexes (GQ), four-stranded structures that are implicated in telomere maintenance and cell immortalization and are targets in cancer therapy. Isolated GQs have been studied in detail, however much less is known about folding in long repeat sequences. Such chains adopt an enormous number of configurations containing various arrangements of GQs and unfolded gaps, leading to a highly frustrated energy landscape. To better understand this phenomenon, we used mutagenesis, thermal melting, and global analysis to determine stability, kinetic, and cooperativity parameters for GQ folding within chains containing 8–12 TTAGGG repeats. We then used these parameters to simulate the folding of 32-repeat chains, more representative of intact telomeres. We found that a combination of folding frustration and negative cooperativity between adjacent GQs increases TTAGGG unfolding by up to 40-fold, providing an abundance of unfolded gaps that are potential binding sites for telomeric proteins. This effect was most pronounced at the chain termini, which could promote telomere extension by telomerase. We conclude that folding frustration is an important and largely overlooked factor controlling the structure of telomeric DNA.     

NMR structural study on the self-trimerization of d(GTTAGG) into a dynamic trimolecular G-quadruplex assembly preferentially in Na+ solution with a moderate K+ tolerance

Vast G-quadruplexes (GQs) are primarily folded by one, two, or four G-rich oligomers, rarely with an exception. Here, we present the first NMR solution structure of a trimolecular GQ (tri-GQ) that is solely assembled by the self-trimerization of d(GTTAGG), preferentially in Na⁺ solution tolerant to an equal amount of K⁺ cation. Eight guanines from three asymmetrically folded strands of d(GTTAGG) are organized into a two-tetrad core, which features a broken G-column and two width-irregular grooves. Fast strand exchanges on a timescale of second at 17°C spontaneously occur between folded tri-GQ and unfolded single-strand of d(GTTAGG) that both species coexist in dynamic equilibrium. Thus, this tri-GQ is not just simply a static assembly but rather a dynamic assembly. Moreover, another minor tetra-GQ that has putatively tetrameric (2+2) antiparallel topology becomes noticeable only at an extremely high strand concentration above 18 mM. The major tri-GQ and minor tetra-GQ are considered to be mutually related, and their reversible interconversion pathways are proposed accordingly. The sequence d(GTTAGG) could be regarded as either a reading frame shifted single repeat of human telomeric DNA or a 1.5 repeat of Bombyx mori telomeric DNA. Overall, our findings provide new insight into GQs and expect more functional applications.     

Conformational polymorphysm of G-rich fragments of DNA Alu-repeats. II. The putative role of G-quadruplex structures in genomic rearrangements

Three evolutionary conserved (G-rich) sites of Alu repeats (PQS2, PQS3, and PQS4) could form in vitro stable inter- and intramolecular G-quadruplexes (GQs). Structures and topologies of these GQs were elucidated using spectral methods. The study of self-association of G-rich Alu fragments performed using a FRET-based method revealed dimeric GQ formation from two distally located sites (PQS2)₂, (PQS3)₂ or PQS2?PQS3 within one extended single stranded DNA. Using DOSY NMR, AFM microscopy and differential CD spectroscopy it has been demonstrated that oligomer PQS4 (folded into a parallel intramolecular GQ) forms stacks of quadruplexes stabilized by stacking interactions of external G-tetrads. Comparative analysis of the properties of various GQs suggests involvement of two universal general mechanisms of GQ-dependent genomic rearrangements: (i) formation of dimeric GQs from fragments of different molecules; (ii) formation of GQ-GQ-stacks from pre-folded intramolecular parallel GQs from different strands. Thus, association of G-rich Alu motifs with sensitivity to double-strand breaks and rearrangements may be attributed not to structural features of G-rich Alu fragments, but also to their high abundance.     

Conformation-sensitive nucleoside analogues as topology-specific fluorescence turn-on probes for DNA and RNA G-quadruplexes

Development of probes that can discriminate G-quadruplex (GQ) structures and indentify efficient GQ binders on the basis of topology and nucleic acid type is highly desired to advance GQ-directed therapeutic strategies. In this context, we describe the development of minimally perturbing and environment-sensitive pyrimidine nucleoside analogues, based on a 5-(benzofuran-2-yl)uracil core, as topology-specific fluorescence turn-on probes for human telomeric DNA and RNA GQs. The pyrimidine residues of one of the loop regions (TTA) of telomeric DNA and RNA GQ oligonucleotide (ON) sequences were replaced with 5-benzofuran-modified 2′-deoxyuridine and uridine analogues. Depending on the position of modification the fluorescent nucleoside analogues distinguish antiparallel, mixed parallel-antiparallel and parallel stranded DNA and RNA GQ topologies from corresponding duplexes with significant enhancement in fluorescence intensity and quantum yield. Further, these GQ sensors enabled the development of a simple fluorescence binding assay to quantify topology- and nucleic acid-specific binding of small molecule ligands to GQ structures. Together, our results demonstrate that these nucleoside analogues are useful GQ probes, which are anticipated to provide new opportunities to study and discover efficient G-quadruplex binders of therapeutic potential.     

Interaction of human telomeric DNA with N-methyl mesoporphyrin IX

The remarkable selectivity of N-methyl mesoporphyrin IX (NMM) for G-quadruplexes (GQs) is long known, however its ability to stabilize and bind GQs has not been investigated in detail. Through the use of circular dichroism, UV-visible spectroscopy and fluorescence resonance energy transfer (FRET) melting assay we have shown that NMM stabilizes human telomeric DNA dAG(3)(TTAG(3))(3) (Tel22) and is selective for its parallel conformation to which it binds in 1:1 stoichiometry with a binding constant of ≈ 1.0 × 10(5)M(-1). NMM does not interact with an antiparallel conformation of Tel22 in sodium buffer and is the second example in the literature, after TOxaPy, of a ligand with an excellent selectivity for a specific GQ structure. NMM's stabilizing ability toward predominantly parallel GQ conformation is universal: it stabilizes a variety of biologically relevant G-rich sequences including telomeres and oncogene promoters. The N-methyl group is integral for selectivity and stabilization, as the unmethylated analogue, mesoporphyrin IX, does not stabilize GQ DNA in FRET melting assays. Finally, NMM induces the isomerization of Tel22 into a structure with increased parallel component in K(+) but not in Na(+) buffer. The ability of NMM to cause structural rearrangement and efficient stabilization of Tel22 may bear biological significance.     

Induction of G-quadruplex DNA structure by Zn(II) 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin

G-quadruplexes (GQ) are formed by the association of guanine-rich stretches of DNA. Certain small molecules can influence kinetics and thermodynamics of this association. Understanding the mechanism of ligand-assisted GQ folding is necessary for the design of more efficient cancer therapeutics. The oligonucleotide d(TAGGG)₂ forms parallel bimolecular GQ in the presence of ≥66 mM K⁺; GQs are not formed under Na⁺, Li⁺ or low K⁺ conditions. The thermodynamic parameters for GQ folding at 60 μM oligonucleotide and 100 mM KCl are ΔH = −35 ± 2 kcal mol⁻¹ and ΔG₃₁₀ = −1.4 kcal mol⁻¹. Quadruplex [d(TAGGG)₂]₂ binds 2-3 K⁺ ions with K_d of 0.5 ± 0.2 mM. Our work addresses the question of whether metal free 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin (TMPyP4) and its Zn(II), Cu(II), and Pt(II) derivatives are capable of facilitating GQ folding of d(TAGGG)₂ from single stranded, or binding to preformed GQ, using UV-vis and circular dichroism (CD) spectroscopies. ZnTMPyP4 is unique among other porphyrins in its ability to induce GQ structure of d(TAGGG)₂, which also requires at least a low amount of potassium. ZnTMPyP4 binds with 2:1 stoichiometry possibly in an end-stacking mode with a ∼10⁶ M⁻¹ binding constant, determined through UV-vis and ITC titrations. This process is entropically driven and has ΔG₂₉₈ of −8.0 kcal mol⁻¹. TMPyP4 binds with 3:1 stoichiometry and K_a of ∼10⁶ M⁻¹. ZnTMPyP4 and TMPyP4 are efficient stabilizers of [d(TAGGG)₂]₂ displaying ΔT_1/2 of 13.5 and 13.8 °C, respectively, at 1:2 GQ to porphyrin ratio; CuTMPyP4 shows a much weaker effect (ΔT_1/2 = 4.7 °C) and PtTMPyP4 is weakly destabilizing (ΔT_1/2 = −2.9 °C). The selectivity of ZnTMPyP4 for GQ versus dsDNA is comparable to that of TMPyP4. The ability of ZnTMPyP4 to bind and stabilize GQ, to induce GQ formation, and speed up its folding may suggest an important biological activity for this molecule.     

Functional characterization of an alkaline exonuclease and single strand annealing protein from the SXT genetic element of Vibrio cholerae

Background

SXT is an integrating conjugative element (ICE) originally isolated from Vibrio cholerae, the bacterial pathogen that causes cholera. It houses multiple antibiotic and heavy metal resistance genes on its ca. 100 kb circular double stranded DNA (dsDNA) genome, and functions as an effective vehicle for the horizontal transfer of resistance genes within susceptible bacterial populations. Here, we characterize the activities of an alkaline exonuclease (S066, SXT-Exo) and single strand annealing protein (S065, SXT-Bet) encoded on the SXT genetic element, which share significant sequence homology with Exo and Bet from bacteriophage lambda, respectively.

Results

SXT-Exo has the ability to degrade both linear dsDNA and single stranded DNA (ssDNA) molecules, but has no detectable endonuclease or nicking activities. Adopting a stable trimeric arrangement in solution, the exonuclease activities of SXT-Exo are optimal at pH 8.2 and essentially require Mn²⁺ or Mg²⁺ ions. Similar to lambda-Exo, SXT-Exo hydrolyzes dsDNA with 5'- to 3'-polarity in a highly processive manner, and digests DNA substrates with 5'-phosphorylated termini significantly more effectively than those lacking 5'-phosphate groups. Notably, the dsDNA exonuclease activities of both SXT-Exo and lambda-Exo are stimulated by the addition of lambda-Bet, SXT-Bet or a single strand DNA binding protein encoded on the SXT genetic element (S064, SXT-Ssb). When co-expressed in E. coli cells, SXT-Bet and SXT-Exo mediate homologous recombination between a PCR-generated dsDNA fragment and the chromosome, analogous to RecET and lambda-Bet/Exo.

Conclusions

The activities of the SXT-Exo protein are consistent with it having the ability to resect the ends of linearized dsDNA molecules, forming partially ssDNA substrates for the partnering SXT-Bet single strand annealing protein. As such, SXT-Exo and SXT-Bet may function together to repair or process SXT genetic elements within infected V. cholerae cells, through facilitating homologous DNA recombination events. The results presented here significantly extend our general understanding of the properties and activities of alkaline exonuclease and single strand annealing proteins of viral/bacteriophage origin, and will assist the rational development of bacterial recombineering systems.     

ATP-dependent G-quadruplex unfolding by Bloom helicase exhibits low processivity

Various helicases and single stranded DNA (ssDNA) binding proteins unfold G-quadruplex (GQ) structures. However, the underlying mechanisms of this activity have only recently come to focus. We report kinetic studies on Bloom (BLM) helicase and human telomeric GQ interactions using single-molecule Förster resonance energy transfer (smFRET). Using partial duplex DNA (pdDNA) constructs with different 5′ ssDNA overhangs, we show that BLM localizes in the vicinity of ssDNA/double-stranded DNA (dsDNA) junction and reels in the ssDNA overhang in an ATP-dependent manner. A comparison of DNA constructs with or without GQ in the overhang shows that GQ unfolding is achieved in 50–70% of reeling attempts under physiological salt and pH conditions. The unsuccessful attempts often result in dissociation of BLM from DNA which slows down the overall BLM activity. BLM-mediated GQ unfolding is typically followed by refolding of the GQ, a pattern that is repeated several times before BLM dissociates from DNA. BLM is significantly less processive compared to the highly efficient GQ destabilizer Pif1 that can repeat GQ unfolding activity hundreds of times before dissociating from DNA. Despite the variations in processivity, our studies point to possible common patterns used by different helicases in minimizing the duration of stable GQ formation.     

Interrogating accessibility of telomeric sequences with FRET-PAINT: evidence for length-dependent telomere compaction

Single-stranded telomeric overhangs are ∼200 nucleotides long and can form tandem G-quadruplex (GQ) structures, which reduce their accessibility to nucleases and proteins that activate DNA damage response. Whether these tandem GQs further stack to form compact superstructures, which may provide better protection for longer telomeres, is not known. We report single-molecule measurements where the accessibility of 24–144 nucleotide long human telomeric DNA molecules is interrogated by a short PNA molecule that is complementary to a single GGGTTA repeat, as implemented in the FRET-PAINT method. Binding of the PNA strand to available GGGTTA sequences results in discrete FRET bursts which were analyzed in terms of their dwell times, binding frequencies, and topographic distributions. The binding frequencies were greater for binding to intermediate regions of telomeric DNA compared to 3′- or 5′-ends, suggesting these regions are more accessible. Significantly, the binding frequency per telomeric repeat monotonically decreased with increasing telomere length. These results are consistent with telomeres forming more compact structures at longer lengths, reducing accessibility of these critical genomic sites.     

Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations

DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.     

Folding and persistence times of intramolecular G-quadruplexes transiently embedded in a DNA duplex

G-quadruplex (G4) DNA structures have emerged as important regulatory elements during DNA metabolic transactions. While many in vitro studies have focused on the kinetics of G4 formation within DNA single-strands, G4 are found in vivo in double-stranded DNA regions, where their formation is challenged by the complementary strand. Since the energy of hybridization of Watson-Crick structures dominates the energy of G4 folding, this competition should play a critical role on G4 persistence. To address this, we designed a single-molecule assay allowing to measure G4 folding and persistence times in the presence of the complementary strand. We quantified both folding and unfolding rates of biologically relevant G4 sequences, such as the cMYC and cKIT oncogene promoters, human telomeres and an avian replication origin. We confirmed that G4s are found much more stable in tested replication origin and promoters than in human telomere repeats. In addition, we characterized how G4 dynamics was affected by G4 ligands and showed that both folding rate and persistence time increased. Our assay opens new perspectives for the measurement of G4 dynamics in double-stranded DNA mimicking a replication fork, which is important to understand their role in DNA replication and gene regulation at a mechanistic level.     

The N-terminal domain of Escherichia coli RecA have multiple functions in promoting homologous recombination

Escherichia coli RecA mediates homologous recombination, a process essential to maintaining genome integrity. In the presence of ATP, RecA proteins bind a single-stranded DNA (ssDNA) to form a RecA-ssDNA presynaptic nucleoprotein filament that captures donor double-stranded DNA (dsDNA), searches for homology, and then catalyzes the strand exchange between ssDNA and dsDNA to produce a new heteroduplex DNA. Based upon a recently reported crystal structure of the RecA-ssDNA nucleoprotein filament, we carried out structural and functional studies of the N-terminal domain (NTD) of the RecA protein. The RecA NTD was thought to be required for monomer-monomer interaction. Here we report that it has two other distinct roles in promoting homologous recombination. It first facilitates the formation of a RecA-ssDNA presynaptic nucleoprotein filament by converting ATP to an ADP-Pi intermediate. Then, once the RecA-ssDNA presynaptic nucleoprotein filament is stably assembled in the presence of ATPγS, the NTD is required to capture donor dsDNA. Our results also suggest that the second function of NTD may be similar to that of Arg243 and Lys245, which were implicated earlier as binding sites of donor dsDNA. A two-step model is proposed to explain how a RecA-ssDNA presynaptic nucleoprotein filament interacts with donor dsDNA.