High-affinity binding at quadruplex-duplex junctions: rather the rule than the exception

Strategy for modeling higher-order G-quadruplex structures recalcitrant to NMR determination

Guanine-rich nucleic acids can form intramolecularly folded four-stranded structures known as G-quadruplexes (G4s). Traditionally, G4 research has focused on short, highly modified DNA or RNA sequences that form well-defined homogeneous compact structures. However, the existence of longer sequences with multiple G4 repeats, from proto-oncogene promoters to telomeres, suggests the potential for more complex higher-order structures with multiple G4 units that might offer selective drug-targeting sites for therapeutic development. These larger structures present significant challenges for structural characterization by traditional high-resolution methods like multi-dimensional NMR and X-ray crystallography due to their molecular complexity. To address this current challenge, we have developed an integrated structural biology (ISB) platform, combining experimental and computational methods to determine self-consistent molecular models of higher-order G4s (xG4s). Here we outline our ISB method using two recent examples from our lab, an extended c-Myc promoter and long human telomere G4 repeats, that highlights the utility and generality of our approach to characterizing biologically relevant xG4s.

Quadruplex-duplex junction in LTR-III: A molecular insight into the complexes with BMH-21, namitecan and doxorubicin

Quadruplex-Duplex (Q-D) junctions are unique structural motifs garnering increasing interest as drug targets, due to their frequent occurrence in genomic sequences. The viral HIV LTR-III sequence was chosen as a Q-D junction model to study the affinity of the selected compounds BMH-21, namitecan (ST-1968), and doxorubicin (DOXO), all containing a planar polycyclic aromatic moiety, linked to either one short aminoalkyl or an aminoglycosyl group. A multidisciplinary approach that combines NMR spectroscopy, molecular modelling, circular dichroism (CD) and fluorescence spectroscopy was employed. The studied ligands induced moderate but clear stabilization to the Q-D junction by interacting with the interfacial tetrad. DOXO was found to be the best Q-D junction binder. Interestingly, the removal of the aminoglycosyl group significantly changed the pattern of the interactions, indicating that highly polar substituents have a stronger affinity with the exposed regions of the Q-D junction, particularly at the level of the interfacial tetrad.

Native Characterization of Noncanonical Nucleic Acid Thermodynamics via Programmable Dynamic DNA Chemistry

Folding thermodynamics, quantitatively described using parameters such as ΔGfold°, ΔHfold°, and ΔSfold°, is essential for characterizing the stability and functionality of noncanonical nucleic acid structures but remains difficult to measure at the molecular level. Leveraging the programmability of dynamic deoxyribonucleic acid (DNA) chemistry, we introduce a DNA-based molecular tool capable of performing a free energy shift assay (FESA) that directly characterizes the thermodynamics of noncanonical DNA structures in their native environments. FESA operates by the rational design of a reference DNA probe that is energetically equivalent to a target noncanonical nucleic acid structure in a series of toehold-exchange reactions, yet is structurally incapable of folding. As a result, a free energy shift (ΔΔGrxn°) is observed when plotting the reaction yield against the free energy of each toehold-exchange. We mathematically demonstrated that ΔGfold°, ΔHfold°, and ΔSfold° of the analyte can be calculated based on ΔΔGrxn°. After validating FESA using six DNA hairpins by comparing the measured ΔGfold°, ΔHfold°, and ΔSfold° values against predictions made by NUPACK software, we adapted FESA to characterize noncanonical nucleic acid structures, encompassing DNA triplexes, G-quadruplexes, and aptamers. This adaptation enabled the successful characterization of the folding thermodynamics for these complex structures under various experimental conditions. The successful development of FESA marks a paradigm shift and a technical advancement in characterizing the thermodynamics of noncanonical DNA structures through molecular tools. It also opens new avenues for probing fundamental chemical and biophysical questions through the lens of molecular engineering and dynamic DNA chemistry.

