G-quadruplex Stability Research Articles

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 ΔGbind 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.

G-quadruplex binding properties of a potent PARP-1 inhibitor derived from 7-azaindole-1-carboxamide

Poly ADP-ribose polymerases (PARP) are key proteins involved in DNA repair, maintenance as well as regulation of programmed cell death. For this reason they are important therapeutic targets for cancer treatment. Recent studies have revealed a close interplay between PARP1 recruitment and G-quadruplex stabilization, showing that PARP enzymes are activated upon treatment with a G4 ligand. In this work the DNA binding properties of a PARP-1 inhibitor derived from 7-azaindole-1-carboxamide, (2-[6-(4-pyrrolidin-1-ylmethyl-phenyl)-pyrrolo[2,3-b]pyridin-1-yl]-acetamide, compound 1) with model duplex and quadruplex DNA oligomers were studied by NMR, CD, fluorescence and molecular modelling. We provide evidence that compound 1 is a strong G-quadruplex binder. In addition we provide molecular details of the interaction of compound 1 with two model G-quadruplex structures: the single repeat of human telomeres, d(TTAGGGT)4, and the c-MYC promoter Pu22 sequence. The formation of defined and strong complexes with G-quadruplex models suggests a dual G4 stabilization/PARP inhibition mechanism of action for compound 1 and provides the molecular bases of its therapeutic potential.

Learning to Model G-Quadruplexes: Current Methods and Perspectives.

G-quadruplexes have raised considerable interest during the past years for the development of therapies against cancer. These noncanonical structures of DNA may be found in telomeres and/or oncogene promoters, and it has been observed that the stabilization of such G-quadruplexes may disturb tumor cell growth. Nevertheless, the mechanisms leading to folding and stabilization of these G-quadruplexes are still not well established, and they are the focus of much current work in this field. In seminal works, stabilization was observed to be produced by cations. However, subsequent studies showed that different kinds of small molecules, from planar and nonplanar organic molecules to square-planar and octahedral metal complexes, may also lead to the stabilization of G-quadruplexes. Thus, the comprehension and rationalization of the interaction of these small molecules with G-quadruplexes are also important topics of current interest in medical applications. To shed light on the questions arising from the literature on the formation of G-quadruplexes, their stabilization, and their interaction with small molecules, synergies between experimental studies and computational works are needed. In this review, we mainly focus on in silico approaches and provide a broad compilation of different leading studies carried out to date by different computational methods. We divide these methods into twomain categories: (a) classical methods, which allow for long-timescale molecular dynamics simulations and the corresponding analysis of dynamical information, and (b) quantum methods (semiempirical, quantum mechanics/molecular mechanics, and density functional theory methods), which allow for the explicit simulation of the electronic structure of the system but, in general, are not capable of being used in long-timescale molecular dynamics simulations and, therefore, give a more static picture of the relevant processes.

A stable uncompleted tetramolecular G-quadruplex formed by d(AGnA) under acidic condition

In this work, the effects of terminal adenines on the formation and stability of tetramolecular G-quadruplexes (G4s) have been studied by electrospray ionization mass spectrometry (ESI-MS), UV, CD and NMR spectroscopy. Several evidences suggested that the sequences d(AGnA) (n = 4 or 5) form stable uncompleted tetramolecular G4 at acidic condition which is different from the canonical one in the neutral condition. In addition, hydrolysis of guanine has also been observed in acidic condition that may occur for unpaired strands rather than in complete G4. Thus, a new G4 topology containing incomplete G-quartet is proposed that is very stable and particularly presents in acidic ammonium ions solution. The information presented in this study provides the new insight on the polymorphism of G4s in acidic environment, which may help understand of the special role of adenines on the formation of G4s.

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.

Single Molecule Investigations of Protein and Small Molecule Interactions with G-Quadruplexes and their Impact on Telomere Structure

G-quadruplex (GQ) structures form in guanine-rich segments of the genome and are particularly more prominent at telomeres and promoters. Telomeric GQ structures (GQs) have emerged as potential targets for small molecule drugs that aim to inhibit telomerase activity in cancer cells, while GQs in promoters are considered to play a role in transcription and translation level gene expression regulation. Due to their high thermal stability, GQs require activity of DNA binding proteins or helicases to be unfolded. Stabilizing small molecules are typically characterized based on their impact on thermal stability of GQs. As a physiologically more relevant approach, we determined the efficacy of several prominent small molecules based on their impact of helicase activity. We measured binding kinetics of individual small molecules and their impact on dynamics and steady state Bloom helicase (BLM) mediated GQ unfolding activities. We also studied Replication Protein A (RPA), GQ, and small molecule interactions while the GQ is maintained under tension by a looped DNA. Finally, we present FRET-PAINT measurements where we studied the accessibility of long telomeric repeats that form multiple tandem GQs and observed length dependent compaction in telomere structure.

