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

Computer Folding of Parallel DNA G-Quadruplex: Hitchhiker's Guide to the Conformational Space

Guanine quadruplexes (GQs) play crucial roles in various biological processes, and understanding their folding pathways provides insight into their stability, dynamics, and functions. This knowledge aids in designing therapeutic strategies, as GQs are potential targets for anticancer drugs and other therapeutics. Although experimental and theoretical techniques have provided valuable insights into different stages of the GQ folding, the structural complexity of GQs poses significant challenges, and our understanding remains incomplete. This study introduces a novel computational protocol for folding an entire GQ from single-strand conformation to its native state. By combining two complementary enhanced sampling techniques, we were able to model folding pathways, encompassing a diverse range of intermediates. Although our investigation of the GQ free energy surface (FES) is focused solely on the folding of the all-anti parallel GQ topology, this protocol has the potential to be adapted for the folding of systems with more complex folding landscapes.

Probing G-quadruplex-ligand binding using DNA intrinsic fluorescence

G-quadruplexes (G4s) are helical four-stranded nucleic acid structures that can form in guanine-rich sequences, which are mostly found in functional cellular regions, such as telomeres, promoters, and DNA replication origins. Great efforts are being made to target these structures towards the development of specific small molecule G4 binders for novel anti-cancer, neurological, and viral therapies. Here, we introduce an optical assay based on quenching of the intrinsic fluorescence of DNA G-quadruplexes for assessing and comparing the G4 binding affinity of various small molecule ligands in solutions. We show that the approach allows direct quantification of ligand binding to distinctive G4 topologies. We believe that this method will facilitate quick and reliable evaluation of small molecule G4 ligands and support their development.

Antiparallel G-Quadruplex Formation Hinders Conversion to a Parallel Topology

G-quadruplexes (G4s) are four-stranded structures formed by guanine-rich sequences. While their structures, properties, and applications have been extensively studied, an understanding of their folding processes remains limited. In this study, we investigated the folding of the sequence d[(G3T2)3G3] in potassium solutions, focusing on the impact of a folding intermediate on the overall folding process. Our results indicate that this sequence eventually folds into a parallel G4 structure, either directly or through an antiparallel conformation intermediate, suggesting the existence of a specific competitive folding process. Detailed kinetic analysis using stopped-flow techniques reveals that the antiparallel conformation forms much faster than the parallel one. This antiparallel G4 slowly converts to the thermodynamically favored parallel topology, thus slowing the overall folding rate. As a result, the formation of the parallel quadruplex via an antiparallel G4 intermediate is slower than the direct process, indicating that this antiparallel conformation negatively impacts the overall folding process in a temperature-dependent manner. Interestingly, sodium was shown to facilitate the conversion from antiparallel to parallel. These analyses highlight the complexity of the G4 folding process, which is crucial for most biological applications.

Early Events in G-quadruplex Folding Captured by Time-Resolved Small-Angle X-Ray Scattering

Time-resolved small-angle X-ray experiments (TR-SAXS) are reported here that capture and quantify a previously unknown rapid collapse of the unfolded oligonucleotide as an early step in G4 folding of hybrid 1 and hybrid 2 telomeric G-quadruplex structures. The rapid collapse, initiated by a pH jump, is characterized by an exponential decrease in the radius of gyration from 20.6 to 12.6 Å. The collapse is monophasic and is complete in less than 600 ms. Additional hand-mixing pH-jump kinetic studies show that slower kinetic steps follow the collapse. The folded and unfolded states at equilibrium were further characterized by SAXS studies and other biophysical tools, to show that G4 unfolding was complete at alkaline pH, but not in LiCl solution as is often claimed. The SAXS Ensemble Optimization Method (EOM) analysis reveals models of the unfolded state as a dynamic ensemble of flexible oligonucleotide chains with a variety of transient hairpin structures. These results suggest a G4 folding pathway in which a rapid collapse, analogous to molten globule formation seen in proteins, is followed by a confined conformational search within the collapsed particle to form the native contacts ultimately found in the stable folded form.

Probing Conical Intersection in the Multipathway Isomerization of CH3Cl Using Coulomb Explosion

