Small molecule binding to a G-hairpin and a G-triplex: a new insight into anticancer drug design targeting G-rich regions

Interconversion between guanine quartets and triads on the Au(111) surface

Based on STM imaging and DFT calculations, we show the real-space experimental evidence of the interconversion between G-quartets and G-triads on the Au(111) surface, and further reveal the relative stabilities of these two elementary motifs, which helps to increase the fundamental understanding of the relationship between G-triplex and G-quadruplex DNA structures.

Stimuli-Responsive and Reversible Nanoassemblies of G-Triplexes

G-triplex (G3) structures formed with three consecutive G-tracts have recently been identified as a new emerging guanine-rich DNA fold. There could likely be a wide range of biological functions for G3s as occurring for G-quadruplex (G4) structures formed with four consecutive G-tracts. However, in comparison to the many reports on G4 nanoassemblies that organize monomers together in a controllable manner, G3-favored nanoassemblies have yet to be explored. In this work, we found that a natural alkaloid of sanguinarine can serve as a dynamic ligand glue to reversibly switch the dimeric nanoassemblies of the thrombin binding aptamer G3 (TBA-G3). The glue planarity was considered to be a crucial factor for realizing this switching. More importantly, external stimuli including pH, sulfite, O₂ and H₂ O₂ can be employed as common regulators to easily modulate the glue's adhesivity for constructing and destructing the G3 nanoassemblies as a result of the ligand converting between isoforms. However, this assembly behavior does not occur with the counterpart TBA-G4. Our work demonstrates that higher-order G3 nanoassemblies can be reversibly operated by manipulating ligand adhesivity. This provides an alternative understanding of the unique behavior of guanine-rich sequences and focuses attention on the G3 fold since the nanoassembly event investigated herein might occur in living cells.

Docking and scoring for nucleic acid-ligand interactions: Principles and current status

Nucleic acid (NA)-ligand interactions have crucial roles in many cellular processes and, thus, are increasingly attracting therapeutic interest in drug discovery. Molecular docking is a valuable tool for studying molecular interactions. However, because NAs differ significantly from proteins in both their physical and chemical properties, traditional docking algorithms and scoring functions for protein-ligand interactions might not be applicable to NA-ligand docking. Therefore, various sampling strategies and scoring functions for NA-ligand interactions have been developed. Here, we review the basic principles and current status of docking algorithms and scoring functions for DNA/RNA-ligand interactions. We also discuss challenges and limitations of current docking and scoring approaches.

A new G-triplex-based strategy for sensitivity enhancement of the detection of endonuclease activity and inhibition

EcoRI is an important biomacromolecule in live cells and protects bacterial cells against foreign DNA. In this work, we developed a simple and convenient G-triplex (G3) (5′-TGGGAAGGGAGGGAATTCCCT-3′)-based colorimetric assay for the rapid and selective detection of EcoRI activity and inhibition. The sequence specifically responds to EcoRI in the presence of K⁺ and hemin to form a G-triplex/hemin complex. Taking advantage of G-triplex, EcoRI activity was investigated under the optimized conditions. The absorption intensity ratio displayed a linear relationship against the concentration of EcoRI in the range 0 to 100 U mL⁻¹, and the detection limit was 5.7 U mL⁻¹. Furthermore, G3 showed good selectivity, and the ability to be used to screen for EcoRI inhibitors, indicating its potential in detection and analysis applications.

A new G-triplex-based probe was developed for detecting EcoRI activity and inhibition. The probe showed good selectivity towards EcoRI. The assay was colorimetric and can be monitored by the naked eye.

Insight into formation propensity of pseudocircular DNA G-hairpins

We recently showed that Saccharomyces cerevisiae telomeric DNA can fold into an unprecedented pseudocircular G-hairpin (PGH) structure. However, the formation of PGHs in the context of extended sequences, which is a prerequisite for their function in vivo and their applications in biotechnology, has not been elucidated. Here, we show that despite its ‘circular’ nature, PGHs tolerate single-stranded (ss) protrusions. High-resolution NMR structure of a novel member of PGH family reveals the atomistic details on a junction between ssDNA and PGH unit. Identification of new sequences capable of folding into one of the two forms of PGH helped in defining minimal sequence requirements for their formation. Our time-resolved NMR data indicate a possibility that PGHs fold via a complex kinetic partitioning mechanism and suggests the existence of K⁺ ion-dependent PGH folding intermediates. The data not only provide an explanation of cation-type-dependent formation of PGHs, but also explain the unusually large hysteresis between PGH melting and annealing noted in our previous study. Our findings have important implications for DNA biology and nanotechnology. Overrepresentation of sequences able to form PGHs in the evolutionary-conserved regions of the human genome implies their functionally important biological role(s).

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.

