Engineering the quadruplex fold: nucleoside conformation determines both folding topology and molecularity in guanine quadruplexes

Defined folding pattern of poly(rG) supports inherent ability to encode biological information

The RNA World hypothesis posits that RNA can represent a primitive life form by reproducing itself and demonstrating catalytic activity. However, this hypothesis is incapable of addressing several major origin-of-life (OoL) questions. A recently described paradox-free alternative OoL hypothesis, the Quadruplex (G4) World, is based on the ability of poly(dG) to fold into a stable architecture with an unambiguous folding pattern using G-tetrads as building elements. Because of the folding pattern of three G-tetrads and single-G loops, dG15 is programmable and has the capability to encode biological information. Here, we address two open questions of the G4 World hypothesis: (1) Does RNA follow the same folding pattern as DNA? (2) How do stable quadruplexes evolve into the present-day system of information transfer, which is based on Watson-Crick base pair complementarity? To address these questions, we systematically studied the thermodynamic and optical properties of both DNA and RNA G15- and G3T (GGGTGGGTGGGTGGG)-derived sequences. Our study revealed that similar to DNA sequences, RNAs adopt quadruplexes with only three G-tetrads. Thus, both poly(dG) and poly(rG) possess inherent ability to fold into 3D quadruplex architecture with strictly defined folding pattern. The study also revealed that despite high stability of both DNA and RNA quadruplexes, they are vulnerable to single-nucleotide substitutions, which drop the thermal stability by ~40°C and can facilitate introduction of the complementarity principle into the G4 World.

Comprehensive analysis of intramolecular G-quadruplex structures: furthering the understanding of their formalism

G-quadruplexes (G4) are helical structures found in guanine-rich DNA or RNA sequences. Generally, their formalism is based on a few dozen structures, which can produce some inconsistencies or incompleteness. Using the website ASC-G4, we analyzed the structures of 333 intramolecular G4s, of all types, which allowed us to clarify some key concepts and present new information. To each of the eight distinguishable topologies corresponds a groove-width signature and a predominant glycosidic configuration (gc) pattern governed by the directions of the strands. The relative orientations of the stacking guanines within the strands, which we quantified and related to their vertical gc successions, determine the twist and tilt of the helices. The latter impact the minimum groove widths, which represent the space available for lateral ligand binding. The G4 four helices have similar twists, even when these twists are irregular, meaning that they have various angles along the strands. Despite its importance, the vertical gc succession has no strict one-to-one relationship with the topology, which explains the discrepancy between some topologies and their corresponding circular dichroism spectra. This study allowed us to introduce the new concept of platypus G4s, which are structures with properties corresponding to several topologies.

High-Resolution Structure of RNA G-Quadruplex Containing Unique Structural Motifs Originating from the 5'-UTR of Human Tyrosine Kinase 2 (TYK2)

Tyrosine kinase 2 (TYK2) is a member of the JAK family of nonreceptor-associated tyrosine kinases together with highly homologous JAK1, JAK2, and JAK3 paralogues. Overexpression of TYK2 is associated with several inflammatory diseases, including severe complications during the COVID-19 infection. Since the downregulation of JAK paralogues could lead to serious health consequences or even death, it is critical to avoid it when designing drugs to suppress TYK2. To achieve the required specificity only for TYK2, researchers have recently selectively targeted TYK2 mRNA by developing antisense oligonucleotides. In this work, we expand the target space of TYK2 mRNA by showing that the mRNA adopts tetra-helical noncanonical structures called G-quadruplexes. We identified a TYKwt RNA oligonucleotide from the 5'-UTR of TYK2 mRNA, which adopts multiple different parallel G-quadruplexes that exist at equilibrium. Using NMR spectroscopy, we showed that some of the G-quadruplexes adopt unique structural motifs, mainly due to the formation of a stable GA bulge. Using guanine to uridine substitutions, we prepared the oligonucleotide TYK3_U6, which serves as an excellent model for the bulged G-quadruplexes formed by the TYKwt oligonucleotide. NMR structural analysis, including data on the residual coupling constants (RDC) of the loop regions, unveiled that the studied three-quartet parallel G-quadruplex contains many unusual structural features such as a G(U)A bulge, a guanine residue in the syn conformation, A and U residues stacked on the top G-quartet, and a well-defined adenine from a three-residue long propeller loop oriented in the groove, all of which could be valuable targets for future drug design.

The Epigenomic Features and Potential Functions of PEG- and PDS-Favorable DNA G-Quadruplexes in Rice

A G-quadruplex (G4) is a typical non-B DNA structure and involved in various DNA-templated events in eukaryotic genomes. PEG and PDS chemicals have been widely applied for promoting the folding of in vivo or in vitro G4s. However, how PEG and PDS preferentially affect a subset of G4 formation genome-wide is still largely unknown. We here conducted a BG4-based IP-seq in vitro under K++PEG or K++PDS conditions in the rice genome. We found that PEG-favored IP-G4s+ have distinct sequence features, distinct genomic distributions and distinct associations with TEGs, non-TEGs and subtypes of TEs compared to PDS-favored ones. Strikingly, PEG-specific IP-G4s+ are associated with euchromatin with less enrichment levels of DNA methylation but with more enriched active histone marks, while PDS-specific IP-G4s+ are associated with heterochromatin with higher enrichment levels of DNA methylation and repressive marks. Moreover, we found that genes with PEG-specific IP-G4s+ are more expressed than those with PDS-specific IP-G4s+, suggesting that PEG/PDS-specific IP-G4s+ alone or coordinating with epigenetic marks are involved in the regulation of the differential expression of related genes, therefore functioning in distinct biological processes. Thus, our study provides new insights into differential impacts of PEG and PDS on G4 formation, thereby advancing our understanding of G4 biology.

Resolving a guanine-quadruplex structure in the SARS-CoV-2 genome through circular dichroism and multiscale molecular modeling

