Folding Landscape of a Parallel G-Quadruplex

Distinguishing G-Quadruplexes Stabilizer and Chaperone for c-MYC Promoter G-Quadruplexes through Single-Molecule Manipulation

G-quadruplex (G4) selective stabilizing ligands can regulate c-MYC gene expression, but the kinetic basis remains unclear. Determining the effects of ligands on c-MYC promoter G4s' folding/unfolding kinetics is challenging due to the polymorphic nature of G4s and the high energy barrier to unfold c-MYC promoter G4s. Here, we used single-molecule magnetic tweezers to manipulate a duplex hairpin containing a c-MYC promoter sequence to mimic the transiently denatured duplex during transcription. We measured the effects of six commonly used G4s binding ligands on the competition between quadruplex and duplex structures, as well as the folding/unfolding kinetics of G4s. Our results revealed two distinct roles for G4s selective stabilization: CX-5461 is mainly acting as c-MYC G4s stabilizer, reducing the unfolding rate (ku) of c-MYC G4s, whereas PDS and 360A also act as G4s chaperone, accelerating the folding rates (kf) of c-MYC G4s. qRT-PCR results obtained from CA46 and Raji cell lines demonstrated that G4s stabilizing ligands can downregulate c-MYC expression, while G4s stabilizer CX-5461 exhibited the strongest c-MYC gene suppression. These results shed light on the potential of manipulating G4s' folding/unfolding kinetics by ligands for precise regulation of promoter G4-associated biological activities.

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.

Stability of Human Telomeric G-Quadruplexes Complexed with Photosensitive Ligands and Irradiated with Visible Light

Guanine-rich DNA sequences can fold into non-canonical nucleic acid structures called G-quadruplexes (G4s). These nanostructures have strong implications in many fields, from medical science to bottom-up nanotechnologies. As a result, ligands interacting with G4s have attracted great attention as candidates in medical therapies, molecular probe applications, and biosensing. In recent years, the use of G4-ligand complexes as photopharmacological targets has shown significant promise for developing novel therapeutic strategies and nanodevices. Here, we studied the possibility of manipulating the secondary structure of a human telomeric G4 sequence through the interaction with two photosensitive ligands, DTE and TMPyP4, whose response to visible light is different. The effect of these two ligands on G4 thermal unfolding was also considered, revealing the occurrence of peculiar multi-step melting pathways and the different attitudes of the two molecules on the quadruplex stabilization.

Structure of a 28.5 kDa duplex-embedded G-quadruplex system resolved to 7.4 Å resolution with cryo-EM

Genomic regions with high guanine content can fold into non-B form DNA four-stranded structures known as G-quadruplexes (G4s). Extensive in vivo investigations have revealed that promoter G4s are transcriptional regulators. Little structural information exists for these G4s embedded within duplexes, their presumed genomic environment. Here, we report the 7.4 Å resolution structure and dynamics of a 28.5 kDa duplex-G4-duplex (DGD) model system using cryo-EM, molecular dynamics, and small-angle X-ray scattering (SAXS) studies. The DGD cryo-EM refined model features a 53° bend induced by a stacked duplex-G4 interaction at the 5' G-tetrad interface with a persistently unstacked 3' duplex. The surrogate complement poly dT loop preferably stacks onto the 3' G-tetrad interface resulting in occlusion of both 5' and 3' tetrad interfaces. Structural analysis shows that the DGD model is quantifiably more druggable than the monomeric G4 structure alone and represents a new structural drug target. Our results illustrate how the integration of cryo-EM, MD, and SAXS can reveal complementary detailed static and dynamic structural information on DNA G4 systems.

Variation in G-quadruplex sequence and topology differentially impacts human DNA polymerase fidelity

