Kinetic partitioning modulates human telomere DNA G-quadruplex structural polymorphism

Balancing act: BRCA2's elaborate management of telomere replication through control of G-quadruplex dynamicity

In billion years of evolution, eukaryotes preserved the chromosome ends with arrays of guanine repeats surrounded by thymines and adenines, which can form stacks of four-stranded planar structure known as G-quadruplex (G4). The rationale behind the evolutionary conservation of the G4 structure at the telomere remained elusive. Our recent study has shed light on this matter by revealing that telomere G4 undergoes oscillation between at least two distinct folded conformations. Additionally, tumor suppressor BRCA2 exhibits a unique mode of interaction with telomere G4. To elaborate, BRCA2 directly interacts with G-triplex (G3)-derived intermediates that form during the interconversion of the two different G4 states. In doing so, BRCA2 remodels the G4, facilitating the restart of stalled replication forks. In this review, we succinctly summarize the findings regarding the dynamicity of telomeric G4, emphasize its importance in maintaining telomere replication homeostasis, and the physiological consequences of losing G4 dynamicity at the telomere.

Exploring the interaction between a fluorescent Ag(I)-biscarbene complex and non-canonical DNA structures: a multi-technique investigation

Silver compounds are mainly studied as antimicrobial agents, but they also have anticancer properties, with the latter, in some cases, being better than their gold counterparts. Herein, we analyse the first example of a new Ag(I)-biscarbene that can bind non-canonical structures of DNA, more precisely G-quadruplexes (G4), with different binding signatures depending on the type of G4. Moreover, we show that this Ag-based carbene binds the i-motif DNA structure. Alternatively, its Au(I) counterpart, which was investigated for comparison, stabilises mitochondrial G4. Theoretical in silico studies elucidated the details of different binding modes depending on the geometry of G4. The two complexes showed increased cytotoxic activity compared to cisplatin, overcoming its resistance in ovarian cancer. The binding of these new drug candidates with other relevant biosubstrates was studied to afford a more complete picture of their possible targets. In particular, the Ag(I) complex preferentially binds DNA structures over RNA structures, with higher binding constants for the non-canonical nucleic acids with respect to natural calf thymus DNA. Regarding possible protein targets, its interaction with the albumin model protein BSA was also tested.

Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations

Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.

Molecular dynamics simulations reveal the parallel stranded d(GGGA)3GGG DNA quadruplex folds via multiple paths from a coil-like ensemble

G-quadruplexes (G4s) are non-canonical nucleic acid structures that fold through complex processes. Characterization of the G4 folding landscape may help to elucidate biological roles of G4s but is challenging both experimentally and computationally. Here, we achieved complete folding of a three-quartet parallel DNA G4 with (GGGA)3GGG sequence using all-atom explicit-solvent enhanced-sampling molecular dynamics (MD) simulations. The simulations suggested early formation of guanine stacks in the G-tracts, which behave as semi-rigid blocks in the folding process. The folding continues via the formation of a collapsed compact coil-like ensemble. Structuring of the G4 from the coil then proceeds via various cross-like, hairpin, slip-stranded and two-quartet ensembles and can bypass the G-triplex structure. Folding of the parallel G4 does not appear to involve any salient intermediates and is a multi-pathway process. We also carried out an extended set of simulations of parallel G-hairpins. While parallel G-hairpins are extremely unstable when isolated, they are more stable inside the coil structure. On the methodology side, we show that the AMBER DNA force field predicts the folded G4 to be less stable than the unfolded ensemble, uncovering substantial force-field issues. Overall, we provide unique atomistic insights into the folding landscape of parallel-stranded G4 but also reveal limitations of current state-of-the-art MD techniques.

Computational Probing of the Folding Mechanism of Human Telomeric G-Quadruplex DNA

The human telomeric (htel) sequences in the terminal regions of human telomeres form diverse G-quadruplex (GQ) structures. Despite much experimental efforts to elucidate the folding pathways of htel GQ, no comprehensive model of htel GQ folding has been presented. Here, we describe folding pathways of the htel GQ determined by state-of-the-art enhanced sampling molecular dynamics simulation at the all-atom level. Briefly, GQ folding is initiated by the formation of a single-hairpin and then followed by the formation of double-hairpins, which then branch via distinct folding pathways to produce different GQ topologies (antiparallel chair, antiparallel basket, hybrids 1 and 2, and parallel propeller). In addition to these double-hairpin states, three-triad and two-tetrad structures in antiparallel backbone alignment serve as key intermediates that connect the GQ folding and transition between two different GQs.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all-anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.

Inhibited complete folding of consecutive human telomeric G-quadruplexes

