Diagonal Loops Research Articles

Hindered intermolecular stacking of anti-parallel telomeric G-quadruplexes.

Telomeric G-quadruplexes (G4s) are non-canonical DNA structures composed of TTAGGG repeats. They are extensively studied both as biomolecules key for genome stability and as promising building blocks and functional elements in synthetic biology and nanotechnology. This is why it is extremely important to understand how the interaction between G4s is affected by their topology. We used small-angle x-ray scattering to investigate the end-to-end stacking of antiparallel telomeric G-quadruplexes formed by the sequence AG3(T2AG3)3. To represent the experimental data, we developed a highly efficient coarse-grained fitting tool, which successfully described the samples as an equilibrium mixture of monomeric and dimeric G4 species. Our findings indicate that the antiparallel topology prevents the formation of long multimeric structures under self-crowding conditions, unlike the hybrid/parallel structures formed by the same DNA sequence. This result supports the idea that the stacking of monomeric G-quadruplexes is strongly affected by the presence of diagonal loops.

Entangled World of DNA Quadruplex Folds.

DNA quadruplexes participate in many biological functions. It takes up a variety of folds based on the sequence and environment. Here, a meticulous analysis of experimentally determined 437 quadruplex structures (433 PDBs) deposited in the PDB is carried out. The analysis reveals the modular representation of the quadruplex folds. Forty-eight unique quadruplex motifs (whose diversity arises out of the propeller, bulge, diagonal, and lateral loops that connect the quartets) are identified, leading to simple to complex inter/intramolecular quadruplex folds. The two-layered structural motifs are further classified into 33 continuous and 15 discontinuous motifs. While the continuous motifs can directly be extended to a quadruplex fold, the discontinuous motif requires an additional loop(s) to complete a fold, as illustrated here with examples. Similarly, higher-order quadruplex folds can also be represented by continuous or discontinuous motifs or their combinations. Such a modular representation of the quadruplex folds may assist in custom engineering of quadruplexes, designing motif-based drugs, and the prediction of the quadruplex structure. Furthermore, it could facilitate understanding of the role of quadruplexes in biological functions and diseases.

Fluorescence of Bimolecular Guanine Quadruplexes: From Femtoseconds to Nanoseconds

The paper deals with the fluorescence of guanine quadruplexes (G4) formed by association of two DNA strands d(GGGGTTTTGGGG) in the presence of K+ cations, noted as OXY/K+ in reference to the protozoon Oxytricha nova, whose telomere contains TTTTGGGG repeats. They were studied by steady-state and time-resolved techniques, time-correlated single photon counting, and fluorescence upconversion. The maximum of the OXY/K+ fluorescence spectrum is located at 334 nm, and the quantum yield is 5.8 × 10-4. About 75% of the photons are emitted before 100 ps and stem from ππ* states, possibly with a small contribution of charge transfer. Time-resolved fluorescence anisotropy measurements indicate that ultrafast (<330 fs) excitation transfer, due to internal conversion among exciton states, is more efficient in OXY/K+ compared to previously studied G4 structures. This is attributed to the arrangement of the peripheral thymines in two diagonal loops with restricted mobility, facilitating the interaction among them and with guanines. Thymines should also be responsible for a weak intensity excimer/exciplex emission band, peaking at 445 nm. Finally, the longest living fluorescence component (∼2.1 ns) is observed at the blue side of the spectrum. So far, high-energy long-lived emitting states had been reported only for double-stranded structures but not for G4.

Efficacy of prestressed SMA diagonal loops in seismic retrofitting of non-seismically detailed RC beam-column joints