Trizbenzim, Cu-Trizbenzim and Zn-Trizbenzim as G-Quadruplex Inducing and Stabilizing Compounds on Human Telomeric Sequence and their Anticancer Properties.

The benzimidazole and their derivatives have rich biological relevance with respect to available natural amino acids and their role in protein folding and quaternary conformations. Thus the ligand trizbenzim and their Cu(II) and Zn(II) metal complexes were prepared to induce G-quadruplex conformation even under no-salt conditions with remarkable anticancer activities. The ligand N,N',N''-Tris-(1H-benzoimidazol-2-ylmethyl)-[1,3,5]triazine-2,4,6-triamine ( trizbenzim) and its Cu and Zn complexes (Cu-trizbenzim, Zn-trizbenzim) were synthesized and characterized by IR, NMR, and MALDI-TOF techniques. The pure ligand and its complexes interacted with human telomere DNA sequence d(TTAGGG), HTelo8and HTelo20and the interactions were followed by circular dichroism spectroscopy, FID assay, and molecular docking techniques. The compounds were tested for anticancer activity towards selected cell lines. All the three compounds stabilized the HTelo8 and HTelo20 in parallel and antiparallel G-quadruplex conformations with salt conditions. Under no-salt conditions, the compounds induce and stabilize the G-quadruplex conformation in antiparallel topology selectively. The pure ligand, Cu-trizbenzim, and Zn-trizbenzim were involved in partial or classical intercalation and some backbone interactions on the strand. The FID assay using thiazole orange intercalator supports the proposed intercalation mode of binding for the three compounds, especially for the pure ligand and the Cu-complex. The MOE docking experiments using X-ray and NMR derived G-quadruplex models with the title compounds extensively support the G-quadruplex induction and stabilization of the telomere sequence by these compounds. The guanines bases involved in the G-tetrad formation interact well with the triazine and the benzimidazole part of the ligand through strong π-π interactions. The primary mode of binding is described as end stacking and intercalation of the compounds to the G-quadruplex structures. The Cu-trizbenzim exhibited more anticancer property in comparison to the pure ligand and the Zntrizbenzim complex. The IC50 values were in the nanomolar range from 50 to 150nM in concentration. This novel self-induction of G-quadruplex is novel without the presence of alkali metal ions.

Stabilization of telomeric G-quadruplex by ligand binding increases susceptibility to S1 nuclease.

The extent of thermodynamic stabilization of telomeric G-quadruplex (G4) by isomers of G4 ligand L2H2-6OTD, a telomestatin analog, is inversely correlated with susceptibility to S1 nuclease. L2H2-6OTD facilitated the S1 nuclease activities through the base flipping in G4, unlike the conventional role of G4 ligands which inhibit the protein binding to DNA/RNA upon ligand interactions.

Interactions of a DNA G-quadruplex with TMAO and urea: a molecular dynamics study on co-solute compensation mechanisms.

We study the individual and combined influence of TMAO and urea on a basket-type DNA G-quadruplex by means of atomistic molecular dynamics (MD) simulations. In combination with the Kirkwood-Buff theory of solutions, we propose a simple mechanism to elucidate the impact of TMAO and urea on the G-quadruplex. Our results reveal the importance of the molecular accumulation around the DNA in terms of stabilizing or destabilizing effects. The results for mixtures show only a weak interaction between both co-solutes, which highlights the additivity of contributions. Despite the fact, that TMAO can to some extent compensate the adverse impact of urea on the G-quadruplex structure, the destabilizing influence is not completely eliminated. This observation opens the door for further research on selective stabilization of DNA G-quadruplexes by modulating the concentrations of TMAO and urea in solution.

Mechanism of recognition of parallel G-quadruplexes by DEAH/RHAU helicase DHX36 explored by molecular dynamics simulations

