Explaining the varied glycosidic conformational, G-tract length and sequence preferences for anti-parallel G-quadruplexes

Flanking Effect on the Structure and Stability of Human Telomeric G-Quadruplex in Varying Salt Concentrations

The stability of the human telomere G-quadruplex (G4) is directly linked to cancer disease. The human telomere is mostly associated with the flanking nucleobases, which can affect the stability of G4. Hence, in this study, the effect of the flanking nucleobases in the context of their chemical nature, number, and position on the structure and stability of G4 has been investigated in varying concentrations of KCl mimicking the normal and cancer KCl microenvironments. The addition of flanking nucleobases does not alter the G4 topology. However, the presence of merely a single flanking nucleobase destabilizes the telomeric G4. This destabilizing effect is more prominent for thymine than adenine flanking nucleobase, probably due to the formation of the intermolecular G4 topology by thymine. Interestingly, the change in the stability of the telomeric G4 in the presence of thymine flanking nucleobase is sensitive to the concentration of KCl relevant to the normal and cancerous microenvironments, in contrast to adenine. Flanking nucleobases have a greater impact at the 5' end compared to the 3' end, particularly noticeable in KCl concentrations resembling the normal microenvironment rather than the cancerous one. These findings indicate that the effect of the flanking nucleobases on telomeric G4 is different in the KCl salt relevant to normal and cancerous microenvironments. This study may be helpful in attaining molecular-level insight into the role of G4 in telomeric length regulation under normal and cancerous KCl salt conditions.

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.

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.

Copper(II) coordination compounds based on bis-hydrazones of 2,6-diacetylpyridine: synthesis, structure, and cytotoxic activity

By reacting a series of 2,6-diacetylpyridine bis-hydrazones containing pyrimidine (H2L1), benzimidazole (H2L2) and phthalazine (H2L3) heterocyclic fragments with copper(II) chloride and bromide, a variety of pentacoordinated complexes of the composition [Cu(H2L1)X]X, [Cu(HL2)X] and [Cu(HL3)X], where X = Cl-, Br-, are formed. The properties and structure of the compounds were studied by means of NMR, IR, UV-vis, ESR, and X-ray absorption spectroscopy, cyclic voltammetry and X-Ray single crystal diffraction methods. It was shown that complexes of the cationic type [Cu(H2L1)X]X have an asymmetric structure with a distorted square-pyramidal geometry of the coordination unit. The coordination polyhedron of metal chelates [Cu(HL2)X] and [Cu(HL3)X] is an almost ideal square pyramid. Investigations of the cytotoxic activity of the obtained compounds in vitro on human hepatocellular carcinoma (HepG2) and non-tumor human lung fibroblast (MRC-5) cell lines demonstrated that complexes show higher activity compared with the well-known anticancer agent cisplatin. In addition, metal chelates [Cu(H2L1)Cl]Cl, [Cu(HL2)Cl], [Cu(HL2)Br] and [Cu(HL3)Cl] were less toxic to non-tumor cells MRC-5. A study of the binding of complexes to bovine serum albumin (BSA) protein using fluorescence spectroscopy showed that copper complexes are strongly bound to BSA. To study the mechanism of interaction of the complexes with the DNA of cancer cells, molecular dynamics simulation of the compound [Cu(HL3)Cl] was carried out. It was shown that the complex enters into π-stacking interactions predominantly with adenine and thymine bases.

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.

Affinity Chromatography-Based Assays for the Screening of Potential Ligands Selective for G-Quadruplex Structures

DNA G‐quadruplexes (G4s) are key structures for the development of targeted anticancer therapies. In this context, ligands selectively interacting with G4s can represent valuable anticancer drugs. Aiming at speeding up the identification of G4‐targeting synthetic or natural compounds, we developed an affinity chromatography‐based assay, named G‐quadruplex on Oligo Affinity Support (G4‐OAS), by synthesizing G4‐forming sequences on commercially available polystyrene OAS. Then, due to unspecific binding of several hydrophobic ligands on nude OAS, we moved to Controlled Pore Glass (CPG). We thus conceived an ad hoc functionalized, universal support on which both the on‐support elongation and deprotection of the G4‐forming oligonucleotides can be performed, along with the successive affinity chromatography‐based assay, renamed as G‐quadruplex on Controlled Pore Glass (G4‐CPG) assay. Here we describe these assays and their applications to the screening of several libraries of chemically different putative G4 ligands. Finally, ongoing studies and outlook of our G4‐CPG assay are reported.