In this study, a detailed experimental and numerical investigation is carried out to study the efficacy of prestressed SMA diagonal compression loops in seismic retrofitting of non-seismically detailed RC beam-column joints. Prestressed diagonal compression loops exert active confinement to the joint core which is found to perform better than the conventional passive confinement in terms of improving the shear strength, ductility, and energy dissipation capacity. In this study, active confinement of the joint core is achieved by applying in-plane diagonal compression forces using externally fastened prestressed NiTiNb shape memory alloy (SMA) diagonal loops. In the experimental part, two different techniques are investigated to prestress the SMA loops, which are: 1) rapid heat-activated prestressing technique utilizing the shape memory feature of SMAs; and 2) conventional prestressing of SMA loops by mechanical means. In the finite element analysis, a wide range of prestress values (post-tension force) is studied, keeping the basic arrangement same. The efficacy of the retrofitting technique is evaluated in terms of enhancement in strength, ductility, energy dissipation capacity, damage reduction in the specimens, and the ease of application. The results from this study suggest that depending on the applied confinement level, the retrofitting scheme can increase the ultimate strength in the range of 20–30% and the energy dissipation capacity in the range of 60–70% however, no significant enhancement in ductility of the retrofitted specimen may be achieved using this approach. The study also highlights the complexity involved in controlling the material behaviour of NiTiNb SMA wires and the risks associated with the application of SMAs in practical problems.

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.

Probing of G-Quadruplex Structures via Ligand-Sensitized Photochemical Reactions in BrU-Substituted DNA

We studied photochemical reactions of BrU-substituted G-quadruplex (G4) DNA substrates with two pyrene-substituted polyazamacrocyclic ligands, M-1PY and M-2PY. Both ligands bind to and stabilize G4-DNA structures without altering their folding topology, as demonstrated by FRET-melting experiments, fluorimetric titrations and CD spectroscopy. Notably, the bis-pyrene derivative (M-2PY) behaves as a significantly more affine and selective G4 ligand, compared with its mono-pyrene counterpart (M-1PY) and control compounds. Upon short UVA irradiation (365 nm) both ligands, in particular M-2PY, efficiently sensitize photoreactions at BrU residues incorporated in G4 structures and give rise to two kinds of photoproducts, namely DNA strand cleavage and covalent ligand–DNA photoadducts. Remarkably, the photoinduced strand cleavage is observed exclusively with G4 structures presenting BrU residues in lateral or diagonal loops, but not with parallel G4-DNA structures presenting only propeller loops. In contrast, the formation of fluorescent photoadducts is observed with all BrU-substituted G4-DNA substrates, with M-2PY giving significantly higher yields (up to 27%) than M-1PY. Both ligand-sensitized photoreactions are specific to BrU-modified G4-DNA structures with respect to double-stranded or stem-loop substrates. Thus, ligand-sensitized photoreactions with BrU-substituted G4-DNA may be exploited (i) as a photochemical probe, allowing “photofootprinting” of G4 folding topologies in vitro and (ii) for covalent trapping of G4 structures as photoadducts with pyrene-substituted ligands.

Structural Dynamics of Lateral and Diagonal Loops of Human Telomeric G-Quadruplexes in Extended MD Simulations.

The NMR solution structures of human telomeric (Htel) G-quadruplexes (GQs) are characterized by the presence of two lateral loops complemented by either diagonal or propeller loops. Bases of a given loop can establish interactions within the loop as well as with other loops and the flanking bases. This can lead to a formation of base alignments above and below the GQ stems. These base alignments are known to affect the loop structures and relative stabilities of different Htel GQ folds. We have carried out a total of 217 μs of classical (unbiased) molecular dynamics (MD) simulations starting from the available solution structures of Htel GQs to characterize structural dynamics of the lateral and diagonal loops, using several recent AMBER DNA force-field variants. As the loops are involved in diverse stacking and H-bonding interactions, their dynamics is slow, and extended sampling is required to capture different conformations. Nevertheless, although the simulations are far from being quantitatively converged, the data suggest that multiple 10 μs-scale simulations can provide a quite good assessment of the loop conformational space as described by the force field. The simulations indicate that the lateral loops may sample multiple coexisting conformations, which should be considered when comparing simulations with the NMR models as the latter include ensemble averaging. The adenine-thymine Watson-Crick arrangement was the most stable base pairing in the simulations. Adenine-adenine and thymine-thymine base pairs were also sampled but were less stable. The data suggest that the description of lateral and diagonal GQ loops in contemporary MD simulations is considerably more realistic than the description of propeller loops, though definitely not flawless.

Quadruplexes in 'Dicty': crystal structure of a four-quartet G-quadruplex formed by G-rich motif found in the Dictyostelium discoideum genome.