G-quadruplexes (G4s), a type of non-B DNA, play important roles in a wide range of molecular processes, including replication, transcription, and translation. Genome integrity relies on efficient and accurate DNA synthesis, and is compromised by various stressors, to which non-B DNA structures such as G4s can be particularly vulnerable. However, the impact of G4 structures on DNA polymerase fidelity is largely unknown. Using an in vitro forward mutation assay, we investigated the fidelity of human DNA polymerases delta (δ4, four-subunit), eta (η), and kappa (κ) during synthesis of G4 motifs representing those in the human genome. The motifs differ in sequence, topology, and stability, features that may affect DNA polymerase errors. Polymerase error rate hierarchy (δ4 < κ < η) is largely maintained during G4 synthesis. Importantly, we observed unique polymerase error signatures during synthesis of VEGF G4 motifs, stable G4s which form parallel topologies. These statistically significant errors occurred within, immediately flanking, and encompassing the G4 motif. For pol δ4, the errors were deletions, insertions and complex errors within the G4 or encompassing the G4 motif and surrounding sequence. For pol η, the errors occurred in 3' sequences flanking the G4 motif. For pol κ, the errors were frameshift mutations within G-tracts of the G4. Because these error signatures were not observed during synthesis of an antiparallel G4 and, to a lesser extent, a hybrid G4, we suggest that G4 topology and/or stability could influence polymerase fidelity. Using in silico analyses, we show that most polymerase errors are predicted to have minimal effects on predicted G4 stability. Our results provide a unique view of G4s not previously elucidated, showing that G4 motif heterogeneity differentially influences polymerase fidelity within the motif and flanking sequences. Thus, our study advances the understanding of how DNA polymerase errors contribute to G4 mutagenesis.

Distribution of Conformational States Adopted by DNA from the Promoter Regions of the VEGF and Bcl-2 Oncogenes

We employed a previously described procedure, based on circular dichroism (CD) spectroscopy, to quantify the distribution of conformational states adopted by equimolar mixtures of complementary G-rich and C-rich DNA strands from the promoter regions of the VEGF and Bcl-2 oncogenes. Spectra were recorded at different pHs, concentrations of KCl, and temperatures. The temperature dependences of the fractional populations of the duplex, G-quadruplex, i-motif, and coiled conformations of each promoter were then analyzed within the framework of a thermodynamic model to obtain the enthalpy and melting temperature of each folded-to-unfolded transition involved in the equilibrium. A comparison of the conformational data on the VEGF and Bcl-2 DNA with similar results on the c-MYC DNA, which we reported previously, reveals that the distribution of conformational states depends on the specific DNA sequence and is modulated by environmental factors. Under the physiological conditions of room temperature, neutral pH, and elevated concentrations of potassium ions, the duplex conformation coexists with the G-quadruplex conformation in proportions that depend on the sequence. This observed conformational diversity has biological implications, and it further supports our previously proposed thermodynamic hypothesis of gene regulation. In that hypothesis, a specific distribution of duplex and tetraplex conformations in a promoter region is fine-tuned to maintain the healthy level of gene expression. Any deviation from a healthy distribution of conformational states may result in pathology stemming from up- or downregulation of the gene.

Biophysical characterisation of the Bcl-x pre-mRNA and binding specificity of the ellipticine derivative GQC-05: Implication for alternative splicing regulation

The BCL2L1 gene expresses two isoforms of Bcl-x protein via the use of either of two alternative 5′ splice sites (5′ss) in exon 2. These proteins have antagonistic actions, Bcl-X_L being anti-apoptotic and Bcl-X_S pro-apoptotic. In a number of cancers the Bcl-X_L isoform is over-expressed, resulting in cancer cell survival and growth, so switching splicing to the X_s isoform could have therapeutic benefits. We have previously proposed that a putative G-quadruplex (G4) exists downstream of the X_S 5′ss and shown that the ellipticine derivative GQC-05, a previously identified DNA G4-specific ligand, induces an increase in the X_S/X_L ratio both in vitro and in cells. Here, we demonstrate that this G4 forms in vitro and that the structure is stabilised in the presence of GQC-05. We also show that GQC-05 binds RNA non-specifically in buffer conditions, but selectively to the Bcl-x G4 in the presence of nuclear extract, highlighting the limitations of biophysical measurements taken outside of a functional environment. We also demonstrate that GQC-05 is able to shift the equilibrium between competing G4 and duplex structures towards the G4 conformation, leading to an increase in accessibility of the X_S 5′ss, supporting our previous model on the mechanism of action of GQC-05.

