Encoding canonical DNA quadruplex structure

Unlocking the potential of protein-derived peptides to target G-quadruplex DNA: from recognition to anticancer activity

Noncanonical nucleic acid structures, particularly G-quadruplexes, have garnered significant attention as potential therapeutic targets in cancer treatment. Here, the recognition of G-quadruplex DNA by peptides derived from the Rap1 protein is explored, with the aim of developing novel peptide-based G-quadruplex ligands with enhanced selectivity and anticancer activity. Biophysical techniques were employed to assess the interaction of a peptide derived from the G-quadruplex-binding domain of the protein with various biologically relevant G-quadruplex structures. Through alanine scanning mutagenesis, key amino acids crucial for G-quadruplex recognition were identified, leading to the discovery of two peptides with improved G-quadruplex-binding properties. However, despite their in vitro efficacy, these peptides showed limited cell penetration and anticancer activity. To overcome this challenge, cell-penetrating peptide (CPP)-conjugated derivatives were designed, some of which exhibited significant cytotoxic effects on cancer cells. Interestingly, selected CPP-conjugated peptides exerted potent anticancer activity across various tumour types via a G-quadruplex-dependent mechanism. These findings underscore the potential of peptide-based G-quadruplex ligands in cancer therapy and pave the way for the development of novel therapeutic strategies targeting these DNA structures.

Structural Basis for Parallel G-Quadruplex Recognition by an Ankyrin Protein

G-Quadruplex (G4) structures formed by guanine-rich DNA and RNA sequences are implicated in various biological processes. Understanding the mechanisms by which proteins recognize G4 structures is crucial for elucidating their functional roles. Here we present the X-ray crystal structure of an ankyrin protein bound to a parallel G4 structure. Our findings reveal a new specific recognition mode in which a bundle of α-helices and loops of the ankyrin form a flat surface to stack on the G-tetrad core. The protein employs a combination of hydrogen bonds and hydrophobic contacts to interact with the G4, and electrostatic interaction is used to enhance the binding affinity. This binding mechanism provides valuable insights into understanding G4 recognition by proteins.

Supramolecular DNA-based catalysis in organic solvents

The distinct folding accompanied by its polymorphic character renders DNA G-quadruplexes promising biomolecular building blocks to construct novel DNA-based and supramolecular assemblies. However, the highly polar nature of DNA limits the use of G-quadruplexes to water as a solvent. In addition, the archetypical G-quadruplex fold needs to be stabilized by metal-cations, which is usually a potassium ion. Here, we show that a noncovalent PEGylation process enabled by electrostatic interactions allows the first metal-free G-quadruplexes in organic solvents. Strikingly, incorporation of an iron-containing porphyrin renders the self-assembled metal-free G-quadruplex catalytically active in organic solvents. Hence, these "supraG4zymes" enable DNA-based catalysis in organic media. The results will allow the broad utilization of DNA G-quadruplexes in nonaqueous environments.

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.

Why Are Left-Handed G-Quadruplexes Scarce?

G-quadruplexes (G4s) are nucleic acid structures crucial for the regulation of gene expression and genome maintenance. While they hold promise as nanodevice components, achieving desired G4 folds requires understanding the interplay between stability and structural properties, like helicity. Although right-handed G4 structures dominate the experimental data, the molecular basis for this preference over left-handed helicity is unclear. To address this, we employ all-atom molecular dynamics simulations and quantum chemical methods. Our results reveal that right-handed G4s exhibit greater thermodynamic and kinetic stability as a result of favorable sugar-phosphate backbone conformations in guanine tracts. Moreover, while hydrogen-bonding patterns influence helicity-specific G4 loop conformations, they minimally affect stability differences. We also elucidate the strong correlation between helicity and the strand progression direction, essential for G4 structures. These findings deepen our understanding of G4s, providing molecular-level insights into their structural and energetic preferences, which could inform the design of novel nanodevices.

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

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

Showcasing Different G-Quadruplex Folds of a G-Rich Sequence: Between Rule-Based Prediction and Butterfly Effect

In better understanding the interactions of G-quadruplexes in a cellular or noncellular environment, a reliable sequence-based prediction of their three-dimensional fold would be extremely useful, yet is often limited by their remarkable structural diversity. A G-rich sequence related to a promoter sequence of the PDGFR-β nuclease hypersensitivity element (NHE) comprises a G3-G3-G2-G4-G3 pattern of five G-runs with two to four G residues. Although the predominant formation of three-layered canonical G-quadruplexes with uninterrupted G-columns can be expected, minimal base substitutions in a non-G-tract domain were shown to guide folding into either a basket-type antiparallel quadruplex, a parallel-stranded quadruplex with an interrupted G-column, a quadruplex with a V-shaped loop, or a (3+1) hybrid quadruplex. A 3D NMR structure for each of the different folds was determined. Supported by thermodynamic profiling on additional sequence variants, formed topologies were rationalized by the identification and assessment of specific critical interactions of loop and overhang residues, giving valuable insights into their contribution to favor a particular conformer. The variability of such tertiary interactions, together with only small differences in quadruplex free energies, emphasizes current limits for a reliable sequence-dependent prediction of favored topologies from sequences with multiple irregularly positioned G-tracts.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations

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

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.

