Strand directionality affects cation binding and movement within tetramolecular G-quadruplexes

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

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

Insight into Tetramolecular DNA G-Quadruplexes Associated with ALS and FTLD: Cation Interactions and Formation of Higher-Ordered Structure

The G4C2 hexanucleotide repeat expansion in the c9orf72 gene is a major genetic cause of familial amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD), with the formation of G-quadruplexes directly linked to the development of these diseases. Cations play a crucial role in the formation and structure of G-quadruplexes. In this study, we investigated the impact of biologically relevant potassium ions on G-quadruplex structures and utilized 15N-labeled ammonium cations as a substitute for K+ ions to gain further insights into cation binding and exchange dynamics. Through nuclear magnetic resonance spectroscopy and molecular dynamics simulations, we demonstrate that the single d(G4C2) repeat, in the presence of 15NH4+ ions, adopts a tetramolecular G-quadruplex with an all-syn quartet at the 5'-end. The movement of 15NH4+ ions through the central channel of the G-quadruplex, as well as to the bulk solution, is governed by the vacant cation binding site, in addition to the all-syn quartet at the 5'-end. Furthermore, the addition of K+ ions to G-quadruplexes folded in the presence of 15NH4+ ions induces stacking of G-quadruplexes via their 5'-end G-quartets, leading to the formation of stable higher-ordered species.

8-Oxoguanine Forms Quartets with a Large Central Cavity

Oxidation of a guanine nucleotide in DNA yields an 8-oxoguanine nucleotide (^oxoG) and is a mutagenic event in the genome. Due to different arrangements of hydrogen-bond donors and acceptors, ^oxoG can affect the secondary structure of nucleic acids. We have investigated base pairing preferences of ^oxoG in the core of a tetrahelical G-quadruplex structure, adopted by analogues of d(TG₄T). Using spectroscopic methods, we have shown that G-quartets can be fully substituted with ^oxoG nucleobases to form an ^oxoG-quartet with a revamped hydrogen-bonding scheme. While an ^oxoG-quartet can be incorporated into the G-quadruplex core without distorting the phosphodiester backbone, larger dimensions of the central cavity change the cation localization and exchange properties.

Unprecedented hour-long residence time of a cation in a left-handed G-quadruplex

Cations are critical for the folding and assembly of nucleic acids. In G-quadruplex structures, cations can bind between stacked G-tetrads and coordinate with negatively charged guanine carbonyl oxygens. They usually exchange between binding sites and with the bulk in solution with time constants ranging from sub-millisecond to seconds. Here we report the first observation of extremely long-lived K+ and NH4 + ions, with an exchange time constant on the order of an hour, when coordinated at the center of a left-handed G-quadruplex DNA. A single-base mutation, that switched one half of the structure from left- to right-handed conformation resulting in a right-left hybrid G-quadruplex, was shown to remove this long-lived behaviour of the central cation.

Overlapping but distinct: a new model for G-quadruplex biochemical specificity

G-quadruplexes are noncanonical nucleic acid structures formed by stacked guanine tetrads. They are capable of a range of functions and thought to play widespread biological roles. This diversity raises an important question: what determines the biochemical specificity of G-quadruplex structures? The answer is particularly important from the perspective of biological regulation because genomes can contain hundreds of thousands of G-quadruplexes with a range of functions. Here we analyze the specificity of each sequence in a 496-member library of variants of a reference G-quadruplex with respect to five functions. Our analysis shows that the sequence requirements of G-quadruplexes with these functions are different from one another, with some mutations altering biochemical specificity by orders of magnitude. Mutations in tetrads have larger effects than mutations in loops, and changes in specificity are correlated with changes in multimeric state. To complement our biochemical data we determined the solution structure of a monomeric G-quadruplex from the library. The stacked and accessible tetrads rationalize why monomers tend to promote a model peroxidase reaction and generate fluorescence. Our experiments support a model in which the sequence requirements of G-quadruplexes with different functions are overlapping but distinct. This has implications for biological regulation, bioinformatics, and drug design.

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

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

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

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

DNA and RNA telomeric G-quadruplexes: what topology features can be inferred from ion mobility mass spectrometry?