Noncanonical DNA structures, termed G-quadruplexes, are present in human genomic DNA and are important elements in many DNA metabolic processes. Multiple sites in the human genome have G-rich DNA stretches able to support formation of several consecutive G-quadruplexes. One of those sites is the telomeric overhang region that has multiple repeats of TTAGGG and is tightly associated with both cancer and aging. We investigated the folding of consecutive G-quadruplexes in both potassium- and sodium-containing solutions using single-molecule FRET spectroscopy, circular dichroism, thermal melting and molecular dynamics simulations. Our observations show coexistence of partially and fully folded DNA, the latter consisting of consecutive G-quadruplexes. Following the folding process over hours in sodium-containing buffers revealed fast G-quadruplex folding but slow establishment of thermodynamic equilibrium. We find that full consecutive G-quadruplex formation is inhibited by the many DNA structures randomly nucleating on the DNA, some of which are off-path conformations that need to unfold to allow full folding. Our study allows describing consecutive G-quadruplex formation in both nonequilibrium and equilibrium conditions by a unified picture, where, due to the many possible DNA conformations, full folding with consecutive G-quadruplexes as beads on a string is not necessarily achieved.

Dynamic interaction of BRCA2 with telomeric G-quadruplexes underlies telomere replication homeostasis

BRCA2-deficient cells precipitate telomere shortening upon collapse of stalled replication forks. Here, we report that the dynamic interaction between BRCA2 and telomeric G-quadruplex (G4), the non-canonical four-stranded secondary structure, underlies telomere replication homeostasis. We find that the OB-folds of BRCA2 binds to telomeric G4, which can be an obstacle during replication. We further demonstrate that BRCA2 associates with G-triplex (G3)-derived intermediates, which are likely to form during direct interconversion between parallel and non-parallel G4. Intriguingly, BRCA2 binding to G3 intermediates promoted RAD51 recruitment to the telomere G4. Furthermore, MRE11 resected G4-telomere, which was inhibited by BRCA2. Pathogenic mutations at the OB-folds abrogated the binding with telomere G4, indicating that the way BRCA2 associates with telomere is innate to its tumor suppressor activity. Collectively, we propose that BRCA2 binding to telomeric G4 remodels it and allows RAD51-mediated restart of the G4-driven replication fork stalling, simultaneously preventing MRE11-mediated breakdown of telomere.

G-quadruplex (G4) can be formed in telomeric DNA. Here the authors show that BRCA2 interacts with telomere G4 structure generated during telomere replication, protecting telomere from nuclease attack.

Studying the Dynamics of a Complex G-Quadruplex System: Insights into the Comparison of MD and NMR Data

Molecular dynamics (MD) simulations are coming of age in the study of nucleic acids, including specific tertiary structures such as G-quadruplexes. While being precious for providing structural and dynamic information inaccessible to experiments at the atomistic level of resolution, MD simulations in this field may still be limited by several factors. These include the force fields used, different models for ion parameters, ionic strengths, and water models. We address various aspects of this problem by analyzing and comparing microsecond-long atomistic simulations of the G-quadruplex structure formed by the human immunodeficiency virus long terminal repeat (HIV LTR)-III sequence for which nuclear magnetic resonance (NMR) structures are available. The system is studied in different conditions, systematically varying the ionic strengths, ion numbers, and water models. We comparatively analyze the dynamic behavior of the G-quadruplex motif in various conditions and assess the ability of each simulation to satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints and structural parameters. The conditions taking into account K⁺-ions to neutralize the system charge, mimicking the intracellular ionic strength, and using the four-atom water model are found to be the best in reproducing the experimental NMR constraints and data. Our analysis also reveals that in all of the simulated environments residues belonging to the duplex moiety of HIV LTR-III exhibit the highest flexibility.

Thermally Induced Transitions of d(G4T4G3) Quadruplexes Can Be Described as Kinetically Driven Processes

DNA sequences that are rich in guanines and can form four-stranded structures are called G-quadruplexes. Due to the growing evidence that they may play an important role in several key biological processes, the G-quadruplexes have captured the interest of several researchers. G-quadruplexes may form in the presence of different metal cations as polymorphic structures formed in kinetically governed processes. Here we investigate a complex polymorphism of d(G4T4G3) quadruplexes at different K+ concentrations. We show that population size of different d(G4T4G3) quadruplex conformations can be manipulated by cooling rate and/or K+ concentration. We use a kinetic model to describe data obtained from DSC, CD and UV spectroscopy and PAGE experiments. Our model is able to describe the observed thermally induced conformational transitions of d(G4T4G3) quadruplexes at different K+ concentrations.

Folding dynamics of polymorphic G-quadruplex structures