Iso-FRET: an isothermal competition assay to analyze quadruplex formation in vitro

Algorithms have been widely used to predict G-quadruplexes (G4s)-prone sequences. However, an experimental validation of these predictions is generally required. We previously reported a high-throughput technique to evidence G4 formation in vitro called FRET-MC. This method, while convenient and reproducible, has one known weakness: its inability to pin point G4 motifs of low thermal stability. As such quadruplexes may still be biologically relevant if formed at physiological temperature, we wanted to develop an independent assay to overcome this limitation. To this aim, we introduced an isothermal version of the competition assay, called iso-FRET, based on a duplex-quadruplex competition and a well-characterized bis-quinolinium G4 ligand, PhenDC3. G4-forming competitors act as decoys for PhenDC3, lowering its ability to stabilize the G4-forming motif reporter oligonucleotide conjugated to a fluorescence quencher (37Q). The decrease in available G4 ligand concentration restores the ability of 37Q to hybridize to its FAM-labeled short complementary C-rich strand (F22), leading to a decrease in fluorescence signal. In contrast, when no G4-forming competitor is present, PhenDC3 remains available to stabilize the 37Q quadruplex, preventing the formation of the F22 + 37Q complex. Iso-FRET was first applied to a reference panel of 70 sequences, and then used to investigate 23 different viral sequences.

Solvent Vibrations as a Proxy of the Telomere G-Quadruplex Rearrangements across Thermal Unfolding

G-quadruplexes (G4s) are noncanonical forms of DNA involved in many key genome functions. Here, we exploited UV Resonance Raman scattering to simultaneously explore the vibrational behavior of a human telomeric G4 (Tel22) and its aqueous solvent as the biomolecule underwent thermal melting. We found that the OH stretching band, related to the local hydrogen-bonded network of a water molecule, was in strict relation with the vibrational features of the G4 structure as a function of temperature. In particular, the modifications to the tetrahedral ordering of the water network were strongly coupled to the DNA rearrangements, showing changes in temperature that mirrored the multi-step melting process of Tel22. The comparison between circular dichroism and Raman results supported this view. The present findings provide novel insights into the impact of the molecular environment on G4 conformation. Improving current knowledge on the solvent structural properties will also contribute to a better understanding of the role played by water arrangement in the complexation of G4s with ligands.

Multimeric G-quadruplexes: A review on their biological roles and targeting

In human cells, nucleic acids adopt several non-canonical structures that regulate key cellular processes. Among them, G-quadruplexes (G4s) are stable structures that form in guanine-rich regions in vitro and in cells. G4 folded/unfolded state shapes numerous cellular processes, including genome replication, transcription, and translation. Moreover, G4 folding is involved in genomic instability. G4s have been described to multimerize, forming high-order structures in both DNA and/or RNA strands. Multimeric G4s can be formed by adjacent intramolecular G4s joined by stacking interactions or connected by short loops. Multimeric G4s can also originate from the assembly of guanines embedded on independent DNA or RNA strands. Notably, crucial regions of the human genome, such as the 3'-terminal overhang of the telomeric DNA as well as the open reading frame of genes involved in the preservation of neuron viability in the human central and peripheral nervous system are prone to form multimeric G4s. The biological importance of such structures has been recently described, with multimeric G4s playing potentially protective or deleterious effects in the pathogenic cascade of various diseases. Here, we portray the multifaceted scenario of multimeric G4s, in terms of structural properties, biological roles, and targeting strategies.

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.

Multistep mechanism of G-quadruplex resolution during DNA replication

