DNA-RNA hybrid G-quadruplex tends to form near the 3′ end of telomere overhang

Telomeric repeat-containing RNA (TERRA) has been suggested to participate in telomere maintenance. TERRA consisting of UUAGGG repeats is capable of forming an intermolecular G-quadruplex (GQ) with single-stranded TTAGGG-repeat DNA in the telomere 3′ overhang. To explore the structural features and potential functions of this DNA-RNA hybrid GQ (HGQ), we used single-molecule FRET to study the folding patterns of DNA with four to seven telomeric tandem repeats annealed with a short RNA consisting of two or five telomeric repeats. Our data highlight that RNA prefers to form DNA-RNA HGQ near the 3′ end of telomeric DNA. Furthermore, the unfolding of secondary structures by a complementary C-rich sequence was observed for DNA GQ but not for DNA-RNA HGQ, which demonstrated the enhanced stability of the telomere 3′ end via hybridization with RNA. These conformational and physical properties of telomeric DNA-RNA HGQ suggest that TERRA might limit access to the 3′ end of the telomeric DNA overhang, which is known to be critical for the interaction with telomerase and other telomere-associated proteins.     

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.     

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.     

Mechanical unfolding of human telomere G-quadruplex DNA probed by integrated fluorescence and magnetic tweezers spectroscopy

Single-molecule techniques facilitate analysis of mechanical transitions within nucleic acids and proteins. Here, we describe an integrated fluorescence and magnetic tweezers instrument that permits detection of nanometer-scale DNA structural rearrangements together with the application of a wide range of stretching forces to individual DNA molecules. We have analyzed the force-dependent equilibrium and rate constants for telomere DNA G-quadruplex (GQ) folding and unfolding, and have determined the location of the transition state barrier along the well-defined DNA-stretching reaction coordinate. Our results reveal the mechanical unfolding pathway of the telomere DNA GQ is characterized by a short distance (<1 nm) to the transition state for the unfolding reaction. This mechanical unfolding response reflects a critical contribution of long-range interactions to the global stability of the GQ fold, and suggests that telomere-associated proteins need only disrupt a few base pairs to destabilize GQ structures. Comparison of the GQ unfolded state with a single-stranded polyT DNA revealed the unfolded GQ exhibits a compacted non-native conformation reminiscent of the protein molten globule. We expect the capacity to interrogate macromolecular structural transitions with high spatial resolution under conditions of low forces will have broad application in analyses of nucleic acid and protein folding.     

Direct unfolding of RuvA-HJ complex at the single-molecule level

The repair of double-stranded DNA breaks via homologous recombination involves a four-way cross-strand intermediate known as Holliday junction (HJ), which is recognized, processed, and resolved by a specific set of proteins. RuvA, a prokaryotic HJ-binding protein, is known to stabilize the square-planar conformation of the HJ, which is otherwise a short-lived intermediate. Despite much progress being made regarding the molecular mechanism of RuvA-HJ interactions, the mechanochemical aspect of this protein-HJ complex is yet to be investigated. Here, we employed an optical-tweezers-based, single-molecule manipulation assay to detect the formation of RuvA-HJ complex and determined its mechanical and thermodynamic properties in a manner that would be impossible with traditional ensemble techniques. We found that the binding of RuvA increases the unfolding force (F_unfold) of the HJ by ～2-fold. Compared with the ΔG_unfold of the HJ alone (54 ± 13 kcal/mol), the increased free energy of the RuvA-HJ complex (101 ± 20 kcal/mol) demonstrates that the RuvA protein stabilizes HJs. Interestingly, the protein remains bound to the mechanically melted HJ, facilitating its refolding at an unusually high force when the stretched DNA molecule is relaxed. These results suggest that the RuvA protein not only stabilizes the HJs but also induces refolding of the HJs. The single-molecule platform that we employed here for studying the RuvA-HJ interaction is broadly applicable to study other HJ-binding proteins involved in the critical DNA repair process.     

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.     

Single-Molecule TPM Studies on the Conversion of Human Telomeric DNA