Guanine-rich DNA has the potential to fold into non-canonical G-quadruplex (G4) structures. Analysis of the genome of the social amoeba Dictyostelium discoideum indicates a low number of sequences with G4-forming potential (249–1055). Therefore, D. discoideum is a perfect model organism to investigate the relationship between the presence of G4s and their biological functions. As a first step in this investigation, we crystallized the dGGGGGAGGGGTACAGGGGTACAGGGG sequence from the putative promoter region of two divergent genes in D. discoideum. According to the crystal structure, this sequence folds into a four-quartet intramolecular antiparallel G4 with two lateral and one diagonal loops. The G-quadruplex core is further stabilized by a G-C Watson–Crick base pair and a A–T–A triad and displays high thermal stability (Tm > 90°C at 100 mM KCl). Biophysical characterization of the native sequence and loop mutants suggests that the DNA adopts the same structure in solution and in crystalline form, and that loop interactions are important for the G4 stability but not for its folding. Four-tetrad G4 structures are sparse. Thus, our work advances understanding of the structural diversity of G-quadruplexes and yields coordinates for in silico drug screening programs and G4 predictive tools.

G-Quadruplex Secondary Structure Obtained from Circular Dichroism Spectroscopy.

A curated library of circular dichroism spectra of 23 G-quadruplexes of known structure was built and analyzed. The goal of this study was to use this reference library to develop an algorithm to derive quantitative estimates of the secondary structure content of quadruplexes from their experimental CD spectra. Principal component analysis and singular value decomposition were used to characterize the reference spectral library. CD spectra were successfully fit to obtain estimates of the amounts of base steps in anti-anti, syn-anti or anti-syn conformations, in diagonal or lateral loops, or in other conformations. The results show that CD spectra of nucleic acids can be analyzed to obtain quantitative structural information about secondary structure content in an analogous way to methods used to analyze protein CD spectra.

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

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

Selective and Cooperative Ligand Binding to Antiparallel Human Telomeric DNA G-Quadruplexes.

The quest for ligands that specifically bind to particular G-quadruplex nucleic acid structures is particularly important to conceive molecules with specific effects on gene expression or telomere maintenance, or conceive structure-specific molecular probes. Using electrospray mass spectrometry in native conditions, we reveal a highly cooperative and selective 2:1 binding of Cu(II) -tolylterpyridine complexes to human telomeric G-quadruplexes. Circular dichroism and comparisons of affinities for different sequences reveal a marked preference for antiparallel structures with diagonal loops and/or wide-medium-narrow-medium groove-width order. The cooperativity is attributed to conformational changes in the polymorphic telomeric G-quadruplex sequences, which convert preferably into an antiparallel three-quartet topology upon binding of two ligands.

Interactions of two cytotoxic organoruthenium(II) complexes with G-quadruplex

Two previously isolated cytotoxic complexes [(η(6)-p-cymene)Ru(κ(2)-CF3COCHCOC5H3O)L](n+) (L=Cl (1); n=0, pta (2) (pta=1,3,5-triaza-7-phosphaadamantane); n=1) were investigated for their selectivity and ability to interact with DNA G-quadruplex adopted by d[G3ATG3ACACAG4ACG3] (3) whose topology exhibits diagonal, edge-type and double-chain reversal loops. Structural changes were followed using high-resolution NMR techniques in the presence of 1 and 2. Results showed weak interaction between the organoruthenium complexes and G-quadruplex. Moreover, no significant changes in thermal stability were confirmed by a UV-melting assay for both 1 and 2. These findings emphasize that anticancer activity of Ru(II) complexes may not be correlated with binding to nucleic acid like G-quadruplex.

Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations.

DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.

Molecular docking and dynamics simulations on the interaction of cationic porphyrin–anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin–anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.

Single-molecule analysis of thymine dimer-containing G-quadruplexes formed from the human telomere sequence.