G-quadruplex (or G4) structures form in guanine-rich DNA sequences and threaten genome stability when not properly resolved. G4 unwinding occurs during S phase via an unknown mechanism. Using Xenopus egg extracts, we define a three-step G4 unwinding mechanism that acts during DNA replication. First, the replicative helicase composed of Cdc45, MCM2-7 and GINS (CMG) stalls at a leading strand G4 structure. Second, the DEAH-box helicase 36 (DHX36) mediates bypass of the CMG past the intact G4 structure, allowing approach of the leading strand to the G4. Third, G4 structure unwinding by the Fanconi anemia complementation group J helicase (FANCJ) enables DNA polymerase to synthesize past the G4 motif. A G4 on the lagging strand template does not stall CMG but still requires DNA replication for unwinding. DHX36 and FANCJ have partially redundant roles, conferring pathway robustness. This previously unknown genome maintenance pathway promotes faithful G4 replication, thereby avoiding genome instability.

Light-induced in situ chemical activation of a fluorescent probe for monitoring intracellular G-quadruplex structures

Light-activated functional materials capable of remote control over duplex and G-quadruplex (G4) nucleic acids formation at the cellular level are still very rare. Herein, we report on the photoinduced macrocyclisation of a helicenoid quinoline derivative of binaphthol that selectively provides easy access to an unprecedented class of extended heteroaromatic structures with remarkable photophysical and DNA/RNA binding properties. Thus, while the native bisquinoline precursor shows no DNA binding activity, the new in situ photochemically generated probe features high association constants to DNA and RNA G4s. The latter inhibits DNA synthesis by selectively stabilizing G4 structures associated with oncogenic promoters and telomere repeat units. Finally, the light sensitive compound is capable of in cellulo photoconversion, localizes primarily in the G4-rich sites of cancer cells, competes with a well-known G4 binder and shows a clear nuclear co-localization with the quadruplex specific antibody BG4. This work provides a benchmark for the future design and development of a brand-new generation of light-activated target-selective G4-binders.

Porphyrin Binding and Irradiation Promote G-Quadruplex DNA Dimeric Structure

Nucleic acid sequences rich in guanines can organize into noncanonical DNA G-quadruplexes (G4s) of variable size. The design of small molecules stabilizing the structure of G4s is a rapidly growing area for the development of novel anticancer therapeutic strategies and bottom-up nanotechnologies. Among a multitude of binders, porphyrins are very attractive due to their light activation that can make them valuable conformational regulators of G4s. Here, a structure-based strategy, integrating complementary probes, is employed to study the interaction between TMPyP4 porphyrin and a 22-base human telomeric sequence (Tel22) before and after irradiation with blue light. Porphyrin binding is discovered to promote Tel22 dimerization, while light irradiation of the Tel22-TMPyP4 complex controls dimer fraction. Such a change in quaternary structure is found to be strictly correlated with modifications at the secondary structure level, thus providing an unprecedented link between the degree of dimerization and the underlying conformational changes in G4s.

A Minimalistic Coumarin Turn-On Probe for Selective Recognition of Parallel G-Quadruplex DNA Structures

G-quadruplex (G4) DNA structures are widespread in the human genome and are implicated in biologically important processes such as telomere maintenance, gene regulation, and DNA replication. Guanine-rich sequences with potential to form G4 structures are prevalent in the promoter regions of oncogenes, and G4 sites are now considered as attractive targets for anticancer therapies. However, there are very few reports of small “druglike” optical G4 reporters that are easily accessible through one-step synthesis and that are capable of discriminating between different G4 topologies. Here, we present a small water-soluble light-up fluorescent probe that features a minimalistic amidinocoumarin-based molecular scaffold that selectively targets parallel G4 structures over antiparallel and non-G4 structures. We showed that this biocompatible ligand is able to selectively stabilize the G4 template resulting in slower DNA synthesis. By tracking individual DNA molecules, we demonstrated that the G4-stabilizing ligand perturbs DNA replication in cancer cells, resulting in decreased cell viability. Moreover, the fast-cellular entry of the probe enabled detection of nucleolar G4 structures in living cells. Finally, insights gained from the structure–activity relationships of the probe suggest the basis for the recognition of parallel G4s, opening up new avenues for the design of new biocompatible G4-specific small molecules for G4-driven theranostic applications.