Because of high stability and slow unfolding rates of G-quadruplexes (G4), cells have evolved specialized helicases that disrupt these non-canonical DNA and RNA structures in an ATP-dependent manner. One example is DHX36, a DEAH-box helicase, which participates in gene expression and replication by recognizing and unwinding parallel G4s. Here, we studied the molecular basis for the high affinity and specificity of DHX36 for parallel-type G4s using all-atom molecular dynamics simulations. By computing binding free energies, we found that the two main G4-interacting subdomains of DHX36, DSM and OB, separately exhibit high G4 affinity but they act cooperatively to recognize two distinctive features of parallel G4s: the exposed planar face of a guanine tetrad and the unique backbone conformation of a continuous guanine tract, respectively. Our results also show that DSM-mediated interactions are the main contributor to the binding free energy and rely on making extensive van der Waals contacts between the GXXXG motifs and hydrophobic residues of DSM and a flat guanine plane. Accordingly, the sterically more accessible 5′-G-tetrad allows for more favorable van der Waals and hydrophobic interactions which leads to the preferential binding of DSM to the 5′-side. In contrast to DSM, OB binds to G4 mostly through polar interactions by flexibly adapting to the 5′-terminal guanine tract to form a number of strong hydrogen bonds with the backbone phosphate groups. We also identified a third DHX36/G4 interaction site formed by the flexible loop missing in the crystal structure.

The combined action of cations and anions of ionic liquids modulates the formation and stability of G-quadruplex DNA.

G-Quadruplex (Gq) formation and stabilization by any molecule is an essential requirement for its application in therapy, especially in oncology. Metal cations have shown higher propensity of the formation of the Gq structure and its stabilization. In this study, the role of both cations and anions of ionic liquids (ILs) on the Gq formation of human telomere (hTeloG) and its stability was investigated using spectroscopic and molecular dynamics simulation techniques. Irrespective of the nature of anions of ILs, tetramethylguanidinium (TMG) cations associated with different anions can form an antiparallel Gq structure in hTeloG. However, the propensity of the formation of an antiparallel Gq structure and its stability depend on the chain length of anions of ILs. Gq is significantly less stable in ILs having longer hydrocarbon chain anions compared to the short chain anions suggesting that the hydrophobicity of the anion plays a critical role in the stability and formation of the Gq structure by ILs. The data indicate that longer hydrocarbon chain anions of ILs preferably interact in the loop region of Gq through hydrophobic interaction which enhances the overall binding of the cation of ILs with Gq causing a decrease in the stacking energy between the G-quartets as well as Hoogsteen hydrogen bonds between the guanine bases leading to the destabilization of the antiparallel Gq structure.

New 2,9-disubstituted-1,10-phenanthroline derivatives with anticancer activity by selective targeting of telomeric G-quadruplex DNA

Fifteen new 1,10-phenanthrolines disubstituted at positions 2 and 9 via amide bonds with different heterocycles have been designed and synthesized as G-quadruplex DNA stabilizers. Ten compounds were evaluated for the in vitro anticancer activity against 60 human tumor cell lines panel, four of them showing a very good inhibitory activity on several cell lines. To assess the ability of the most active compounds to interact with G-quadruplex DNA (G4-DNA), circular dichroism experiments were performed. The potency of the compounds to stabilize the G4-DNA has been shown from the thermal denaturation experiments. The mechanism of compounds binding to DNA and to G4-DNA was theoretically investigated by molecular docking studies.The experimental results demonstrated excellent capacity of the two compounds bearing two pyridin-3-yl residues (methylated and non-methylated) to act as selective G-quadruplex binders with promising anticancer activity.

Ion Binding Properties and Dynamics of the bcl-2 G-Quadruplex Using a Polarizable Force Field.

G-quadruplexes (GQs) are topologically diverse, highly thermostable noncanonical nucleic acid structures that form in guanine-rich sequences in DNA and RNA. GQs are implicated in transcriptional and translational regulation and genome maintenance, and deleterious alterations to their structures contribute to diseases such as cancer. The expression of the B-cell lymphoma 2 (Bcl-2) antiapoptotic protein, for example, is under transcriptional control of a GQ in the promoter of the bcl-2 gene. Modulation of the bcl-2 GQ by small molecules is of interest for chemotherapeutic development but doing so requires knowledge of the factors driving GQ folding and stabilization. To develop a greater understanding of the electrostatic properties of the bcl-2 promoter GQ, we performed molecular dynamics simulations using the Drude-2017 polarizable force field and compared relevant outcomes to the nonpolarizable CHARMM36 force field. Our simulation outcomes highlight the importance of dipole-dipole interactions in the bcl-2 GQ, particularly during the recruitment of a bulk K+ ion to the solvent-exposed face of the tetrad stem. We also predict and characterize an "electronegative pocket" at the tetrad-long loop junction that induces local backbone conformational change and may induce local conformational changes at cellular concentrations of K+. These outcomes suggest that moieties within the bcl-2 GQ can be targeted by small molecules to modulate bcl-2 GQ stability.