The genome of SARS-CoV-2 coronavirus is made up of a single-stranded RNA fragment that can assume a specific secondary structure, whose stability can influence the virus's ability to reproduce. Recent studies have identified putative guanine quadruplex sequences in SARS-CoV-2 genome fragments that are involved in coding for both structural and non-structural proteins. In this contribution, we focus on a specific G-rich sequence referred to as RG-2, which codes for the non-structural protein 10 (Nsp10) and assumes a guanine-quadruplex (G4) arrangement. We provide the secondary structure of RG-2 G4 at atomistic resolution by molecular modeling and simulation, validated by the superposition of experimental and calculated electronic circular dichroism spectra. Through both experimental and simulation approaches, we have demonstrated that pyridostatin (PDS), a widely recognized G4 binder, can bind to and stabilize RG-2 G4 more strongly than RG-1, another G4 forming sequence that was previously proposed as a potential target for antiviral drug candidates. Overall, this study highlights RG-2 as a valuable target to inhibit the translation and replication of SARS-CoV-2, paving the way towards original therapeutic approaches against emerging RNA viruses.

Pervasive transcriptome interactions of protein-targeted drugs

The off-target toxicity of drugs targeted to proteins imparts substantial health and economic costs. Proteome interaction studies can reveal off-target effects with unintended proteins; however, little attention has been paid to intracellular RNAs as potential off-targets that may contribute to toxicity. To begin to assess this, we developed a reactivity-based RNA profiling methodology and applied it to uncover transcriptome interactions of a set of Food and Drug Administration-approved small-molecule drugs in vivo. We show that these protein-targeted drugs pervasively interact with the human transcriptome and can exert unintended biological effects on RNA functions. In addition, we show that many off-target interactions occur at RNA loci associated with protein binding and structural changes, allowing us to generate hypotheses to infer the biological consequences of RNA off-target binding. The results suggest that rigorous characterization of drugs' transcriptome interactions may help assess target specificity and potentially avoid toxicity and clinical failures.

In Silico Identification of Potential Quadruplex Forming Sequences in LncRNAs of Cervical Cancer

Long non-coding RNAs (lncRNAs) have emerged as auxiliary regulators of gene expression influencing tumor microenvironment, metastasis and radio-resistance in cancer. The presence of lncRNA in extracellular fluids makes them promising diagnostic markers. LncRNAs deploy higher-order structures to facilitate a complex range of functions. Among such structures, G-quadruplexes (G4s) can be detected or targeted by small molecular probes to drive theranostic applications. The in vitro identification of G4 formation in lncRNAs can be a tedious and expensive proposition. Bioinformatics-driven strategies can provide comprehensive and economic alternatives in conjunction with suitable experimental validation. We propose a pipeline to identify G4-forming sequences, protein partners and biological functions associated with dysregulated lncRNAs in cervical cancer. We identified 17 lncRNA clusters which possess transcripts that can fold into a G4 structure. We confirmed in vitro G4 formation in the four biologically active isoforms of SNHG20, MEG3, CRNDE and LINP1 by Circular Dichroism spectroscopy and Thioflavin-T-assisted fluorescence spectroscopy and reverse-transcriptase stop assay. Gene expression data demonstrated that these four lncRNAs can be potential prognostic biomarkers of cervical cancer. Two approaches were employed for identifying G4 specific protein partners for these lncRNAs and FMR2 was a potential interacting partner for all four clusters. We report a detailed investigation of G4 formation in lncRNAs that are dysregulated in cervical cancer. LncRNAs MEG3, CRNDE, LINP1 and SNHG20 are shown to influence cervical cancer progression and we report G4 specific protein partners for these lncRNAs. The protein partners and G4s predicted in lncRNAs can be exploited for theranostic objectives.

Regulation of thrombin activity by ligand-induced topological alteration in a thrombin-binding aptamer

Thrombin-binding aptamer (TBA), which forms a G-quadruplex (G4) structure with anti-parallel topology, interacts with thrombin to inhibit its enzymatic activity. Here we show that the G4-topology-altering ligand L2H2-2M2EA-6LCO (6LCO) changes the anti-parallel topology of TBA G4 to the parallel topology, thereby abrogating the thrombin-inhibitory activity of TBA. This finding suggests that G4 ligands that alter topology may be promising drug candidates for diseases involving G4-binding proteins.

Unlike RNA-TBA (rTBA), iso-rTBA, the 2'-5'-linked RNA-thrombin-binding aptamer, is a functional equivalent of TBA

An antiparallel, functional RNA G-quadruplex of the 2'-5'-linked thrombin-binding aptamer (iso-rTBA) is reported for the first time. It can inhibit clotting and is remarkably stable to nuclease-degradation, besides having high thermal stability. It is thus, a superior candidate to TBA, rTBA or isoTBA, for further development as an anticoagulant.

Homopurine guanine-rich sequences in complex with N-methyl mesoporphyrin IX form parallel G-quadruplex dimers and display a unique symmetry tetrad

DNA can fold into G-quadruplexes (GQs), non-canonical secondary structures formed by π-π stacking of G-tetrads. GQs are important in many biological processes, which makes them promising therapeutic targets. We identified a 42-nucleotide long, purine-only G-rich sequence from human genome, which contains eight G-stretches connected by A and AAAA loops. We divided this sequence into five unique segments, four guanine stretches each, named GA1-5. In order to investigate the role of adenines in GQ structure formation, we performed biophysical and X-ray crystallographic studies of GA1-5 and their complexes with a highly selective GQ ligand, N-methyl mesoporphyrin IX (NMM). Our data indicate that all variants form parallel GQs whose stability depends on the number of flexible AAAA loops. GA1-3 bind NMM with 1:1 stoichiometry. The Ka for GA1 and GA3 is modest, ∼0.3 μM -1, and that for GA2 is significantly higher, ∼1.2 μM -1. NMM stabilizes GA1-3 by 14.6, 13.1, and 7.0 °C, respectively, at 2 equivalents. We determined X-ray crystal structures of GA1-NMM (1.98 Å resolution) and GA3-NMM (2.01 Å). The structures confirm the parallel topology of GQs with all adenines forming loops and display NMM binding at the 3' G-tetrad. Both complexes dimerize through the 5' interface. We observe two novel structural features: 1) a 'symmetry tetrad' at the dimer interface, which is formed by two guanines from each GQ monomer and 2) a NMM dimer in GA1-NMM. Our structural work confirms great flexibility of adenines as structural elements in GQ formation and contributes greatly to our understanding of the structural diversity of GQs and their modes of interaction with small molecule ligands.

G-quadruplexes in bacteria: insights into the regulatory roles and interacting proteins of non-canonical nucleic acid structures