ASC-G4, an algorithm to calculate advanced structural characteristics of G-quadruplexes

ASC-G4 is an algorithm for the calculation of the advanced structural characteristics of G-quadruplexes (G4). It allows the unambiguous determination of the intramolecular G4 topology, based on the oriented strand numbering. It also resolves the ambiguity in the determination of the guanine glycosidic configuration. With this algorithm, we showed that the use of the C3' or C5' atoms to calculate the groove width in G4 is more appropriate than the P atoms and that the groove width does not always reflect the space available within the groove. For the latter, the minimum groove width is more appropriate. The application of ASC-G4 to 207 G4 structures guided the choices made for the calculations. A website based on ASC-G4 (http://tiny.cc/ASC-G4) was created, where the user uploads his G4 structure and gets its topology, the types of its loops and their lengths, the presence of snapbacks and bulges, the distribution of guanines in the tetrads and strands, the glycosidic configuration of these guanines, their rise, the groove widths, the minimum groove widths, the tilt and twist angles, the backbone dihedral angles, etc. It also provides a large number of atom-atom and atom-plane distances that are relevant to evaluating the quality of the structure.

Polymorphism and Ligand Binding Modulate Fast Dynamics of Human Telomeric G-Quadruplexes

Telomeric G-quadruplexes (G4s) are promising targets in the design and development of anticancer drugs. Their actual topology depends on several factors, resulting in structural polymorphism. In this study, we investigate how the fast dynamics of the telomeric sequence AG3(TTAG3)3 (Tel22) depends on the conformation. By using Fourier transform Infrared spectroscopy, we show that, in the hydrated powder state, Tel22 adopts parallel and mixed antiparallel/parallel topologies in the presence of K+ and Na+ ions, respectively. These conformational differences are reflected in the reduced mobility of Tel22 in Na+ environment in the sub-nanosecond timescale, as probed by elastic incoherent neutron scattering. These findings are consistent with the G4 antiparallel conformation being more stable than the parallel one, possibly due to the presence of ordered hydration water networks. In addition, we study the effect of Tel22 complexation with BRACO19 ligand. Despite the quite similar conformation in the complexed and uncomplexed state, the fast dynamics of Tel22-BRACO19 is enhanced compared to that of Tel22 alone, independently of the ions. We ascribe this effect to the preferential binding of water molecules to Tel22 against the ligand. The present results suggest that the effect of polymorphism and complexation on the G4 fast dynamics is mediated by hydration water.

Unusual Paradigm for DNA-DNA Recognition and Binding: "Socket-Plug" Complementarity

DNA is the key informational polymer in biology by virtue of its precisely defined self-assembling properties. Watson-Crick complementarity, which underlies DNA's self-assembly, is required not only in biology but has also proved powerful in the field of nanoscience, where it has been utilized to assemble complex 2D and 3D architectures and nanodevices built from the DNA double-helix. Aside from Watson-Crick base-pairing, however, DNA also participates in alternative base pairing schemes, giving rise to DNA triplexes and G-quadruplexes. Herein, we describe "sticky-ended" DNA triplex-quadruplex composites that specifically recognize and bind to each other using a wholly different logic, "socket-plug" complementarity, a shape-sensing fitting of guanine "prongs" into guanine-lacking "cavities." A remarkable property of this kind of complementarity is the key role played in it by specific counter-cations: thus, exclusive "self" socket-plug recognition occurs over "other" in sodium salt solutions while precisely the reverse occurs in potassium salt solutions. We have used gel electrophoresis, Förster resonance energy transfer, alkylation protection, and structural modeling to study this remarkable fundamental property of DNA, that we anticipate will find wide practical application.

Role of fast dynamics in the complexation of G-quadruplexes with small molecules

G-quadruplexes (G4s) formed by the human telomeric sequence AG₃ (TTAG₃)₃ (Tel22) play a key role in cancer and aging. We combined elastic incoherent neutron scattering (EINS) and quasielastic incoherent neutron scattering (QENS) to characterize the internal dynamics of Tel22 G4s and to assess how it is affected by complexation with two standard ligands, Berberine and BRACO19. We show that the interaction with the two ligands induces an increase of the overall mobility of Tel22 as quantified by the mean squared displacements (MSD) of hydrogen atoms. At the same time, the complexes display a lower stiffness than G4 alone. Two different types of motion characterize the G4 nanosecond timescale dynamics. Upon complexation, an increasing fraction of G4 atomic groups participate in this fast dynamics, along with an increase in the relevant characteristic length scales. We suggest that the entropic contribution to the conformational free energy of these motions might be crucial for the complexation mechanisms.

Nanocellulose Composites as Smart Devices With Chassis, Light-Directed DNA Storage, Engineered Electronic Properties, and Chip Integration

The rapid development of green and sustainable materials opens up new possibilities in the field of applied research. Such materials include nanocellulose composites that can integrate many components into composites and provide a good chassis for smart devices. In our study, we evaluate four approaches for turning a nanocellulose composite into an information storage or processing device: 1) nanocellulose can be a suitable carrier material and protect information stored in DNA. 2) Nucleotide-processing enzymes (polymerase and exonuclease) can be controlled by light after fusing them with light-gating domains; nucleotide substrate specificity can be changed by mutation or pH change (read-in and read-out of the information). 3) Semiconductors and electronic capabilities can be achieved: we show that nanocellulose is rendered electronic by iodine treatment replacing silicon including microstructures. Nanocellulose semiconductor properties are measured, and the resulting potential including single-electron transistors (SET) and their properties are modeled. Electric current can also be transported by DNA through G-quadruplex DNA molecules; these as well as classical silicon semiconductors can easily be integrated into the nanocellulose composite. 4) To elaborate upon miniaturization and integration for a smart nanocellulose chip device, we demonstrate pH-sensitive dyes in nanocellulose, nanopore creation, and kinase micropatterning on bacterial membranes as well as digital PCR micro-wells. Future application potential includes nano-3D printing and fast molecular processors (e.g., SETs) integrated with DNA storage and conventional electronics. This would also lead to environment-friendly nanocellulose chips for information processing as well as smart nanocellulose composites for biomedical applications and nano-factories.