Ubiquitous ultrafast isomerization is paramount in photoexcited molecules, in which non-adiabatic coupling among multiple electronic states can occur. We use the pump-probe Coulomb explosion imaging method to study the isomerization of CH3Cl molecules. We find that the isomerization under our strong field pump-probe scheme proceeds along multiple pathways, which are encoded in several distinct branches of the time-resolved kinetic energy release spectra for the CH2++HCl+ Coulomb explosion channel. Apart from the isomerized dissociative pathway in neutral and cationic excited states, the pump laser can also induce coherent vibrational dynamics in two coupled intermediate states and set up the initial conditions for the two concurrently proceeding isomerization pathways. The isomerization of CH3Cl provides an intriguing example of a chemical reaction consisting of multiple pathways and non-adiabatic dynamics.

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.

A TeZla micromixer for interrogating the early and broad folding landscape of G-quadruplex via multistage velocity descending

Biomacromolecular folding kinetics involves fast folding events and broad timescales. Current techniques face limitations in either the required time resolution or the observation window. In this study, we developed the TeZla micromixer, integrating Tesla and Zigzag microstructures with a multistage velocity descending strategy. TeZla achieves a significant short mixing dead time (40 µs) and a wide time window covering four orders of magnitude (up to 300 ms). Using this unique micromixer, we explored the folding landscape of c-Myc G4 and its noncanonical-G4 derivatives with different loop lengths or G-vacancy sites. Our findings revealed that c-Myc can bypass folding intermediates and directly adopt a G4 structure in the cation-deficient buffer. Moreover, we found that the loop length and specific G-vacancy site could affect the folding pathway and significantly slow down the folding rates. These results were also cross-validated with real-time NMR and circular dichroism. In conclusion, TeZla represents a versatile tool for studying biomolecular folding kinetics, and our findings may ultimately contribute to the design of drugs targeting G4 structures.

A General Framework to Interpret Hydrogen-Deuterium Exchange Native Mass Spectrometry of G-Quadruplex DNA

G-quadruplexes (G4s) are secondary structures formed by guanine-rich oligonucleotides involved in various biological processes. However, characterizing G4s is challenging, because of their structural polymorphism. Here, we establish how hydrogen-deuterium exchange native mass spectrometry (HDX/MS) can help to characterize the G4 structures and dynamics in solution. We correlated the time range of G4 exchange to the number of guanines involved in the inner and outer tetrads. We also established relationships among exchange rates, numbers of tetrads and bound cations, and stability. The use of HDX/native MS allows for the determination of tetrads formed and assessment of G4 stability at a constant temperature. A key finding is that stable G4s exchange through local fluctuations (EX2 exchange), whereas less stable G4s also undergo exchange through partial or complete unfolding (EX1 exchange). Deconvolution of the bimodal isotope distributions resulting from EX1 exchange provides valuable insight into the kinetics of folding and unfolding processes and allows one to detect and characterize transiently unfolded intermediates, even if scarcely populated. HDX/native MS thus represents a powerful tool for a more comprehensive exploration of the folding landscapes of G4s.

Dynamics of G-Quadruplex Formation under Molecular Crowding

G-quadruplex (G4) structures assemble from guanine-rich DNA sequences and are believed to regulate several key cellular processes. G4 formation and conformational interconversions are well-established to occur dynamically in vitro. However, a clear understanding of G4 formation dynamics in cells as well as under conditions mimicking the cellular environment is missing. To fill this gap, we have investigated the G4 dynamics in molecularly crowded solutions, thus replicating the effect of the excluded volume present in cells. The results show that the volume exclusion exerted by large crowding agents accelerates the rate of G4 formation by at least an order of magnitude, leading to significant G4 stabilization. Extrapolation from our experimental data predicts crowding-induced G4 stabilization by more than 3 kcal/mol, under crowding levels found in the cellular environment. Such effects are likely to be important for G4-driven regulatory functions.

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.

Detecting secondary structure formation with FRET-PAINT

We present single molecule studies demonstrating the capabilities of the FRET-PAINT method to detect secondary structures that would be challenging to detect with alternative methods, particularly single molecule FRET (smFRET). Instead of relying on the change in end-to-end separation as in smFRET, we use the change in accessibility to a small probe as the criterion for secondary structure formation and relative stability. As a model system, we study G-triplex formation by human telomeric repeat sequences in different structural contexts.

Kinetic and Thermodynamic Control of G-Quadruplex Polymorphism by Na+ and K+ Cations

