Discovery of a structural-element specific G-quadruplex "light-up" probe.

The development of a fluorescent probe capable of detecting and distinguishing the wide diversity of G-quadruplex structures is particularly challenging. Herein, we report a novel BODIPY-based fluorescent sensor (GQR) that shows unprecedented selectivity to parallel-stranded G-quadruplexes with exposed ends and four medium grooves. Mechanistic studies suggest that GQR associates with G-quadruplex grooves close to the end of the tetrad core, which may explain the dye's specificity to only a subset of parallel structures. This specific recognition favours the disaggregation of GQR in aqueous solutions thereby recovering the inherent fluorescence of the dye. Due to its unique features, GQR represents a valuable tool for basic biological research and the rapid discovery of novel, specific ligands that target similar structural features of G-quadruplexes.

Nucleic acid sequences rich in guanines (G) have a propensity to arrange into G-quadruplexes (G4)1. These non-canonical structures have been implicated in transcription regulation and telomere maintenance; hence, they are targeted for cancer treatment and detection23. Conversely, engineered G4 oligonucleotides can act as drug delivery vehicles456, materials for nanodevices78, support catalysts9 or even drugs for cancer, HIV and other diseases10. For these reasons, there is strong impetus to develop probes to improve our understanding of G4 and their ligand binding characteristics11.

Studies suggest that G4 structures are highly polymorphic and exhibit distinct differences between each other1213. These myriad of topologies arises from the combination of several well-defined structural elements such as nucleic acid type (DNA or RNA), molecularity (monomer, dimer or tetramer), strand orientation (parallel, antiparallel or hybrid), loops (orientation, sequence and length), grooves type (narrow, wide and/or medium), end-cap morphology and number of G-tetrads114. Efficient probing of these motifs provides opportunities for discrimination and is also a prerequisite to the discovery and study of new G4 specific ligands1516. We were thus motivated to develop fluorescent probes that are specific for these fundamental elements.

Several fluorescent probes capable of discerning G4 from duplex DNA have been developed. These emissive probes are generally planar structures which achieve selectivity through end-stacking with G-tetrads171819. In contrast, structural specific G4 probes are considerably rarer. Moreover, to achieve higher specificity in distinguishing different G4s, simultaneous association with multiple structural elements is desired19. Reported probes are conceived from ligands which have been found to interact with the grooves of G4 a priori202122. Therefore, such rational design strategies can have limited efficiency for discovering new G4 probes with novel selectivity21. Conversely, employing the diversity-oriented fluorescent library approach (DOFLA) is a more effective strategy232425; screening of diversity oriented fluorescent libraries (DOFL) can accelerate the identification of fluorescent probes with the desired qualities without former insight on the structural guidelines needed.

Results

Sensor discovery

A collection of 5000 potential fluorescent sensors were thus gathered and an unbiased, high throughput screening was performed to uncover primary hits responding to G4 (Fig. S1). Additional modifications to improve the quantum yield of the hit after binding led to the discovery of a sensitive and highly selective fluorescent sensor for G4 (GQR, Fig. 1). In the presence of 93del – an interlocked, dimeric, parallel-stranded G4 with 4 medium grooves26 – GQR displayed up to 30-fold increase in fluorescence at 597 nm and a 12 nm bathochromic shift in emission maximum (Fig. 2). A Job plot analysis of the GQR-93del complex revealed a binding stoichiometry of 1:1 (Fig. S2); the dissociation constant (Kd) was thus determined as 25.18 ± 0.02 μM (Fig. 2, inset and Table S1).

Figure 1

The structures of GQR and its analogues.

Figure 2

GQR reacts with 93del to give a turn-on fluorescence response.

Fluorescence spectra of GQR (10 μM) upon incubation with serial dilutions of 93del (from 0–80 μM) in buffer (20 mM K2HPO4/KH2PO4, 100 mM KCl, pH 7.0). λex: 360 nm. ΦF (without 93del) = 0.014, ΦF (in 40 μM 93del) = 0.28. Inset: photographic image of GQR (10 μM) in the presence of 93del (from 0–40 μM) and data plot of fluorescent emission intensity upon addition of 93del (from 0–80 μM). F0 and Fmax are the fluorescent maximum intensities of the GQR in the absence and presence of 93del respectively. Kd = 25.18 ± 0.02 μM (one site specific binding model). Values are represented as means (n = 3). Measurements were taken at room temperature (RT).

Selectivity of GQR