G-quadruplexes (G4), found in numerous places within the human genome, are involved in essential processes of cell regulation. Chromosomal DNA G4s are involved for example, in replication and transcription as first steps of gene expression. Hence, they influence a plethora of downstream processes. G4s possess an intricate structure that differs from canonical B-form DNA. Identical DNA G4 sequences can adopt multiple long-lived conformations, a phenomenon known as G4 polymorphism. A detailed understanding of the molecular mechanisms that drive G4 folding is essential to understand their ambivalent regulatory roles. Disentangling the inherent dynamic and polymorphic nature of G4 structures thus is key to unravel their biological functions and make them amenable as molecular targets in novel therapeutic approaches. We here review recent experimental approaches to monitor G4 folding and discuss structural aspects for possible folding pathways. Substantial progress in the understanding of G4 folding within the recent years now allows drawing comprehensive models of the complex folding energy landscape of G4s that we herein evaluate based on computational and experimental evidence.

Mechanical diversity and folding intermediates of parallel-stranded G-quadruplexes with a bulge

A significant number of sequences in the human genome form noncanonical G-quadruplexes (G4s) with bulges or a guanine vacancy. Here, we systematically characterized the mechanical stability of parallel-stranded G4s with a one to seven nucleotides bulge at various positions. Our results show that G4-forming sequences with a bulge form multiple conformations, including fully-folded G4 with high mechanical stability (unfolding forces > 40 pN), partially-folded intermediates (unfolding forces < 40 pN). The folding probability and folded populations strongly depend on the positions and lengths of the bulge. By combining a single-molecule unfolding assay, dimethyl sulfate (DMS) footprinting, and a guanine-peptide conjugate that selectively stabilizes guanine-vacancy-bearing G-quadruplexes (GVBQs), we identified that GVBQs are the major intermediates of G4s with a bulge near the 5′ or 3′ ends. The existence of multiple structures may induce different regulatory functions in many biological processes. This study also demonstrates a new strategy for selectively stabilizing the intermediates of bulged G4s to modulate their functions.

Folding Dynamics of DNA G-Quadruplexes Probed by Millisecond Temperature Jump Circular Dichroism