Maintenance of the telomeres is key to chromosome integrity and cell proliferation. The G-quadruplex structures formed by telomeric DNA and RNA (TTAGGG and UUAGGG repeats, respectively) are key to this process. However, because these sequences are particularly polymorphic, solving high-resolution structures is not always possible, and there is a need for new methodologies to characterize the multiple structures coexisting in solution. In this context, we evaluated whether ion mobility spectrometry coupled to native mass spectrometry could help separate and assign the G-quadruplex topologies. We explored the circular dichroism spectra, multimer formation, cation binding, and ion mobility spectra of several 4-repeat and 8-repeat telomeric DNA and RNA sequences, both in NH₄⁺ and in K⁺. In 1 mM K⁺ and 100 mM trimethylammonium acetate, all RNAs fold intramolecularly (no multimer). In 8-repeat sequences, the subunits are not independent: in DNA the first subunit disfavors the folding of the second one, whereas in RNA the two subunits fold cooperatively via cation-mediated stacking. Ion mobility spectrometry shows that gas-phase structures keep a memory of - but are not identical to - the solution ones. At the native charge states, the loops can rearrange in a variety of ways (unless they are constrained by pre-formed hydrogen bonds), thereby wrapping the core and masking the strand arrangements. Our study highlights that, to progress towards structural assignment from IM-MS experiments, deeper understanding of the solution-to-gas-phase rearrangement mechanisms is warranted.

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

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

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.

Neutral and Negatively Charged Phosphate Modifications Altering Thermal Stability, Kinetics of Formation and Monovalent Ion Dependence of DNA G-Quadruplexes

The effect of phosphate group modifications on formation and properties of G-quadruplexes (G4s) has not been investigated in detail. Here, we evaluated the structural, thermodynamic and kinetic properties of the parallel G-quadruplexes formed by oligodeoxynucleotides d(G₄ T), d(TG₄ T) and d(TG₅ T), in which all phosphates were replaced with N-methanesulfonyl (mesyl) phosphoramidate or phosphoryl guanidine groups resulting in either negatively charged or neutral DNA sequences, respectively. We established that all modified sequences were able to form G-quadruplexes of parallel topology; however, the presence of modifications led to a decrease in thermal stability relative to unmodified G4s. In contrast to negatively charged G4s, assembly of neutral G4 DNA species was faster in the presence of sodium ions than potassium ions, and was independent of the salt concentration used. Formation of mixed G4s composed of both native and neutral G-rich strands has been detected using native gel electrophoresis, size-exclusion chromatography and ESI-MS. In summary, our results indicate that the phosphate modifications studied are compatible with G-quadruplex formation, which could be used for the design of biologically active compounds.

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.

Loop Length Affects Syn-Anti Conformational Rearrangements in Parallel G-Quadruplexes

A G-quadruplex forming sequence from the MYC promoter region was modified with syn-favoring 8-bromo-2'-deoxyguanosine residues. Depending on the number and position of modifications in the intramolecular parallel G-quadruplex, substitutions with the bromoguanosine analogue at the 5'-tetrad induce conformational rearrangements with concerted all-anti to all-syn transitions for all residues of the modified G-quartet. No unfavorable steric interactions of the C8-substituents in the medium grooves are apparent in the high-resolution structure as determined for a tetrasubstituted MYC quadruplex that exclusively forms the all-syn isomer. In contrast, considerable steric clashes with 5'-phosphate oxygen atoms for those analogues that follow a less flexible 1-nucleotide loop in the native all-anti conformation seem to constitute the major driving force for the tetrad inversion and allow for the rational design of appropriately substituted sequences. Correlations found between the population of species subjected to a tetrad flip and melting temperatures indicate that more effective conformational transitions are compromised by lower thermal stabilities of the modified parallel quadruplexes.

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.

Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes

G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]₄ GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]₄ simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na⁺ ions move inside the GQs in a concerted manner, while larger relocations of the K⁺ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K⁺ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.

G-quadruplex dynamics

G-quadruplexes (GQs) are four-stranded nucleic acid secondary structures formed by guanosine (G)-rich DNA and RNA sequences. It is becoming increasingly clear that cellular processes including gene expression and mRNA translation are regulated by GQs. GQ structures have been extensively characterized, however little attention to date has been paid to their conformational dynamics, despite the fact that many biological GQ sequences populate multiple structures of similar free energies, leading to an ensemble of exchanging conformations. The impact of these dynamics on biological function is currently not well understood. Recently, structural dynamics have been demonstrated to entropically stabilize GQ ensembles, potentially modulating gene expression. Transient, low-populated states in GQ ensembles may additionally regulate nucleic acid interactions and function. This review will underscore the interplay of GQ dynamics and biological function, focusing on several dynamic processes for biological GQs and the characterization of GQ dynamics by nuclear magnetic resonance (NMR) spectroscopy in conjunction with other biophysical techniques. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Unexpected Position-Dependent Effects of Ribose G-Quartets in G-Quadruplexes