Mechanistic Insights into Multiple-step Transport of Mitochondrial ADP/ATP Carrier

The ADP/ATP carrier (AAC) is crucial for mitochondrial functions by importing ADP and exporting ATP across the inner mitochondrial membrane. However, the mechanism of highly specific ADP recognition and transport by AAC remains largely elusive. In this work, spontaneous ADP binding process to the ground c-state AAC was investigated through rigorous molecular dynamics simulations of over 31 microseconds in total. With improved simulation strategy, we have successfully identified a highly specific ADP binding site in the upper region of the cavity, and this site exhibits selectivity for ADP over ATP based on free-energy calculations. Sequence analyses on adenine nucleotide transporters also suggest that this subgroup uses the upper region of the cavity, rather than the previously proposed central binding site located at the bottom of the cavity to discriminate their substrates. Identification of the new site unveils the unusually high substrate specificity of AAC and explains the dependence of transport on the flexibility between anti and syn glycosidic conformers of ADP. Moreover, this new site together with the central site supports early biochemical findings. In light of these early findings, our simulations described a multi-step model in which the carrier uses different sites for substrate attraction, recognition and conformational transition. These results provide new insights into the transport mechanism of AAC and other adenine nucleotide transporters.

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.

Natural compounds from Juncus plants interacting with telomeric and oncogene G-quadruplex structures as potential anticancer agents

Aiming at discovering novel, putative anticancer drugs featuring low-to-null side effects, natural compounds isolated from Juncaceae were studied here for their ability to target G-quadruplex structures originating from cancer-related telomeric and oncogene DNA sequences. Particularly, various dihydrophenanthrene, benzocoumarin and dihydrodibenzoxepin derivatives were firstly screened by the affinity chromatography-based G4-CPG assay, and the compound with the highest affinity and selectivity for G-quadruplexes (named J10) was selected for further studies. Fluorescence spectroscopy and circular dichroism experiments corroborated its capability to selectively recognize and stabilize G-quadruplexes over duplex DNA, also showing a preference for parallel G-quadruplexes. Molecular docking proved that the selective G-quadruplex interactions over duplex interactions could be due to the ability of J10 to bind to the grooves of the telomeric and oncogene G-quadruplex structures. Finally, biological assays demonstrated that J10 induces significant antiproliferative effects on human leukemia cells, with no relevant effects on healthy human fibroblasts. Interestingly, J10 exerts its antiproliferative action on tumor cells by activating the apoptotic pathway.

Combining Electrospray Mass Spectrometry (ESI-MS) and Computational Techniques in the Assessment of G-Quadruplex Ligands: A Hybrid Approach to Optimize Hit Discovery

Guanine-rich sequences forming G-quadruplexes (GQs) are present in several genomes, ranging from viral to human. Given their peculiar localization, the induction of GQ formation or GQ stabilization with small molecules represents a strategy for interfering with crucial biological functions. Investigating the recognition event at the molecular level, with the aim of fully understanding the triggered pharmacological effects, is challenging. Native electrospray ionization mass spectrometry (ESI-MS) is being optimized to study these noncovalent assemblies. Quantitative parameters retrieved from ESI-MS studies, such as binding affinity, the equilibrium binding constant, and sequence selectivity, will be overviewed. Computational experiments supporting the ESI-MS investigation and boosting its efficiency in the search for GQ ligands will also be discussed with practical examples. The combination of ESI-MS and in silico techniques in a hybrid high-throughput-screening workflow represents a valuable tool for the medicinal chemist, providing data on the quantitative and structural aspects of ligand–GQ interactions.

A world beyond double-helical nucleic acids: the structural diversity of tetra-stranded G-quadruplexes

Nucleic acids can adopt various secondary structures including double-, triple-, and tetra-stranded helices that differ by the specific hydrogen bond mediated pairing pattern between their nucleobase constituents. Whereas double-helical DNA relies on Watson–Crick base pairing to play a prominent role in storing genetic information, G-quadruplexes are tetra-stranded structures that are formed by the association of guanine bases from G-rich DNA and RNA sequences. During the last few decades, G-quadruplexes have attracted considerable interest after the realization that they form and exert regulatory functions in vivo. In addition, quadruplex architectures have also been recognized as versatile and powerful tools in a growing number of technological applications. To appreciate the astonishing structural diversity of these tetra-stranded structures and to give some insight into basic interactions that govern their folding, this article gives an overview of quadruplex structures and rules associated with the formation of different topologies. A brief discussion will also focus on nonconventional quadruplexes as well as on general principles when targeting quadruplexes with ligands.