G-quadruplexes (G4s) are highly stable, non-canonical DNA or RNA structures that can form in guanine-rich stretches of nucleic acids. G4-forming sequences have been found in all domains of life, and proteins that bind and/or resolve G4s have been discovered in both bacterial and eukaryotic organisms. G4s regulate a variety of cellular processes through inhibitory or stimulatory roles that depend upon their positions within genomes or transcripts. These include potential roles as impediments to genome replication, transcription, and translation or, in other contexts, as activators of genome stability, transcription, and recombination. This duality suggests that G4 sequences can aid cellular processes but that their presence can also be problematic. Despite their documented importance in bacterial species, G4s remain understudied in bacteria relative to eukaryotes. In this review, we highlight the roles of bacterial G4s by discussing their prevalence in bacterial genomes, the proteins that bind and unwind G4s in bacteria, and the processes regulated by bacterial G4s. We identify limitations in our current understanding of the functions of G4s in bacteria and describe new avenues for studying these remarkable nucleic acid structures.

A Chimeric DNA/RNA Antiparallel Quadruplex with Improved Stability

Nucleic acid quadruplexes are proposed to play a role in the regulation of gene expression, are often present in aptamers selected for specific binding functions and have potential applications in medicine and biotechnology. Therefore, understanding their structure and thermodynamic properties and designing highly stable quadruplexes is desirable for a variety of applications. Here, we evaluate DNA→RNA substitutions in the context of a monomolecular, antiparallel quadruplex, the thrombin-binding aptamer (TBA, GGTTGGTGTGGTTGG) in the presence of either K+ or Sr2+ . TBA predominantly folds into a chair-type configuration containing two G-tetrads, with G residues in both syn and anti conformation. All chimeras with DNA→RNA substitutions (G→g) at G residues requiring the syn conformation demonstrated strong destabilization. In contrast, G→g substitutions at Gs with anti conformation increased stability without affecting the monomolecular chair-type topology. None of the DNA→RNA substitutions in loop positions affected the quadruplex topology; however, these substitutions varied widely in their stabilizing or destabilizing effects in an unpredictable manner. This analysis allowed us to design a chimeric DNA/RNA TBA construct that demonstrated substantially improved stability relative to the all-DNA construct. These results have implications for a variety of quadruplex-based applications including for the design of dynamic nanomachines.

Action and function of helicases on RNA G-quadruplexes

•

Methodological progresses and piling evidence prove the rG4 biology in vivo.

•

rG4s step in virtually every aspect of RNA biology.

•

Helicases unwinding of rG4s is a fine regulatory layer to the downstream processes and general cell homeostasis.

•

The current knowledge is however limited to a few cell lines.

•

The regulation of helicases themselves is delineating as a important question.

•

Non-helicase rG4-processing proteins likely play a role.

The nucleic acid structure called G-quadruplex (G4) is currently discussed to function in nucleic acid-based mechanisms that influence several cellular processes. They can modulate the cellular machinery either positively or negatively, both at the DNA and RNA level. The majority of what we know about G4 biology comes from DNA G4 (dG4) research. RNA G4s (rG4), on the other hand, are gaining interest as researchers become more aware of their role in several aspects of cellular homeostasis. In either case, the correct regulation of G4 structures within cells is essential and demands specialized proteins able to resolve them. Small changes in the formation and unfolding of G4 structures can have severe consequences for the cells that could even stimulate genome instability, apoptosis or proliferation. Helicases are the most relevant negative G4 regulators, which prevent and unfold G4 formation within cells during different pathways. Yet, and despite their importance only a handful of rG4 unwinding helicases have been identified and characterized thus far. This review addresses the current knowledge on rG4s-processing helicases with a focus on methodological approaches. An example of a non-helicase rG4s-unwinding protein is also briefly described.

Selective Binding and Redox-Activity on Parallel G-Quadruplexes by Pegylated Naphthalene Diimide-Copper Complexes

G-quadruplexes (G4s) are higher-order supramolecular structures, biologically important in the regulation of many key processes. Among all, the recent discoveries relating to RNA-G4s, including their potential involvement as antiviral targets against COVID-19, have triggered the ever-increasing need to develop selective molecules able to interact with parallel G4s. Naphthalene diimides (NDIs) are widely exploited as G4 ligands, being able to induce and strongly stabilize these structures. Sometimes, a reversible NDI-G4 interaction is also associated with an irreversible one, due to the cleavage and/or modification of G4s by functional-NDIs. This is the case of NDI-Cu-DETA, a copper(II) complex able to cleave G4s in the closest proximity to the target binding site. Herein, we present two original Cu(II)-NDI complexes, inspired by NDI-Cu-DETA, differently functionalized with 2-(2-aminoethoxy)ethanol side-chains, to selectively drive redox-catalyzed activity towards parallel G4s. The selective interaction toward parallel G4 topology, controlled by the presence of 2-(2-aminoethoxy)ethanol side chains, was already firmly demonstrated by us using core-extended NDIs. In the present study, the presence of protonable moieties and the copper(II) cavity, increases the binding affinity and specificity of these two NDIs for a telomeric RNA-G4. Once defined the copper coordination relationship and binding constants by competition titrations, ability in G4 stabilization, and ROS-induced cleavage were analyzed. The propensity in the stabilization of parallel topology was highlighted for both of the new compounds HP2Cu and PE2Cu. The results obtained are particularly promising, paving the way for the development of new selective functional ligands for binding and destructuring parallel G4s.

Beyond small molecules: targeting G-quadruplex structures with oligonucleotides and their analogues

G-Quadruplexes (G4s) are widely studied secondary DNA/RNA structures, naturally occurring when G-rich sequences are present. The strategic localization of G4s in genome areas of crucial importance, such as proto-oncogenes and telomeres, entails fundamental implications in terms of gene expression regulation and other important biological processes. Although thousands of small molecules capable to induce G4 stabilization have been reported over the past 20 years, approaches based on the hybridization of a synthetic probe, allowing sequence-specific G4-recognition and targeting are still rather limited. In this review, after introducing important general notions about G4s, we aim to list, explain and critically analyse in more detail the principal approaches available to target G4s by using oligonucleotides and synthetic analogues such as Locked Nucleic Acids (LNAs) and Peptide Nucleic Acids (PNAs), reporting on the most relevant examples described in literature to date.

RNA G-quadruplexes (rG4s): genomics and biological functions