Solution-State Structure of a Long-Loop G-Quadruplex Formed Within Promoters of Plasmodium falciparum B var Genes

We report the high-resolution NMR solution-state structure of an intramolecular G-quadruplex with a diagonal loop of ten nucleotides. The G-quadruplex is formed by a 34-nt DNA sequence, d[CAG3T2A2G3TATA2CT3AG4T2AG3T2], named UpsB-Q-1. This sequence is found within promoters of the var genes of Plasmodium falciparum, which play a key role in malaria pathogenesis and evasion of the immune system. The [3+1]-hybrid G-quadruplex formed under physiologically relevant conditions exhibits a unique equilibrium between two structures, both stabilized by base stacking and non-canonical hydrogen bonding. Unique equilibrium of the two closely related 3D structures originates from a North-South repuckering of deoxyribose moiety of residue T27 in the lateral loop. Besides the 12 guanines involved in three G-quartets, most residues in loop regions are involved in interactions at both G-quartet-loop interfaces.

Structural Differences at Quadruplex-Duplex Interfaces Enable Ligand-Induced Topological Transitions

Quadruplex-duplex (QD) junctions, which represent unique structural motifs of both biological and technological significance, have been shown to constitute high-affinity binding sites for various ligands. A QD hybrid construct based on a human telomeric sequence, which harbors a duplex stem-loop in place of a short lateral loop, is structurally characterized by NMR. It folds into two major species with a (3+1) hybrid and a chair-type (2+2) antiparallel quadruplex domain coexisting in a K+ buffer solution. The antiparallel species is stabilized by an unusual capping structure involving a thymine and protonated adenine base AH+ of the lateral loop facing the hairpin duplex to form a T·AH+·G·C quartet with the interfacial G·C base pair at neutral pH. Addition and binding of Phen-DC3 to the QD hybrid mixture by its partial intercalation at corresponding QD junctions leads to a topological transition with exclusive formation of the (3+1) hybrid fold. In agreement with the available experimental data, such an unprecedented discrimination of QD junctions by a ligand can be rationalized following an induced fit mechanism.

Structural elucidation of HIV-1 G-quadruplexes in a cellular environment and their ligand binding using responsive 19F-labeled nucleoside probes

Understanding the structure and recognition of highly conserved regulatory segments of the integrated viral DNA genome that forms unique topologies can greatly aid in devising novel therapeutic strategies to counter chronic infections. In this study, we configured a probe system using highly environment-sensitive nucleoside analogs, 5-fluoro-2'-deoxyuridine (FdU) and 5-fluorobenzofuran-2'-deoxyuridine (FBFdU), to investigate the structural polymorphism of HIV-1 long terminal repeat (LTR) G-quadruplexes (GQs) by fluorescence and 19F NMR. FdU and FBFdU, serving as hairpin and GQ sensors, produced distinct spectral signatures for different GQ topologies adopted by LTR G-rich oligonucleotides. Importantly, systematic 19F NMR analysis in Xenopus laevis oocytes gave unprecedented information on the structure adopted by the LTR G-rich region in the cellular environment. The results indicate that it forms a unique GQ-hairpin hybrid architecture, a potent hotspot for selective targeting. Furthermore, structural models generated using MD simulations provided insights on how the probe system senses different GQs. Using the responsiveness of the probes and Taq DNA polymerase stop assay, we monitored GQ- and hairpin-specific ligand interactions and their synergistic inhibitory effect on the replication process. Our findings suggest that targeting GQ and hairpin motifs simultaneously using bimodal ligands could be a new strategy to selectively block the viral replication.

A pH-Responsive Topological Switch Based on a DNA Quadruplex-Duplex Hybrid