G-quadruplex structure of the C. elegans telomeric repeat: a two tetrads basket type conformation stabilized by a non-canonical C-T base-pair

The Caenorhabditis elegans model has greatly contributed to the understanding of the role of G-quadruplexes in genomic instability. The GGCTTA repeats of the C. elegans telomeres resemble the GGGTTA repeats of the human telomeres. However, the comparison of telomeric sequences (Homo sapiens, Tetrahymena, Oxytricha, Bombyx mori and Giardia) revealed that small changes in these repeats can drastically change the topology of the folded G-quadruplex. In the present work we determined the structure adopted by the C. elegans telomeric sequence d[GG(CTTAGG)3]. The investigated C. elegans telomeric sequence is shown to fold into an intramolecular two G-tetrads basket type G-quadruplex structure that includes a C-T base pair in the diagonal loop. This work sheds light on the telomeric structure of the widely used C. elegans animal model.

Guiding the folding of G-quadruplexes through loop residue interactions

A G-rich sequence was designed to allow folding into either a stable parallel or hybrid-type topology. With the parent sequence featuring coexisting species, various related sequences with single and double mutations and with a shortened central propeller loop affected the topological equilibrium. Two simple modifications, likewise introduced separately to all sequences, were employed to lock folds into one of the topologies without noticeable structural alterations. The unique combination of sequence mutations, high-resolution NMR structural information, and the thermodynamic stability for both topological competitors identified critical loop residue interactions. In contrast to first loop residues, which are mostly disordered and exposed to solvent in both propeller and lateral loops bridging a narrow groove, the last loop residue in a lateral three-nucleotide loop is engaged in stabilizing stacking interactions. The propensity of single-nucleotide loops to favor all-parallel topologies by enforcing a propeller-like conformation of an additional longer loop is shown to result from their preference in linking two outer tetrads of the same tetrad polarity. Taken together, the present studies contribute to a better structural and thermodynamic understanding of delicate loop interactions in genomic and artificially designed quadruplexes, e.g. when employed as therapeutics or in other biotechnological applications.

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

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

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.

Detection of SARS-CoV-2 RNA Using a DNA Aptamer Mimic of Green Fluorescent Protein

RNA detection is important in diverse diagnostic and analytical applications. RNAs can be rapidly detected using molecular beacons, which fluoresce upon hybridizing to a target RNA but require oligonucleotides with complex fluorescent dye and quencher conjugations. Here, we describe a simplified method for rapid fluorescence detection of a target RNA using simple unmodified DNA oligonucleotides. To detect RNA, we developed Lettuce, a fluorogenic DNA aptamer that binds and activates the fluorescence of DFHBI-1T, an otherwise nonfluorescent molecule that resembles the chromophore found in green fluorescent protein. Lettuce was selected from a randomized DNA library based on binding to DFHBI-agarose. We further show that Lettuce can be split into two separate oligonucleotide components, which are nonfluorescent on their own but become fluorescent when their proximity is induced by a target RNA. We designed several pairs of split Lettuce fragments that contain an additional 15-20 nucleotides that are complementary to adjacent regions of the SARS-CoV-2 RNA, resulting in Lettuce fluorescence only in the presence of the viral RNA. Overall, these studies describe a simplified RNA detection approach using fully unmodified DNA oligonucleotides that reconstitute the Lettuce aptamer templated by RNA.