Kinetic Target-Guided Synthesis of Small-Molecule G-Quadruplex Stabilizers.

The formation of a G‐quadruplex motif in the promoter region of the c‐MYC protooncogene prevents its expression. Accordingly, G‐quadruplex stabilization by a suitable ligand may be a viable approach for anticancer therapy. In our study, we used the 4‐(4‐methylpiperazin‐1‐yl)aniline molecule, previously identified as a fragment library screen hit, as a template for the SAR‐guided design of a new small library of clickable fragments and subjected them to click reactions, including kinetic target‐guided synthesis in the presence of a G‐quadruplex forming oligonucleotide Pu24. We tested the clickable fragments and products of click reactions for their G‐quadruplex stabilizing activity and determined their mode of binding to the c‐MYC G‐quadruplex by NMR spectroscopy. The enhanced stabilizing potency of click products in biology assays (FRET, Polymerase extension assay) matched the increased yields of in situ click reactions. In conclusion, we identified the newly synthesized click products of bis‐amino derivatives of 4‐(4‐methylpiperazin‐1‐yl)aniline as potent stabilizers of c‐MYC G‐quadruplex, and their further evolution may lead to the development of an efficient tool for cancer treatment.

DNA G-Quadruplex and i-Motif Structure Formation Is Interdependent in Human Cells.

Guanine- and cytosine-rich nucleic acid sequences have the potential to form secondary structures such as G-quadruplexes and i-motifs, respectively. We show that stabilization of G-quadruplexes using small molecules destabilizes the i-motifs, and vice versa, indicating these gene regulatory controllers are interdependent in human cells. This has important implications as these structures are predominately considered as isolated structural targets for therapy, but their interdependency highlights the interplay of both structures as an important gene regulatory switch.

Biophysical and X-ray structural studies of the (GGGTT)3GGG G-quadruplex in complex with N-methyl mesoporphyrin IX.

The G-quadruplex (GQ) is a well-studied non-canonical DNA structure formed by G-rich sequences found at telomeres and gene promoters. Biological studies suggest that GQs may play roles in regulating gene expression, DNA replication, and DNA repair. Small molecule ligands were shown to alter GQ structure and stability and thereby serve as novel therapies, particularly against cancer. In this work, we investigate the interaction of a G-rich sequence, 5'-GGGTTGGGTTGGGTTGGG-3' (T1), with a water-soluble porphyrin, N-methyl mesoporphyrin IX (NMM) via biophysical and X-ray crystallographic studies. UV-vis and fluorescence titrations, as well as a Job plot, revealed a 1:1 binding stoichiometry with an impressively tight binding constant of 30-50 μM-1 and ΔG298 of -10.3 kcal/mol. Eight extended variants of T1 (named T2 -T9) were fully characterized and T7 was identified as a suitable candidate for crystallographic studies. We solved the crystal structures of the T1- and T7-NMM complexes at 2.39 and 2.34 Å resolution, respectively. Both complexes form a 5'-5' dimer of parallel GQs capped by NMM at the 3' G-quartet, supporting the 1:1 binding stoichiometry. Our work provides invaluable details about GQ-ligand binding interactions and informs the design of novel anticancer drugs that selectively recognize specific GQs and modulate their stability for therapeutic purposes.

Lighting Up the Thrombin-Binding Aptamer G-Quadruplex with an Internal Cyanine-Indole-Quinolinium Nucleobase Surrogate. Direct Fluorescent Intensity Readout for Thrombin Binding without Topology Switching.

Fluorescent nucleobases represent an important class of molecular reporters of nucleic acid interactions. In this work, the advantages of utilizing a noncanonical fluorescent nucleobase surrogate for monitoring thrombin binding by the 15-mer thrombin binding aptamer (TBA) is presented. TBA folds into an antiparallel G-quadruplex (GQ) with loop thymidine (T) residues interacting directly with the protein in the thrombin-TBA complex. In the free GQ, T3 is solvent-exposed and does not form canonical base-pairs within the antiparallel GQ motif. Upon thrombin binding, T3 interacts directly with a hydrophobic protein binding pocket. Replacing T3 with a cyanine-indole-quinolinium (4QI) hemicyanine dye tethered to an acyclic 1,2-propanediol linker is shown to have minimal impact on GQ stability and structure with the internal 4QI displaying a 40-fold increase in emission intensity at 586 nm (excitation 508 nm) compared to the free dye in solution. Molecular dynamics (MD) simulations demonstrate that the 4QI label π-stacks with T4 and T13 within the antiparallel GQ fold, which is supported by strong energy transfer (ET) fluorescence from the GQ (donor) to the 4QI label (acceptor). Thrombin binding to 4QI-TBA diminishes π-stacking interactions between 4QI and the GQ structure to cause a turn-off emission intensity response with an apparent dissociation constant (Kd) of 650 nM and a limit of detection (LoD) of 150 nM. These features highlight the utility of internal noncanonical fluorescent surrogates for monitoring protein binding by GQ-folding aptamers in the absence of DNA topology switching.