G-quadruplexes play important roles in cellular regulatory functions, but despite significant experimental and theoretical efforts, their folding mechanisms remain poorly understood. In this context, we developed a T-jump experiment to access the thermal denaturation and renaturation dynamics of short intramolecular G-quadruplexes in vitro, on the time scale of a few hundred milliseconds. With this new setup, we compared the thermal denaturation and renaturation kinetics of three antiparallel topologies made of the human telomeric sequences d[(5'-GGG(TTAGGG)₃-3']/Na⁺ and d[5'-AGGG(TTAGGG)₃-3']/Na⁺ and the thrombin-binding aptamer sequence d[5'-GGTTGGTGTGGTTGG-3']/K⁺, with those of the parallel topology made of the human CEB25 minisatellite d[5'-AAGGGTGGGTGTAAGTGTGGGTGGGT-3']/Na⁺. In all cases, exponential kinetics of the order of several hundred milliseconds were observed. Measurements performed for different initial temperatures revealed distinct denaturation and renaturation dynamics, ruling out a simple two-state mechanism. The parallel topology, in which all guanines adopt an anti conformation, displays much slower dynamics than antiparallel topologies associated with very low activation barriers. This behavior can be explained by the constrained conformational space due to the presence of the single-base propeller loops that likely hinders the movement of the coiled DNA strand and reduces the contribution of the entropy during the renaturation process at high temperatures.

Single-Molecular Förster Resonance Energy Transfer Measurement on Structures and Interactions of Biomolecules

Single-molecule Förster resonance energy transfer (smFRET) inherits the strategy of measurement from the effective "spectroscopic ruler" FRET and can be utilized to observe molecular behaviors with relatively high throughput at nanometer scale. The simplicity in principle and configuration of smFRET make it easy to apply and couple with other technologies to comprehensively understand single-molecule dynamics in various application scenarios. Despite its widespread application, smFRET is continuously developing and novel studies based on the advanced platforms have been done. Here, we summarize some representative examples of smFRET research of recent years to exhibit the versatility and note typical strategies to further improve the performance of smFRET measurement on different biomolecules.

Insights into G-Quadruplex-Hemin Dynamics Using Atomistic Simulations: Implications for Reactivity and Folding

Guanine quadruplex nucleic acids (G4s) are involved in key biological processes such as replication or transcription. Beyond their biological relevance, G4s find applications as biotechnological tools since they readily bind hemin and enhance its peroxidase activity, creating a G4-DNAzyme. The biocatalytic properties of G4-DNAzymes have been thoroughly studied and used for biosensing purposes. Despite hundreds of applications and massive experimental efforts, the atomistic details of the reaction mechanism remain unclear. To help select between the different hypotheses currently under investigation, we use extended explicit-solvent molecular dynamics (MD) simulations to scrutinize the G4/hemin interaction. We find that besides the dominant conformation in which hemin is stacked atop the external G-quartets, hemin can also transiently bind to the loops and be brought to the external G-quartets through diverse delivery mechanisms. The simulations do not support the catalytic mechanism relying on a wobbling guanine. Similarly, the catalytic role of the iron-bound water molecule is not in line with our results; however, given the simulation limitations, this observation should be considered with some caution. The simulations rather suggest tentative mechanisms in which the external G-quartet itself could be responsible for the unique H₂O₂-promoted biocatalytic properties of the G4/hemin complexes. Once stacked atop a terminal G-quartet, hemin rotates about its vertical axis while readily sampling shifted geometries where the iron transiently contacts oxygen atoms of the adjacent G-quartet. This dynamics is not apparent from the ensemble-averaged structure. We also visualize transient interactions between the stacked hemin and the G4 loops. Finally, we investigated interactions between hemin and on-pathway folding intermediates of the parallel-stranded G4 fold. The simulations suggest that hemin drives the folding of parallel-stranded G4s from slip-stranded intermediates, acting as a G4 chaperone. Limitations of the MD technique are briefly discussed.

Dynamics and Barrier of Movements of Sodium and Potassium Ions Across the Oxytricha nova G-Quadruplex Core

G-quadruplexes (GQs) are highly stable noncanonical forms of nucleic acids that are present in important genomic regions. The central core of the GQ is lined up by four closely spaced carbonyl groups from the G-quartets, and the resulting electrostatic repulsion is neutralized by the coordinating cations. In spite of several reports on GQ structure and cation-GQ interactions, the atomic- to molecular-level understanding of the ion dynamics and ion exchange in the GQ core is quite poor. Here, we attempt to elucidate the mechanism of Na⁺ and K⁺ binding to the GQ core and trace the exchange of these ions with the ions in bulk by means of all-atomic molecular dynamics (MD) simulations. One of the most studied GQs, Oxytricha nova telomeric G-quadruplex (OxyGQ), is taken as the representative GQ. Subsequently, umbrella sampling MD simulations were performed to elucidate the energetics of ion translocation from one end to the other end of the GQ central core. Our study highlights the importance of ion hydration for the uptake and correct positioning of the cations in the core. The free-energy landscape of ion transport has shown favorable in-plane binding of Na⁺ ions with GQ quartets, which matches very well with the crystal structure. The binding of K⁺ ions, on the other hand, was out-of-plane and its translocation required a larger barrier to cross.

G-Quadruplexes at Telomeres: Friend or Foe?

Telomeres are DNA-protein complexes that cap and protect the ends of linear chromosomes. In almost all species, telomeric DNA has a G/C strand bias, and the short tandem repeats of the G-rich strand have the capacity to form into secondary structures in vitro, such as four-stranded G-quadruplexes. This has long prompted speculation that G-quadruplexes play a positive role in telomere biology, resulting in selection for G-rich tandem telomere repeats during evolution. There is some evidence that G-quadruplexes at telomeres may play a protective capping role, at least in yeast, and that they may positively affect telomere maintenance by either the enzyme telomerase or by recombination-based mechanisms. On the other hand, G-quadruplex formation in telomeric DNA, as elsewhere in the genome, can form an impediment to DNA replication and a source of genome instability. This review summarizes recent evidence for the in vivo existence of G-quadruplexes at telomeres, with a focus on human telomeres, and highlights some of the many unanswered questions regarding the location, form, and functions of these structures.

A mechanism for the extension and unfolding of parallel telomeric G-quadruplexes by human telomerase at single-molecule resolution

Telomeric G-quadruplexes (G4) were long believed to form a protective structure at telomeres, preventing their extension by the ribonucleoprotein telomerase. Contrary to this belief, we have previously demonstrated that parallel-stranded conformations of telomeric G4 can be extended by human and ciliate telomerase. However, a mechanistic understanding of the interaction of telomerase with structured DNA remained elusive. Here, we use single-molecule fluorescence resonance energy transfer (smFRET) microscopy and bulk-phase enzymology to propose a mechanism for the resolution and extension of parallel G4 by telomerase. Binding is initiated by the RNA template of telomerase interacting with the G-quadruplex; nucleotide addition then proceeds to the end of the RNA template. It is only through the large conformational change of translocation following synthesis that the G-quadruplex structure is completely unfolded to a linear product. Surprisingly, parallel G4 stabilization with either small molecule ligands or by chemical modification does not always inhibit G4 unfolding and extension by telomerase. These data reveal that telomerase is a parallel G-quadruplex resolvase.

Quantifying the impact of small molecule ligands on G-quadruplex stability against Bloom helicase

G-quadruplex (GQ) stabilizing small molecule (SM) ligands have been used to stabilize human telomeric GQ (hGQ) to inhibit telomerase activity, or non-telomeric GQs to manipulate gene expression at transcription or translation level. GQs are known to inhibit DNA replication unless destabilized by helicases, such as Bloom helicase (BLM). Even though the impact of SM ligands on thermal stability of GQs is commonly used to characterize their efficacy, how these ligands influence helicase-mediated GQ unfolding is not well understood. Three prominent SM ligands (an oxazole telomestatin derivative, pyridostatin, and PhenDC₃), which thermally stabilize hGQ at different levels, were utilized in this study. How these ligands influence BLM-mediated hGQ unfolding was investigated using two independent single-molecule approaches. While the frequency of dynamic hGQ unfolding events was used as the metric in the first approach, the second approach was based on quantifying the cumulative unfolding activity as a function of time. All three SM ligands inhibited BLM activity at similar levels, 2–3 fold, in both approaches. Our observations suggest that the impact of SM ligands on GQ thermal stability is not an ideal predictor for their inhibition of helicase-mediated unfolding, which is physiologically more relevant.

Parallel G-triplexes and G-hairpins as potential transitory ensembles in the folding of parallel-stranded DNA G-Quadruplexes

Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.

The Presence and Localization of G-Quadruplex Forming Sequences in the Domain of Bacteria

The role of local DNA structures in the regulation of basic cellular processes is an emerging field of research. Amongst local non-B DNA structures, the significance of G-quadruplexes was demonstrated in the last decade, and their presence and functional relevance has been demonstrated in many genomes, including humans. In this study, we analyzed the presence and locations of G-quadruplex-forming sequences by G4Hunter in all complete bacterial genomes available in the NCBI database. G-quadruplex-forming sequences were identified in all species, however the frequency differed significantly across evolutionary groups. The highest frequency of G-quadruplex forming sequences was detected in the subgroup Deinococcus-Thermus, and the lowest frequency in Thermotogae. G-quadruplex forming sequences are non-randomly distributed and are favored in various evolutionary groups. G-quadruplex-forming sequences are enriched in ncRNA segments followed by mRNAs. Analyses of surrounding sequences showed G-quadruplex-forming sequences around tRNA and regulatory sequences. These data point to the unique and non-random localization of G-quadruplex-forming sequences in bacterial genomes.

Extreme mechanical diversity of human telomeric DNA revealed by fluorescence-force spectroscopy

G-quadruplexes (GQs) can adopt diverse structures and are functionally implicated in transcription, replication, translation, and maintenance of telomere. Their conformational diversity under physiological levels of mechanical stress, however, is poorly understood. We used single-molecule fluorescence-force spectroscopy that combines fluorescence resonance energy transfer with optical tweezers to measure human telomeric sequences under tension. Abrupt GQ unfolding with K⁺ in solution occurred at as many as four discrete levels of force. Added to an ultrastable state and a gradually unfolding state, there were six mechanically distinct structures. Extreme mechanical diversity was also observed with Na⁺, although GQs were mechanically weaker. Our ability to detect small conformational changes at low forces enabled the determination of refolding forces of about 2 pN. Refolding was rapid and stochastically redistributed molecules to mechanically distinct states. A single guanine-to-thymine substitution mutant required much higher ion concentrations to display GQ-like unfolding and refolded via intermediates, contrary to the wild type. Contradicting an earlier proposal, truncation to three hexanucleotide repeats resulted in a single-stranded DNA-like mechanical behavior under all conditions, indicating that at least four repeats are required to form mechanically stable structures.

Structural dynamics of propeller loop: towards folding of RNA G-quadruplex

We have carried out an extended set of standard and enhanced-sampling MD simulations (for a cumulative simulation time of 620 μs) with the aim to study folding landscapes of the rGGGUUAGGG and rGGGAGGG parallel G-hairpins (PH) with propeller loop. We identify folding and unfolding pathways of the PH, which is bridged with the unfolded state via an ensemble of cross-like structures (CS) possessing mutually tilted or perpendicular G-strands interacting via guanine-guanine H-bonding. The oligonucleotides reach the PH conformation from the unfolded state via a conformational diffusion through the folding landscape, i.e. as a series of rearrangements of the H-bond interactions starting from compacted anti-parallel hairpin-like structures. Although isolated PHs do not appear to be thermodynamically stable we suggest that CS and PH-types of structures are sufficiently populated during RNA guanine quadruplex (GQ) folding within the context of complete GQ-forming sequences. These structures may participate in compact coil-like ensembles that involve all four G-strands and already some bound ions. Such ensembles can then rearrange into the fully folded parallel GQs via conformational diffusion. We propose that the basic atomistic folding mechanism of propeller loops suggested in this work may be common for their formation in RNA and DNA GQs.

Insilico direct folding of thrombin-binding aptamer G-quadruplex at all-atom level

The reversible folding of the thrombin-binding DNA aptamer G-quadruplexes (GQs) (TBA-15) starting from fully unfolded states was demonstrated using a prolonged time scale (10-12 μs) parallel tempering metadynamics (PTMetaD) simulation method in conjunction with a modified version of the AMBER bsc1 force field. For unbiased descriptions of the folding free energy landscape of TBA-15, this force field was minimally modified. From this direct folding simulation using the modified bsc1 force field, reasonably converged free energy landscapes were obtained in K+-rich aqueous solution (150 mM), providing detailed atomistic pictures of GQ folding mechanisms for TBA-15. This study found that the TBA folding occurred via multiple folding pathways with two major free energy barriers of 13 and 15 kcal/mol in the presence of several intermediate states of G-triplex variants. The early formation of these intermediates was associated with a single K+ ion capturing. Interestingly, these intermediate states appear to undergo facile transitions among themselves through relatively small energy barriers.

A Comparative Study of G-Quadruplex Unfolding and DNA Reeling Activities of Human RECQ5 Helicase

RECQ5 is one of five members of the RecQ family of helicases in humans, which include RECQ1, Bloom (BLM), Werner (WRN), RECQ4, and RECQ5. Both WRN and BLM have been shown to resolve G-quadruplex (GQ) structures. Deficiencies in unfolding GQ are known to result in DNA breaks and genomic instability, which are prominent in Werner and Bloom syndromes. RECQ5 is significant in suppressing sister chromatid exchanges during homologous recombination but its GQ unfolding activity are not known. We performed single-molecule studies under different salt (50-150 mM KCl or NaCl) and ATP concentrations on different GQ constructs including human telomeric GQ (with different overhangs and polarities) and GQ formed by thrombin-binding aptamer to investigate this activity. These studies demonstrated a RECQ5-mediated GQ unfolding activity that was an order of magnitude weaker than BLM and WRN. On the other hand, BLM and RECQ5 demonstrated similar single-stranded DNA (ssDNA) reeling activities that were not coupled to GQ unfolding. These results demonstrate overlap in function between these RecQ helicases; however, the relatively weak GQ destabilization activity of RECQ5 compared to BLM and WRN suggests that RECQ5 is not primarily associated with GQ destabilization, but could substitute for the more efficient helicases under conditions where their activity is compromised. In addition, these results implicate a more general role for helicase-promoted ssDNA reeling activity such as removal of proteins at the replication fork, whereas the association of ssDNA reeling with GQ destabilization is more helicase-specific.

Single-Molecule Studies of Telomeres and Telomerase

Telomeres are specialized chromatin structures that protect chromosome ends from dangerous processing events. In most tissues, telomeres shorten with each round of cell division, placing a finite limit on cell growth. In rapidly dividing cells, including the majority of human cancers, cells bypass this growth limit through telomerase-catalyzed maintenance of telomere length. The dynamic properties of telomeres and telomerase render them difficult to study using ensemble biochemical and structural techniques. This review describes single-molecule approaches to studying how individual components of telomeres and telomerase contribute to function. Single-molecule methods provide a window into the complex nature of telomeres and telomerase by permitting researchers to directly visualize and manipulate the individual protein, DNA, and RNA molecules required for telomere function. The work reviewed in this article highlights how single-molecule techniques have been utilized to investigate the function of telomeres and telomerase.

Structure, Properties, and Biological Relevance of the DNA and RNA G-Quadruplexes: Overview 50 Years after Their Discovery

G-quadruplexes (G4s), which are known to have important roles in regulation of key biological processes in both normal and pathological cells, are the most actively studied non-canonical structures of nucleic acids. In this review, we summarize the results of studies published in recent years that change significantly scientific views on various aspects of our understanding of quadruplexes. Modern notions on the polymorphism of DNA quadruplexes, on factors affecting thermodynamics and kinetics of G4 folding–unfolding, on structural organization of multiquadruplex systems, and on conformational features of RNA G4s and hybrid DNA–RNA G4s are discussed. Here we report the data on location of G4 sequence motifs in the genomes of eukaryotes, bacteria, and viruses, characterize G4-specific small-molecule ligands and proteins, as well as the mechanisms of their interactions with quadruplexes. New information on the structure and stability of G4s in telomeric DNA and oncogene promoters is discussed as well as proof being provided on the occurrence of G-quadruplexes in cells. Prominence is given to novel experimental techniques (single molecule manipulations, optical and magnetic tweezers, original chemical approaches, G4 detection in situ, in-cell NMR spectroscopy) that facilitate breakthroughs in the investigation of the structure and functions of G-quadruplexes.

A practical guide to studying G-quadruplex structures using single-molecule FRET

In this article, we summarize the knowledge and best practices learned from bulk and single-molecule measurements to address some of the frequently experienced difficulties in single-molecule Förster resonance energy transfer (smFRET) measurements on G-quadruplex (GQ) structures. The number of studies that use smFRET to investigate the structure, function, dynamics, and interactions of GQ structures has grown significantly in the last few years, with new applications already in sight. However, a number of challenges need to be overcome before reliable and reproducible smFRET data can be obtained in measurements that include GQ. The annealing and storage conditions, the location of fluorophores on the DNA construct, and the ionic conditions of the experiment are some of the factors that are of critical importance for the outcome of measurements, and many of these manifest themselves in unique ways in smFRET assays. By reviewing these aspects and providing a summary of best practices, we aim to provide a practical guide that will help in successfully designing and performing smFRET studies on GQ structures.

Folding of guanine quadruplex molecules-funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies

BACKGROUND

Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding.

SCOPE OF REVIEW

We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles.

MAJOR CONCLUSIONS

The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables.

GENERAL SIGNIFICANCE

We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

A direct view of the complex multi-pathway folding of telomeric G-quadruplexes

G-quadruplexes (G4s) are DNA secondary structures that are capable of forming and function in vivo. The propensity of G4s to exhibit extreme polymorphism and complex dynamics is likely to influence their cellular function, yet a clear microscopic picture of their folding process is lacking. Here we employed single-molecule FRET microscopy to obtain a direct view of the folding and underlying conformational dynamics of G4s formed by the human telomeric sequence in potassium containing solutions. Our experiments allowed detecting several folded states that are populated in the course of G4 folding and determining their folding energetics and timescales. Combining the single-molecule data with molecular dynamics simulations enabled obtaining a structural description of the experimentally observed folded states. Our work thus provides a comprehensive thermodynamic and kinetic description of the folding of G4s that proceeds through a complex multi-route pathway, involving several marginally stable conformational states.

Folding and misfolding pathways of G-quadruplex DNA

G-quadruplexes adopt various folding topologies, but information on their folding pathways remains scarce. Here, we used electrospray mass spectrometry to detect and quantify the specifically bound potassium ions, and circular dichroism to characterize the stacking topology of each ensemble. For human telomeric (hTel) sequences containing the d((GGGTTA)₃GGG) core, K⁺ binding affinity and cooperativity strongly depends on the chosen construct. The shortest sequences bind only one K⁺ at low KCl concentration, and this 2-quartet G-quadruplex is antiparallel. Flanking bases increase the K⁺ binding cooperativity. To decipher the folding pathways, we investigated the kinetics of K⁺ binding to telomeric (hybrid) and c-myc (parallel) G-quadruplexes. G-quadruplexes fold via branched pathways with multiple parallel reactions. Up to six states (one ensemble without K⁺, two ensembles with 1-K⁺ and three ensembles with 2-K⁺) are separated based on their formation rates and ion mobility spectrometry. All G-quadruplexes first form long-lived misfolded structures (off-pathway compared to the most stable structures) containing one K⁺ and two quartets in an antiparallel stacking arrangement. The results highlight the particular ruggedness of G-quadruplex nucleic acid folding landscapes. Misfolded structures can play important roles for designing artificial G-quadruplex based structures, and for conformational selection by ligands or proteins in a biological context.

Folding Kinetics of Single Human Telomeric G-Quadruplex Affected by Cisplatin

G-Quadruplex DNA structure has been proven to be a binding target for small molecular organic compounds and hence regarded as a promising pharmacological target. Cisplatin is a widely used chemotherapy drug that targets duplex DNA and was recently shown to react also with G-quadruplex, implying that cisplatin actually may also target G-quadruplex. In this work, we employed magnetic tweezers to investigate the influence of cisplatin on the folding kinetics of single human telomeric G-quadruplex. It was revealed that cisplatin and G-quadruplex interact in two different and competitive ways that depend on cisplatin concentration.

G-quadruplex recognition and remodeling by the FANCJ helicase

Guanine rich nucleic acid sequences can form G-quadruplex (G4) structures that interfere with DNA replication, repair and RNA transcription. The human FANCJ helicase contributes to maintaining genomic integrity by promoting DNA replication through G4-forming DNA regions. Here, we combined single-molecule and ensemble biochemical analysis to show that FANCJ possesses a G4-specific recognition site. Through this interaction, FANCJ targets G4-containing DNA where its helicase and G4-binding activities enable repeated rounds of stepwise G4-unfolding and refolding. In contrast to other G4-remodeling enzymes, FANCJ partially stabilizes the G-quadruplex. This would preserve the substrate for the REV1 translesion DNA synthesis polymerase to incorporate cytosine across from a replication-stalling G-quadruplex. The residues responsible for G-quadruplex recognition also participate in interaction with MLH1 mismatch-repair protein, suggesting that the FANCJ activity supporting replication and its participation in DNA interstrand crosslink repair and/or heteroduplex rejection are mutually exclusive. Our findings not only describe the mechanism by which FANCJ recognizes G-quadruplexes and mediates their stepwise unfolding, but also explain how FANCJ chooses between supporting DNA repair versus promoting DNA replication through G-rich sequences.

The lighthouse at the end of the chromosome

Fluorescence microscopy can be used to assess the dynamic localization and intensity of single entities in vitro or in living cells. It has been applied with aplomb to many different cellular processes and has significantly enlightened our understanding of the heterogeneity and complexity of biological systems. Recently, high-resolution fluorescence microscopy has been brought to bear on telomeres, leading to new insights into telomere spatial organization and accessibility, and into the mechanistic nuances of telomere elongation. We provide a snapshot of some of these recent advances with a focus on mammalian systems, and show how three-dimensional, time-lapse microscopy and single-molecule fluorescence shine a new light on the end of the chromosome.

Folding dynamics and conformational heterogeneity of human telomeric G-quadruplex structures in Na+ solutions by single molecule FRET microscopy

G-quadruplex structures can occur throughout the genome, including at telomeres. They are involved in cellular regulation and are potential drug targets. Human telomeric G-quadruplex structures can fold into a number of different conformations and show large conformational diversity. To elucidate the different G-quadruplex conformations and their dynamics, we investigated telomeric G-quadruplex folding using single molecule FRET microscopy in conditions where it was previously believed to yield low structural heterogeneity. We observed four FRET states in Na⁺ buffers: an unfolded state and three G-quadruplex related states that can interconvert between each other. Several of these states were almost equally populated at low to medium salt concentrations. These observations appear surprising as previous studies reported primarily one G-quadruplex conformation in Na⁺ buffers. Our results permit, through the analysis of the dynamics of the different observed states, the identification of a more stable G-quadruplex conformation and two transient G-quadruplex states. Importantly these results offer a unique view into G-quadruplex topological heterogeneity and conformational dynamics.

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.

Construction of anti-parallel G-quadruplexes through sequential templated click

Biologically relevant anti-parallel DNA G-quadruplexes were constrained and stabilised onto addressable cyclopeptide scaffolds through sequential oxime and CuAAc reactions.

Targeting unimolecular G-quadruplex nucleic acids: a new paradigm for the drug discovery?

INTRODUCTION

G-quadruplexes (G4s) are targets of great interest because of their roles in crucial biological processes, such as aging and cancer. G4s are based on the formation of G-quartets, stabilised by Hoogsteen-type hydrogen bonds and by interaction with cations between the tetrads. These biologically relevant conformations were first discovered in eukaryotic chromosomal telomeric DNA, but have also been found in the proximal location of promoters in a number of human genes. Therefore, the extensive analysis of an intriguing target could move towards the rational drug design of new selective anticancer agents.

AREAS COVERED

The authors review G4 structural characterisation, with detailed insight related to the polymorphism issue. The authors describe the topologically distinct G4 structural forms and the factors involved in their interconversion mechanisms, such as the sequence of the oligonucleotides, the strand stoichiometry and orientation, the syn-anti conformation of the guanine glycosidic bonds and the G4 loop types and the environmental factors. Furthermore, the authors report several studies related to folding and unfolding kinetic profiles in order to understand the conformational view of monomolecular G4 formations.

EXPERT OPINION

G4 unimolecular nucleic acids can be considered as valid targets for the rational drug development of novel anticancer agents. Structural biology represents an essential link between the biology and medicinal chemistry knowledge in this field. In silico methods have already been demonstrated to be useful, especially if well integrated with biophysical tests. If this proves successful, the G4-targeting paradigm could also be extended to drug discovery beyond neoplastic pathologies.

Dynamics and stability of polymorphic human telomeric G-quadruplex under tension

As critical DNA structures capping the human chromosome ends, the stability and structural polymorphism of human telomeric G-quadruplex (G4) have drawn increasing attention in recent years. This work characterizes the equilibrium transitions of single-molecule telomeric G4 at physiological K⁺ concentration. We report three folded states of telomeric G4 with markedly different lifetime and mechanical stability. Our results show that the kinetically favored folding pathway is through a short-lived intermediate state to a longer-lived state. By examining the force dependence of transition rates, the force-dependent transition free energy landscape for this pathway is determined. In addition, an ultra-long-lived form of telomeric G4 structure with a much stronger mechanical stability is identified.

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.