On the thermodynamics of folding of an i-motif DNA in solution under favorable conditions

Under slightly acidic conditions, cytosine-rich DNA sequences can form non-canonical secondary structures called i-motifs, which occur as four stretches of cytosine repeats form hemi-protonated C·C+ base pairs. The growing interest in the i-motif structures as important components in functional DNA-based nanotechnology or as potential targets of anticancer drugs, increases the need for a deep understanding of the energetics of their structural transitions. Here, a combination of spectroscopic and calorimetric techniques is used to unravel the thermodynamics of folding of an i-motif DNA under favorable conditions. The results give new insights into the energetic aspects of i-motifs and show that thermodynamic and thermal stability are related but not identical properties of such DNA structures.

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.

The Folding Landscapes of Human Telomeric RNA and DNA G-Quadruplexes are Markedly Different

We investigated the folding kinetics of G‐quadruplex (G4) structures by comparing the K⁺‐induced folding of an RNA G4 derived from the human telomeric repeat‐containing RNA (TERRA25) with a sequence homologous DNA G4 (wtTel25) using CD spectroscopy and real‐time NMR spectroscopy. While DNA G4 folding is biphasic, reveals kinetic partitioning and involves kinetically favoured off‐pathway intermediates, RNA G4 folding is faster and monophasic. The differences in kinetics are correlated to the differences in the folded conformations of RNA vs. DNA G4s, in particular with regard to the conformation around the glycosidic torsion angle χ that uniformly adopts anti conformations for RNA G4s and both, syn and anti conformation for DNA G4s. Modified DNA G4s with ¹⁹F bound to C2′ in arabino configuration adopt exclusively anti conformations for χ. These fluoro‐modified DNA (antiTel25) reveal faster folding kinetics and monomorphic conformations similar to RNA G4s, suggesting the correlation between folding kinetics and pathways with differences in χ angle preferences in DNA and RNA, respectively.

Unraveling the Kinetics of Spare-Tire DNA G-Quadruplex Folding

The folding of DNA G-quadruplexes (G4) is essential to regulate expression of oncogenes and involves polymorphic long-lived intermediate states. G4 formation requires four G-tracts, but human gene-promoters often contain multiple G-tracts that act as spare-tires. These additional G-tracts are highly conserved and add multiple layers of functional complexity, as they are crucial to maintain G4 function after oxidative damage. Herein, we unravel the folding dynamics of the G4 sequence containing five G-tracts from cMYC, the major proliferation-driving oncogene. We devise a general method to induce folding at constant experimental conditions using a photochemical trapping strategy. Our data dissect the individual kinetics and thermodynamics of the spare-tire mechanism of cMYC-G4.

Insight into formation propensity of pseudocircular DNA G-hairpins

We recently showed that Saccharomyces cerevisiae telomeric DNA can fold into an unprecedented pseudocircular G-hairpin (PGH) structure. However, the formation of PGHs in the context of extended sequences, which is a prerequisite for their function in vivo and their applications in biotechnology, has not been elucidated. Here, we show that despite its ‘circular’ nature, PGHs tolerate single-stranded (ss) protrusions. High-resolution NMR structure of a novel member of PGH family reveals the atomistic details on a junction between ssDNA and PGH unit. Identification of new sequences capable of folding into one of the two forms of PGH helped in defining minimal sequence requirements for their formation. Our time-resolved NMR data indicate a possibility that PGHs fold via a complex kinetic partitioning mechanism and suggests the existence of K⁺ ion-dependent PGH folding intermediates. The data not only provide an explanation of cation-type-dependent formation of PGHs, but also explain the unusually large hysteresis between PGH melting and annealing noted in our previous study. Our findings have important implications for DNA biology and nanotechnology. Overrepresentation of sequences able to form PGHs in the evolutionary-conserved regions of the human genome implies their functionally important biological role(s).