To understand the role of ribose G-quartets and how they affect the properties of G-quadruplex structures, we studied three systems in which one, two, three, or four deoxyribose G-quartets were substituted with ribose G-quartets. These systems were a parallel DNA intramolecular G-quadruplex, d(TTGGGTGGGTTGGGTGGGTT), and two tetramolecular G-quadruplexes, d(TGGGT) and d(TGGGGT). Thermal denaturation experiments revealed that ribose G-quartets have position-dependent and cumulative effects on G-quadruplex stability. An unexpected destabilization was observed when rG quartets were presented at the 5'-end of the G stack. This observation challenges the general belief that RNA residues stabilize G-quadruplexes. Furthermore, in contrast to past proposals, hydration is not the main factor determining the stability of our RNA/DNA chimeric G-quadruplexes. Interestingly, the presence of rG residues in a central G-quartet facilitated the formation of additional tetramolecular G-quadruplex topologies showing positive circular dichroism signals at 295 nm. 2D NMR analysis of the tetramolecular TGgGGT (lowercase letter indicates ribose) indicates that Gs in the 5'-most G-quartet adopt the syn conformation. These analyses highlight several new aspects of the role of ribose G-quartets on G-quadruplex structure and stability, and demonstrate that the positions of ribose residues are critical for tuning G-quadruplex properties.

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

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

Observation of a Dynamic G-Tetrad Flip in Intramolecular G-Quadruplexes

A MYC sequence forming an intramolecular G-quadruplex with a parallel topology was modified by the incorporation of 8-bromoguanosine (^BrG) analogues in one of its outer G-tetrads. The propensity of the ^BrG analogues to adopt a syn glycosidic torsion angle results in an exceptional monomolecular quadruplex conformation featuring a complete flip of one tetrad while keeping a parallel orientation of all G-tracts as shown by circular dichroism and nuclear magnetic resonance spectroscopic studies. When substituting three of the four G-tetrad residues with ^BrG analogues, two coexisting quadruplex conformational isomers with an all-syn and all-anti outer G-quartet are approximately equally populated in solution. A dynamic interconversion of the two quadruplexes with an exchange rate (k_ex) of 0.2 s^-1 is demonstrated through the observation of exchange crosspeaks in rotating frame Overhauser effect spectroscopy and nuclear Overhauser effect spectroscopy experiments at 50 °C. The kinetic properties suggest disruption of the corresponding outer G-tetrad but not of the whole quadruplex core during the tetrad flip. Conformational syn-anti isomers with homopolar and heteropolar stacking interactions are nearly isoenergetic with a transition enthalpy of 18.2 kJ/mol in favor of the all-syn isomer.

Backbone modified TBA analogues endowed with antiproliferative activity

BACKGROUND

The thrombin binding aptamer (TBA) is endowed with antiproliferative properties but its potential development is counteracted by the concomitant anticoagulant activity.

METHODS

Five oligonucleotides (ODNs) based on TBA sequence (GGTTGGTGTGGTTGG) and containing l-residues or both l-residues and inversion of polarity sites have been investigated by NMR and CD techniques for their ability to form G-quadruplex structures. Furthermore, their anticoagulant (PT assay) and antiproliferative properties (MTT assay), and their resistance in fetal bovine serum have been tested.

RESULTS

CD and NMR data suggest that the investigated ODNs are able to form right- and left-handed G-quadruplex structures. All ODNs do not retain the anticoagulant activity characteristic of TBA but are endowed with a significant antiproliferative activity against two cancerous cell lines. Their resistance in biological environment after six days is variable, depending on the ODN.

CONCLUSIONS

A comparison between results and literature data suggests that the antiproliferative activity of the TBA analogues investigated could depends on two factors: a) biological pathways and targets different from those already identified or proposed for other antiproliferative G-quadruplex aptamers, and b) the contribution of the guanine-based degradation products.

GENERAL SIGNIFICANCE

Modified TBA analogues containing l-residues and inversion of polarity sites lose the anticoagulant activity but gain antiproliferative properties against two cancer cell lines. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

Metal Cations in G-Quadruplex Folding and Stability

This review is focused on the structural and physicochemical aspects of metal cation coordination to G-Quadruplexes (GQ) and their effects on GQ stability and conformation. G-quadruplex structures are non-canonical secondary structures formed by both DNA and RNA. G-quadruplexes regulate a wide range of important biochemical processes. Besides the sequence requirements, the coordination of monovalent cations in the GQ is essential for its formation and determines the stability and polymorphism of GQ structures. The nature, location, and dynamics of the cation coordination and their impact on the overall GQ stability are dependent on several factors such as the ionic radii, hydration energy, and the bonding strength to the O6 of guanines. The intracellular monovalent cation concentration and the localized ion concentrations determine the formation of GQs and can potentially dictate their regulatory roles. A wide range of biochemical and biophysical studies on an array of GQ enabling sequences have generated at a minimum the knowledge base that allows us to often predict the stability of GQs in the presence of the physiologically relevant metal ions, however, prediction of conformation of such GQs is still out of the realm.