Ruthenium Polypyridyl Complex Bound to a Unimolecular Chair-Form G-Quadruplex

The DNA G-quadruplex is known for forming a range of topologies and for the observed lability of the assembly, consistent with its transient formation in live cells. The stabilization of a particular topology by a small molecule is of great importance for therapeutic applications. Here, we show that the ruthenium complex Λ-[Ru(phen)₂(qdppz)]²⁺ displays enantiospecific G-quadruplex binding. It crystallized in 1:1 stoichiometry with a modified human telomeric G-quadruplex sequence, GGGTTAGGGTTAGGGTTTGGG (htel21T₁₈), in an antiparallel chair topology, the first structurally characterized example of ligand binding to this topology. The lambda complex is bound in an intercalation cavity created by a terminal G-quartet and the central narrow lateral loop formed by T₁₀–T₁₁–A₁₂. The two remaining wide lateral loops are linked through a third K⁺ ion at the other end of the G-quartet stack, which also coordinates three thymine residues. In a comparative ligand-binding study, we showed, using a Klenow fragment assay, that this complex is the strongest observed inhibitor of replication, both using the native human telomeric sequence and the modified sequence used in this work.

ONQUADRO: a database of experimentally determined quadruplex structures

ONQUADRO is an advanced database system that supports the study of the structures of canonical and non-canonical quadruplexes. It combines a relational database that collects comprehensive information on tetrads, quadruplexes, and G4-helices; programs to compute structure parameters and visualise the data; scripts for statistical analysis; automatic updates and newsletter modules; and a web application that provides a user interface. The database is a self-updating resource, with new information arriving once a week. The preliminary data are downloaded from the Protein Data Bank, processed, annotated, and completed. As of August 2021, ONQUADRO contains 1,661 tetrads, 518 quadruplexes, and 30 G4-helices found in 467 experimentally determined 3D structures of nucleic acids. Users can view and download their description: sequence, secondary structure (dot-bracket, classical diagram, arc diagram), tertiary structure (ball-and-stick, surface or vdw-ball model, layer diagram), planarity, twist, rise, chi angle (value and type), loop characteristics, strand directionality, metal ions, ONZ, and Webba da Silva classification (the latter by loop topology and tetrad combination), origin structure ID, assembly ID, experimental method, and molecule type. The database is freely available at https://onquadro.cs.put.poznan.pl/. It can be used on both desktop computers and mobile devices.

OUP accepted manuscript

G-quadruplexes (GQs) are non-canonical DNA structures composed of stacks of stabilized G-tetrads. GQs play an important role in a variety of biological processes and may form at telomeres and oncogene promoters among other genomic locations. Here, we investigate nine variants of telomeric DNA from Tetrahymena thermophila with the repeat (TTGGGG)_n. Biophysical data indicate that the sequences fold into stable four-tetrad GQs which adopt multiple conformations according to native PAGE. Excitingly, we solved the crystal structure of two variants, TET25 and TET26. The two variants differ by the presence of a 3′-T yet adopt different GQ conformations. TET25 forms a hybrid [3 + 1] GQ and exhibits a rare 5′-top snapback feature. Consequently, TET25 contains four loops: three lateral (TT, TT, and GTT) and one propeller (TT). TET26 folds into a parallel GQ with three TT propeller loops. To the best of our knowledge, TET25 and TET26 are the first reported hybrid and parallel four-tetrad unimolecular GQ structures. The results presented here expand the repertoire of available GQ structures and provide insight into the intricacy and plasticity of the 3D architecture adopted by telomeric repeats from T. thermophila and GQs in general.

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.

Chemical Modulation of DNA Replication along G-Quadruplex Based on Topology-Dependent Ligand Binding

Ligands that bind to and stabilize guanine-quadruplex (G4) structures to regulate DNA replication have therapeutic potential for cancer and neurodegenerative diseases. Because there are several G4 topologies, ligands that bind to their specific types may have the ability to preferentially regulate the replication of only certain genes. Here, we demonstrated that binding ligands stalled the replication of template DNA at G4, depending on different topologies. For example, naphthalene diimide derivatives bound to the G-quartet of G4 with an additional interaction between the ligand and the loop region of a hybrid G4 type from human telomeres, which efficiently repressed the replication of the G4. Thus, these inhibitory effects were not only stability-dependent but also topology-selective based on the manner in which G4 structures interacted with G4 ligands. Our original method, referred to as a quantitative study of topology-dependent replication (QSTR), was developed to evaluate correlations between replication rate and G4 stability. QSTR enabled the systematic categorization of ligands based on topology-dependent binding. It also demonstrated accuracy in determining quantitatively how G4 ligands control the intermediate state of replication and the kinetics of G4 unwinding. Hence, the QSTR index would facilitate the design of new drugs capable of controlling the topology-dependent regulation of gene expression.