Human telomere contains guanine-rich (G-rich) tandem repeats of single-stranded DNA sequences at its 3′ tail. The G-rich sequences can be folded into various secondary structures, termed G-quadruplexes (G4s), by Hoogsteen basepairing in the presence of monovalent cations (such as Na⁺ and K⁺). We developed a single-molecule tethered particle motion (TPM) method to investigate the unfolding process of G4s in the human telomeric sequence AGGG(TTAGGG)₃ in real time. The TPM method monitors the DNA tether length change caused by formation of the G4, thus allowing the unfolding process and structural conversion to be monitored at the single-molecule level. In the presence of its antisense sequence, the folded G4 structure can be disrupted and converted to the unfolded conformation, with apparent unfolding time constants of 82 s and 3152 s. We also observed that the stability of the G4 is greatly affected by different monovalent cations. The folding equilibrium constant of G4 is strongly dependent on the salt concentration, ranging from 1.75 at 5 mM Na⁺ to 3.40 at 15 mM Na⁺. Earlier spectral studies of Na⁺- and K⁺-folded states suggested that the spectral conversion between these two different folded structures may go through a structurally unfolded intermediate state. However, our single-molecule TPM experiments did not detect any totally unfolded intermediate within our experimental resolution when sodium-folded G4 DNA molecules were titrated with high-concentration, excess potassium ions. This observation suggests that a totally unfolding pathway is likely not the major pathway for spectral conversion on the timescale of minutes, and that interconversion among folded states can be achieved by the loop rearrangement. This study also demonstrates that TPM experiments can be used to study conformational changes in single-stranded DNA molecules.     

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.     

Dynamics of TRF1 organizing a single human telomere

Chromosome stability is primarily determined by telomere length. TRF1 is the core subunit of shelterin that plays a critical role in telomere organization and replication. However, the dynamics of TRF1 in scenarios of telomere-processing activities remain elusive. Using single-molecule magnetic tweezers, we here investigated the dynamics of TRF1 upon organizing a human telomere and the protein-DNA interactions at a moving telomeric fork. We first developed a method to obtain telomeres from human cells for directly measuring the telomere length by single-molecule force spectroscopy. Next, we examined the compaction and decompaction of a telomere by TRF1 dimers. TRF1 dissociates from a compacted telomere with heterogenous loops in ∼20 s. We also found a negative correlation between the number of telomeric loops and loop sizes. We further characterized the dynamics of TRF1 at a telomeric DNA fork. With binding energies of 11 k_BT, TRF1 can modulate the forward and backward steps of DNA fork movements by 2–9 s at a critical force of F_1/2, temporarily maintaining the telomeric fork open. Our results shed light on the mechanisms of how TRF1 organizes human telomeres and facilitates the efficient replication of telomeric DNA. Our work will help future research on the chemical biology of telomeres and shelterin-targeted drug discovery.     

Single-Molecule Studies of the Im7 Folding Landscape

Under appropriate conditions, the four-helical Im7 (immunity protein 7) folds from an ensemble of unfolded conformers to a highly compact native state via an on-pathway intermediate. Here, we investigate the unfolded, intermediate, and native states populated during folding using diffusion single-pair fluorescence resonance energy transfer by measuring the efficiency of energy transfer (or proximity or P ratio) between pairs of fluorophores introduced into the side chains of cysteine residues placed in the center of helices 1 and 4, 1 and 3, or 2 and 4. We show that while the native states of each variant give rise to a single narrow distribution with high P values, the distributions of the intermediates trapped at equilibrium (denoted I^eqm) are fitted by two Gaussian distributions. Modulation of the folding conditions from those that stabilize the intermediate to those that destabilize the intermediate enabled the distribution of lower P value to be assigned to the population of the unfolded ensemble in equilibrium with the intermediate state. The reduced stability of the I^eqm variants allowed analysis of the effect of denaturant concentration on the compaction and breadth of the unfolded state ensemble to be quantified from 0 to 6 M urea. Significant compaction is observed as the concentration of urea is decreased in both the presence and absence of sodium sulfate, as previously reported for a variety of proteins. In the presence of Na₂SO₄ in 0 M urea, the P value of the unfolded state ensemble approaches that of the native state. Concurrent with compaction, the ensemble displays increased peak width of P values, possibly reflecting a reduction in the rate of conformational exchange among iso-energetic unfolded, but compact conformations. The results provide new insights into the initial stages of folding of Im7 and suggest that the unfolded state is highly conformationally constrained at the outset of folding.     

Assembly and activation of a kinase ribozyme

RNA activities can be regulated by modulating the relative energies of all conformations in a folding landscape; however, it is often unknown precisely how peripheral elements perturb the overall landscape in the absence of discrete alternative folds (inactive ensemble). This work explores the effects of sequence and secondary structure in governing kinase ribozyme activity. Kin.46 catalyzes thiophosphoryl transfer from ATPγS onto the 5′ hydroxyl of polynucleotide substrates, and is regulated 10,000-fold by annealing an effector oligonucleotide to form activator helix P4. Transfer kinetics for an extensive series of ribozyme variants identified several dispensable internal single-stranded segments, in addition to a potential pseudoknot at the active site between segments J1/4 and J3/2 that is partially supported by compensatory rescue. Standard allosteric mechanisms were ruled out, such as formation of discrete repressive structures or docking P4 into the rest of the ribozyme via backbone 2′ hydroxyls. Instead, P4 serves both to complete an important structural element (100-fold contribution to the reaction relative to a P4-deleted variant) and to mitigate nonspecific, inhibitory effects of the single-stranded tail (an additional 100-fold contribution to the apparent rate constant, k_obs). Thermodynamic activation parameters ΔH^‡ and ΔS^‡, calculated from the temperature dependence of k_obs, varied with tail length and sequence. Inhibitory effects of the unpaired tail are largely enthalpic for short tails and are both enthalpic and entropic for longer tails. These results refine the structural view of this kinase ribozyme and highlight the importance of nonspecific ensemble effects in conformational regulation by peripheral elements.     