A guanine-rich oligonucleotide based on a human telomeric sequence but with the first three-nucleotide intervening stretch replaced by a putative 15-nucleotide hairpin-forming sequence shows a pH-dependent folding into different quadruplex-duplex hybrids in a potassium containing buffer. At slightly acidic pH, the quadruplex domain adopts a chair-type conformation. Upon increasing the pH, a transition with a midpoint close to neutral pH to a major and minor (3+1) hybrid topology with either a coaxially stacked or orthogonally oriented duplex stem-loop occurs. NMR-derived high-resolution structures reveal that an adenine protonation is prerequisite for the formation of a non-canonical base quartet, capping the outer G-tetrad at the quadruplex-duplex interface and stabilizing the antiparallel chair conformation in an acidic environment. Being directly associated with interactions at the quadruplex-duplex interface, this unique pH-dependent topological transition is fully reversible. Coupled with a conformation-sensitive optical readout demonstrated as a proof of concept using the fluorescent dye thiazole orange, the present quadruplex-duplex hybrid architecture represents a potentially valuable pH-sensing system responsive in a physiological pH range of 7±1.

Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junction

G-quadruplexes (GQs) play key regulatory roles within the human genome and have also been identified to play similar roles in other eukaryotes, bacteria, archaea, and viruses. Human immunodeficiency virus 1, the etiological agent of acquired immunodeficiency syndrome, can form two GQs in its long terminal repeat (LTR) promoter region, each of which act to regulate viral gene expression in opposing manners. The major LTR GQ, called LTR-III, is a distinct hybrid GQ containing a 12-nucleotide duplex loop attached to the quadruplex motif. The resulting quadruplex:duplex junction (QDJ) has been hypothesized to serve as a selective drug targeting site. To better understand the dynamics of this QDJ, we performed conventional and enhanced-sampling molecular dynamics simulations using the Drude-2017 force field. We observed unbiased and reversible formation of additional base pairs in the QDJ, between Ade4:Thy14 and Gua3:Thy14. Both base pairs were electrostatically favored, but geometric constraints within the junction may drive the formation of, and preference for, the Ade4:Thy14 base pair. Finally, we demonstrated that the base pairs are separated only by small energy barriers that may enable transitions between both base-paired states. Together, these simulations provide new insights into the dynamics, electrostatics, and thermodynamics of the LTR-III QDJ.

Steric hindrance and structural flexibility shape the functional properties of a guanine-rich oligonucleotide

Ligand/protein molecular recognition involves a dynamic process, whereby both partners require a degree of structural plasticity to regulate the binding/unbinding event. Here, we present the characterization of the interaction between a highly dynamic G-rich oligonucleotide, M08s-1, and its target protein, human α-thrombin. M08s-1 is the most active anticoagulant aptamer selected thus far. Circular dichroism and gel electrophoresis analyses indicate that both intramolecular and intermolecular G-quadruplex structures are populated in solution. The presence of thrombin stabilises the antiparallel intramolecular chair-like G-quadruplex conformation, that provides by far the main contribution to the biological activity of the aptamer. The crystal structure of the thrombin-oligonucleotide complex reveals that M08s-1 adopts a kinked structural organization formed by a G-quadruplex domain and a long duplex module, linked by a stretch of five purine bases. The quadruplex motif hooks the exosite I region of thrombin and the duplex region is folded towards the surface of the protein. This structural feature, which has never been observed in other anti-exosite I aptamers with a shorter duplex motif, hinders the approach of a protein substrate to the active site region and may well explain the significant increase in the anticoagulant activity of M08s-1 compared to the other anti-exosite I aptamers.

Probing juxtaposed G-quadruplex and hairpin motifs using a responsive nucleoside probe: a unique scaffold for chemotherapy