Graphic abstract

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.

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.

NMR structural study on the self-trimerization of d(GTTAGG) into a dynamic trimolecular G-quadruplex assembly preferentially in Na+ solution with a moderate K+ tolerance

Vast G-quadruplexes (GQs) are primarily folded by one, two, or four G-rich oligomers, rarely with an exception. Here, we present the first NMR solution structure of a trimolecular GQ (tri-GQ) that is solely assembled by the self-trimerization of d(GTTAGG), preferentially in Na⁺ solution tolerant to an equal amount of K⁺ cation. Eight guanines from three asymmetrically folded strands of d(GTTAGG) are organized into a two-tetrad core, which features a broken G-column and two width-irregular grooves. Fast strand exchanges on a timescale of second at 17°C spontaneously occur between folded tri-GQ and unfolded single-strand of d(GTTAGG) that both species coexist in dynamic equilibrium. Thus, this tri-GQ is not just simply a static assembly but rather a dynamic assembly. Moreover, another minor tetra-GQ that has putatively tetrameric (2+2) antiparallel topology becomes noticeable only at an extremely high strand concentration above 18 mM. The major tri-GQ and minor tetra-GQ are considered to be mutually related, and their reversible interconversion pathways are proposed accordingly. The sequence d(GTTAGG) could be regarded as either a reading frame shifted single repeat of human telomeric DNA or a 1.5 repeat of Bombyx mori telomeric DNA. Overall, our findings provide new insight into GQs and expect more functional applications.

Oxidative lesions modulate G-quadruplex stability and structure in the human BCL2 promoter

Misregulation of BCL2 expression has been observed with many diseases and is associated with cellular exposure to reactive oxygen species. A region upstream of the P1 promoter in the human BCL2 gene plays a major role in regulating transcription. This G/C-rich region is highly polymorphic and capable of forming G-quadruplex structures. Herein we report that an oxidative event simulated with an 8-oxo-7,8-dihydroguanine (^oxoG) substitution within a long G-tract results in a reduction of structural polymorphism. Surprisingly, ^oxoG within a 25-nt construct boosts thermal stability of the resulting G-quadruplex. This is achieved by distinct hydrogen bonding properties of ^oxoG, which facilitate formation of an antiparallel basket-type G-quadruplex with a three G-quartet core and a G·^oxoG·C base triad. While ^oxoG has previously been considered detrimental for G-quadruplex formation, its stabilizing effect within a promoter described in this study suggests a potential novel regulatory role of oxidative stress in general and specifically in BCL2 gene transcription.

An atomistic investigation on the interaction of distamycin A and its derivative with the telomeric G-Quadruplex as anticancer agents by molecular dynamics simulation

Human telomerase that activates within cancer cells has a telomeric sequence at the 3' end. Each factor that stabilizes the G-quadruplex in guanine-rich telomeric sequences can inhibit the regular telomerase activity. Therefore, the telomeric G-quadruplex is known as a promising target in cancer treatment. In this work, we studied the binding of positively charged distamycin A and its uncharged derivative to the G-quadruplex in a solution environment by Molecular Dynamics (MD) simulation. The binding mechanism and subtle conformational changes were investigated as a result of the ligand attachment. Moreover, binding free energy and clustering analysis describe the stability and flexibility of G-quadruplexes upon ligand binding. Structural analyses displayed that the favorable binding of both ligands imposes significant stability and rigidity in G-quadruplex conformation compared to free G-quadruplex, especially charged distamycin. Hydration pattern and ion distribution were different for free G-quadruplex and both of the ligand complexes. Energy decomposition reveals the electrostatic effect on the stability of G-quadruplex. The radial distribution function displayed the solvent shell and ion moving away from the groove. The hydrogen bond played an essential role in the binding of both ligands, especially for the charged derivative. van der Waals interaction is the only factor that is more important in binding uncharged distamycin into G-quadruplex than the charged one. The calculated ΔG_bind showed the stability of both ligands within grooves and good agreement with the experimental binding free energy data. Finally, the results suggest that ligand modification improves the binding mode toward stabilizing G-quadruplexes.

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.

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.

Targeting Telomeres: Molecular Dynamics and Free Energy Simulation of Gold-Carbene Binding to DNA