Human Telomerase (hTERT) is Directly Regulated by the Non-Telomeric Telomere-Binding Factor TRF2

Human telomerase reverse transcriptase hTERT remains suppressed in most normal adult somatic cells. Resulting erosion of telomeres leads eventually to replicative senescence. Reactivation of hTERT maintains telomeres and triggers progression of >90% human cancers. However any direct causal link between telomeres and telomerase regulation remains unclear. Here we show that the telomere-repeat-binding-factor TRF2 binds to the hTERT promoter and recruits the polycomb-repressor EZH2/PRC2 complex. This is causal for H3K27 trimethylation at the hTERT promoter and represses hTERT in cancer as well as normal cells. TRF2 interacts with hTERT promoter G-quadruplexes. Two highly recurrent hTERT promoter mutations found in many cancers, including ~83% glioblastoma-multiforme, that are known to destabilize hTERT promoter G-quadruplexes, obliterated TRF2 binding in patient-derived primary glioblastoma-multiforme cells. Ligand-induced G-quadruplex stabilization restored TRF2 binding, H3K27-trimethylation and hTERT re-suppression. These results uncover a previously unknown mechanism of hTERT regulation through a telomeric factor, implicating telomere-telomerase molecular links that might be important in neoplastic transformation, ageing and regenerative therapy.

Effect of interaction between loop bases and ions on stability of G-quadruplex DNA* *Project supported by the National Natural Science Foundation of China (Grant Nos. 11705160 and 11647074).

G-quadruplexes (GQs) are guanine-rich, non-canonical nucleic acid structures that play fundamental roles in biological processes. The topology of GQs is associated with the sequences and lengths of DNA, the types of linking loops, and the associated metal cations. However, our understanding on the basic physical properties of the formation process and the stability of GQs is rather limited. In this work, we employed ab initio, molecular dynamics (MD), and steered MD (SMD) simulations to study the interaction between loop bases and ions, and the effect on the stability of G-quadruplex DNA, the Drude oscillator model was used in MD and SMD simulations as a computationally efficient manner method for modeling electronic polarization in DNA ion solutions. We observed that the binding energy between DNA bases and ions (K+/Na+) is about the base stacking free energies indicates that there will be a competition among the binding of M+-base, H-bonds between bases, and the base-stacking while ions were bound in loop of GQs. Our SMD simulations indicated that the side loop inclined to form the base stacking while the loop sequence was Thy or Ade, and the cross-link loop upon the G-tetrads was not easy to form the base stacking. The base stacking side loop complex K+ was found to have a good stabilization synergy. Although a stronger interaction was observed to exist between Cyt and K+, such an interaction was unable to promote the stability of the loop with the sequence Cyt.

Structural basis for stabilization of human telomeric G-quadruplex [d-(TTAGGGT)]4 by anticancer drug epirubicin

Anthracycline anticancer drugs show multiple strategies of action on gene functioning by regulation of telomerase enzyme by apoptotic factors, e.g. ceramide level, p53 activity, bcl-2 protein levels, besides inhibiting DNA/RNA synthesis and topoisomerase-II action. We report binding of epirubicin with G-quadruplex (G4) DNA, [d-(TTAGGGT)]4, comprising human telomeric DNA sequence TTAGGG, using 1H and 31P NMR spectroscopy. Diffusion ordered spectroscopy, sequence selective changes in chemical shift (~0.33 ppm) and line broadening in DNA signals suggest formation of a well-defined complex. Presence of sequential nuclear Overhauser enhancements at all base quartet steps and absence of large downfield shifts in 31P resonances preclude intercalative mode of interaction. Restrained molecular dynamics simulations using AMBER force field incorporating intermolecular drug to DNA interproton distances, involving ring D protons of epirubicin depict external binding close to T1-T2-A3 and G6pT7 sites. Binding induced thermal stabilization of G4 DNA (~36 °C), obtained from imino protons and differential scanning calorimetry, is likely to come in the way of telomerase association with telomeres. The findings pave the way for drug-designing with modifications at ring D and daunosamine sugar.