Two coexisting pseudo-mirror heteromolecular telomeric G-quadruplexes in opposite loop progressions differentially recognized by a low equivalent of Thioflavin T

The final 3′-terminal residue of the telomeric DNA G-overhang is inherently less precise. Here, we describe how alteration of the last 3′-terminal base affects the mutual recognition between two different G-rich oligomers of human telomeric DNA in the formation of heteromolecular G-quadruplexes (hetero-GQs). Associations between three- and single-repeat fragments of human telomeric DNA, target d(GGGTTAGGGTTAGGG) and probe d(TAGGGT), in Na⁺ solution yield two coexisting forms of (3 + 1) hybrid hetero-GQs: the kinetically favourable LLP-form (left loop progression) and the thermodynamically controlled RLP-form (right loop progression). However, only the adoption of a single LLP-form has been previously reported between the same probe d(TAGGGT) and a target variant d(GGGTTAGGGTTAGGGT) having one extra 3′-end thymine. Moreover, the flanking base alterations of short G-rich probe variants also significantly affect the loop progressions of hetero-GQs. Although seemingly two pseudo-mirror counter partners, the RLP-form exhibits a preference over the LLP-form to be recognized by a low equivalent of fluorescence dye thioflavin T (ThT). To a greater extent, ThT preferentially binds to RLP hetero-GQ than with the corresponding telomeric DNA duplex context or several other representative unimolecular GQs.

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

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

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

The beginning and the end: flanking nucleotides induce a parallel G-quadruplex topology

Genomic sequences susceptible to form G-quadruplexes (G4s) are always flanked by other nucleotides, but G4 formation in vitro is generally studied with short synthetic DNA or RNA oligonucleotides, for which bases adjacent to the G4 core are often omitted. Herein, we systematically studied the effects of flanking nucleotides on structural polymorphism of 371 different oligodeoxynucleotides that adopt intramolecular G4 structures. We found out that the addition of nucleotides favors the formation of a parallel fold, defined as the 'flanking effect' in this work. This 'flanking effect' was more pronounced when nucleotides were added at the 5'-end, and depended on loop arrangement. NMR experiments and molecular dynamics simulations revealed that flanking sequences at the 5'-end abolish a strong syn-specific hydrogen bond commonly found in non-parallel conformations, thus favoring a parallel topology. These analyses pave a new way for more accurate prediction of DNA G4 folding in a physiological context.

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.

Expanding the Topological Landscape by a G-Column Flip of a Parallel G-Quadruplex

Canonical G‐quadruplexes can adopt a variety of different topologies depending on the arrangement of propeller, lateral, or diagonal loops connecting the four G‐columns. A novel intramolecular G‐quadruplex structure is derived through inversion of the last G‐tract of a three‐layered parallel fold, associated with the transition of a single propeller into a lateral loop. The resulting (3+1) hybrid fold features three syn⋅anti⋅anti⋅anti G‐tetrads with a 3’‐terminal all‐syn G‐column. Although the ability of forming a duplex stem‐loop between G‐tracts seems beneficial for a propeller‐to‐lateral loop rearrangement, unmodified G‐rich sequences resist folding into the new (3+1) topology. However, refolding can be driven by the incorporation of syn‐favoring guanosine analogues into positions of the fourth G‐stretch. The presented hybrid‐type G‐quadruplex structure as determined by NMR spectroscopy may provide for an additional scaffold in quadruplex‐based technologies.

Structural motifs and intramolecular interactions in non-canonical G-quadruplexes

Guanine(G)-rich DNA or RNA sequences can assemble or intramolecularly fold into G-quadruplexes formed through the stacking of planar G·G·G·G tetrads in the presence of monovalent cations. These secondary nucleic acid structures have convincingly been shown to also exist within a cellular environment exerting important regulatory functions in physiological processes. For identifying nucleic acid segments prone to quadruplex formation, a putative quadruplex sequence motif encompassing closely spaced tracts of three or more guanosines is frequently employed for bioinformatic search algorithms. Depending on the number and type of intervening residues as well as on solution conditions, such sequences may fold into various canonical G4 topologies with continuous G-columns. On the other hand, a growing number of sequences capable of quadruplex formation feature G-deficient guanine tracts, escaping the conservative consensus motif. By folding into non-canonical quadruplex structures, they adopt unique topologies depending on their specific sequence context. These include G-columns with only two guanines, bulges, snapback loops, D- and V-shaped loops as well as interlocked structures. This review focuses on G-quadruplex species carrying such distinct structural motifs. It evaluates characteristic features of their non-conventional scaffold and highlights principles of stabilizing interactions that also allow for their folding into stable G-quadruplex structures.