To investigate the uniqueness of GQR, we examined its selectivity towards various G4 oligonucleotides with different structural elements (Table 1). Like 93del, J19, T95, T95 and T95-2T form parallel-stranded G4 with exposed ends and four medium grooves (Fig. S3)26272829. When incubated with GQR, all four sequences elicited a fluorescence increase from the dye (Fig. 3). Specifically, the enhancement was most pronounced with 93del. In contrast, both c-kit1 and Pu24T are parallel-stranded but they possess additional snapback motifs that cap the G-tetrad core3031. Interestingly, both sequences returned negligible fluorescence when mixed with GQR (Fig. 3 and Fig. S4). On the contrary, Oxy and HT adopt non-parallel-stranded conformations (Fig. S3)3233; GQR similarly remained quenched when incubated with either G4 (Fig. 3 and Fig. S4). Likewise, GQR remained virtually non-fluorescent in the presence of other conventional nucleic acids (Fig. 3 and Fig. S4). Previous reports of parallel-stranded selective G4 sensors generally do not display selectivity for additional motifs within the parallel-stranded structures202234. To the best of our knowledge, GQR is the first fluorescent dye that is able to discern among various types of parallel-stranded G4.

Table 1

Oligonucleotide sequences used in this work

Name	Sequence (5′ → 3′)	Structural Elements	Origin
93del	GGGG TGGG AGGA GGGT	parallel, interlocked dimeric	Aptamer, HIV-1 integrase
J19	GIGT GGGT GGGT GGGT	parallel, dimeric	Aptamer, HIV-1 integrase
T95	GGGT GGGT GGGT GGGT	parallel, dimeric	Aptamer, HIV-1 integrase
T95-2T	TT GGGT GGGT GGGT GGGT	parallel, monomeric	Aptamer, HIV-1 integrase
c-kit1	AGGG AGGG CGCT GGGA GGAG GG	snap-back parallel	c-kit oncogene promoter
Pu24T	TGAG GGTG GTGA GGGT GGGG AAGG	snap-back parallel	c-myc oncogene promoter
Oxy	GGGG TTTT GGGG	antiparallel	Oxytricha nova telomere
HT	TT GGGTTA GGGTTA GGGTTA GGGA	hybrid	Human telomere
dsDNA	Sigma-D8515	genomic dsDNA	Calf thymus
ssDNA	Sigma-D8899	genomic ssDNA	Calf thymus
RNA	Sigma-R7250	single-stranded	Calf liver

Figure 3

Selectivity of GQR to oligonucleotides.

(a) Fluorescence titration of GQR with various G-quadruplexes, DNA and RNA. Conditions: GQR (10 μM), nucleic acid (from 0 to 40 μM), buffer (20 mM K2HPO4/KH2PO4, 100 mM KCl, pH 7.0), RT. λex: 360 nm, λem: 600 nm. (b) Photographic image of GQR (10 μM) mixed with various nucleic acids (40 μM). Irradiation with a hand-held UV lamp at 365 nm.

Mechanism studies

Disaggregation-induced emission of GQR

Following experiments were aimed at understanding the sensing mechanism of GQR. First, we examined the photophysical properties of GQR. GQR exhibited strong fluorescence emission in organic solvents, but only negligible emission in buffer (Table 2, entries 13–16 vs. entry 1). Further measurements of the absorption spectra showed significant peak broadening and a bathochromic shifts in buffer (Fig. 4a). With higher concentrations of GQR, more pronounced red-shifts were observed suggesting that this phenomenon may be a result of dye aggregation (Fig. S5). Transmission electron microscope and dynamic light scattering analysis revealed the existence of GQR-aggregates with sizes related to the dye concentration (Fig. S6 and S7). Together, these results indicate that the low fluorescence of GQR in buffer is a result of aggregation-caused quenching (ACQ)35 thereby affirming the potential of GQR to behave as fluorescent turn-on sensors for G4223637. To confirm this potential, the absorption spectra of GQR with 93del were studied. Indeed, when GQR was mixed with increasing concentrations of the G4, the spectral shape gradually shifted towards that similar to organic solvents (disaggregated state) (Fig. 4b).

Table 2

Fluorescence quantum yields (ΦF) of respective probes in various environments

Entry	Probe	Solvent	Analytea	ΦFb
1	GQR	Bufferc	-	0.014
2	GQR	Buffer	93del	0.28
3	GQR	Buffer	J19	0.22
4	GQR	Buffer	T95	0.19
5	GQR	Buffer	T95-2T	0.14
6	GQR	Buffer	c-kit1	0.009
7	GQR	Buffer	Pu24T	0.005
8	GQR	Buffer	Oxy	0.027
9	GQR	Buffer	HT	0.007
10	GQR	Buffer	dsDNA	0.011
11	GQR	Buffer	ssDNA	0.008
12	GQR	Buffer	RNA	0.016
13	GQR	DMSO	-	0.55
14	GQR	MeOH	-	0.57
15	GQR	PEG-400	-	0.46
16	GQR	Toluene	-	0.30
17	2	DMSO	-	0.30
18	3	DMSO	-	0.36
19	4	DMSO	-	0.39
20	5	DMSO	-	0.21
21	2	Buffer	-	0.013
22	3	Buffer	-	0.030
23	4	Buffer	-	0.33
24	5	Buffer	-	0.24
25	2	Buffer	93del	0.20
26	3	Buffer	93del	0.21
27	4	Buffer	93del	0.34
28	5	Buffer	93del	0.23