G-quadruplexes (G4s) are non-classical DNA or RNA secondary structures that have been first observed decades ago. Over the years, these four-stranded structural motifs have been demonstrated to have significant regulatory roles in diverse biological processes, but challenges remain in detecting them globally and reliably. Compared to DNA G4s (dG4s), the study of RNA G4s (rG4s) has received less attention until recently. In this review, we will summarize the innovative high-throughput methods recently developed to detect rG4s on a transcriptome-wide scale, highlight the many novel and important functions of rG4 being discovered in vivo across the tree of life, and discuss the key biological questions to be addressed in the near future.

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.

HIV-1 genomic RNA U3 region forms a stable quadruplex-hairpin structure

The U3 promoter region of the HIV-1 long terminal repeat (LTR) has previously been shown to fold into a series of dynamic G-quadruplex structures. Among the G-quadruplexes identified in the LTR sequence, LTR-III was shown to be the most stable in vitro. NMR studies of this 28-nucleotide (nt) DNA revealed a unique quadruplex-hairpin structure. Whether the hairpin forms in RNA element is unknown and the role of the hairpin in the structure and stability of quadruplexes has not been characterized. Here, we used optical and thermodynamic studies to address these questions. The wild-type LTR-III RNA formed a monomolecular quadruplex with a parallel topology using only propeller loops, including the hairpin loop element. By comparison to the WT and variant RNAs, LTR-III DNA structures were more heterogeneous and less stable. Increased stability of the RNA suggests that the RNA quadruplex-hairpin structure may be a more attractive therapeutic target than the analogous DNA element.

Proximal Single-Stranded RNA Destabilizes Human Telomerase RNA G-Quadruplex and Induces Its Distinct Conformers

Single-stranded guanine-rich RNA sequences have a propensity to fold into compact G-quadruplexes (RG4s). The conformational transitions of these molecules provide an important way to regulate their biological functions. Here, we examined the stability and conformation of an RG4-forming sequence identified near the end of human telomerase RNA. We found that a proximal single-stranded (ss) RNA significantly impairs RG4 stability at physiological K⁺ concentrations, resulting in a reduced RG4 rupture force of ∼ 24.4 pN and easier accessibility of the G-rich sequence. The destabilizing effect requires a minimum of six nucleotides of ssRNA and is effective at either end of RG4. Remarkably, this RG4-forming sequence, under the influence of such a proximal ssRNA, exhibits interconversions between at least three less stable RG4 conformers that might represent potential intermediates along its folding/unfolding pathway. This work provides insights into the stability and folding dynamics of RG4 that are essential for understanding its biological functions.

Recent progress in non-native nucleic acid modifications

While Nature harnesses RNA and DNA to store, read and write genetic information, the inherent programmability, synthetic accessibility and wide functionality of these nucleic acids make them attractive tools for use in a vast array of applications.

Improving Thermodynamic Stability and Anticoagulant Activity of a Thrombin Binding Aptamer by Incorporation of 8-trifluoromethyl-2'-deoxyguanosine

In this study, we incorporated 8-trifluoromethyl-2'-deoxyguanosine (^FG) into a thrombin binding aptamer (TBA). Circular dichroism, nuclear magnetic resonance (NMR), electrophoresis, and prothrombin time (PT) assay were performed to investigate the structure, thermodynamic stability, biological stability, and anticoagulant activity of the ^FG-modified TBA sequences. We found that the replacement of ^FG into TBA sequences led to a remarkable improvement in the melting temperature up to 30 °C compared with the native sequence. The trifluoromethyl group allowed us to investigate the TBA G-quadruplex structure by ¹⁹F NMR spectroscopy. Furthermore, PT assays showed that the modified sequences can significantly improve the anticoagulant activity in comparison with the native TBA. Finally, we demonstrated that the trifluoromethyl-modified TBA sequence could function as an anticoagulant reagent in live rats. Our results strongly suggested that ^FG is a powerful nucleoside derivative to increase the thermodynamic stability and anticoagulant activity of TBA.

Assembly and Characterization of RNA/DNA Hetero-G-Quadruplexes

Transient association of guanine-rich RNA and DNA in the form of hetero-G-quadruplexes (RDQs) has emerged as an important mechanism for regulating genome transcription and replication but relatively little is known about the structure and biophysical properties of RDQs compared with DNA and RNA homo-G-quadruplexes. Herein, we report the assembly and characterization of three RDQs based on sequence motifs found in human telomeres and mitochondrial nucleic acids. Stable RDQs were assembled using a duplex scaffold, which prevented segregation of the DNA and RNA strands into separate homo-GQs. Each of the RDQs exhibited UV melting temperatures above 50 °C in 100 mM KCl and predominantly parallel morphologies, evidently driven by the RNA component. The fluorogenic dye thioflavin T binds to each RDQ with low micromolar K_D values, similar to its binding to RNA and DNA homo-GQs. These results establish a method for assembling RDQs that should be amenable to screening compounds and libraries to identify selective RDQ-binding small molecules, oligonucleotides, and proteins.

Conformation of G-quadruplex Controlled by Click Reaction

G-quadruplexes are non-canonical four stranded secondary structures possessing great biological importance. Controlling G-quadruplex conformation for further regulating biological processes is both exciting and challenging. In this study, we described a method for regulating G-quadruplex conformation by click chemistry for the first time. 8-ethynyl-2'-deoxyguanosine was synthesized and incorporated into a 12-nt telomere DNA sequence. Such a sequence, at first, formed mixed parallel/anti-parallel G-quadruplexes, while it changed to anti-parallel after reaction with azidobenzene. Meanwhile, the click reaction can give the sequence intense fluorescence.

G-quadruplex-based aptamers targeting human thrombin: Discovery, chemical modifications and antithrombotic effects

First studies on thrombin-inhibiting DNA aptamers were reported in 1992, and since then a large number of anticoagulant aptamers has been discovered. TBA - also named HD1, a 15-mer G-quadruplex (G4)-forming oligonucleotide - is the best characterized thrombin binding aptamer, able to specifically recognize the protein exosite I, thus inhibiting the conversion of soluble fibrinogen into insoluble fibrin strands. Unmodified nucleic acid-based aptamers, in general, and TBA in particular, exhibit limited pharmacokinetic properties and are rapidly degraded in vivo by nucleases. In order to improve the biological performance of aptamers, a widely investigated strategy is the introduction of chemical modifications in their backbone at the level of the nucleobases, sugar moieties or phosphodiester linkages. Besides TBA, also other thrombin binding aptamers, able to adopt a well-defined G4 structure, e.g. mixed duplex/quadruplex sequences, as well as homo- and hetero-bivalent constructs, have been identified and optimized. Considering the growing need of new efficient anticoagulant agents associated with the strong therapeutic potential of these thrombin inhibitors, the research on thrombin binding aptamers is still a very hot and intriguing field. Herein, we comprehensively described the state-of-the-art knowledge on the DNA-based aptamers targeting thrombin, especially focusing on the optimized analogues obtained by chemically modifying the oligonucleotide backbone, and their biological performances in therapeutic applications.