G-Quadruplexes (G4s) are ubiquitous nucleic acid folding motifs that exhibit structural diversity that is dependent on cationic conditions. In this work, we exploit temperature-controlled single-molecule fluorescence resonance energy transfer (smFRET) to elucidate the kinetic and thermodynamic mechanisms by which monovalent cations (K+ and Na+) impact folding topologies for a simple G-quadruplex sequence (5'-GGG-(TAAGGG)3-3') with a three-state folding equilibrium. Kinetic measurements indicate that Na+ and K+ influence G4 formation in two distinctly different ways: the presence of Na+ modestly enhances an antiparallel G4 topology through an induced fit (IF) mechanism with a low affinity (Kd = 228 ± 26 mM), while K+ drives G4 into a parallel/hybrid topology via a conformational selection (CS) mechanism with much higher affinity (Kd = 1.9 ± 0.2 mM). Additionally, temperature-dependent studies of folding rate constants and equilibrium ratios reveal distinctly different thermodynamic driving forces behind G4 binding to K+ (ΔH°bind > 0, ΔS°bind > 0) versus Na+ (ΔH°bind < 0, ΔS°bind < 0), which further illuminates the diversity of the possible pathways for monovalent facilitation of G-quadruplex folding.

An overview of single-molecule techniques and applications in the study of nucleic acid structure and function

Nucleic acids are an indispensable component in all known life forms. The biological processes are regulated by Nucleic acids, which associate to form special high-order structures. since the high-level structures of nucleic acids are related to gene expression in cancer cells or viruses, it is very likely to become a potential drug target. Traditional biochemical methods are limited to distinguish the conformational distribution and dynamic transition process of single nucleic acid structure. The ligands based on the intermediate and transition states between different conformations are not designed by traditional biochemical methods. The single-molecule techniques enable real-time observation of the individual nucleic acid behavior due to its high resolution. Here, we introduce the application of single-molecule techniques in the study of small molecules to recognize nucleic acid structures, such as single-molecule FRET, magnetic tweezers, optical tweezers and atomic force microscopy. At the same time, we also introduce the specific advantages of single-molecule technology compared with traditional biochemical methods and some problems arisen in current research.

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.

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.

Insights into the free energy landscape and salt-controlled mechanism of the conformational conversions between human telomeric G-quadruplex structures

G-quadruplexes have become attractive drug targets in cancer therapy. However, due to the polymorphism of G-quadruplex structures, it is difficult to experimentally verify the relevant structures of multiple intermediates and transition states in dynamic equilibrium. Hence, understanding the mechanism by which structural conversions of G-quadruplexes occur is still challenging. We conducted targeted molecular dynamics simulation with umbrella sampling to investigate how salt affects the conformational conversion of human telomeric G-quadruplex. Our results explore a unique view into the structures and energy barrier of the intermediates and transition states in the interconversion process. The pathway of G-quadruplex conformational interconversion was mapped out by a free energy landscape, consisting of branched parallel pathways with multiple energy basins. We propose a salt-controlled mechanism that as the salt concentration increases, the conformational conversion mechanism switches from multi-pathway folding to sequential folding pathways. The hybrid-I and hybrid-II structures are intermediates in the basket-propeller transformation. In high-salt solutions, the conformational conversion upon K⁺ binding is more feasible than upon Na⁺ binding. The free energy barrier for conformational conversions ranges from 1.6 to 4.6 kcal/mol. Our work will be beneficial in developing anticancer agents.

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.

Thermodynamic Stability of G-Quadruplexes: Impact of Sequence and Environment

G-quadruplexes have attracted growing interest in recent years due to their occurrence in vivo and their possible biological functions. In addition to being promising targets for drug design, these four-stranded nucleic acid structures have also been recognized as versatile tools for various technological applications. Whereas a large number of studies have yielded insight into their remarkable structural diversity, our current knowledge on G-quadruplex stabilities as a function of sequence and environmental factors only gradually emerges with an expanding collection of thermodynamic data. This minireview provides an overview of general rules that may be used to better evaluate quadruplex thermodynamic stabilities but also discusses present challenges in predicting most stable folds for a given sequence and environment.

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.

Ligand Binding to Dynamically Populated G-Quadruplex DNA

