Interactions of GW-2974 and SCH-442416 with human telemoric G-quadruplex DNA were evaluated using UV-absorption, fluorescence emission, fluorescence quenching and circular dichroism spectroscopies. Binding parameters that include binding constants, binding stoichiometries, binding modes, melting temperatures and selectivity towards G-quadruplex over dsDNA and ct-DNA were investigated. The above data were supported with molecular docking.
3.1. Solubility of GW-2974 and SCH-442416
GW-2974 and SCH-442416 are insoluble in water and soluble in DMSO and ethylene glycol. Solutions of GW-2974 and SCH-442416 in ethylene glycol (100%) were found stable over more than four months at 4 °C. Stabilities of GW-2974 and SCH-442416 in 10% ethylene glycol - Tris-KCl buffer were tested using UV-Vis absorption spectrophotometry. Solutions showed no significant changes in absorption over four hours. This time was considered enough to run all experiments and we adopted 10% ethylene glycol in all our measurements.
Effect of ethylene glycol on structural conformation of telomeric G-quadruplex DNA was also tested (Figure 1). Additions of ethylene glycol to G-quadruplex (up to 20%.) in Tris-KCl buffer showed no change in the shape and intensity of G-quadruplex’ CD spectrum. These results indicated that 20% ethylene glycol is safe and do not change the conformation or denaturize telomeric G-quadruplex DNA. The results are consistent with Bonner and Klibanov's report stated that up to 90% ethylene glycol is safe and non-denaturizing for DNA [75]. Therefore, freshly prepared - 10% ethylene glycol buffer solutions were used throughout this study to improve the solubility of GW-2974 and SCH-442416.
3.2. Spectrophotometric Titration
UV-Vis Absorption
Figure 2 shows the UV-Vis absorption spectra of GW-2974 and SCH-442416 in Tris-KCl buffer pH 7.4. GW-2974 spectrum shows strong bands at 243 nm (ε = 6100) and 320 nm (ε = 5500) and a weak band at 421 nm (ε = 2600) (Fig. 2a). The spectrum of SCH-442416 shows a strong absorption band at 267 nm (ε = 7700) (Fig. 2b).
Changing the pH of GW-2974 and SCH-442416 solutions (10-5 M) indicated no significant spectral changes in the pH range 4.0–10.0. Consequently, the two compounds were considered chemically stable within this pH range and can be safely used for further investigations on their interactions with G-quadruplex DNA at the physiological pH of 7.4.
Figure 2a also shows the effect of successive additions of human telomere G-quadruplex DNA (1.18x10-4 M) to GW-2974 (5x10−5 M) in Tris-KCl buffer pH 7.4. The absorption bands at 421 and 320 nm continuously decreased (hypochromicity) without forming isosbestic points or being shifted. On contrary, the intensity of absorption band at 243 nm increased due to the overlap with the DNA absorption band at 260 nm. Figure 2b shows that additions of human telomeric G-quadruplex DNA to SCH-442416 decreased the intensity of the absorption band at 350 nm and increased the band at 267 nm. The latter can be attributed to interference with the DNA's absorption band at 260 nm (Fig. 2b).
The above hypochromic effects (≥ 50%) obtained upon addition of DNA to GW-2974 or SCH-442416 suggested intercalation binding modes for both compounds with telomere G-quadruplex DNA.
Hypochromicity along with the binding stichometry and binding constants were used to suggest an intercalation binding via stacking on G-quartet and outside binding mode [63]. Further evidences will be obtained using fluorescence and CD spectroscopic measurements.
Fluorescence
A confirmation for binding interactions of GW-2974 and SCH-442416 to human telomeric G-quadruplex DNA were obtained by fluorescence titrations. Figure 3 shows the fluorescence spectra of GW-2974 and SCH-442416. GW-2974 gave a fluorescence band at 410 nm when exited at 243 nm while SCH-442416 gave a broad band centered at 420 nm when excited at 267 nm. The latter band seems to be composed of two overlapping bands at 408 and 423 nm.
Sequential additions of G-quadruplex DNA to GW-2974 resulted in a continuous decrease in the fluorescence intensity at 410 nm. A 65.0% hypochromicity was observed after adding 20 folds of G-quadruplex DNA (Fig. 3a). The band at 484 nm is attributed to Rayleigh scattering, normally shown at double the excitation wavelength. Similarly, additions of G-quadruplex DNA to SCH-442416 continuous decreased the fluorescence emission at 420 nm (Fig. 3b). At 20 folds of added DNA, the band was slightly red shifted and reduced its intensity by 91.4% (hypochromic effect).
Reductions in fluorescence emissions of GW-2974 and SCH-442416 by ≥ 65 – 91. %, indicated that both compounds bound to human telomeric G-quadruplex DNA through intercalation binding modes [63].
Fluorescence Quenching
Additional evidence on intercalation of GW-2974 and SCH-442416 with telomere G-quadruplex DNA was obtained from fluorescence quenching measurements. Fluorescein labeled telomeric G-quadruplex DNA (5’-Fl-AGGG(TTAGGG)3-3') gave a strong fluorescence band at 518 nm when excited at 494 nm. Sequential additions of GW-2974 or SCH-442416 (10-4 M) to human telomeric Fl-G-quadruplex DNA (2x10-6 M) quenched its fluorescence emission at 518 nm by more than 90% (Fig. 4). The extent of quenching depends on the number of molecules bound per DNA molecule.
Decreases of fluorescence emission at 518 nm confirmed intercalations of both compounds in proximity of the fluorescein flag moiety. Since fluorescein molecule is connected to the 5’ prime, one may infer that both GW-2974 and SCH-442416 bind between the G-quartets adjacent to the TTA loop cavity.
These results are consistent with previous findings from UV-Vis absorption and fluorescence measurements and confirm intercalation binding modes between each of GW-2974 and SCH-442416 with G-quadruplex DNA.
Circular Dichroism Titration
CD spectroscopy is the tool of choice for studying conformational changes in DNA upon interactions with drugs. Structural conformation of human telomeric G-quadruplex DNA in solution depends on its sequence, length, and environment. It forms antiparallel structure in Na+ solution whereas in K+ solution, a hybrid mixed with parallel / antiparallel structure is formed [76–77].
Figure 5 shows the CD spectrum of the twenty-two bases - human telomere G-quadruplex DNA in Tris-KCl buffer pH 7.4. The negative band centered at 235 nm, the positive shoulder at 253 nm and the positive band at 293 nm confirmed the formation of a hybrid parallel-antiparallel structure in K+ solution [67–68].
Figure 6 shows the effect of GW-2974 and SCH-442416 additions to telomeric G-quadruplex DNA. Gradual decreases in the intensity of the bands at 235, 253 and 293 nm with no change in the spectral shapes or bands' positions indicated no change in the hybrid parallel antiparallel conformation of the telomeric G-quadruplex DNA at the end of titration. Changes in CD spectra of G-quadruplex DNA during titration have been correlated with the binding mode of the ligand involved. A continuous decrease in CD intensity was correlated with intercalation-binding mode while increase was correlated with groove-binding mode [78].
Thus, decreases in CD intensities upon sequential additions of GW-2974 or SCH-442416 could be correlated with intercalation binding modes by π-π stacking between the G-quartets' faces of G-quadruplex (Fig. 6). These results are also consistent with previous results obtained from absorption and fluorescence titrations.
3.3. Binding Stoichiometry
Stoichiometry of GW-2974 and SCH-442416 interactions with human telomere G-quadruplex DNA was estimated using the molar ratio method based on fluorescence measurements. Figure 7 shows the molar ratio curves obtained by adding successive increments of telomeric G-quadruplex (1.33x10-4 M) to 2.0 ml of 5x10-6 M of GW-2974 or SCH-442416. Plots of fluorescence intensity at 410 and 420 nm versus molar ratio gave the molar stoichiometric ratios of GW-2974 and SCH-442416 to G-quadruplex DNA (Figs. 7a and 7b). A 2:1 molar ratio was obtained.These ratios indicated that two molecules of GW-2974 or SCH-442416 bound per each G-quadruplex molecule. This ratio is consistent with qurecitin to human telomeric G-quadruplex DNA molar ratio.
3.4. Binding Affinity
Several modes have been proposed for binding G-quadruplexe DNA to ligands. These included intercalative binding between adjacent G-quartets, face π-π stacking on external G-quartets or with the loops through and a weaker external binding
Decreases in the intensity of the GW-2974 and SCH-442416 fluorescence soret bands upon additions of G-quadruplex favored their intercalations by face π-π stacking in the ratio 2:1 (GW-2974 or SCH-442416 to G-quadruplex DNA) (Fig. 7).
Binding constants of the two compounds to G-quadruplex DNA were estimated using Scatchard plots. Figure 8 shows non-linear downward concave curves of r/Cf versus r. The plots suggested that more than one type of dependent binding sites exist on G-quadruplex molecule. Dependent sites can synergistically or antagonistically affect each other (neighbor exclusion effect). In the former case, the first bound ligand encourages the next binding ligand while in the latter the first bound ligand suppress the next binding [79–80].
Resolving the nonlinear Scatchard plots in Fig. 7 gave two intersecting straight lines each represents one type of binding site. The slopes of these lines gave binding constants of 1.3x108 and 1.72x106 M-1 for GW-2974 and 1.55x107 and 3.74x105 M-1 for SCH-442416. The number of binding sites ranged between 0.7-1.27 for GW-2974 and 0.4-1.5 for SCH-442416. (Figs. 8a and 8c, Table 1). These results apparently show a type of dependent binding sites.
Table 1
Binding constants (K) and number of binding sites (n) per G-quadruplex DNA molecule for GW-2974 and SCH-442416.
G-quadruplex complex
|
K1 (M-1)
|
K2 (M-1)
|
n1
|
n2
|
Graphical Method
|
GW2974-G-quadruplex
|
1.3x108
|
1.72x106
|
0.7
|
1.27
|
SCH442416-G-quadruplex
|
1.55x107
|
3.74x105
|
0.4
|
1.5
|
Nonlinear Curve Fitting
|
GW2974-G-quadruplex
|
2.65x106
|
--
|
1.15
|
--
|
SCH442416-G-quadruplex
|
7.0x105
|
--
|
1.88
|
--
|
Applying the Scatchard modified equation (5), A non-linear r versus Cf curves were obtained (Figs. 8b and 8d). Nonlinear fitting of the two curves resulted in an average binding constant of 2.65x106 and 7.0x105 M-1 with number of binding sites of 1.15 and 1.88 for GW-2974 and SCH-442416, respectively (Table 1).
These results reflected high stabilization effects for both compounds on human telomeric G-quadruplex DNA. However, GW-2974 showed higher binding affinities than SCH-442416. The number of binding sites are also consistent with the molar ratio analysis and published values (Section 3.3).
3.5. Melting Temperature
Figure 9 shows the melting temperature curves for G-quadruplex, ct-DNA and their complexes with GW-2974 and SCH-442416, based on CD measurements. Melting temperature (Tm) was calculated as the midway CD signal on the melting curve and represents the temperature at which DNA is half unfolded to its single strands.
The melting curve of human telomere G-quadruplex DNA gave a Tm value of 65.5°C while its complexes with GW-2974 and SCH-442416 gave Tm values of 75.4 and 75.1°C (Figs. 9a and 9b) (Table 2). These results revealed stabilization effect for G-quadruplex DNA by ΔTm = 9.9 and 9.6°C respectively. These results indicated a higher stabilizing effect for GW2974 compared to SCH-442416. The results are also consistent with the binding constants shown in Table 1. On the other hand, Fig. 9c shows a melting temperature of 88.0 oC for ct-DNA. Complexations with GW-2974 and SCH-442416 stabilized the ct-DNA by ΔTm of 2.1 and 2.4°C, respectively.
Table 2
Melting temperature of human telomere G-quardruplex DNA and its GW-2974 and SCH-442416 complexes at different [drug]/[DNA] ratios.
Ligand
|
[Ligand]/[G-quadruplex] molar ratio
|
1:1
|
3:1
|
5:1
|
Tm °C
|
ΔTm °C
|
Tm °C
|
ΔTm °C
|
Tm °C
|
ΔTm °C
|
G-quadruplex DNA
|
65.5
|
0.0
|
65.5
|
0.0
|
65.5
|
0.0
|
GW-2974-G-quadruplex
|
75.4
|
9.9
|
69.7
|
4.2
|
69.0
|
3.5
|
SCH-442416-G-quadruplex
|
75.1
|
9.6
|
69.3
|
3.8
|
68.2
|
2.7
|
ct-DNA
|
88.0
|
0.0
|
--
|
--
|
--
|
--
|
GW-2974-ct-DNA
|
90.1
|
2.1
|
--
|
--
|
--
|
--
|
SCH-442416-ct-DNA
|
90.4
|
2.4
|
--
|
--
|
--
|
--
|
These results indicated that GW-2974 and SCH-442416 have ≥ 4 four folds stronger stabilization effects on human telomeric G-quadruplex compared to duplex ct-DNA (4.71 times for GW2974 and 4.00 times for SCH442416). The results also indicated that both compounds have preferential and selective binding to G-quadruplex over duplex ct-DNA.
Figure 9 also shows the effect of GW-2974 and SCH-442416 concentrations on stabilizing human telomeric G-quadruplex DNA. Complexes with molar ratios of 1:1, 3:1 and 5:1 (ligand:G-quadruplex) were measured for their melting curves. Increasing the molar ratio of GW-2974 from 1:1 to 5:1 resulted in decreasing ΔTm values from 9.9 to 3.5 °C. Similar effect was observed for SCH-442416 where ΔTm decreased from 9.6 to 2.7 °C upon changing the ratio from 1:1 to 5:1 (Table 2). These results demonstrated that increasing the ligand’s concentrations above 1:1 molar ratio has destabilized the human telomeric G-quadruplex DNA.
These results suggest that the two investigated compounds act as stabilizers for G-quadruplex DNA. They also teach us that they have higher affinity and selectivity towards human telomeric G-quadruplex DNA over duplex ct-DNA [81].
3.6. Selectivity of binding towards G-quadruplex DNA
Binding selectivity of GW-2974 and SCH-442416 towards telomeric G-quadruplex DNA over duplex DNA was estimated using telomeric dsDNA and ct-DNA as interfering species. A 10.0, 50.0 and 100.0 folds of duplex DNA were added to complexes formed between Fl-G-quadruplex (10-10 M) and either GW-2974 or SCH-442416. The mixtures were scanned for their fluorescence intensity at 518 nm as described in section 2.4.5.
Figure 10 shows changes in fluorescence intensity of Fl-G-quadruplex DNA-GW-2794 upon adding different folds of telomere dsDNA (10a and 10c) and ct-DNA (10b and 10d). Selectivity coefficient was calculated by dividing the change in fluorescence intensity by the intensity of Fl-G-quadruplex complex. Coefficients having values of (1.1-2.2)x10-2 and (2.2-4.4)x10-2 were obtained. These results indicated preferential binding selectivity for binding GW-2974 and SCH-442416 to bind with human telomere G-quadruplex DNA over duplex dsDNA and ct-DNA.
Thus, the two compounds may selectively bind towards G-quadruplex DNA over dsDNA predominant in human cells. These results are also consistent with the results obtained from melting temperature experiments (section 3.5).
3.7. Molecular Docking
Several crystal structures for human telomere G-quadruplex DNA were obtained from the protein data bank and tried in this work. The crystal structure PDB:2MS6 was found the most suitable since its co-crystallized ligand; quercetin has almost a similar size to that of GW-2974 and SCH-442416 molecules. The presence of two co-crystalized ligand molecules bound to a human telomeric G-quadruplex DNA indicated the presence of two binding sites on the crystal structure of DNA. Therefore, our compounds were docked into these two pockets as per the conditions given in the experimental part.
Table 3 shows the docking scores of GW-2974 and SCH-442416 along with the co-crystallized ligand quercetin. The two compounds gave favorable exothermic binding energies ranged from -6.07 to -8.14 kcal/mol. GW-2974 has the best docking score (-8.14 kcal/mol) which is comparable to quercetin (-9.03 kcal/mol) when docked into the same pocket (Site 1, Table 1). On site 2, GW-2974 and SCH-442416 gave comparable binding scores of -6.17 and -6.37 kcal/mol which are inferior to the score obtained by the co-crystallized ligand: quercetin (-9.27 kcal/mol).
Table 3
Docking scores of SCH-442416 and GW-2974 along with the co-crystallized ligand into the binding site of Telomeric G-quadruplex DNA.
Ligand
|
Docking score (kcal/mol)
|
Site 1
|
Site 2
|
GW2974
|
-8.14
|
-6.17
|
SCH-442416
|
-6.07
|
-6.37
|
The co-crystallized ligand; quercetin
|
-9.03
|
-9.72
|
Figure 11 shows that GW-2974 was able to slide in-between the nitrogen bases surrounding site 1 and make hydrogen bonding with DG20, associated with multiple stacking interactions to DG6 and DG27. The Fig. 11 also shows that SCH-442416 was able to slide in-between the nitrogen bases comprising the binding site 2, making the important interaction; hydrogen bonding with DT16.
To sum up, our docking results are -in line with the experimental data- suggest that GW-2974 and SCH-442416 are favorably and exothermally bind to human telomeric G-quadruplex DNA. The two compounds conveniently bind on the quercetin’s binding sites with comparable binding scores.