Paucity of efficient probes and small molecule ligands that can distinguish different G-quadruplex (GQ) topologies poses challenges not only in understanding their basic structure but also in targeting an individual GQ form from others. Alternatively, G-rich sequences that harbour unique chimeric structural motifs (e.g., GQ-duplex or GQ-hairpin junctions) are perceived as new therapeutic hotspots. In this context, the epidermal growth factor receptor (EGFR) gene, implicated in many cancers, contains a 30 nucleotide G-rich segment in the promoter region, which adopts in vitro two unique architectures each composed of a GQ topology (parallel and hybrid-type) juxtaposed with a hairpin domain. Here, we report the use of a novel dual-app probe, C5-trifluoromethyl benzofuran-modified 2'-deoxyuridine (TFBF-dU), in the systematic analysis of EGFR GQs and their interaction with small molecules by fluorescence and 19F NMR techniques. Notably, distinct fluorescence and 19F NMR signals exhibited by the probe enabled the quantification of the relative population of random, parallel and hybrid-type GQ structures under different conditions, which could not be obtained by conventional CD and 1H NMR techniques. Using the fluorescence component, we quantified ligand binding properties of GQs, whereas the 19F label enabled the assessment of ligand-induced changes in GQ dynamics. Studies also revealed that mutations in the hairpin domain affected GQ formation and stability, which was further functionally verified in polymerase stop assay. We anticipate that these findings and useful properties of the nucleoside probe could be utilized in designing and evaluating binders that jointly target both GQ and hairpin domains for enhanced selectivity and druggability.

Taq-Polymerase Stop Assay to Determine Target Selectivity of G4 Ligands in Native Promoter Sequences of MYC, TERT, and KIT Oncogenes

Computational and high-throughput experimental methods predict thousands of potential quadruplex sequences (PQSs) in the human genome. Often these PQSs contain more than four G-runs, which introduce additional uncertainty into the conformational polymorphism of the G4 DNA. G4-specific ligands, which are currently being actively developed as potential anticancer agents or tools for studying G4 structures in genomes, may preferentially bind to specific G4 structures over the others that can be potentially formed in the extended G-rich genomic region. We propose a simple technique that identifies the sequences that tend to form G4 in the presence of potassium ions or a specific ligand. Thermostable DNA Taq-polymerase stop assay can detect the preferential position of the G4 -ligand binging within a long PQS-rich genomic DNA fragment. This technique was tested for four G4 binders PDS, PhenDC3, Braco-19, and TMPyP4 at three promoter sequences of MYC, KIT, and TERT that contain several PQSs each. We demonstrate that the intensity of polymerase pausing reveals the preferential binding of a ligand to particular G4 structures within the promoter. However, the strength of the polymerase stop at a specific site does not always correlate with the ligand-induced thermodynamic stabilization of the corresponding G4 structure.

The Potent G-Quadruplex-Binding Compound QN-302 Downregulates S100P Gene Expression in Cells and in an In Vivo Model of Pancreatic Cancer

The naphthalene diimide compound QN-302, designed to bind to G-quadruplex DNA sequences within the promoter regions of cancer-related genes, has high anti-proliferative activity in pancreatic cancer cell lines and anti-tumor activity in several experimental models for the disease. We show here that QN-302 also causes downregulation of the expression of the S100P gene and the S100P protein in cells and in vivo. This protein is well established as being involved in key proliferation and motility pathways in several human cancers and has been identified as a potential biomarker in pancreatic cancer. The S100P gene contains 60 putative quadruplex-forming sequences, one of which is in the promoter region, 48 nucleotides upstream from the transcription start site. We report biophysical and molecular modeling studies showing that this sequence forms a highly stable G-quadruplex in vitro, which is further stabilized by QN-302. We also report transcriptome analyses showing that S100P expression is highly upregulated in tissues from human pancreatic cancer tumors, compared to normal pancreas material. The extent of upregulation is dependent on the degree of differentiation of tumor cells, with the most poorly differentiated, from more advanced disease, having the highest level of S100P expression. The experimental drug QN-302 is currently in pre-IND development (as of Q1 2023), and its ability to downregulate S100P protein expression supports a role for this protein as a marker of therapeutic response in pancreatic cancer. These results are also consistent with the hypothesis that the S100P promoter G-quadruplex is a potential therapeutic target in pancreatic cancer at the transcriptional level for QN-302.