An enzyme-free probe based on G-triplex assisted by silver nanocluster pairs for sensitive detection of microRNA-21

A sensitive ratiometric fluorescence probe based on hybridization chain reaction (HCR) was constructed for sensitive detection of miRNA-21 by using G-triplex and silver nanocluster pairs (AgNC pairs) as an enzyme-free and label-free signal output group. miRNA-21 was used as the primer for the hybridization chain reaction of molecular beacon 1 (MB1) containing the locked G-triplex sequence and molecular beacon 2 (MB2) with intact AgNC pairs at the 5' and 3' end activation. The double-stranded product was obtained along with the opening of the G-triplex and the separation of the AgNC pairs. A detection limit of 67 pM and a linear detection range of 0.1-300 nM were obtained for miRNA-21 determination. The proposed strategy enabled the monitoring of miRNA-21 levels in at least three cell lines, indicating that it provided new ideas for detecting miRNA in real samples. MB1 and MB2 contained the locked G-triplex sequence and silver nanocluster pairs (AgNC pairs), respectively. In the presence of target, the hybridization chain reaction (HCR) between MB1 and MB2 was initiated. At the same time, the locked G-triplex was released and combined to the dye thioflavin T (THT) to increase fluorescence, while the separation of the AgNC pairs caused the fluorescence to decrease. The double-stranded (ds) DNA product was generated to form a ratiometric signal to be detected.

A pH-triggered G-triplex switch with K+ tolerance

A pH-triggered G-triplex (G3) switch is demonstrated to operate in K⁺ using a planar ligand enabling reversible iminium-alkanolamine conversion as the G3 structuring-destructuring initiator.

A comprehensive evaluation of a typical plant telomeric G-quadruplex (G4) DNA reveals the dynamics of G4 formation, rearrangement, and unfolding

Telomeres are specific nucleoprotein structures that are located at the ends of linear eukaryotic chromosomes and play crucial roles in genomic stability. Telomere DNA consists of simple repeats of a short G-rich sequence: TTAGGG in mammals and TTTAGGG in most plants. In recent years, the mammalian telomeric G-rich repeats have been shown to form G-quadruplex (G4) structures, which are crucial for modulating telomere functions. Surprisingly, even though plant telomeres are essential for plant growth, development, and environmental adaptions, only few reports exist on plant telomeric G4 DNA (pTG4). Here, using bulk and single-molecule assays, including CD spectroscopy, and single-molecule FRET approaches, we comprehensively characterized the structure and dynamics of a typical plant telomeric sequence, d[GGG(TTTAGGG)3]. We found that this sequence can fold into mixed G4s in potassium, including parallel and antiparallel structures. We also directly detected intermediate dynamic transitions, including G-hairpin, parallel G-triplex, and antiparallel G-triplex structures. Moreover, we observed that pTG4 is unfolded by the AtRecQ2 helicase but not by AtRecQ3. The results of our work shed light on our understanding about the existence, topological structures, stability, intermediates, unwinding, and functions of pTG4.

Parallel G-triplexes and G-hairpins as potential transitory ensembles in the folding of parallel-stranded DNA G-Quadruplexes

Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.

Scaffold stabilization of a G-triplex and study of its interactions with G-quadruplex targeting ligands

G-triplex nucleic acid structures (G3) have been conjectured to form in vivo but little is known about their physiological functions. The identification of ligands capable of specific binding to G3 structures is therefore highly appealing but remains elusive. Here we report on the assembly of a DNA conjugate which folds into a stable G3 structure. The structural mimic was used to probe the interactions between a G3 ligand and first-in-class G4 ligands, revealing signification binding promiscuity.

Studying the Potassium-Induced G-Quadruplex DNA Folding Process Using Microscale Thermophoresis

Guanine (G) quadruplexes (G4s) can be formed by G-rich sequences when stabilized by the binding of cations (typically K⁺ or Na⁺) and play an essential role in replication, recombination, transcription, and telomere maintenance. Understanding of the G4 folding process is crucial for determining their cellular functions. However, G4-K⁺ interactions and folding pathways are still not well understood. By using human telomeric G4 (hTG4) as an example, two binding states corresponding to two K⁺ cations binding to hTG4 were distinguished clearly and fitted precisely. The basic binding parameters during G4-K⁺ interactions were measured and calculated by taking advantage of microscale thermophoresis (MST), which monitors the changes in charge and size at the same time. The G-hairpin and G-triplex have been suggested as intermediates during G4 folding and unfolding. We further analyzed the equilibrium dissociation constants of 10 possible folding intermediates using MST; thus, the energetically favorable folding/unfolding pathways were proposed. The results might not only shed new light on G4-K⁺ interactions and G4 folding pathways but also provide an example for experimentally studying DNA-ion interactions.