The human telomere plays crucial roles in maintaining genome stability. In the presence of suitable cations, the repetitive 5′-TTAGGG-3′ human telomere sequence can fold into G-quadruplexes that adopt the hybrid, basket, or propeller fold. The telomere sequence is hypersensitive to UV-induced thymine dimer (T=T) formation, yet it does not cause telomere shortening. In this work, the potential structural disruption and thermodynamic stability of the T=T-containing natural telomere sequences were studied to understand why this damage is tolerated in telomeres. First, established methods, such as thermal melting measurements, electrophoretic mobility shift assays, and circular dichroism spectroscopy, were utilized to determine the effects of the damage on these structures. Second, a single-molecule ion channel recording technique using α-hemolysin (α-HL) was employed to examine further the structural differences between the damaged sequences. It was observed that the damage caused slightly lower thermal stabilities and subtle changes in the circular dichroism spectra for hybrid and basket folds. The α-HL experiments determined that T=Ts disrupt double-chain reversal loop formation but are tolerated in edgewise and diagonal loops. The largest change was observed for the T=T-containing natural telomere sequence when the propeller fold (all double-chain reversal loops) was studied. On the basis of the α-HL experiments, it was determined that a triplexlike structure exists under conditions that favor a propeller structure. The biological significance of these observations is discussed.

Triplex intermediates in folding of human telomeric quadruplexes probed by microsecond-scale molecular dynamics simulations.

We have carried out extended set of μs-scale explicit solvent MD simulations of all possible G-triplexes which can participate in folding pathways of the human telomeric quadruplex. Our study accumulates almost 60μs of simulation data, which is by about three orders of magnitude larger sampling compared to the earlier simulations of human telomeric G-DNA triplexes. Starting structures were obtained from experimental quadruplex structures by deleting either the first or the last strand. The life-times of antiparallel triplexes with lateral and diagonal loops are at least on μs-scale, which should be sufficient to contribute to the folding pathways. However, the triplex states may involve structures with various local deviations from the ideal triplexes, such as strand tilting and various alternative and incomplete triads. The simulations reveal easy rearrangements between lateral and diagonal loop triplex topologies. Propeller loops of antiparallel triplexes may to certain extent interfere with the G-triplexes but these structures are still viable candidates to participate in the folding. In contrast, all-parallel all-anti triplexes are very unstable and are unlikely to contribute to the folding. Although our simulations demonstrate that antiparallel G-triplexes, if folded, would have life-times sufficient to participate in the quadruplex folding, the results do not rule out the possibility that the G-triplexes are out-competed by other structures not included in our study. Among them, numerous possible misfolded structures containing guanine quartets can act as off-path intermediates with longer life-times than the triplexes. Besides analyzing the structural dynamics of a diverse set of G-DNA triplexes, we also provide a brief discussion of the limitations of the simulation methodology, which is necessary for proper understanding of the simulation data.

Solution-state structure of an intramolecular G-quadruplex with propeller, diagonal and edgewise loops

We herein report on the formation and high-resolution NMR solution-state structure determination of a G-quadruplex adopted by d[G3ATG3ACACAG4ACG3] comprised of four G-tracts with the third one consisting of four guanines that are intervened with non-G streches of different lengths. A single intramolecular antiparallel (3+1) G-quadruplex exhibits three stacked G-quartets connected with propeller, diagonal and edgewise loops of different lengths. The propeller and edgewise loops are well structured, whereas the longer diagonal loop is more flexible. To the best of our knowledge, this is the first high-resolution G-quadruplex structure where all of the three main loop types are present.

Insight into G-DNA Structural Polymorphism and Folding from Sequence and Loop Connectivity through Free Energy Analysis

The lengths of G-tracts and their connecting loop sequences determine G-quadruplex folding and stability. Complete understanding of the sequence–structure relationships remains elusive. Here, single-loop G-quadruplexes were investigated using explicit solvent molecular dynamics (MD) simulations to characterize the effect of loop length, loop sequence, and G-tract length on the folding topologies and stability of G-quadruplexes. Eight loop types, including different variants of lateral, diagonal, and propeller loops, and six different loop sequences [d0 (i.e., no intervening residues in the loop), dT, dT2, dT3, dTTA, and dT4] were considered through MD simulation and free energy analysis. In most cases the free energetic estimates agree well with the experimental observations. The work also provides new insight into G-quadruplex folding and stability. This includes reporting the observed instability of the left propeller loop, which extends the rules for G-quadruplex folding. We also suggest a plausible explanation why human telomere sequences predominantly form hybrid-I and hybrid-II type structures in K+ solution. Overall, our calculation results indicate that short loops generally are less stable than longer loops, and we hypothesize that the extreme stability of sequences with very short loops could possibly derive from the formation of parallel multimers. The results suggest that free energy differences, estimated from MD and free energy analysis with current force fields and simulation protocols, are able to complement experiment and to help dissect and explain loop sequence, loop length, and G-tract length and orientation influences on G-quadruplex structure.