Parallel reaction pathways accelerate folding of a guanine quadruplex

G-quadruplexes (G4s) are four-stranded, guanine-rich nucleic acid structures that can influence a variety of biological processes such as the transcription and translation of genes and DNA replication. In many cases, a single G4-forming nucleic acid sequence can adopt multiple different folded conformations that interconvert on biologically relevant timescales, entropically stabilizing the folded state. The coexistence of different folded conformations also suggests that there are multiple pathways leading from the unfolded to the folded state ensembles, potentially modulating the folding rate and biological activity. We have developed an experimental method for quantifying the contributions of individual pathways to the folding of conformationally heterogeneous G4s that is based on mutagenesis, thermal hysteresis kinetic experiments and global analysis, and validated our results using photocaged kinetic NMR experiments. We studied the regulatory Pu22 G4 from the c-myc oncogene promoter, which adopts at least four distinct folded isomers. We found that the presence of four parallel pathways leads to a 2.5-fold acceleration in folding; that is, the effective folding rate from the unfolded to folded ensembles is 2.5 times as large as the rate constant for the fastest individual pathway. Since many G4 sequences can adopt many more than four isomers, folding accelerations of more than an order of magnitude are possible via this mechanism.

High Mechanical Stability and Slow Unfolding Rates Are Prevalent in Parallel-Stranded DNA G-Quadruplexes

Guanine-rich repeat sequences are known to adopt diverse G-quadruplex (G4) topologies. Determining the unfolding rates of individual G4 species is challenging due to the coexistence of multiple G4 conformations in a solution. Here, using single-molecule magnetic tweezers, we systematically measured the unfolding force distributions of 4 oncogene promoter G4s, 12 model sequences with two 1-nucleotide (nt) thymine loops that predominantly adopt parallel-stranded G4 structures, and 6 sequences forming multiple G4 structures. All parallel-stranded G4s reveal an unfolding force peak at 40-60 pN, which is associated with extremely slow unfolding rates on the order of 10^-5-10^-7 s^-1. In contrast, nonparallel G4s and partially folded intermediate states reveal an unfolding force peak <40 pN. These results suggest a strong correlation between the parallel-stranded G4s folding topology and the slow unfolding rates and provide important insights into the mechanism that govern the stability and the transition kinetics of G4s.

Diversity of Parallel Guanine Quadruplexes Induced by Guanine Substitutions

Recently, we reported an inhibitory effect of guanine substitutions on the conformational switch from antiparallel to parallel quadruplexes (G4) induced by dehydrating agents. As a possible cause, we proposed a difference in the sensitivity of parallel and antiparallel quadruplexes to the guanine substitutions in the resulting thermodynamic stability. Reports on the influence of guanine substitutions on the biophysical properties of intramolecular parallel quadruplexes are rare. Moreover, such reports are often complicated by the multimerisation tendencies of parallel quadruplexes. To address this incomplete knowledge, we employed circular dichroism spectroscopy (CD), both as stopped-flow-assisted fast kinetics measurements and end-point measurements, accompanied by thermodynamic analyses, based on UV absorption melting profiles, and electrophoretic methods. We showed that parallel quadruplexes are significantly more sensitive towards guanine substitutions than antiparallel ones. Furthermore, guanine-substituted variants, which in principle might correspond to native genomic sequences, distinctly differ in their biophysical properties, indicating that the four guanines in each tetrad of parallel quadruplexes are not equal. In addition, we were able to distinguish by CD an intramolecular G4 from intermolecular ones resulting from multimerisation mediated by terminal tetrad association, but not from intermolecular G4s formed due to inter-strand Hoogsteen hydrogen bond formation. In conclusion, our study indicates significant variability in parallel quadruplex structures, otherwise disregarded without detailed experimental analysis.

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.