G-quadruplexes offer a conserved structural motif for NONO recruitment to NEAT1 architectural lncRNA

The long non-coding RNA NEAT1 serves as a scaffold for the assembly of paraspeckles, membraneless nuclear organelles involved in gene regulation. Paraspeckle assembly requires NEAT1 recruitment of the RNA-binding protein NONO, however the NEAT1 elements responsible for recruitment are unknown. Herein we present evidence that previously unrecognized structural features of NEAT1 serve an important role in these interactions. Led by the initial observation that NONO preferentially binds the G-quadruplex conformation of G-rich C9orf72 repeat RNA, we find that G-quadruplex motifs are abundant and conserved features of NEAT1. Furthermore, we determine that NONO binds NEAT1 G-quadruplexes with structural specificity and provide evidence that G-quadruplex motifs mediate NONO-NEAT1 association, with NONO binding sites on NEAT1 corresponding largely to G-quadruplex motifs, and treatment with a G-quadruplex-disrupting small molecule causing dissociation of native NONO-NEAT1 complexes. Together, these findings position G-quadruplexes as a primary candidate for the NONO-recruiting elements of NEAT1 and provide a framework for further investigation into the role of G-quadruplexes in paraspeckle formation and function.

Epigenetic Modulation of Chromatin States and Gene Expression by G-Quadruplex Structures

G-quadruplexes are four-stranded helical nucleic acid structures formed by guanine-rich sequences. A considerable number of studies have revealed that these noncanonical structural motifs are widespread throughout the genome and transcriptome of numerous organisms, including humans. In particular, G-quadruplexes occupy strategic locations in genomic DNA and both coding and noncoding RNA molecules, being involved in many essential cellular and organismal functions. In this review, we first outline the fundamental structural features of G-quadruplexes and then focus on the concept that these DNA and RNA structures convey a distinctive layer of epigenetic information that is critical for the complex regulation, either positive or negative, of biological activities in different contexts. In this framework, we summarize and discuss the proposed mechanisms underlying the functions of G-quadruplexes and their interacting factors. Furthermore, we give special emphasis to the interplay between G-quadruplex formation/disruption and other epigenetic marks, including biochemical modifications of DNA bases and histones, nucleosome positioning, and three-dimensional organization of chromatin. Finally, epigenetic roles of RNA G-quadruplexes in post-transcriptional regulation of gene expression are also discussed. Undoubtedly, the issues addressed in this review take on particular importance in the field of comparative epigenetics, as well as in translational research.

A synthetic Pur-based peptide binds and alters G-quadruplex secondary structure present in the expanded RNA repeat of C9orf72 ALS/FTD

Increased Pur-alpha (Pura) protein levels in animal models alleviate certain cellular symptoms of the disease spectrum amyotrophic lateral sclerosis/frontotemporal dementia (ALS/FTD). Pura is a member of the Pur family of evolutionarily conserved guanine-rich polynucleotide binding proteins containing a repeated signature PUR domain of 60-80 amino acids. Here we have employed a synthetic peptide, TZIP, similar to a Pur domain, but with sequence alterations based on a consensus of evolutionarily conserved Pur family binding domains and having an added transporter sequence. A major familial form of ALS/FTD, C9orf72 (C9), is due to a hexanucleotide repeat expansion (HRE) of (GGGGCC), a Pur binding element. We show by circular dichroism that RNA oligonucleotides containing this purine-rich sequence consist largely of parallel G-quadruplexes. TZIP peptide binds this repeat sequence in both DNA and RNA. It binds the RNA element, including the G-quadruplexes, with a high degree of specificity versus a random oligonucleotide. In addition, TZIP binds both linear and G-quadruplex repeat RNA to form higher order G-quadruplex secondary structures. This change in conformational form by Pur-based peptide represents a new mechanism for regulating G quadruplex secondary structure within the C9 repeat. TZIP modulation of C9 RNA structural configuration may alter interaction of the complex with other proteins. This Pur-based mechanism provides new targets for therapy, and it may help to explain Pura alleviation of certain cellular pathological aspects of ALS/FTD.

Nanomechanics and co-transcriptional folding of Spinach and Mango

Recent advances in fluorogen-binding “light-up” RNA aptamers have enabled protein-free detection of RNA in cells. Detailed biophysical characterization of folding of G-Quadruplex (GQ)-based light-up aptamers such as Spinach, Mango and Corn is still lacking despite the potential implications on their folding and function. In this work we employ single-molecule fluorescence-force spectroscopy to examine mechanical responses of Spinach2, iMangoIII and MangoIV. Spinach2 unfolds in four discrete steps as force is increased to 7 pN and refolds in reciprocal steps upon force relaxation. In contrast, GQ-core unfolding in iMangoIII and MangoIV occurs in one discrete step at forces >10 pN and refolding occurred at lower forces showing hysteresis. Co-transcriptional folding using superhelicases shows reduced misfolding propensity and allowed a folding pathway different from refolding. Under physiologically relevant pico-Newton levels of force, these aptamers may unfold in vivo and subsequently misfold. Understanding of the dynamics of RNA aptamers will aid engineering of improved fluorogenic modules for cellular applications.

Light-up aptamers are widely used for fluorescence visualization of non-coding RNA in vivo. Here the authors employ single-molecule fluorescence-force spectroscopy to characterize the mechanical responses of the G-Quadruplex based light-up aptamers Spinach2, iMangoIII and MangoIV, which is of interest for the development of improved fluorogenic modules for imaging applications.

MOV10L1 Binds RNA G-Quadruplex in a Structure-Specific Manner and Resolves It More Efficiently Than MOV10