DNA sequences in regulatory regions and in telomers at the ends of chromosomes frequently contain tandem repeats of guanine nucleotides that can form stacked structures stabilized by Hoogsten pairing and centrally bound monovalent cations. The replication and elongation of telomeres requires the disruption of these G-quadruplex structures. Hence, drug molecules such as gold (Au)-carbene that stabilize G-quadruplexes may also interfere with the elongation of telomeres and, in turn, could be used to control cell replication and growth. To better understand the molecular mechanism of Au-carbene binding to G-quadruplexes, we employed molecular dynamics simulations and free energy simulations. Whereas very restricted mobility of two Au-carbene ligands was found upon binding as a doublet to one side of the G-quadruplex, much larger translational and orientational mobility was observed for a single Au-carbene binding at the second G-quadruplex surface. Comparative simulations on duplex DNA in the presence of Au-carbene ligands indicates a preference for the minor groove and weaker unspecific and more salt-dependent binding than to the G-quadruplex surface. Analysis of energetic contributions reveals a dominance of nonpolar and van der Waals interactions to drive binding. The simulations can also be helpful for proposing possible modifications that could improve Au-carbene affinity and specificity for G-quadruplex binding.

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.

Structure of two G-quadruplexes in equilibrium in the KRAS promoter

KRAS is one of the most mutated oncogenes and still considered an undruggable target. An alternative strategy would consist in targeting its gene rather than the protein, specifically the formation of G-quadruplexes (G4) in its promoter. G4 are secondary structures implicated in biological processes, which can be formed among G-rich DNA (or RNA) sequences. Here we have studied the major conformations of the commonly known KRAS 32R, or simply 32R, a 32 residue sequence within the KRAS Nuclease Hypersensitive Element (NHE) region. We have determined the structure of the two major stable conformers that 32R can adopt and which display slow equilibrium (>ms) with each other. By using different biophysical methods, we found that the nucleotides G9, G25, G28 and G32 are particularly implicated in the exchange between these two conformations. We also showed that a triad at the 3' end further stabilizes one of the G4 conformations, while the second conformer remains more flexible and less stable.

Transcription-facilitating histone chaperons interact with genomic and synthetic G4 structures

Affinity for G-quadruplex (G4) structures may be a common feature of transcription-facilitating histone chaperons (HCs). This assumption is based on previous unmatched studies of HCs FACT, nucleolin (NCL), BRD3, and ATRX. We verified this assumption and considered its implications for the therapeutic applications of synthetic (exogenous) G4s and the biological significance of genomic G4s. First, we questioned whether exogenous G4s that recognize cell-surface NCL and could trap other HCs in the nucleus are usable as anticancer agents. We performed in vitro binding assays and selected leading multi-targeted G4s. They exhibited minor effects on cell viability. The presumed NCL-regulated intracellular transport of G4s was inefficient or insufficient for tumor-specific G4 delivery. Next, to clarify whether G4s in the human genome could recruit HCs, we compared available HC ChIP-seq data with G4-seq/G4-ChIP-seq data. Several G4s, including the well-known c-Myc quadruplex structure, were found to be colocalized with HC occupancy sites in cancer cell lines. As evidenced by our molecular modeling data, c-Myc G4 might interfere with the HC function of BRD3 but is unlikely to prevent the BRD3-driven assembly of the chromatin remodeling complex. The c-Myc case illustrates the intricate role of genomic G4s in chromatin remodeling, nucleosome remodeling, and transcription.

A DNA G-quadruplex/i-motif hybrid

DNA can form many structures beyond the canonical Watson–Crick double helix. It is now clear that noncanonical structures are present in genomic DNA and have biological functions. G-rich G-quadruplexes and C-rich i-motifs are the most well-characterized noncanonical DNA motifs that have been detected in vivo with either proscribed or postulated biological roles. Because of their independent sequence requirements, these structures have largely been considered distinct types of quadruplexes. Here, we describe the crystal structure of the DNA oligonucleotide, d(CCAGGCTGCAA), that self-associates to form a quadruplex structure containing two central antiparallel G-tetrads and six i-motif C–C⁺ base pairs. Solution studies suggest a robust structural motif capable of assembling as a tetramer of individual strands or as a dimer when composed of tandem repeats. This hybrid structure highlights the growing structural diversity of DNA and suggests that biological systems may harbor many functionally important non-duplex structures.

Fine-Tuning of the AMBER RNA Force Field with a New Term Adjusting Interactions of Terminal Nucleotides