A novel minimal motif for left-handed G-quadruplex formation

A recent study on the left-handed G-quadruplex (LHG4) DNA revealed a 12-nt minimal motif GTGGTGGTGGTG with the ability to independently form an LHG4 and to drive an adjacent sequence to LHG4 formation. Here we have identified a second LHG4-forming motif, GGTGGTGGTGTG, and determined the X-ray crystal structure of an LHG4 involving this motif. Our structural analysis indicated the role of split guanines and single thymine loops in promoting LHG4 formation.

Concurrent formation of H- and J-aggregates of dyes with chiralities individually determined by G-quadruplex handedness

DNA templated dye assemblies pave an easy way to regulate the optical properties of molecular aggregates. G-quadruplexes (G4s) provide versatile DNA platforms for the dye assemblies since their foldings can be easily tuned by cation ions and sequences. In this work, we found that the G4 handedness can be used to control the aggregate chirality of a dye of 3,3'-diethylthiacarbocyanine (DiSC₂(3)). The left-handed and right-handed G4s can template the concurrent formation of the J- and H-aggregates of DiSC₂(3) with emergence of the featured absorption spectra. However, the chiral J-aggregate of DiSC₂(3) can be formed only on the left-handed G4s, while the chiral H-aggregate is otherwise grown only on the right-handed G4s, as confirmed by the induced circular dichroism (ICD) spectra with the characteristic splitting bands. Additionally, these G4s even at tens of nM level are efficient to produce these chiral aggregates, demonstrating the high sensitivity of G4s in creating these optically active dye assemblies. The possible growth sites of the aggregates are proposed by the sequence length-dependent assemblies. Our work will provide a new way to control the chiral assemblies of dye aggregates via the G4 handedness.

DNA G-quadruplexes for native mass spectrometry in potassium: a database of validated structures in electrospray-compatible conditions