MOV10L1 and its paralog MOV10 are evolutionally conserved RNA helicases involved in distinct RNA regulatory pathways. The testis-specific MOV10L1 is essential for spermatogenesis and PIWI-interacting RNAs biogenesis, whereas MOV10 is ubiquitous and multifunctional. Although both proteins have been implied to correlate with RNA G-quadruplex (RG4) in vivo, their capabilities in binding and resolving RG4 and their respective biological significance remain unclear. Herein, we comprehensively characterize and compare the activities of these two helicases on various nucleic acid substrates in vitro, with a focus on RG4 structure. We find that both MOV10L1 and MOV10 are able to resolve RG4, with MOV10L1 being more efficient in that. In contrast to MOV10, MOV10L1 prefers to bind to a junction between single-stranded RNA and RG4, which is mediated by both its N and C termini. Furthermore, we show that RG4 unwinding by MOV10L1 facilitates the cleavage of this specific RNA structure by an endonuclease.

•

Both MOV10L1 and MOV10 can resolve RG4 structure in an ATP-dependent manner

•

MOV10L1 unwinds RG4 more efficiently than MOV10

•

MOV10L1 preferentially binds to an ssRNA-RG4 junction

•

RG4 unwinding by MOV10L1 facilitates its endonucleolytic cleavage

Manipulating DNA G-Quadruplex Structures by Using Guanosine Analogues

The ability to control the folding topology of DNA G-quadruplexes allows for rational design of quadruplex-based scaffolds for potential use in various therapeutic and technological applications. By exploiting the distinct conformational properties of some base- and sugar-modified guanosine surrogates, conformational transitions can be induced through their judicious incorporation at specific sites in the quadruplex core. Changes may involve tetrad polarity inversions with conservation of the global fold or complete refolding to new topologies. Reliable predictions relating to low-energy conformers formed upon specific chemical perturbations of the system and the rational design of modified sequences suffer from our still limited understanding of the subtle interplay of various favorable and unfavorable interactions within a particular quadruplex scaffold. However, aided by an increasing number of systematic substitution experiments and high-resolution structures of modified quadruplex variants, critical interactions, in addition to glycosidic bond angle propensities, are starting to emerge as important contributors to modification-driven quadruplex refolding.

Stability Factors of the Parallel Quadruplexes: DNA Versus RNA

One of the most stable quadruplexes is formed by the G3T sequence (GGGTGGGTGGGTGGG) that folds into a parallel quadruplex with three G-tetrads and chain-reversal T-loops. For example, in 1 mM K⁺, it unfolds at 75 °C and at physiological conditions, it unfolds above 100 °C. The RNA analogue, ggguggguggguggg (g3u), which employs exactly same folding topology, demonstrates even higher thermal stability. Here, we performed melting experiments of G3T, g3u, and more than 30 chimeric constructs (G3T with RNA nucleotides at certain positions). Although the g3u quadruplex is 13 °C more stable than G3T, majority of G → g (DNA-for-RNA) substitutions destabilize G3T. Only three G → g and loop T → u substitutions stabilize the structure. However, stabilization effects of these six substitutions overcome destabilization of other nine G → g, resulting in higher stability of all-RNA g3u. The present work clearly indicates that the stacking interactions are more favorable in parallel DNA quadruplexes, whereas the chain-reversal loops play an important role in higher stability of RNA quadruplexes. In addition, we have shown that the 5'-end of RNA quadruplexes represents a more favorable target for stacking interactions than the 3'-end. Based on the current study, rational design of the quadruplexes for particular biotechnological applications and drugs, targeting the quadruplexes, may be envisaged.

α-2′-Deoxyguanosine can switch DNA G-quadruplex topologies from antiparallel to parallel

α-2′-Deoxyguanosine (α-dG) converts antiparallel, dimeric G-quadruplex DNA into a parallel, tetramolecular complex.

G-quadruplexes formation in the 5'UTRs of mRNAs associated with colorectal cancer pathways

RNA G-quadruplexes (rG4) are stable non-canonical secondary structures composed of G-rich sequences. Many rG4 structures located in the 5'UTRs of mRNAs act as translation repressors due to their high stability which is thought to impede ribosomal scanning. That said, it is not known if these are mRNA-specific examples, or if they are indicative of a global expression regulation mechanism of the mRNAs involved in a common pathway based on structure folding recognition. Gene-ontology analysis of mRNAs bearing a predicted rG4 motif in their 5'UTRs revealed an enrichment for mRNAs associated with the colorectal cancer pathway. Bioinformatic tools for rG4 prediction, and experimental in vitro validations were used to confirm and compare the folding of the predicted rG4s of the mRNAs associated with dysregulated pathways in colorectal cancer. The rG4 folding was confirmed for the first time for 9 mRNAs. A repressive effect of 3 rG4 candidates on the expression of a reporter gene was also measured in colorectal cancer cell lines. This work highlights the fact that rG4 prediction is not yet accurate, and that experimental characterization is still essential in order to identify the precise rG4 folding sequences and the possible common features shared between the rG4 overrepresented in important biological pathways.

Structuring polarity-inverted TBA to G-quadruplex for selective recognition of planarity of natural isoquinoline alkaloids

Efficient structuring of DNA by small molecules is very crucial in developing DNA-based novel switches with an ideal performance. In this work, we found that inverting only the polarity of the 3' terminal guanine of the thrombin-binding aptamer (3iTBA) totally eradicates the original TBA G-quadruplex (G4) structure in K+. The unstructured 3iTBA can be further refolded upon specifically interacting with small molecules of natural isoquinoline alkaloids (IAs) due to their fruitful binding patterns with variant nucleic acid structures. We identified that 3iTBA can serve as a topology selector for planar IAs. Nitidine (NIT), owing to the planar aromatic ring and coplanar substituents, is the most efficient to restructure the 3iTBA random coil toward the anti-parallel G4 conformation. However, common metal ions can't realize this structuring. The topology selector competency of 3iTBA toward IAs' planarity can be visualized using gold nanoparticles (AuNPs) as the chromogenic readout. Our work expands the G4 repertoire by exploring the polarity inversion regulation and provides a new approach to switch nucleic acid structures toward a small molecule structure-sensitive sensor.

Encoding canonical DNA quadruplex structure

The main challenge in DNA quadruplex design is to encode a three-dimensional structure into the primary sequence, despite its multiple, repetitive guanine segments. We identify and detail structural elements describing all 14 feasible canonical quadruplex scaffolds and demonstrate their use in control of design. This work outlines a new roadmap for implementation of targeted design of quadruplexes for material, biotechnological, and therapeutic applications.