Determination of RNA structural-dynamic properties is challenging for experimental methods. Thus, atomistic molecular dynamics (MD) simulations represent a helpful technique complementary to experiments. However, contemporary MD methods still suffer from limitations of force fields (ffs), including imbalances in the nonbonded ff terms. We have recently demonstrated that some improvement of state-of-the-art AMBER RNA ff can be achieved by adding a new term for H-bonding called gHBfix, which increases tuning flexibility and reduces risk of side-effects. Still, the first gHBfix version did not fully correct simulations of short RNA tetranucleotides (TNs). TNs are key benchmark systems due to availability of unique NMR data, although giving too much weight on improving TN simulations can easily lead to overfitting to A-form RNA. Here we combine the gHBfix version with another term called tHBfix, which separately treats H-bond interactions formed by terminal nucleotides. This allows to refine simulations of RNA TNs without affecting simulations of other RNAs. The approach is in line with adopted strategy of current RNA ffs, where the terminal nucleotides possess different parameters for terminal atoms than the internal nucleotides. Combination of gHBfix with tHBfix significantly improves the behavior of RNA TNs during well-converged enhanced-sampling simulations using replica exchange with solute tempering. TNs mostly populate canonical A-form like states while spurious intercalated structures are largely suppressed. Still, simulations of r(AAAA) and r(UUUU) TNs show some residual discrepancies with primary NMR data which suggests that future tuning of some other ff terms might be useful. Nevertheless, the tHBfix has a clear potential to improve modeling of key biochemical processes, where interactions of RNA single stranded ends are involved.

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.

Understanding The Role of Reline, a Natural DES, on Temperature-Induced Conformational Changes of C-Kit G-Quadruplex DNA: A Molecular Dynamics Study

The noncanonical guanine-rich DNAs have drawn particular attention to the scientific world due to their controllable diverse and polymorphic structures. Apart from biological and medical significance, G-quadruplex DNAs are widely used in various fields such as nanotechnology, nanomachine, biosensors, and biocatalyst. So far, the applications of the G-quadruplex DNA are mainly limited in the water medium. Recently, a new generation of solvent named deep eutectic solvent (DES) has become very popular and has been widely used as a reaction medium of biocatalytic reactions and long-term storage medium for nucleic acids, even at high temperature. Hence, it is essential to understand the role of DES on temperature-induced conformational changes of a G-quadruplex DNA. In this research work, we have explored the temperature-mediated conformational dynamics of c-kit oncogene promoter G-quadruplex DNA in reline medium in the temperature range of 300-500 K, using a total of 10 μs unbiased all-atom molecular dynamics simulation. Here, from RMSD, RMSF, R_g and principal component analyses, we notice that the c-kit G-quadruplex DNA is stable up to 450 K in reline medium. However, it unfolds in water medium at 450 K. It is found that the hydrogen bonding interactions between c-kit G-quadruplex DNA and reline play a key role in the stabilization of the G-quadruplex DNA even at high temperature. Furthermore, in this work we have observed a very interesting and distinctive phenomenon of the central cation of the G-quadruplex DNA. Its position was seen to fluctuate between the two tetrad cores, that is, the region between tetrad-1 and tetrad-2 and that between tetrad-2 and tetrad-3 and vice versa at 450 and 500 K in reline medium which is absent in water medium at 450 K. Moreover, the rate of its oscillation is increased when temperature is increased.

Revealing conformational dynamics of 2'-O-methyl-RNA guanine modified G-quadruplex by replica exchange molecular dynamics