aAnalyte concentration: 40 μM.

bValues are represented as means (n = 3). Measurements were taken at room temperature (RT).

cBuffer: 20 mM KH2PO4/K2HPO4, 100 mM KCl, pH 7.0.

Figure 4

Absorption spectra of GQR.

(a) GQR (10 μM) in various organic solvents. (b) GQR (10 μM) upon incubation with serial dilutions of 93del (from 0–40 μM) in buffer. Values are represented as means (n = 3). Measurements were taken at RT.

Environmental sensitivity of GQR

To investigate the contribution of environment polarity to the emission of GQR, we examined its emission spectra in various toluene-methanol mixtures. We observed that the emission maximum of GQR red-shifted from 582 to 590 nm as solvent polarity decreases from methanol to toluene (Fig. S8). These observations suggests that the bathochromic shift in the emission of GQR when incubated with 93del is a consequence of the G4 providing a hydrophobic binding pocket for GQR thereby affirming the interaction between the two entities.

Binding mode studies

We were further interested in the interaction of GQR with its respective G4s at the molecular level. To this end, four additional analogues of GQR were studied in greater detail (Fig. 1). Like GQR, compounds 2, 3, 4 and 5 shared the same absorbance spectra with GQR in DMSO (Fig. S9) and displayed high quantum yields (Table 2, entries 17–20). However, in buffer, only 2 and 3 showed aggregation-induced absorption spectra (Fig. S9) and quenched quantum yields (Table 2, entries 21–24) analogous to GQR. Similarly, when incubated with 93del, only 2 and 3 similarly exhibited a distinct enhancement in fluorescence albeit to a lesser extent than GQR (Table 2, entries 25–28 vs. 6). These results validate the importance of disaggregation and the carboxamide group of GQR in achieving effective fluorescence probing of G4.

Molecular docking was subsequently performed to understand the binding mode3839. As depicted in Fig. 5a, binding simulation suggests that GQR associates with the groove region of the 93del (PDB ID: 1Y8D) core close to the ends. In contrast, a simulation of an end-stacking mode revealed a dramatically weaker affinity and higher energy compared to a groove binding mode (Table S2). This result is consistent with the general lack of peak shifts in the 1H imino proton NMR spectra of 93del and T95-2T upon addition of GQR – an observation that contrasts that of end-stacking ligands (Fig. S10)151718. A closer examination of the binding model suggests plausible hydrogen bonds with the carboxamide group which may account for the enhanced fluorescence of GQR compared to 2 and 3 (Fig. 5b).

Figure 5

Molecular docking results of GQR binding to the 93del.

Illustration of groove-binding mode of GQR with 93del residues. Carbon atoms are colored in white for GQR, pink and green 93del; oxygen and nitrogen atoms are colored in red and blue respectively. (b) Plausible hydrogen bonds between GQR and 93del residues G10 and A12 (black dashed lines). Carbon atoms are colored in white for GQR, pink and green for 93del; oxygen and nitrogen atoms are colored in red and blue respectively.

Discussion

Common BODIPY dyes have a tendency to aggregate which results in self-quenching4041. This ACQ effect often limits the label-to-analyte ratio and narrows the practical applications42. However, BODIPY can recover its high fluorescence in response to a target molecule through a recognition-induced disassembly of the aggregates. In this study, GQR-aggregates exhibit remarkable fluorescence enhancement in the presence 93del. However, water-soluble analogues show no such response under the same conditions. The discovery of GQR as an effective fluorescent sensor demonstrates that the conventional ACQ problem can, in reality, be a general mechanism for probe development4344.

The structural polymorphism of G-quadruplexes has been supported by NMR and X-ray crystallography. Nonetheless, a simple method to distinguish different G-quadruplex structures conveniently and sensitively is highly desirable. The use of specific fluorescent probes, such as GQR, is one such strategy. Fluorescent probes can be potentially employed for mapping and tracking G-quadruplexes-based structure motifs1745 and can also be used for rapid discovery of novel, specific anticancer drugs that recognize similar G-quadruplex structural elements to the probe.