Role of Alkali Metal Ions in G-Quadruplex Nucleic Acid Structure and Stability

G-quadruplexes are guanine-rich nucleic acids that fold by forming successive quartets of guanines (the G-tetrads), stabilized by intra-quartet hydrogen bonds, inter-quartet stacking, and cation coordination. This specific although highly polymorphic type of secondary structure deviates significantly from the classical B-DNA duplex. G-quadruplexes are detectable in human cells and are strongly suspected to be involved in a number of biological processes at the DNA and RNA levels. The vast structural polymorphism exhibited by G-quadruplexes, together with their putative biological relevance, makes them attractive therapeutic targets compared to canonical duplex DNA. This chapter focuses on the essential and specific coordination of alkali metal cations by G-quadruplex nucleic acids, and most notably on studies highlighting cation-dependent dissimilarities in their stability, structure, formation, and interconversion. Section 1 surveys G-quadruplex structures and their interactions with alkali metal ions while Section 2 presents analytical methods used to study G-quadruplexes. The influence of alkali cations on the stability, structure, and kinetics of formation of G-quadruplex structures of quadruplexes will be discussed in Sections 3 and 4. Section 5 focuses on the cation-induced interconversion of G-quadruplex structures. In Sections 3 to 5, we will particularly emphasize the comparisons between cations, most often K(+) and Na(+) because of their prevalence in the literature and in cells.

Unusual Chair-Like G-Quadruplex Structures: Heterochiral TBA Analogues Containing Inversion of Polarity Sites

Heterochiral oligodeoxynucleotides based on the thrombin binding aptamer sequence, namely, 5′gg3′-3′TT5′-5′ggtgtgg3′-3′TT5′-5′gg3′ (H1), 5′gg3′-3′TT5′-5′gg3′-3′TGT5′-5′gg3′-3′TT5′-5′gg3′ (H2), and 5′gGTTGgtgtgGTTGg3′ (H3), where lower case letters indicate L-residues, have been investigated in their ability to fold in G-quadruplex structures through a combination of gel electrophoresis, circular dichroism, and UV spectroscopy techniques. InH1andH2inversions of polarity sites have been introduced to control the strand direction in the loop regions. Collected data suggest that all modified sequences are able to fold in chair-like G-quadruplexes mimicking the originalTBAstructure.

Solid-state NMR studies of nucleic acid components

Recent applications of solid-state NMR spectroscopy to studies of nucleic acids and their components.

Tuning supramolecular G-quadruplexes with mono- and divalent cations

Supramolecular G-quadruplexes (SGQs) are formed via the cation promoted self-assembly of guanine derivatives into stacks of planar hydrogen-bonded tetramers. Here, we present results on the formation of SGQs made from the 8-(m-acetylphenyl)-2'-deoxyguanosine (mAGi) derivative in the presence of various mono- and divalent cations. NMR and HR ESI-MS data indicate that varying the cation can efficiently tune the molecularity, the fidelity and stability (thermal and kinetic) of the resulting SGQs. The results show that, parallel to the previously reported potassium-templated hexadecamer (mAGi₁₆·3K⁺), Na⁺, Rb⁺ and [Formula: see text] also promote the formation of similar supramolecules with high fidelity and molecularity. In contrast, the divalent cations Pb²⁺, Sr²⁺ and Ba²⁺ template the formation of octamers (mAGi₈), with the latter two inducing higher thermal stabilities. Molecular dynamics simulations for the hexadecamers containing monovalent cations enabled critical insights that help explain the experimental observations.

A straightforward modification in the thrombin binding aptamer improving the stability, affinity to thrombin and nuclease resistance

Degradation of nucleic acids in biological environments is the major drawback of the therapeutic use of aptamers. Among the approaches used to circumvent this negative aspect, the introduction of 3'-3' inversion of polarity sites at the sequence 3'-end has successfully been proposed. However, the introduction of inversion of polarity at the ends of the sequence has never been exploited for G-quadruplex forming aptamers. In this communication we describe CD, UV, electrophoretic and biochemical investigations concerning thrombin binding aptamer analogues containing one or two inversions of polarity sites at the oligonucleotide ends. Data indicate that, in some cases, this straightforward chemical modification is able to improve, at the same time, the thermal stability, affinity to thrombin and nuclease resistance in biological environments, thus suggesting its general application as a post-SELEX modification also for other therapeutically promising aptamers adopting G-quadruplex structures.