Several small-molecule ligands specifically bind and stabilize G-quadruplex (G4) nucleic acid structures, which are considered to be promising therapeutic targets. G4s are polymorphic structures of varying stability, and their formation is dynamic. Here, we investigate the mechanisms of ligand binding to dynamically populated human telomere G4 DNA by using the bisquinolinium based ligand Phen-DC3 and a combination of single-molecule FRET microscopy, ensemble FRET and CD spectroscopies. Different cations are used to tune G4 polymorphism and folding dynamics. We find that ligand binding occurs to pre-folded G4 structures and that Phen-DC3 also induces G4 formation in unfolded single strands. Following ligand binding to dynamically populated G4s, the DNA undergoes pronounced conformational redistributions that do not involve direct ligand-induced G4 conformational interconversion. On the contrary, the redistribution is driven by ligand-induced G4 folding and trapping of dynamically populated short-lived conformation states. Thus, ligand-induced stabilization does not necessarily require the initial presence of stably folded G4s.

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.

Competition of Ligands and the 18-mer Binding Domain of the RHAU Helicase for G-Quadruplexes: Orthosteric or Allosteric Binding Mechanism?

Stabilizing the DNA and RNA structures known as G-quadruplexes (G4s) using specific ligands is a strategy that has been proposed to fight cancer. However, although G-quadruplex:ligand (G4:L) interactions have often been investigated, whether or not ligands are able to disrupt G-quadruplex:protein (G4:P) interactions remains poorly studied. In this study, using native mass spectrometry, we have investigated ternary G4:L:P complexes formed by G4s, some of the highest affinity ligands, and the binding domain of the RHAU helicase. Our results suggest that RHAU binds not only preferentially to parallel G4s, but also to free external G-quartets. We also found that, depending on the G4, ligands could prevent the binding of the peptide, either by direct competition for the binding sites (orthosteric inhibition) or by inducing conformational changes (allosteric inhibition). Notably, the ligand Cu-ttpy (ttpy=4'-tolyl-2,2':6',2''-terpyridine) induced a conformational change that increased the binding of the peptide. This study illustrates that it is important to not only characterize drug-target interactions, but also how the binding to other partners is affected.

Insights into the Kinetic Partitioning Folding Dynamics of the Human Telomeric G-Quadruplex from Molecular Simulations and Machine Learning

The human telomeric DNA G-quadruplex follows a kinetic partitioning folding mechanism. The underlying folding landscape potentially has many minima separated by high free-energy barriers. However, using current theoretical models to characterize this complex folding landscape has remained a challenging problem. In this study, by developing a hybrid atomistic structure-based model that merges structural information on the hybrid-1, hybrid-2, and chair-type G-quadruplex topologies, we investigated a kinetic partitioning folding process of human telomeric DNA involving three native folds. The model was validated as it reproduced the experimental observation that the hybrid-1 conformation is the major fold and the hybrid-2 conformation is kinetically more accessible. A three-step mechanism was revealed for the formation of the hybrid-1 conformation, while a two-step mechanism was demonstrated for the formation of hybrid-2 and chair-type conformations. Likewise, a class of state in which structures adopted inappropriate combinations of syn/anti guanine nucleotides was found to greatly slow down the folding process. In addition, by employing the XGBoost machine learning algorithm, three interatom distances and six dihedral angles were identified as essential internal coordinates to represent the low-dimensional folding landscape. The strategy of coupling the multibasin model and the machine learning algorithm may be useful to investigate the conformational dynamics of other multistate biomolecules.

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 comprehensive evaluation of a typical plant telomeric G-quadruplex (G4) DNA reveals the dynamics of G4 formation, rearrangement, and unfolding

Telomeres are specific nucleoprotein structures that are located at the ends of linear eukaryotic chromosomes and play crucial roles in genomic stability. Telomere DNA consists of simple repeats of a short G-rich sequence: TTAGGG in mammals and TTTAGGG in most plants. In recent years, the mammalian telomeric G-rich repeats have been shown to form G-quadruplex (G4) structures, which are crucial for modulating telomere functions. Surprisingly, even though plant telomeres are essential for plant growth, development, and environmental adaptions, only few reports exist on plant telomeric G4 DNA (pTG4). Here, using bulk and single-molecule assays, including CD spectroscopy, and single-molecule FRET approaches, we comprehensively characterized the structure and dynamics of a typical plant telomeric sequence, d[GGG(TTTAGGG)3]. We found that this sequence can fold into mixed G4s in potassium, including parallel and antiparallel structures. We also directly detected intermediate dynamic transitions, including G-hairpin, parallel G-triplex, and antiparallel G-triplex structures. Moreover, we observed that pTG4 is unfolded by the AtRecQ2 helicase but not by AtRecQ3. The results of our work shed light on our understanding about the existence, topological structures, stability, intermediates, unwinding, and functions of pTG4.

Deciphering the pH-dependence of ground- and excited-state equilibria of thienoguanine

The thienoguanine nucleobase (^thG_b) is an isomorphic fluorescent analogue of guanine. In aqueous buffer at neutral pH, ^thG_b exists as a mixture of two ground-state H1 and H3 keto-amino tautomers with distinct absorption and emission spectra and high quantum yield. In this work, we performed the first systematic photophysical characterization of ^thG_b as a function of pH (2 to 12). Steady-state and time-resolved fluorescence spectroscopies, supplemented with theoretical calculations, enabled us to identify three additional ^thG_b forms, resulting from pH-dependent ground-state and excited-state reactions. Moreover, a thorough analysis allowed us to retrieve their individual absorption and emission spectra as well as the equilibrium constants which govern their interconversion. From these data, the complete photoluminescence pathway of ^thG_b in aqueous solution and its dependence as a function of pH was deduced. As the identified forms differ by their spectra and fluorescence lifetime, ^thG_b could be used as a probe for sensing local pH changes under acidic conditions.

The Potential of Telomeric G-quadruplexes Containing Modified Oligoguanosine Overhangs in Activation of Bacterial Phagocytosis and Leukotriene Synthesis in Human Neutrophils

Human neutrophils are the first line of defense against bacterial and viral infections. They eliminate pathogens through phagocytosis, which activate the 5-lipoxygenase (5-LOX) pathway resulting in synthesis of leukotrienes. Using HPLC analysis, flow cytometry, and other biochemical methods, we studied the effect of synthetic oligodeoxyribonucleotides (ODNs) able to fold into G-quadruplex structures on the main functions of neutrophils. Designed ODNs contained four human telomere TTAGGG repeats (G4) including those with phosphorothioate oligoguanosines attached to the end(s) of G-quadruplex core. Just modified analogues of G4 was shown to more actively than parent ODN penetrate into cells, improve phagocytosis of Salmonella typhimurium bacteria, affect 5-LOX activation, the cytosol calcium ion level, and the oxidative status of neutrophils. As evident from CD and UV spectroscopy data, the presence of oligoguanosines flanking G4 sequence leads to dramatic changes in G-quadruplex topology. While G4 folds into a single antiparallel structure, two main folded forms have been identified in solutions of modified ODNs: antiparallel and dominant, more stable parallel. Thus, both the secondary structure of ODNs and their ability to penetrate into the cytoplasm of cells are important for the activation of neutrophil cellular effects. Our results offer new clues for understanding the role of G-quadruplex ligands in regulation of integral cellular processes and for creating the antimicrobial agents of a new generation.

Studying biomolecular folding and binding using temperature-jump mass spectrometry

Characterizing folding and complex formation of biomolecules provides a view into their thermodynamics, kinetics and folding pathways. Deciphering kinetic intermediates is particularly important because they can often be targeted by drugs. The key advantage of native mass spectrometry over conventional methods that monitor a single observable is its ability to identify and quantify coexisting species. Here, we show the design of a temperature-jump electrospray source for mass spectrometry that allows one to perform fast kinetics experiments (0.16-32 s) at different temperatures (10-90 °C). The setup allows recording of both folding and unfolding kinetics by using temperature jumps from high to low, and low to high, temperatures. Six biological systems, ranging from peptides to proteins to DNA complexes, exemplify the use of this device. Using temperature-dependent experiments, the folding and unfolding of a DNA triplex are studied, providing detailed information on its thermodynamics and kinetics.

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 crystal structure of an antiparallel chair-type G-quadruplex formed by Bromo-substituted human telomeric DNA

Human telomeric guanine-rich DNA, which could adopt different G-quadruplex structures, plays important roles in protecting the cell from recombination and degradation. Although many of these structures were determined, the chair-type G-quadruplex structure remains elusive. Here, we present a crystal structure of the G-quadruplex composed of the human telomeric sequence d[GGGTTAGG8GTTAGGGTTAGG20G] with two dG to 8Br-dG substitutions at positions 8 and 20 with syn conformation in the K+ solution. It forms a novel three-layer chair-type G-quadruplex with two linking trinucleotide loops. Particularly, T5 and T17 are coplanar with two water molecules stacking on the G-tetrad layer in a sandwich-like mode through a coordinating K+ ion and an A6•A18 base pair. While a twisted Hoogsteen A12•T10 base pair caps on the top of G-tetrad core. The three linking TTA loops are edgewise and each DNA strand has two antiparallel adjacent strands. Our findings contribute to a deeper understanding and highlight the unique roles of loop and water molecule in the folding of the G-quadruplex.

Sneaking Out for Happy Hour: Yeast-Based Approaches to Explore and Modulate Immune Response and Immune Evasion

Many pathogens (virus, bacteria, fungi, or parasites) have developed a wide variety of mechanisms to evade their host immune system. The budding yeast Saccharomyces cerevisiae has successfully been used to decipher some of these immune evasion strategies. This includes the cis-acting mechanism that limits the expression of the oncogenic Epstein–Barr virus (EBV)-encoded EBNA1 and thus of antigenic peptides derived from this essential but highly antigenic viral protein. Studies based on budding yeast have also revealed the molecular bases of epigenetic switching or recombination underlying the silencing of all except one members of extended families of genes that encode closely related and highly antigenic surface proteins. This mechanism is exploited by several parasites (that include pathogens such as Plasmodium, Trypanosoma, Candida, or Pneumocystis) to alternate their surface antigens, thereby evading the immune system. Yeast can itself be a pathogen, and pathogenic fungi such as Candida albicans, which is phylogenetically very close to S. cerevisiae, have developed stealthiness strategies that include changes in their cell wall composition, or epitope-masking, to control production or exposure of highly antigenic but essential polysaccharides in their cell wall. Finally, due to the high antigenicity of its cell wall, yeast has been opportunistically exploited to create adjuvants and vectors for vaccination.

Charge-dipole interactions in G-quadruplex thrombin-binding aptamer

DNAs form various structures through hydrogen-bonding, base-stacking and electrostatic interactions. Although these noncovalent interactions are known to be cooperative in stabilizing a G-quadruplex (G4) structure of DNA, we find from all-atom molecular dynamics simulations that the electrostatic charge-dipole interaction is competitive with both hydrogen-bonding and base-stacking interactions. For the thrombin-binding aptamer (TBA) forming a chair-type antiparallel G4 structure, we have examined effects of an intercalating metal ion [K+, Sr2+, Mn+: an ion having a charge of n+ (n = 1-4) with the ionic radius of K+] on structural properties and noncovalent interactions. When K+ in the TBA·K+ complex is replaced with Sr2+, guanine dipoles in the two G-tetrads are realigned toward the central metal ion, thereby distorting the planar G4 geometry. Replacing K+ with Sr2+ significantly enhances the charge-dipole interaction but substantially reduces the number of hydrogen bonds in the G-tetrads. In the case of TBA·Mn+ complexes, as the charge n increases, the charge-dipole interaction increases but both of the hydrogen-bonding and base-stacking interactions decrease. These results suggest that the charge-dipole interaction realigning guanine dipoles in the G-tetrads is not cooperative but competitive with both hydrogen-bonding and base-stacking interactions favoring the planar G-tetrad geometry. Obviously, the charge state of an intercalating metal ion is as important as the ionic radius in forming a stable G4 structure. Thus, a delicate balance between these competing noncovalent interactions makes the chair-type antiparallel G4 structure of TBA selective for intercalating metal ions.

Excited-State Dynamics of Thienoguanosine, an Isomorphic Highly Fluorescent Analogue of Guanosine

Thienoguanosine (^th G) is an isomorphic analogue of guanosine with promising potentialities as fluorescent DNA label. As a free probe in protic solvents, ^th G exists in two tautomeric forms, identified as the H1, being the only one observed in nonprotic solvents, and H3 keto-amino tautomers. We herein investigate the photophysics of ^th G in solvents of different polarity, from water to dioxane, by combining time-resolved fluorescence with PCM/TD-DFT and CASSCF calculations. Fluorescence lifetimes of 14.5-20.5 and 7-13 ns were observed for the H1 and H3 tautomers, respectively, in the tested solvents. In methanol and ethanol, an additional fluorescent decay lifetime (≈3 ns) at the blue emission side (λ≈430 nm) as well as a 0.5 ns component with negative amplitude at the red edge of the spectrum, typical of an excited-state reaction, were observed. Our computational analysis explains the solvent effects observed on the tautomeric equilibrium. The main radiative and nonradiative deactivation routes have been mapped by PCM/TD-DFT calculations in solution and CASSCF in the gas phase. The most easily accessible conical intersection, involving an out-of plane motion of the sulfur atom in the five-membered ring of ^th G, is separated by a sizeable energy barrier (≥0.4 eV) from the minimum of the spectroscopic state, which explains the large experimental fluorescence quantum yield.

Fast-Folding Pathways of the Thrombin-Binding Aptamer G-Quadruplex Revealed by a Markov State Model

G-quadruplex structures participate in many important cellular processes. For a better understanding of their functions, knowledge of the mechanism by which they fold into the functional native structures is necessary. In this work, we studied the folding process of the thrombin-binding aptamer G-quadruplex. Enabled by a computational paradigm that couples an advanced sampling method and a Markov state model, four folding intermediates were identified, including an antiparallel G-hairpin, two G-triplex structures, and a double-hairpin conformation. Likewise, a misfolded structure with a nonnative distribution of syn/anti guanines was also observed. Based on these states, a transition path analysis revealed three fast-folding pathways, along which the thrombin-binding aptamer would fold to the native state directly, with no evidence of potential nonnative competing conformations. The results also showed that the TGT-loop plays an important role in the folding process. The findings of this research may provide general insight about the folding of other G-quadruplex structures.

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.

DNA Quadruple Helices in Nanotechnology

DNA has played an early and powerful role in the development of bottom-up nanotechnologies, not least because of DNA's precise, predictable, and controllable properties of assembly on the nanometer scale. Watson-Crick complementarity has been used to build complex 2D and 3D architectures and design a number of nanometer-scale systems for molecular computing, transport, motors, and biosensing applications. Most of such devices are built with classical B-DNA helices and involve classical A-T/U and G-C base pairs. However, in addition to the above components underlying the iconic double helix, a number of alternative pairing schemes of nucleobases are known. This review focuses on two of these noncanonical classes of DNA helices: G-quadruplexes and the i-motif. The unique properties of these two classes of DNA helix have been utilized toward some remarkable constructions and applications: G-wires; nanostructures such as DNA origami; reconfigurable structures and nanodevices; the formation and utilization of hemin-utilizing DNAzymes, capable of generating varied outputs from biosensing nanostructures; composite nanostructures made up of DNA as well as inorganic materials; and the construction of nanocarriers that show promise for the therapeutics of diseases.

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.

Environmentally Sensitive Fluorescent Nucleoside Analogues for Surveying Dynamic Interconversions of Nucleic Acid Structures

Nucleic acids are characterized by a variety of dynamically interconverting structures that play a major role in transcriptional and translational regulation as well as recombination and repair. To monitor these interconversions, Förster resonance energy transfer (FRET)-based techniques can be used, but require two fluorophores that are typically large and can alter the DNA/RNA structure and protein binding. Additionally, events that do not alter the donor/acceptor distance and/or angular relationship are frequently left undetected. A more benign approach relies on fluorescent nucleobases that can substitute their native counterparts with minimal perturbation, such as the recently developed 2-thienyl-3-hydroxychromone (3HCnt) and thienoguanosine (^th G). To demonstrate the potency of 3HCnt and ^th G in deciphering interconversion mechanisms, we used the conversion of the (-)DNA copy of the HIV-1 primer binding site (-)PBS stem-loop into (+)/(-)PBS duplex, as a model system. When incorporated into the (-)PBS loop, the two probes were found to be highly sensitive to the individual steps both in the absence and the presence of a nucleic acid chaperone, providing the first complete mechanistic description of this critical process in HIV-1 replication. The combination of the two distinct probes appears to be instrumental for characterizing structural transitions of nucleic acids under various stimuli.

Probing DNA Hybridization Equilibrium by Cationic Conjugated Polymer for Highly Selective Detection and Imaging of Single-Nucleotide Mutation

Hybridization-based probes emerge as a promising tool for nucleic acid target detection and imaging. However, the single-nucleotide selectivity is still challenging because the specificity of hybridization reaction is typically low at room temperature. We disclose an effective and simple method for highly selective detection and in situ imaging of single-nucleotide mutation (SNM) by taking the advantages of the specific hybridization of short duplex and the signal amplifying effect of cationic conjugated polymer (CCP). Excellent discrimination of the nucleic acid strands only differing by single nucleotide was achieved enabling the sensitive detection of SNM at the abundance as low as 0.1%. Single-molecule fluorescence resonance energy transfer (smFRET) study reveals that the presence of CCP enhances the perfect matched duplex and the mismatched duplex to a different extent, which can be an explanation for the high single-nucleotide selectivity. Due to the simple design of the probe and the stable brightness of CCP, highly selective mRNA in situ imaging was achieved in fixed cells. Melanoma cell line A375 with BRAF V600E point mutation exhibits higher FRET efficiency than liver cancer cell line HegG2 that was not reported having the mutation at this point.

Single-Stranded Condensation Stochastically Blocks G-Quadruplex Assembly in Human Telomeric RNA

TERRA is an RNA molecule transcribed from human subtelomeric regions toward chromosome ends potentially involved in regulation of heterochromatin stability, semiconservative replication, and telomerase inhibition, among others. TERRA contains tandem repeats of the sequence GGGUUA, with a strong tendency to fold into a four-stranded arrangement known as a parallel G-quadruplex. Here, we demonstrate by using single-molecule force spectroscopy that this potential is limited by the inherent capacity of RNA to self-associate randomly and further condense into entropically more favorable structures. We stretched RNA constructions with more than four and less than eight hexanucleotide repeats, thus unable to form several G-quadruplexes in tandem, flanked by non-G-rich overhangs of random sequence by optical tweezers on a one by one basis. We found that condensed RNA stochastically blocks G-quadruplex folding pathways with a near 20% probability, a behavior that is not found in DNA analogous molecules.

G-quadruplexes in Prague: A Bohemian Rhapsody

The Sixth International Meeting on Quadruplex Nucleic Acids was held at the Hotel Internationale in Prague, Czech Republic from 31 May - 3 June 2017. A vibrant interdisciplinary community of over 300 scientists gathered to share their newest results in this exciting field and exchange ideas for further investigations.

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.

Involvement of G-triplex and G-hairpin in the multi-pathway folding of human telomeric G-quadruplex

G-quadruplex (G4) can be formed by G-rich DNA sequences that are widely distributed throughout the human genome. Although G-triplex and G-hairpin have been proposed as G4 folding intermediates, their formation still requires further investigation by experiments. Here, we employed single-molecule FRET to characterize the folding dynamics of G4 from human telomeric sequence. First, we observed four states during G4 folding initially assigned to be anti-parallel G4, G-triplex, G-hairpin and unfolded ssDNA. Then we constructed putative intra-strand G-triplex, G-hairpin structures and confirmed their existences in both NaCl and KCl. Further studies revealed those structures are going through dynamic transitions between different states and show relatively weak dependence on cations, unlike G4. Based on those results and molecular dynamics simulations, we proposed a multi-pathway folding mechanism for human telomeric G4. The present work may shed new light on our current understanding about the existence and stability of G4 intermediate states.

Kinetic evidence for interaction of TMPyP4 with two different G-quadruplex conformations of human telomeric DNA

BACKGROUND

Stabilization of G-quadruplex helices by small ligands has attracted growing attention because they inhibit the activity of the enzyme telomerase, which is overexpressed in >80% cancer cells. TMPyP4, one of the most studied G-quadruplex ligands, is used as a model to show that the ligands can exhibit different binding features with different conformations of a human telomeric specific sequence.

METHODS

UV-Vis, FRET melting Assay, Isothermal Titration Calorimetry, Time-resolved Fluorescence lifetime, T-Jump and Molecular Dynamics.

RESULTS

TMPyP4 yields two different complexes with two Tel22 telomeric conformations in the presence of Na⁺ or K⁺. T-Jump kinetic experiments show that the rates of formation and dissociation of these complexes in the ms time scale differ by one order of magnitude. MD simulations reveal that, in K⁺ buffer, "hybrid 1" conformation yields kinetic constants on interaction with TMPyP4 one order lower than "hybrid 2". The binding involves π-π stacking with external loop bases.

CONCLUSIONS

For the first time we show that for a particular buffer TMPyP4 interacts in a kinetically different way with the two Tel22 conformations even if the complexes formed are thermodynamically indistinguishable.

GENERAL SIGNIFICANCE

G-quadruplexes, endowed with technological applications and potential impact on regulation mechanisms, define a new research field. The possibility of building different conformations from same sequence is a complex issue that confers G-quadruplexes very interesting features. The obtaining of reliable kinetic data constitutes an efficient tool to determine reaction mechanisms between conformations and small molecules.