RPA-Mediated Unfolding of Systematically Varying G-Quadruplex Structures

G-quadruplex (GQ) is a noncanonical nucleic acid structure that is formed by guanine rich sequences. Unless it is destabilized by proteins such as replication protein A (RPA), GQ could interfere with DNA metabolic functions, such as replication or repair. We studied RPA-mediated GQ unfolding using single-molecule FRET on two groups of GQ structures that have different loop lengths and different numbers of G-tetrad layers. We observed a linear increase in the steady-state stability of the GQ against RPA-mediated unfolding with increasing number of layers or decreasing loop length. The stability demonstrated by different GQ structures varied by at least three orders of magnitude. Those with shorter loops (less than three nucleotides long) or a greater number of layers (more than three layers) maintained a significant folded population even at physiological RPA concentration (≈1 μM), raising the possibility of physiological viability of such GQ structures. Finally, we measured the transition time between the start and end of the RPA-mediated GQ unfolding process to be 0.35 ± 0.10 s for all GQ constructs we studied, despite significant differences in their steady-state stabilities. We propose a two-step RPA-mediated GQ unfolding mechanism that is consistent with our observations.     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

Structural Characteristic of the Initial Unfolded State on Refolding Determines Catalytic Efficiency of the Folded Protein in Presence of Osmolytes

Osmolytes are low molecular weight organic molecules accumulated by organisms to assist proper protein folding, and to provide protection to the structural integrity of proteins under denaturing stress conditions. It is known that osmolyte-induced protein folding is brought by unfavorable interaction of osmolytes with the denatured/unfolded states. The interaction of osmolyte with the native state does not significantly contribute to the osmolyte-induced protein folding. We have therefore investigated if different denatured states of a protein (generated by different denaturing agents) interact differently with the osmolytes to induce protein folding. We observed that osmolyte-assisted refolding of protein obtained from heat-induced denatured state produces native molecules with higher enzyme activity than those initiated from GdmCl- or urea-induced denatured state indicating that the structural property of the initial denatured state during refolding by osmolytes determines the catalytic efficiency of the folded protein molecule. These conclusions have been reached from the systematic measurements of enzymatic kinetic parameters (K _m and k _cat), thermodynamic stability (T _m and ΔH _m) and secondary and tertiary structures of the folded native proteins obtained from refolding of various denatured states (due to heat-, urea- and GdmCl-induced denaturation) of RNase-A in the presence of various osmolytes.     

Secondary structure polymorphism in Oxytricha nova telomeric DNA

Tandem repeats of the telomeric DNA sequence d(T₄G₄) of Oxytricha nova are capable of forming unusually stable secondary structures incorporating Hoogsteen hydrogen bonding interactions. The biological significance of such DNA structures is supported by evidence of specific recognition of telomere end-binding proteins in the crystal state. To further characterize structural polymorphism of Oxytricha telomeric DNAs, we have obtained and interpreted Raman, ultraviolet resonance Raman (UVRR) and circular dichroism (CD) spectra of the tandem repeats d(G₄T₄G₄) (Oxy1.5), d(T₄G₄)₂ (Oxy2) and dT₆(T₄G₄)₂ (T6Oxy2) and related non-telomeric isomers in aqueous salt solutions. Raman markers of Oxy1.5 identify both C2′-endo/anti and C2′-endo/syn conformations of the deoxyguanosine residues and Hoogsteen hydrogen bonded guanine quartets, consistent with the quadruplex fold determined previously by solution NMR spectroscopy. Raman, UVRR and CD signatures and Raman dynamic measurements, to monitor imino NH→ND exchanges, show that the Oxy1.5 antiparallel quadruplex fold is distinct from the hairpin structures of Oxy2 and T6Oxy2, single-stranded structures of d(TG)₈ and dT₆(TG)₈ and previously reported quadruplex structures of d(T₄G₄)₄ (Oxy4) and dG₁₂. Spectral markers of the telomeric and telomere-related DNA structures are tabulated and novel Raman and UVRR indicators of thymidine and deoxyguanosine conformations are identified. The results will be useful for probing structures of Oxytricha telomeric repeats in complexes with telomere end-binding proteins.     

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.     

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.     

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.