More than one non-canonical phosphodiester bond in the G-tract: formation of unusual parallel G-quadruplex structures

In this article, we report an investigation, based on NMR and CD spectroscopic and electrophoretic techniques, of 5'TGGGGT3' analogues containing two or three 3'-3' or 5'-5' inversion sites in the G-run, namely 5'TG3'-3'G5'-5'GGT3' (Q350), 5'TG3'-3'GG5'-5'GT3' (Q305), 5'TGG3'-3'G5'-5'GT3' (Q035), 5'TG3'-3'G5'-5'G3'-3'GT5' (Q353) and 3'TG5'-5'G3'-3'G5'-5'GT3' (Q535). Although the sequences investigated contain either no or only one natural 3'-5' linkage in the G-tract, all modified oligodeoxyribonucleotides (ODNs) have been shown to form stable tetramolecular quadruplex structures. The ability of the 3'-3' or 5'-5' inversion sites to affect the glycosidic conformation of guanosines and, consequently, base stacking, has also been investigated. The results of this study allow us to propose some generalizations concerning strand arrangements and the glycosidic conformational preference of residues adjacent to inverted polarity sites. These rules could be of general interest in the design of modified quadruplex structures, in view of their application as G-wires and modified aptamers.

Effects of six-membered carbohydrate rings on structure, stability, and kinetics of G-quadruplexes

We have evaluated the conformational, thermal, and kinetic properties of d(TGGGGT) analogues with one or five of the ribose nucleotides replaced with the carbohydrate residues hexitol nucleic acid (HNA), cyclohexenyl nucleic acid (CeNA), or altritol nucleic acid (ANA). All of the modified oligonucleotides formed G-quadruplexes, but substitution with the six-membered rings resulted in a mixture of G-quadruplex structures. UV and CD melting analyses showed that the structure formed by d(TGGGGT) modified with HNA was stabilized whereas that modified with CeNA was destabilized, relative to the structure formed by the unmodified oligonucleotide. Substitution at the fourth base of the G-tract with ANA resulted in a greater stabilization effect than substitution at the first G residue; substitution with five ANA residues resulted in significant stabilization of the G-quadruplex. A single substitution with CeNA at the first base of the G-tract or five substitutions with HNA resulted in striking deceleration or acceleration of G-quadruplex formation, respectively. Our results shed light on the effect of the sugar moiety on the properties of G-quadruplex structures.

Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations

Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.

A novel equilibrium relating to the helix handedness in G-quadruplexes formed by heterochiral oligonucleotides with an inversion of polarity site

Investigations of heterochiral oligodeoxynucleotides 5'-TD1GD2GD3-3'-3'-GL3GL2TL1-5' (L33) and 3'-TD1GD2GD3-5'-5'-GL3GL2TL1-3' (L55) forming quadruplex structures are reported. Data indicate the presence of enantiomeric left- and right-handed quadruplex helices. In the case of L55, NMR experiments point to an unusual equilibrium between them.

Relative stability of different DNA guanine quadruplex stem topologies derived using large-scale quantum-chemical computations

We provide theoretical predictions of the intrinsic stability of different arrangements of guanine quadruplex (G-DNA) stems. Most computational studies of nucleic acids have applied Molecular Mechanics (MM) approaches using simple pairwise-additive force fields. The principle limitation of such calculations is the highly approximate nature of the force fields. In this study, we for the first time apply accurate QM computations (DFT-D3 with large atomic orbital basis sets) to essentially complete DNA building blocks, seven different folds of the cation-stabilized two-quartet G-DNA stem, each having more than 250 atoms. The solvent effects are approximated by COSMO continuum solvent. We reveal sizable differences between MM and QM descriptions of relative energies of different G-DNA stems, which apparently reflect approximations of the DNA force field. Using the QM energy data, we propose correction to earlier free energy estimates of relative stabilities of different parallel, hybrid, and antiparallel G-stem folds based on classical simulations. The new energy ranking visibly improves the agreement between theory and experiment. We predict the 5'-anti-anti-3' GpG dinucleotide step to be the most stable one, closely followed by the 5'-syn-anti-3' step. The results are in good agreement with known experimental structures of 2-, 3-, and 4-quartet G-DNA stems. Besides providing specific results for G-DNA, our study highlights basic limitations of force field modeling of nucleic acids. Although QM computations have their own limitations, mainly the lack of conformational sampling and the approximate description of the solvent, they can substantially improve the quality of calculations currently relying exclusively on force fields.