Recent Progress of Targeted G-Quadruplex-Preferred Ligands Toward Cancer Therapy

A G-quadruplex (G4) is a well-known nucleic acid secondary structure comprising guanine-rich sequences, and has profound implications for various pharmacological and biological events, including cancers. Therefore, ligands interacting with G4s have attracted great attention as potential anticancer therapies or in molecular probe applications. To date, a large variety of DNA/RNA G4 ligands have been developed by a number of laboratories. As protein-targeting drugs face similar situations, G-quadruplex-interacting drugs displayed low selectivity to the targeted G-quadruplex structure. This low selectivity could cause unexpected effects that are usually reasons to halt the drug development process. In this review, we address the recent research on synthetic G4 DNA-interacting ligands that allow targeting of selected G4s as an approach toward the discovery of highly effective anticancer drugs.

Involvement of G-triplex and G-hairpin in the multi-pathway folding of human telomeric G-quadruplex

G-quadruplex (G4) can be formed by G-rich DNA sequences that are widely distributed throughout the human genome. Although G-triplex and G-hairpin have been proposed as G4 folding intermediates, their formation still requires further investigation by experiments. Here, we employed single-molecule FRET to characterize the folding dynamics of G4 from human telomeric sequence. First, we observed four states during G4 folding initially assigned to be anti-parallel G4, G-triplex, G-hairpin and unfolded ssDNA. Then we constructed putative intra-strand G-triplex, G-hairpin structures and confirmed their existences in both NaCl and KCl. Further studies revealed those structures are going through dynamic transitions between different states and show relatively weak dependence on cations, unlike G4. Based on those results and molecular dynamics simulations, we proposed a multi-pathway folding mechanism for human telomeric G4. The present work may shed new light on our current understanding about the existence and stability of G4 intermediate states.

Escherichia coli DNA polymerase I can disrupt G-quadruplex structures during DNA replication

Non-canonical four-stranded G-quadruplex (G4) DNA structures can form in G-rich sequences that are widely distributed throughout the genome. The presence of G4 structures can impair DNA replication by hindering the progress of replicative polymerases (Pols), and failure to resolve these structures can lead to genetic instability. In the present study, we combined different approaches to address the question of whether and how Escherichia coli Pol I resolves G4 obstacles during DNA replication and/or repair. We found that E. coli Pol I-catalyzed DNA synthesis could be arrested by G4 structures at low protein concentrations and the degree of inhibition was strongly dependent on the stability of the G4 structures. Interestingly, at high protein concentrations, E. coli Pol I was able to overcome some kinds of G4 obstacles without the involvement of other molecules and could achieve complete replication of G4 DNA. Mechanistic studies suggested that multiple Pol I proteins might be implicated in G4 unfolding, and the disruption of G4 structures requires energy derived from dNTP hydrolysis. The present work not only reveals an unrealized function of E. coli Pol I, but also presents a possible mechanism by which G4 structures can be resolved during DNA replication and/or repair in E. coli.

Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li+, Na+, K+, and Rb+ Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential

Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K⁺ and Na⁺, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li⁺, Na⁺, K⁺, or Rb⁺) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K⁺ or two Rb⁺ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.

Structure of a Stable G-Hairpin

In this study, we report the first atomic resolution structure of a stable G-hairpin formed by a natively occurring DNA sequence. An 11-nt long G-rich DNA oligonucleotide, 5'-d(GTGTGGGTGTG)-3', corresponding to the most abundant sequence motif in irregular telomeric DNA from Saccharomyces cerevisiae (yeast), is demonstrated to adopt a novel type of mixed parallel/antiparallel fold-back DNA structure, which is stabilized by dynamic G:G base pairs that transit between N1-carbonyl symmetric and N1-carbonyl, N7-amino base-pairing arrangements. Although the studied sequence first appears to possess a low capacity for base pairing, it forms a thermodynamically stable structure with a rather complex topology that includes a chain reversal arrangement of the backbone in the center of the continuous G-tract and 3'-to-5' stacking of the terminal residues. The structure reveals previously unknown principles of the folding of G-rich oligonucleotides that could be applied to the prediction of natural and/or the design of artificial recognition DNA elements. The structure also demonstrates that the folding landscapes of short DNA single strands is much more complex than previously assumed.