Fluorescence-based tools to probe G-quadruplexes in cell-free and cellular environments

Biophysical and biochemical investigations provide compelling evidence connecting the four-stranded G-quadruplex (GQ) structure with its role in regulating multiple cellular processes. Hence, modulating the function of GQs by using small molecule binders is being actively pursued as a strategy to develop new chemotherapeutic agents. However, sequence diversity and structural polymorphism of GQs have posed immense challenges in terms of understanding what conformation a G-rich sequence adopts inside the cell and how to specifically target a GQ motif amidst several other GQ-forming sequences. In this context, here we review recent developments in the applications of biophysical tools that use fluorescence readout to probe the GQ structure and recognition in cell-free and cellular environments. First, we provide a detailed discussion on the utility of covalently labeled environment-sensitive fluorescent nucleoside analogs in assessing the subtle difference in GQ structures and their ligand binding abilities. Furthermore, a detailed discussion on structure-specific antibodies and small molecule probes used to visualize and confirm the existence of DNA and RNA GQs in cells is provided. We also highlight the open challenges in the study of tetraplexes (GQ and i-motif structures) and how addressing these challenges by developing new tools and techniques will have a profound impact on tetraplex-directed therapeutic strategies.

Monomolecular tetrahelix of polyguanine with a strictly defined folding pattern

The G₃TG₃TG₃TG₃ (G3T) sequence folds into a monomolecular quadruplex with all-parallel G₃ segments connected to each other by chain-reversal loops. The homopolymer consisting of n number of G3T domains directly conjugated to each other folds into an uninterrupted and unusually stable polymer, tetrahelical monomolecular DNA (tmDNA). It was demonstrated that the tmDNA architecture has strong potential in nanotechnologies as highly programmable building material, high affinity coupler and the driving force for endergonic reactions. Here, we explore capability of analogous DNA sequences (i.e., monomolecular quadruplexes with G₂ or G₄ segments) to construct tmDNA architecture. The study demonstrates that tmDNA can have only one building pattern based on a quadruplex domain with three G-tetrads and single-nucleotide loops, G3N (N = G, A, C and T); all other domains demonstrate antiparallel topologies unsuitable for tmDNA. The present study also suggests that polyguanine is capable of tmDNA formation with strictly defined building pattern; G₃ segments connected to each other by chain-reversal G-loops. These findings can have significant impact on (i) DNA nanotechnologies; (ii) structure prediction of G-rich sequences of genome; and (iii) modeling of abiogenesis.

Structural and Dynamical Characterization of DNA and RNA Quadruplexes Obtained from the GGGGCC and GGGCCT Hexanucleotide Repeats Associated with C9FTD/ALS and SCA36 Diseases

A (GGGGCC) hexanucleotide repeat (HR) expansion in the C9ORF72 gene has been considered the major cause behind both frontotemporal dementia and amyotrophic lateral sclerosis, while a (GGGCCT) is associated with spinocerebellar ataxia 36. Recent experiments involving NMR, CD, optical melting and 1D 1H NMR spectroscopy, suggest that the r(GGGGCC) HR can adopt a hairpin structure with G-G mismatches in equilibrium with a G-quadruplex structure. G-Quadruplexes have also been identified for d(GGGGCC). As these experiments lack molecular resolution, we have used molecular dynamics microsecond simulations to obtain a structural characterization of the G-quadruplexes associated with both HRs. All DNA G-quadruplexes, parallel or antiparallel, with or without loops are stable, while only parallel and one antiparallel (stabilized by diagonal loops) RNA G-quadruplexes are stable. It is known that antiparallel G-quadruplexes require alternating guanines to be in a syn conformation that is hindered by the C3'-endo pucker preferred by RNA. Initial RNA antiparallel quadruplexes built with C2'-endo sugars evolve such that the transition (C2'-endo)-to-(C3'-endo) triggers unwinding and buckling of the flat G-tetrads, resulting in the unfolding of the RNA antiparallel quadruplex. Finally, a parallel G-quadruplex stabilizes an adjacent C-tetrad in both DNA and RNA (thus effectively becoming a mixed quadruplex of 5 layers). The C-tetrad is stabilized by the stacking interactions with the preceding G-tetrad, by cyclical hydrogen bonds C(N4)-(O2), and by an ion between the G-tetrad and the C-tetrad. In addition, antiparallel DNA G-quadruplexes also stabilize flat C-layers at the ends of the quadruplexes.

Thermodynamic, Anticoagulant, and Antiproliferative Properties of Thrombin Binding Aptamer Containing Novel UNA Derivative

Thrombin is a serine protease that plays a crucial role in hemostasis, fibrinolysis, cell proliferation, and migration. Thrombin binding aptamer (TBA) is able to inhibit the activity of thrombin molecule via binding to its exosite I. This 15-nt DNA oligonucleotide forms an intramolecular, antiparallel G-quadruplex structure with a chair-like conformation. In this paper, we report on our investigations on the influence of certain modified nucleotide residues on thermodynamic stability, folding topology, and biological properties of TBA variants. In particular, the effect of single incorporation of a novel 4-thiouracil derivative of unlocked nucleic acid (UNA), as well as single incorporation of 4-thiouridine and all four canonical UNAs, was evaluated. The studies presented herein have shown that 4-thiouridine in RNA and UNA series, as well as all four canonical UNAs, can efficiently modulate G-quadruplex thermodynamic and biological stability, and that the effect is strongly position dependent. Interestingly, TBA variants containing the modified nucleotide residues are characterized by unchanged folding topology. Thrombin time assay revealed that incorporation of certain UNA residues may improve G-quadruplex anticoagulant properties. Noteworthy, some TBA variants, characterized by decreased ability to inhibit thrombin activity, possess significant antiproliferative properties reducing the viability of the HeLa cell line even by 95% at 10 μM concentration.

Nonconventional C-H···F Hydrogen Bonds Support a Tetrad Flip in Modified G-Quadruplexes

A G-quadruplex adopting a (3 + 1)-hybrid structure was substituted at its 5'-tetrad by 2'-deoxy-2'-fluoro-arabinoguanosine (^FaraG) analogs. Incorporation of anti-favoring ^FaraG at syn-positions of the 5'-outer tetrad induced a reversal of the tetrad polarity without noticeably compromising the quadruplex stability. This conformational change is shown to be promoted by nonconventional C-H···F hydrogen bonds acting within the anti-^FaraG residues.