Folding topology, structural polymorphism, and dimerization of intramolecular DNA G-quadruplexes with inverted polarity strands and non-natural loops

Synthetically modified DNA G-quadruplexes (GQs) have great potential in the development of designer molecules for a wide range of applications. Identification of the role of various structural elements in the folding and final topology of artificial GQs is necessary to predict their secondary structure. We report here the results of spectroscopic and electrophoretic studies of GQ scaffolds formed by G-rich sequences comprising four G3-tracts of different polarity connected by either a single-nucleotide thymine loop or a non-natural tetraethyleneglycol loop. Depending on G-strand polarities, loop arrangement and the presence of extra 5'-base G-rich oligonucleotides form monomeric, dimeric, or multimeric species of different topologies. In most cases, oligonucleotides were able to fold into stable parallel or hybrid GQs. However, certain specific arrangements of loops and G-tracts resulted in a diverse mixture of low stable structures. Comparative analysis of topology, stability, and structural heterogeneity of different G-rich sequences suggests the important role of loop type and arrangement, G3-tract polarities, and the presence of 5'-capping residues in the outcome of the folding process. The results also imply that the formation of anti-parallel G-hairpin intermediates is a key event in major favourable folding pathways.

A site-specific self-assembled light-up rotor probe for selective recognition and stabilization of c-MYC G-quadruplex DNA

Direct and unambiguous evidence of the formation of G-quadruplexes (G4s) in human cells have shown their implication in several key biological events and has emphasized their role as important targets for small-molecule cancer therapeutics. Here, we report on the first example of a self-assembled molecular-rotor G4-binder able to discriminate between an extensive panel of G4 and non-G4 structures and to selectively light-up (up to 64-fold), bind (nanomolar range), and stabilize the c-MYC promoter G4 DNA. In particular, association with the c-MYC G4 triggers the disassembly of its supramolecular state (disaggregation-induced emission, DIE) and induces geometrical restrictions (motion-induced change in emission, MICE) leading to a significant enhancement of its emission yield. Moreover, this optical reporter is able to selectively stabilize the c-MYC G4 and inhibit DNA synthesis. Finally, by using confocal laser-scanning microscopy (CLSM) we show the ability of this compound to localize primarily in the subnuclear G4-rich compartments of cancer cells. This work provides a benchmark for the future design and development of a new generation of smart sequence-selective supramolecular G4-binders that combine outstanding sensing and stability properties, to be utilized in anti-cancer therapy.

Conformational Preferences of DNA Strands from the Promoter Region of the c-MYC Oncogene

We characterized the conformational preferences of DNA in an equimolar mixture of complementary G-rich and C-rich strands from the promoter region of the c-MYC oncogene. Our CD-based approach presupposes that the CD spectrum of such a mixture is the spectral sum of the constituent duplex, G-quadruplex, i-motif, and coiled conformations. Spectra were acquired over a range of temperatures at different pHs and concentrations of KCl. Each spectrum was unmixed in terms of the predetermined spectra of the constituent conformational states to obtain the corresponding weighting factors for their fractional contributions to the total population of DNA. The temperature dependences of those contributions then were analyzed in concert according to a model based on a thermodynamic representation of the underlying equilibria. Fitted estimates of the melting enthalpy and temperature obtained for the duplex, G-quadruplex, and i-motif imply that the driving force behind dissociation of the duplex and the concomitant formation of tetrahelical structures is the folding of the G-strand into the G-quadruplex. The liberated C-strand adopts the i-motif conformation at acidic pH and exists in the coiled state at neutral pH. The i-motif alone cannot induce dissociation of the duplex even at pH 5.0, at which it is most stable. Under the physiological conditions of neutral pH, elevated potassium, and room temperature, the duplex and G-quadruplex conformations coexist with the C-strand in the coiled state. Taken together, our results suggest a novel, thermodynamically controlled mechanism for the regulation of gene expression.

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.