Folding of guanine quadruplex molecules-funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies

BACKGROUND

Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding.

SCOPE OF REVIEW

We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles.

MAJOR CONCLUSIONS

The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables.

GENERAL SIGNIFICANCE

We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

Titania Composites with 2 D Transition Metal Carbides as Photocatalysts for Hydrogen Production under Visible-Light Irradiation

MXenes, a family of two-dimensional transition-metal carbides, were successfully demonstrated as co-catalysts with rutile TiO2 for visible-light-induced solar hydrogen production from water splitting. The physicochemical properties of Ti3 C2 Tx MXene coupled with TiO2 were investigated by a variety of characterization techniques. The effect of the Ti3 C2 Tx loading on the photocatalytic performance of the TiO2 /Ti3 C2 Tx composites was elucidated. With an optimized Ti3 C2 Tx content of 5 wt %, the TiO2 /Ti3 C2 Tx composite shows a 400 % enhancement in the photocatalytic hydrogen evolution reaction compared with that of pure rutile TiO2 . We also expanded our exploration to other MXenes (Nb2 CTx and Ti2 CTx ) as co-catalysts coupled with TiO2 , and these materials also exhibited enhanced hydrogen production. These results manifest the generality of MXenes as effective co-catalysts for solar hydrogen production.

NALDB: nucleic acid ligand database for small molecules targeting nucleic acid

Nucleic acid ligand database (NALDB) is a unique database that provides detailed information about the experimental data of small molecules that were reported to target several types of nucleic acid structures. NALDB is the first ligand database that contains ligand information for all type of nucleic acid. NALDB contains more than 3500 ligand entries with detailed pharmacokinetic and pharmacodynamic information such as target name, target sequence, ligand 2D/3D structure, SMILES, molecular formula, molecular weight, net-formal charge, AlogP, number of rings, number of hydrogen bond donor and acceptor, potential energy along with their K_i, K_d, IC₅₀ values. All these details at single platform would be helpful for the development and betterment of novel ligands targeting nucleic acids that could serve as a potential target in different diseases including cancers and neurological disorders. With maximum 255 conformers for each ligand entry, our database is a multi-conformer database and can facilitate the virtual screening process. NALDB provides powerful web-based search tools that make database searching efficient and simplified using option for text as well as for structure query. NALDB also provides multi-dimensional advanced search tool which can screen the database molecules on the basis of molecular properties of ligand provided by database users. A 3D structure visualization tool has also been included for 3D structure representation of ligands. NALDB offers an inclusive pharmacological information and the structurally flexible set of small molecules with their three-dimensional conformers that can accelerate the virtual screening and other modeling processes and eventually complement the nucleic acid-based drug discovery research. NALDB can be routinely updated and freely available on bsbe.iiti.ac.in/bsbe/naldb/HOME.php.

Database URL:http://bsbe.iiti.ac.in/bsbe/naldb/HOME.php

G-quadruplexes with (4n - 1) guanines in the G-tetrad core: formation of a G-triad·water complex and implication for small-molecule binding

G-quadruplexes are non-canonical structures of nucleic acids, in which guanine bases form planar G-tetrads (G·G·G·G) that stack on each other in the core of the structure. G-quadruplexes generally contain multiple times of four (4n) guanines in the core. Here, we study the structure of G-quadruplexes with only (4n - 1) guanines in the core. The solution structure of a DNA sequence containing 11 guanines showed the formation of a parallel G-quadruplex involving two G-tetrads and one G-triad with a vacant site. Molecular dynamics simulation established the formation of a stable G-triad·water complex, where water molecules mimic the position of the missing guanine in the vacant site. The concept of forming G-quadruplexes with missing guanines in the core broadens the current definition of G-quadruplex-forming sequences. The potential ability of such structures to bind different metabolites, including guanine, guanosine and GTP, in the vacant site, could have biological implications in regulatory functions. Formation of this unique binding pocket in the G-triad could be used as a specific target in drug design.

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.

Unimolecular antiparallel G-quadruplex folding topology of 2′–5′-isoTBA sequences remains unaltered by loop composition

Unlike 3′–5′-linked TBA, the 2′–5′-linked isoTBA formed only unimolecular antiparallel G-quadruplexes independent of loop length.