Synchrotron radiation circular dichroism of various G‐quadruplex structures

Here we report synchrotron radiation circular dichroism spectra of various G-quadruplexes from 179 to 350 nm, and a number of bands in the vacuum ultraviolet (VUV) are reported for the first time. For a tetramolecular parallel structure, the strongest band in the spectrum is a negative band in the VUV at 182 nm; for a bimolecular antiparallel structure with diagonal loops, a new strong positive band is found at 190 nm; for a bimolecular parallel structure with edgewise loops, a strong positive band at 189 nm is observed; and for a self-folded chair-type structure, the strongest band in the spectrum is a positive band at 187 nm. For the tetramolecular parallel structure, the CD signals at all wavelengths are dominated by contributions from quartets of G bases, and the signal strength is approximately proportional to the number of quartets. Our experiments on well-characterized G-quadruplex structures lead us to question past attributions of CD signals to helix handedness and G quartet polarity. Although differences can be observed in the VUV region for the various quadruplex types, there do not appear to be clear-cut spectral features that can be used to identify specific topological features. It is suggested that this is because a dominant positive band in the VUV seen near 190 nm in all quadruplex structures is due to intrastrand guanine-guanine base stacking. However, our spectra can serve as reference spectra for the G-quadruplex structures investigated and, not least, to benchmark theoretical calculations and empirical models.

Single Stranded Loops of Quadruplex DNA As Key Benchmark for Testing Nucleic Acids Force Fields.

We have carried out a set of explicit solvent molecular dynamics (MD) simulations on two DNA quadruplex (G-DNA) molecules, namely the antiparallel d(G4T4G4)2 dimeric quadruplex with diagonal loops and the parallel-stranded human telomeric monomolecular quadruplex d[AGGG(TTAGGG)3] with three propeller loops. The main purpose of the paper was testing of the capability of the MD simulation technique to describe single-stranded topologies of G-DNA loops, which represent a very challenging task for computational methods. The total amount of conventional and locally enhanced sampling (LES) simulations analyzed in this study exceeds 1.5 μs, while we tested several versions of the AMBER force field (parm99, parmbsc0, and a version with modified glycosidic χ torsion profile) and the CHARMM27 force field. Further, we compared minimal salt and excess salt simulations. Postprocessing MM-PBSA (Molecular Mechanics, Poisson-Boltzmann, Surface Area) free energy calculations are also reported. None of the presently available force fields is accurate enough in describing the G-DNA loops. The imbalance is best seen for the propeller loops, as their experimental structure is lost within a few ns of standard simulations with all force fields. Among them, parmbsc0 provides results that are clearly closest to the experimental target values but still not in full agreement. This confirms that the improvement of the γ torsional profile penalizing the γ trans substates in the parmbsc0 parametrization was a step in the right direction, albeit not sufficient to treat all imbalances. The modified χ parametrization appears to rigidify the studied systems but does not change the ultimate outcome of the present simulations. The structures obtained in simulations with the modified χ profile are predetermined by its combination with either parm99 or parmbsc0. Experimental geometries of diagonal loops of d(G4T4G4)2 are stable in standard simulations on the ∼10 ns time scale but are becoming progressively lost in longer and LES simulations. In addition, the d(G4T4G4)2 quadruplex contains, besides the three genuine binding sites for cations in the channel of its stem, also an ion binding site at each stem-loop junction. This arrangement of five cations in the quadruplex core region is entirely unstable in all 24 simulations that we attempted. Overall, our results confirm that G-DNA loops represent one of the most difficult targets for molecular modeling approaches and should be considered as reference structures in any future studies aiming to develop or tune nucleic acids force fields.