In What Ways Do Synthetic Nucleotides and Natural Base Lesions Alter the Structural Stability of G-Quadruplex Nucleic Acids?

Synthetic analogs of natural nucleotides have long been utilized for structural studies of canonical and noncanonical nucleic acids, including the extensively investigated polymorphic G-quadruplexes (GQs). Dependence on the sequence and nucleotide modifications of the folding landscape of GQs has been reviewed by several recent studies. Here, an overview is compiled on the thermodynamic stability of the modified GQ folds and on how the stereochemical preferences of more than 70 synthetic and natural derivatives of nucleotides substituting for natural ones determine the stability as well as the conformation. Groups of nucleotide analogs only stabilize or only destabilize the GQ, while the majority of analogs alter the GQ stability in both ways. This depends on the preferred syn or anti N-glycosidic linkage of the modified building blocks, the position of substitution, and the folding architecture of the native GQ. Natural base lesions and epigenetic modifications of GQs explored so far also stabilize or destabilize the GQ assemblies. Learning the effect of synthetic nucleotide analogs on the stability of GQs can assist in engineering a required stable GQ topology, and exploring the in vitro action of the single and clustered natural base damage on GQ architectures may provide indications for the cellular events.

Targeting of Polycomb Repressive Complex 2 to RNA by Short Repeats of Consecutive Guanines

Polycomb repressive complex 2 (PRC2) is a histone methyltransferase that trimethylates H3K27, a mark of repressed chromatin. Mammalian PRC2 binds RNA promiscuously, with thousands of target transcripts in vivo. But what does PRC2 recognize in these RNAs? Here we show that purified human PRC2 recognizes G > C,U ≫ A in single-stranded RNA and has a high affinity for folded guanine quadruplex (G4) structures but little binding to duplex RNAs. Importantly, G-tract motifs are significantly enriched among PRC2-binding transcripts in vivo. DNA sequences coding for PRC2-binding RNA motifs are enriched at PRC2-binding sites on chromatin and H3K27me3-modified nucleosomes. Collectively, the abundance of PRC2-binding RNA motifs rationalizes the promiscuous RNA binding of PRC2, and their enrichment at Polycomb target genes provides a means for RNA-mediated regulation.

Mutagenesis Reveals an Unusual Combination of Guanines in RNA G-Quadruplex Formation

The classic G-quadruplex motif consists of a continuous array of 3-4 guanine residues with an intermittent loop size of 1-7 nucleotides (G_3-4N_1-7G_3-4N_1-7G_3-4N_1-7G_3-4). An RNA G-quadruplex is able to attain only one parallel G-quadruplex topology owing to steric constraints. Investigating the possibilities of the formation of RNA G-quadruplexes with a stretch of sequences deviating from this classic motif will add to the overall conformations of RNA G-quadruplexes, bestowing diversity to this structure. Here, we report unusual combinations of guanine residues involved in RNA G-quadruplex formation in the 5' untranslated region (UTR) of the von Willebrand factor (VWF) mRNA using the mutagenesis approach. Different permutations and combinations of guanine residues form G-quadruplexes. Upon investigation, G-quadruplexes in 5' UTR of VWF mRNA are shown to exhibit an inhibitory function.

Antiparallel RNA G-quadruplex Formed by Human Telomere RNA Containing 8-Bromoguanosine

In this study, by combining nuclear magnetic resonance (NMR), circular dichroism (CD), liquid chromatography-electrospray ionization-mass spectrometry (LC-ESI-MS), and gel electrophoresis, we report an unusual topological structure of the RNA G-quadruplex motif formed by human telomere RNA r(UAGGGU) containing 8-bromoguanosine. Results showed that the RNA sequence formed an antiparallel tetramolecular G-quadruplex, in which each pair of diagonal strands run in opposite directions. Furthermore, guanosines were observed both in syn- and anti-conformations. In addition, two of these G-quadruplex subunits were found to be stacking on top of each other, forming a dimeric RNA G-quadruplex. Our findings provide a new insight into the behavior of RNA G-quadruplex structures.

Stable Domain Assembly of a Monomolecular DNA Quadruplex: Implications for DNA-Based Nanoswitches

In the presence of K⁺ ions, the 5′-GGGTGGGTGGGTGGG-3′ (G3T) sequence folds into a monomolecular quadruplex with unusually high thermal stability and unique optical properties. In this study we report that although single G3T molecules unfold and fold rapidly with overlapping melting and refolding curves, G3T multimers (G3T units covalently attached to each other) demonstrate highly reproducible hysteretic behavior. We demonstrate that this behavior necessitates full-length tandem G3T monomers directly conjugated to each other. Any modification of the tandem sequences eliminates the hysteresis. The experimentally measured kinetic parameters and equilibrium transition profiles suggest a highly specific two-state transition in which the folding and unfolding of the first G3T monomer is rate-limiting for both annealing and melting processes. The highly reproducible hysteretic behavior of G3T multimers has the potential to be used in the design of heat-stimulated DNA switches or transistors.

Molecular recognition of a carboxy pyridostatin toward G-quadruplex structures: Why does it prefer RNA?

The pyridostatin (PDS) represents the lead compound of a family of G-quadruplex (G4) stabilizing synthetic small molecules based on a N,N'-bis(quinolinyl)pyridine-2,6-dicarboxamide scaffold. Its mechanism of action involves the induction of telomere dysfunction by competing for binding with telomere-associated proteins, such as human POT1. Recently, through a template-directed "in situ" click chemistry approach, a PDS derivative, the carboxypyridostatin (cPDS), was discovered. It has the peculiarity to exhibit high molecular specificity for RNA over DNA G4, while PDS is a good generic RNA and DNA G4-interacting small molecule. Structural data on the binding modes of these compounds are not available, and the selectivity mode of cPDS toward TERRA G4 is unknown too. Therefore, this work is aimed at rationalizing the selectivity of cPDS versus TERRA G4 by means of molecular dynamics and docking simulations, coupled to better understand the binding mode of these compounds to telomeric G4 structures. The comprehensive analysis of cPDS binding mode and its conformational behavior demonstrates the importance of the ligand conformation properties coupled with a remarkable solvation contribution. This work is expected to provide valuable clues for further rational design of novel and selective TERRA G4 binders.