In summary, we have described the systematic discovery of a novel fluorescent dye – GQR – which specifically “light-up” when bound to parallel G4 with end exposed medium grooves. Specific recognition and favourable binding to G4 induces the disaggregation of GQR which results in recovery of the inherent fluorescence of the dye. Preliminary mechanistic studies suggest a groove-binding mode close to the end of the quadruplex core, thereby accounting for the high specificity. Therefore, GQR represents a promising tool for the future discovery of similar G4 groove binding ligands that can be important for disease studies, G4-based drug deliveries and nanotechnology.

Methods

Probe synthesis

Detailed description of the synthesis of each probe can be found in the Supplementary Information. Each step was characterized by high-resolution mass spectra, 1H and 13C NMR.

Nucleic acids

DNA oligonucleotides were synthesized and structurally verified by NMR spectroscopy as previously reported2627282930313233. dsDNA, ssDNA and RNA were purchased from Sigma Aldrich (product number D8515, D8899 and D7250 respectively) and used as such. The nucleic acids were dissolved in buffer (20 mM K2HPO4/KH2PO4 100 mM KCl, pH 7.0). DNA concentration was expressed in strand molarity using a nearest-neighbor approximation for the absorption coefficients of the unfolded species.

Absorbance

UV/Vis absorption spectra of dyes and nucleic acids in buffer (20 mM K2HPO4/KH2PO4, 100 mM KCl, pH 7.0) were recorded from 200 to 700 nm using a SpectraMax M2 spectrophotometer.

Fluorescence

Fluorescence measurements were carried out on a SpectraMax M2 spectrophotometer in 96-well plates by scanning the emission spectra between 540 and 700 nm (λex = 360 nm). All experiments were repeated three times. Data analysis was performed using Origin 8.0 (OriginLab Corporation, MA).

Quantum yield measurements

Quantum yields were calculated by measuring the integrated emission area of the fluorescent spectra in its respective solvents and comparing to the area measured for Coumarin 1 (reference) (10 μM), ΦF = 0.73 in ethanol (η = 1.36), λex = 360 nm. Quantum yields were calculated using the equation: where F represents the area of fluorescent emission, η is the refractive index of the solvent, and Abs is absorbance at the excitation wavelength. Emission was integrated between 540 and 700 nm.

Dissociation constant measurements

The Kd of GQR to the respective G-quadruplexes were analyzed by Origin 8.0 (OriginLab Corporation, MA) using the following equation for a one site specific binding model: Where y represents the fluorescence fold change of GQR, ymax the fluorescence fold change of GQR when saturated with G-quadruplex and x the concentration of the G-quadruplex in μM.

Transmission electron microscope

GQR (10 μM) was first prepared in buffer and deposited on a thin copper-support film, followed by drying in vacuo. Images of the samples were obtained with JEOL JEM 3010 HRTEM microscope and operated at 100 kV without any contrast agent.

Dynamic light scattering

The dynamic light scattering of different concentration GQR and other compounds were measured at 25°C in buffer using quartz cell. All measurements were performed in triplicate in Zetasizer Nano ZS.

Nuclear magnetic resonance spectroscopy

Titration with GQR was performed on a 600 MHz NMR Bruker spectrometers equipped with a cryoprobe at 25°C. Oligonucleotides (200 μM) were dissolved in buffer (20 mM K2HPO4/KH2PO4, 70 mM KCl, D2O/H2O (1:9), pH 7.0). The spectra were recorded immediately after each addition of GQR. Water suppression was achieved using excitation sculpting.

Molecular modeling

The coordinates of 93del structures were retrieved from the Protein Data Bank (ID code 1Y8D). DNA structures were prepared for docking. The GQR structure was optimized using the Gaussian03 program (B3LYP/6-31G* level). By using Autodock 4.0, docking studies were carried out with the Lamarckian genetic algorithm following the procedure developed for G-quadruplex DNA and ligand docking3839. Two rounds of simulation were performed. In the first round, simulated annealing was used to find a rough binding mode of GQR with 50 runs while keeping all other parameters default. The search space was subsequently reduced and another 200 runs were conducted to get a more precise result. Following the docking studies of GQR with 93del quadruplex DNA, all the figures were rendered using PyMOL v0.99.

Author Contributions

Y.T.C. and A.T.P. conceived the idea and directed the work. L.Z. and J.C.E. designed the experiments and co-wrote the manuscript. L.Z. designed, screened and identified the G-quadruplex sensor and analyzed the photophysical properties and binding mechanism. J.C.E. and J.Y. designed, synthesized and characterized all dye compounds used in this study. K.K.G. designed and synthesized the diversity-oriented fluorescent libraries. W.J.C. synthesized and characterized all G-quadruplex sequences. W.X. performed the DLS analysis. W.Z. performed the molecular docking studies. All authors contributed to data analysis and manuscript writing.