Thrombin-binding DNA aptamer (TBA) can fold into an antiparallel unimolecular G-quadruplex (G4) structure. Different types of modifications lead to various effects on the structure and stability of the G4 structure. Previous study has shown that a modified TBA (mTBA) that 2'-deoxy guanine (dG) at positions 10 and 11 in the TBA sequence were replaced by 2'-O-methyl-RNA guanine (2'OMe-G) can't fold into a well-defined G4 structure. In order to investigate the detailed structural information and probe the instability factors, we successfully employed the replica exchange molecular dynamics (REMD) to characterize the large conformational variations of the mTBA and systemically describe the influences of the 2'OMe-G on the mTBA in terms of conformation variations and the probability distributions of Hoogsteen hydrogen bonds, dihedral, sugar pucker and glycosyl torsion angle. Replacing position 10 with the 2'OMe-G (2'OMe-G10) induced a strong destabilization of the aptamer, while the 2'OMe-G at position 11(2'OMe-G11) was less destabilizing. More importantly, the glycosyl torsion angle and sugar pucker of 2'OMe-G10 were the most critical destabilization factors. These results were in good agreement with the theoretical and experimental results. Moreover, the structure information can be used as guidelines for the further design of modifications on G4 structure.

Effects of Central Loop Length and Metal Ions on the Thermal Stability of G-Quadruplexes

The central loop of G-quadruplex molecular beacons is a key element to sense target DNA or RNA sequences. In this study, circular dichroism spectroscopy (CD), thermal difference spectrum (TDS), non-denatured non-denaturing gel electrophoresis, and thermal stability analysis were used to investigate the effect of the central loop length on G-quadruplex features. Two series of G-quadruplexes, AG₃TTAG₃-(TTA)n-G₃TTAG₃T (n = 1–8) (named TTA series) and AG₃TTTG₃-(TTA)n-G₃TTTG₃T (n = 1–8) (named TTT series) were examined in K⁺ and Na⁺ solutions, respectively. CD and TDS spectral data indicated that TTA series adopted an antiparallel G-quadruplex structure in Na⁺ solution and a hybrid G-quadruplex structure in K⁺ solution respectively. TTT series exhibited a hybrid G-quadruplex structure in both Na⁺ and K⁺ solutions. UV melting curves indicated that the stability of G-quadruplex in both series was reduced by the elongation of central loop. Thermal stability analysis concluded that the G-quadruplex destabilization with long central loop is an entropy-driven process due to more flexible and longer central loops.

G-Quadruplexes in Human Ribosomal RNA

rRNA is the single most abundant polymer in most cells. Mammalian rRNAs are nearly twice as large as those of prokaryotes. Differences in rRNA size are due to expansion segments, which contain extended tentacles in metazoans. Here we show that the terminus of an rRNA tentacle of Homo sapiens contains 10 tandem G-tracts that form highly stable G-quadruplexes in vitro. We characterized rRNA of the H. sapiens large ribosomal subunit by computation, circular dichroism, UV melting, fluorescent probes, nuclease accessibility, electrophoretic mobility shifts, and blotting. We investigated Expansion Segment 7 (ES7), oligomers derived from ES7, intact 28S rRNA, 80S ribosomes, and polysomes. We used mass spectrometry to identify proteins that bind to rRNA G-quadruplexes in cell lysates. These proteins include helicases (DDX3, CNBP, DDX21, DDX17) and heterogeneous nuclear ribonucleoproteins. Finally, by multiple sequence alignments, we observe that G-quadruplex-forming sequences are a general feature of LSU rRNA of Chordata but not, as far as we can tell, of other species. Chordata ribosomes present polymorphic tentacles with the potential to switch between inter- and intramolecular G-quadruplexes. To our knowledge, G-quadruplexes have not been reported previously in ribosomes.

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.

Impact of Oxidative Lesions on the Human Telomeric G-Quadruplex

Telomere attrition is closely associated with cell aging and exposure to reactive oxygen species (ROS). While oxidation products of nucleotides have been studied extensively in the past, the underlying secondary/tertiary structural changes in DNA remain poorly understood. In this work, we systematically probed guanine positions in the human telomeric oligonucleotide sequence (hTel) by substitutions with the major product of ROS, 8-oxo-7,8-dihydroguanine (^oxoG), and evaluated the G-quadruplex forming ability of such oligonucleotides. Due to reduced hydrogen-bonding capability caused by ^oxoG, a loss of G-quadruplex structure was observed for most oligonucleotides containing oxidative lesions. However, some positions in the hTel sequence were found to tolerate substitutions with ^oxoG. Due to ^oxoG’s preference for the syn conformation, distinct responses were observed when replacing guanines with different glycosidic conformations. Accommodation of ^oxoG at sites originally in syn or anti in nonsubstituted hTel G-quadruplex requires a minor structural rearrangement or a major conformational shift, respectively. The system responds by retaining or switching to a fold where ^oxoG is in syn conformation. Most importantly, these G-quadruplex structures are still stable at physiological temperatures and should be considered detrimental in higher-order telomere structures.

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.

Conformational deviation of Thrombin binding G-quadruplex aptamer (TBA) in presence of divalent cation Sr2+: A classical molecular dynamics simulation study

Thrombin binding TBA-G-quadruplex aptamer (TBA) plays a major role in blood coagulation cascade. The 15-mer TBA sequence tends to form four-stranded TBA-G-quadruplex structure. In this research work, a series of explicit solvent classical MD simulations of the TBA is carried out using different salt (SrCl₂) concentrations (0, 50, 100 and 200 mM). Here we have also testified the effect of salt concentration of divalent cation Sr²⁺ on the conformational change of quadruplex DNA. The structural deviations, fluctuations, torsional angles and the affinity of the ion are explored at different salt concentrations. It is found that the conformation of TBA-G-quadruplex at 0 mM and 50 mM salt concentrations, is very much different than the other salt concentrations (100 mM and 200 mM). Also observed are as follows: (i) no exchange of Sr²⁺ ion between inside and outside of the channel, (ii) an enhancement in the Sr²⁺ ion density around the phosphate region of the loop residues as salt concentration increases and (iii) the stacking of T3 and T4 residues of loop-1 that appears up to 50 mM concentration, vanishes as the salt concentration is increased further.

Encoding canonical DNA quadruplex structure

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

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.

Structural switch from a multistranded G-quadruplex to single strands as a consequence of point mutation in the promoter of the human GRIN1 gene

A huge number of G-rich sequences forming quadruplexes are found in the human genome, especially in telomeric regions, UTRs, and the promoter regions of a number of genes. One such gene is GRIN1 encoding the NR1 subunit of the N-methyl-d-aspartate receptor (NMDA). Several lines of reports have implicated that attenuated function of NMDA results in schizophrenia, a genetic disorder characterized by hallucinations, delusions, and psychosis. Involvement of the GRIN1 gene in the pathogenesis of schizophrenia has been extensively analysed. Recent reports have demonstrated that polymorphism in the promoter region of GRIN1 at position -855 (G/C) has a possible association with schizophrenia. The binding site for the NF-κB transcription factor gets altered due to this mutation, resulting in reduced gene expression as well as NMDA activity. By combining gel electrophoresis (PAGE), circular dichroism (CD) and CD melting techniques, the G → C single nucleotide polymorphism (SNP) at the G-rich sequence (d-CTTAGCCCGAGGAG[combining low line]GGGGGTCCCAAGT; GRIN1) was investigated. We report that the GRIN1 sequence can form an octameric/multistranded quadruplex structure with parallel conformation in the presence of K⁺ as well as Na⁺. CD and gel studies are in good correlation in order to detect molecularity and strand conformation. The parallel G-quadruplex species was hypothesized to be octameric in K⁺/Na⁺ salts. The mutated sequence (d-CTTAGCCCGAGGAC[combining low line]GGGGGTCCCAAGT; GRIN1M) remained single stranded under physiological conditions. CD melting studies support the formation of an interstranded G-quadruplex structure by the GRIN1 sequence. Two structural models are propounded for a multistranded parallel G-quadruplex conformation which might be responsible for regulating the gene expression normally underlying memory and learning.

Structure-Property Relationships in CuII -Binding Tetramolecular G-Quadruplex DNA

A series of artificial metal-base tetrads composed of a Cu^II cation coordinating to four pyridines, covalently attached to the ends of tetramolecular G-quadruplex DNA strands [L^A-D d(G₄ )]₄ (L^A-D =ligand derivatives), was systematically studied. Structurally, the square-planar [Cu(pyridine)₄ ] complex behaves analogously to the canonical guanine quartet. Copper coordination to all studied ligand derivatives was found to increase G-quadruplex thermodynamic stability, tolerating a great variety of ligand linker lengths (1-5 atoms) and thus demonstrating the robustness of the chosen ligand design. Only at long linker lengths, the stabilizing effect of copper binding is compensated by the loss of conformational freedom. A previously reported ligand L^E with chiral backbone enables incorporation at any oligonucleotide position. We show that ligand chirality distinctly steers Cu^II -induced G-quadruplex stabilization. 5'-End formation of two metal-base tetrads by tetramolecular G-quadruplex [L^E₂ d(G)₄ ]₄ shows that stabilization in the presence of Cu^II is not additive. All results are based on UV/Vis thermal denaturation, thermal difference, circular dichroism experiments and molecular dynamics simulations.

Molecular mechanisms by which oxidative DNA damage promotes telomerase activity

Telomeres are highly susceptible to oxidative DNA damage, which if left unrepaired can lead to dysregulation of telomere length homeostasis. Here we employed single molecule FRET, single molecule pull-down and biochemical analysis to investigate how the most common oxidative DNA lesions, 8-oxoguanine (8oxoG) and thymine glycol (Tg), regulate the structural properties of telomeric DNA and telomerase extension activity. In contrast to 8oxoG which disrupts the telomeric DNA structure, Tg exhibits substantially reduced perturbation of G-quadruplex folding. As a result, 8oxoG induces high accessibility, whereas Tg retains limited accessibility, of telomeric G-quadruplex DNA to complementary single stranded DNA and to telomere binding protein POT1. Surprisingly, the Tg lesion stimulates telomerase loading and activity to a similar degree as an 8oxoG lesion. We demonstrate that this unexpected stimulation arises from Tg-induced conformational alterations and dynamics in telomeric DNA. Despite impacting structure by different mechanisms, both 8oxoG and Tg enhance telomerase binding and extension activity to the same degree, potentially contributing to oncogenesis.

Efficient Long-Range Hole Transport Through G-Quadruplexes

DNA offers a means of long-range charge transport for biology and electric nanodevices. Here, a series of tetra-stranded G-quadruplexes were assembled within a dendritic DNA architecture to explore oxidative charge transport (hole transport) through the G-quadruplex. Efficient charge transport was achieved over 28 Å upon UV irradiation. Over a longer G-quadruplex bridge, hole transport was escalated to a higher efficiency, which resulted in a higher yield than that of the optimal duplex DNA for charge transport, that is, the adenine tract. Efficient long-range hole transport suggests tetra-stranded G-quadruplexes, instead of an oxidation hotspot, hold better potential as an electron conduit than duplex DNA.

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.

Exploring the Dynamics of Propeller Loops in Human Telomeric DNA Quadruplexes Using Atomistic Simulations

We have carried out a series of extended unbiased molecular dynamics (MD) simulations (up to 10 μs long, ∼162 μs in total) complemented by replica-exchange with the collective variable tempering (RECT) approach for several human telomeric DNA G-quadruplex (GQ) topologies with TTA propeller loops. We used different AMBER DNA force-field variants and also processed simulations by Markov State Model (MSM) analysis. The slow conformational transitions in the propeller loops took place on a scale of a few μs, emphasizing the need for long simulations in studies of GQ dynamics. The propeller loops sampled similar ensembles for all GQ topologies and for all force-field dihedral-potential variants. The outcomes of standard and RECT simulations were consistent and captured similar spectrum of loop conformations. However, the most common crystallographic loop conformation was very unstable with all force-field versions. Although the loss of canonical γ-trans state of the first propeller loop nucleotide could be related to the indispensable bsc0 α/γ dihedral potential, even supporting this particular dihedral by a bias was insufficient to populate the experimentally dominant loop conformation. In conclusion, while our simulations were capable of providing a reasonable albeit not converged sampling of the TTA propeller loop conformational space, the force-field description still remained far from satisfactory.

Weak Supramolecular Interactions Governing Parallel and Antiparallel DNA Quadruplexes: Insights from Large-Scale Quantum Mechanics Analysis of Experimentally Derived Models

The topology and energetics of guanine (G) quadruplexes is governed by supramolecular interactions within their strands. In this work, an extensive quantum mechanical (QM) study has been performed to analyze supramolecular interactions that shape the stems of (4+0) parallel (P) and (2+2) antiparallel (AP) quadruplex systems. The large-scale (≈400 atoms) models of P and AP were constructed from high-quality experimental structures. The results provide evidence that each of the P and AP structures is shaped by a distinct network of supramolecular interactions. Analysis of electron topological characteristics of hydrogen bonds in P and AP systems indicates that the P model benefits from stronger intratetrad hydrogen bonding. For intertetrad stacking interactions, both noncovalent interaction plot and energy decomposition analysis approaches suggest that the stem of the P quadruplex benefits more from stacking than that of the AP stem; the difference in energetic stabilization for the two topologies is about 10 %. Stronger hydrogen-bonding and stacking interactions in the stem of the P quadruplex, relative to those in the AP system, can be an important indicator to explain the experimental observations that guanine-rich oligonucleotides tend to form all-parallel stems with an all-anti orientation of nucleobases. However, in addition to intrinsic stabilization, partial desolvation effects, which affect the energetics and dynamics of the G-quadruplex folding process, call for further investigations.

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.

Quantum mechanical investigation of G-quartet systems of DNA

Minima of the electric field and positions of K⁺ and Na⁺ (zero of the x-coordinate is the center of the cavity).