G-quadruplex DNA structures have become attractive drug targets, and native mass spectrometry can provide detailed characterization of drug binding stoichiometry and affinity, potentially at high throughput. However, the G-quadruplex DNA polymorphism poses problems for interpreting ligand screening assays. In order to establish standardized MS-based screening assays, we studied 28 sequences with documented NMR structures in (usually ∼100 mM) potassium, and report here their circular dichroism (CD), melting temperature (Tm), NMR spectra and electrospray mass spectra in 1 mM KCl/100 mM trimethylammonium acetate. Based on these results, we make a short-list of sequences that adopt the same structure in the MS assay as reported by NMR, and provide recommendations on using them for MS-based assays. We also built an R-based open-source application to build and consult a database, wherein further sequences can be incorporated in the future. The application handles automatically most of the data processing, and allows generating custom figures and reports. The database is included in the g4dbr package (https://github.com/EricLarG4/g4dbr) and can be explored online (https://ericlarg4.github.io/G4_database.html).

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

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

Reversible photo-regulation on the folding/unfolding of telomere G-quadruplexes with solid-state nanopores

The formation of G-quadruplexes (G4) in human telomere and other important biological regions inhibits the replication and transcription of DNA, thereby influencing further cell proliferation. The investigation of G4 formation and unfolding is vital for understanding their modulation in biological processes and life science. Photo regulation is a facile and sensitive approach for monitoring the structures of biomacromolecules and material surface properties. The nanopore-based technique is also prevalent for label-free single-molecule characterization with high accuracy. This study provides a combination of solid-state nanopore technology with light-switch as a platform for the modulation of human telomere G4 formation and splitting under switchable light exposure. The introduction of molecular switch, namely azobenzene moiety at different positions of the DNA sequence influences the formation and stability of G4. Three azobenzenes immobilized on each of the G-quartet plane (hTelo-3azo-p) or four azobenzenes on the same plane (hTelo-4azo-4p) of the human telomere G4 sequence realized the reversible control of G4 folding/unfolding at the temporal scale upon photo regulation, and the formation and splitting of G4 with hTelo-4azo-4p is slower and not thorough compared to that with hTelo-3azo-p due to the coplanar steric hindrance. Moreover, the G4 formation recorded with the combined nanopore and photo-responsive approach was also characterized with fluorescence, and the variation in the fluorescence intensity of the NMM and G4 complex exhibited a different tendency under reverse light irradiation due to the distinct interactions of NMM with the azobenzene-modified G4. Our study demonstrated a controllable and sensitive way for the manipulation of G4 structures, which will be inspiring for the intervention of G4-related cell senescence, cancer diagnosis and drug exploration.

The origin of the high stability of 3'-terminal uridine tetrads: contributions of hydrogen bonding, stacking interactions, and steric factors evaluated using modified oligonucleotide analogs

RNA G-quadruplexes fold almost exclusively into parallel-stranded structures and thus display much less structural diversity than their DNA counterparts. However, also among RNA G-quadruplexes peculiar structural elements can be found which are capable of reshaping the physico-chemical properties of the folded structure. A striking example is provided by a uridine tetrad (U-tetrad) placed on the 3'-terminus of the tetramolecular G-quadruplex. In this context, the U-tetrad adopts a unique conformation involving chain reversal and is responsible for a tremendous stabilization of the G-quadruplex (ΔT_m up to 30°C). In this report, we attempt to rationalize the origin of this stabilizing effect by concurrent structural, thermal stability, and molecular dynamics studies of a series of G-quadruplexes with subtle chemical modifications at their 3'-termini. Our results provide detailed insights into the energetics of the "reversed" U-tetrad motif and the requirements for its formation. They point to the importance of the 2'OH to phosphate hydrogen bond and preferential stacking interactions for the formation propensity and stability of the motif.

Native de novo structural determinations of non-canonical nucleic acid motifs by X-ray crystallography at long wavelengths

Obtaining phase information remains a formidable challenge for nucleic acid structure determination. The introduction of an X-ray synchrotron beamline designed to be tunable to long wavelengths at Diamond Light Source has opened the possibility to native de novo structure determinations by the use of intrinsic scattering elements. This provides opportunities to overcome the limitations of introducing modifying nucleotides, often required to derive phasing information. In this paper, we build on established methods to generate new tools for nucleic acid structure determinations. We report on the use of (i) native intrinsic potassium single-wavelength anomalous dispersion methods (K-SAD), (ii) use of anomalous scattering elements integral to the crystallization buffer (extrinsic cobalt and intrinsic potassium ions), (iii) extrinsic bromine and intrinsic phosphorus SAD to solve complex nucleic acid structures. Using the reported methods we solved the structures of (i) Pseudorabies virus (PRV) RNA G-quadruplex and ligand complex, (ii) PRV DNA G-quadruplex, and (iii) an i-motif of human telomeric sequence. Our results highlight the utility of using intrinsic scattering as a pathway to solve and determine non-canonical nucleic acid motifs and reveal the variability of topology, influence of ligand binding, and glycosidic angle rearrangements seen between RNA and DNA G-quadruplexes of the same sequence.

Polymorphism of human telomeric quadruplexes with drugs: a multi-technique biophysical study

The human telomeric G-quadruplex structural motif of DNA has come to be known as a new and stimulating target for anticancer drug discovery. Small molecules that interact with G-quadruplex structures in a selective way have gained impressive interest in recent years as they may serve as potential therapeutic agents. Here, we show how circular dichroism, UV resonance Raman and small angle X-ray scattering spectroscopies can be effectively combined to provide insights into structural and molecular aspects of the interaction between human telomeric quadruplexes and ligands. This study focuses on the ability of berberine and palmatine to bind with human telomeric quadruplexes and provides analysis of the conformational landscape visited by the relevant complexes upon thermal unfolding. With increasing temperature, both free and bound G-quadruplexes undergo melting through a multi-state process, populating different intermediate states. Despite the structural similarity of the two ligands, valuable distinctive features characterising their interaction with the G-quadruplex emerged from our multi-technique approach.

Topology-based classification of tetrads and quadruplex structures

Motivation

Quadruplexes attract the attention of researchers from many fields of bio-science. Due to a specific structure, these tertiary motifs are involved in various biological processes. They are also promising therapeutic targets in many strategies of drug development, including anticancer and neurological disease treatment. The uniqueness and diversity of their forms cause that quadruplexes show great potential in novel biological applications. The existing approaches for quadruplex analysis are based on sequence or 3D structure features and address canonical motifs only.

Results

In our study, we analyzed tetrads and quadruplexes contained in nucleic acid molecules deposited in Protein Data Bank. Focusing on their secondary structure topology, we adjusted its graphical diagram and proposed new dot-bracket and arc representations. We defined the novel classification of these motifs. It can handle both canonical and non-canonical cases. Based on this new taxonomy, we implemented a method that automatically recognizes the types of tetrads and quadruplexes occurring as unimolecular structures. Finally, we conducted a statistical analysis of these motifs found in experimentally determined nucleic acid structures in relation to the new classification.

Availability and implementation

https://github.com/tzok/eltetrado/

Supplementary information

Supplementary data are available at Bioinformatics online.

How bioinformatics resources work with G4 RNAs

Quadruplexes (G4s) are of interest, which increases with the number of identified G4 structures and knowledge about their biomedical potential. These unique motifs form in many organisms, including humans, where their appearance correlates with various diseases. Scientists store and analyze quadruplexes using recently developed bioinformatic tools—many of them focused on DNA structures. With an expanding collection of G4 RNAs, we check how existing tools deal with them. We review all available bioinformatics resources dedicated to quadruplexes and examine their usefulness in G4 RNA analysis. We distinguish the following subsets of resources: databases, tools to predict putative quadruplex sequences, tools to predict secondary structure with quadruplexes and tools to analyze and visualize quadruplex structures. We share the results obtained from processing specially created RNA datasets with these tools.

Contact: mszachniuk@cs.put.poznan.pl

Supplementary information: Supplementary data are available at Briefings in Bioinformatics online.

Infrared Spectroscopy Reveals the Preferred Motif Size and Local Disorder in Parallel Stranded DNA G-Quadruplexes

Infrared spectroscopy detects the formation of G-quadruplexes in guanine-rich nucleic acid sequences through shifts in the guanine C=O stretch mode. Here, we use ultrafast 2D infrared (IR) spectroscopy and isotope substitution to show that these shifts arise from vibrational delocalization among stacked G-quartets. This provides a direct measure of the sizes of locally ordered motifs in heterogeneous samples with substantial disordered regions. We find that parallel-stranded, potassium-bound DNA G-quadruplexes are limited to five consecutive G-quartets and 3-4 consecutive layers are preferred for longer polyguanine tracts. The resulting potassium-dependent G-quadruplex assembly landscape reflects the polyguanine tract lengths found in genomes, the ionic conditions prevalent in healthy mammalian cells, and the onset of structural disorder in disease states. Our study describes spectral markers that can be used to probe other G-quadruplex structures and provides insight into the fundamental limits of their formation in biological and artificial systems.

Intra-locked G-quadruplex structures formed by irregular DNA G-rich motifs

G-rich DNA sequences with tracts of three or more continuous guanines (G≥3) are known to have high propensity to adopt stable G-quadruplex (G4) structures. Bioinformatic analyses suggest high prevalence of G-rich sequences with short G-tracts (G≤2) in the human genome. However, due to limited structural studies, the folding principles of such sequences remain largely unexplored and hence poorly understood. Here, we present the solution NMR structure of a sequence named AT26 consisting of irregularly spaced G2 tracts and two isolated single guanines. The structure is a four-layered G4 featuring two bi-layered blocks, locked between themselves in an unprecedented fashion making it a stable scaffold. In addition to edgewise and propeller-type loops, AT26 also harbors two V-shaped loops: a 2-nt V-shaped loop spanning two G-tetrad layers and a 0-nt V-shaped loop spanning three G-tetrad layers, which are named as VS- and VR-loop respectively, based on their distinct structural features. The intra-lock motif can be a basis for extending the G-tetrad core and a very stable intra-locked G4 can be formed by a sequence with G-tracts of various lengths including several G2 tracts. Findings from this study will aid in understanding the folding of G4 topologies from sequences containing irregularly spaced multiple short G-tracts.

Stability of Two-Quartet G-Quadruplexes and Their Dimers in Atomistic Simulations

G-quadruplexes (GQs) are four-stranded noncanonical DNA and RNA architectures that can be formed by guanine-rich sequences. The stability of GQs increases with the number of G-quartets, and three G-quartets generally form stable GQs. However, the stability of two-quartet GQs is an open issue. To understand the intrinsic stability of two-quartet GQ stems, we have carried out a series of unbiased molecular dynamics (MD) simulations (505 μs in total) of two- and four-quartet DNA and RNA GQs, with attention paid mainly to parallel-stranded arrangements. We used AMBER DNA parmOL15 and RNA parmOL3 force fields and tested different ion and water models. Two-quartet parallel-stranded DNA GQs unfolded in all the simulations, while the equivalent RNA GQ was stable in most of the simulations. GQs composed of two stacked units of two-quartet GQs were stable for both DNA and RNA. The simulations suggest that a minimum of three quartets are needed to form an intrinsically stable all-anti parallel-stranded DNA GQ. Parallel two-quartet DNA GQ may exist if substantially stabilized by another molecule or structural element, including multimerization. On the other hand, we predict that isolated RNA two-quartet parallel GQs may form, albeit being weakly stable. We also show that ionic parameters and water models should be chosen with caution because some parameter combinations can cause spurious instability of GQ stems. Some in-so-far unnoticed limitations of force-field description of multiple ions inside the GQs are discussed, which compromise the capability of simulations to fully capture the effect of increase in the number of quartets on the GQ stability.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.