Manish
Debnath
a,
Shirsendu
Ghosh
b,
Deepanjan
Panda
a,
Irene
Bessi
c,
Harald
Schwalbe
c,
Kankan
Bhattacharyya
b and
Jyotirmayee
Dash
*a
aDepartment of Organic Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700032, India. E-mail: ocjd@iacs.res.in
bDepartment of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700032, India
cInstitute of Organic Chemistry and Chemical Biology, Goethe University Frankfurt and Center for Biomolecular Magnetic Resonance (BMRZ), Max-von-Laue Strasse 7, 60438, Frankfurt am Main, Germany
First published on 12th February 2016
The changes in structure and dynamics of oncogenic (c-MYC) and telomeric (h-TELO) G-rich DNA sequences due to the binding of a novel carbazole derivative (BTC) are elucidated using single-molecule Förster resonance energy transfer (sm-FRET), fluorescence correlation spectroscopy (FCS) and NMR spectroscopy. In contrast to the previous reports on the binding of ligands to pre-folded G-quadruplexes, this work illustrates how ligand binding changes the conformational equilibria of both unstructured G-rich DNA sequences and K+-folded G-quadruplexes. The results demonstrate that K+ free c-MYC and h-TELO exist as unfolded and partially folded conformations. The binding of BTC shifts the equilibria of both investigated DNA sequences towards the folded G-quadruplex structure, increases the diffusion coefficients and induces faster end-to-end contact formation. BTC recognizes a minor conformation of the c-MYC quadruplex and the two-tetrad basket conformations of the h-TELO quadruplex.
Single-molecule Förster resonance energy transfer (sm-FRET) is a powerful method that provides key information about the structure, population distribution of folded or unfolded species and the end-to-end distance of bio-molecules.19–35 It has been used to investigate the effect of binding of either protein26 or metal ions (K+/Na+)27–29 on the conformational dynamics of G-quadruplexes. However, only a few studies have been conducted to elucidate the influence of small molecules on the conformational dynamics of G-quadruplex DNA.30,31
Fluorescence correlation spectroscopy (FCS) tracks single or several molecules in solution to investigate changes in diffusion coefficients, molecular size and intra-molecular contact dynamics.32–38 FCS has been used to investigate the intra-molecular dynamics of a DNA hairpin tagged with a donor–acceptor FRET pair.36–38 A combination of sm-FRET and FCS has been used by Majima and co-workers to quantitatively analyse the pH-induced intra-molecular folding dynamics of an i-motif DNA.39 We herein illustrate the changes in structure and dynamics in unstructured and K+-folded c-MYC and h-TELO DNA G-quadruplex forming sequences triggered by a G-quadruplex binding ligand bis-triazolylcarbazole40 (BTC) using a combination of sm-FRET and FCS. Interaction of BTC with G-quadruplexes has also been substantiated by FRET melting, circular dichroism, fluorescence lifetime and NMR spectroscopy studies.
Following our recently developed procedure, BTC was prepared using a one-pot Cu(I) catalyzed azide and alkyne cycloaddition in high yields40 (Fig. 1b, and S1, ESI†). FRET melting revealed that BTC exhibits maximum stabilization potentials for K+-folded c-MYC-(A) and h-TELO-(A) (i.e. pre-annealed in K+ buffer) with ΔTm values of 24.0 °C and 38.7 °C, respectively. It is important to note that BTC shows a similar increase in the Tm values of unfolded c-MYC-(A) and h-TELO-(A) sequences (Fig. 1c and Table S1, ESI†). These results indicate that BTC can stabilize c-MYC-(A) and h-TELO-(A) G-quadruplexes even in the absence of K+ ions (Fig. 1d).40 Furthermore, BTC is found to be selective for G-quadruplexes over duplex DNA as it did not significantly alter the Tm value of duplex (ds) DNA (Table S1, ESI†).
Dual labeled sequences of the highest purity are used to exclude the signals from the donor only sample and no appreciable donor bleaching is observed under current experimental conditions (Fig. S2 and S3, ESI†). We observed that the donor–acceptor fluorescence intensities of the dual-labeled G-rich sequences produce anti-correlated fluctuations in the presence and absence of BTC and K+ ions (Fig. 2 and S4, ESI†). The histograms obtained from the time traces were fitted with tri- and single Gaussian distributions (Fig. 2b and S5, ESI†). The contributions of shot noise in each FRET peak are determined (Table S2, ESI†). The analysis of the FRET histogram of K+-free c-MYC-(A) shows two major peaks with FRET efficiencies (εFRET) centered at ∼0.4 (30%) and 0.6 (68%) (Table S3, ESI†). The distances between donor and acceptor (RDA) dyes corresponding to each FRET state (0.4 and 0.6) calculated using eqn S3 (ESI†), are 58.4 Å and 51.7 Å, respectively. The FRET histogram of the K+-folded c-MYC-(A) shows only one major peak at ∼0.8 (97%). Upon addition of BTC (1 equiv.), the single peak at ∼0.8 is preserved in the histogram (Fig. 2b). The average RDA determined for the K+-folded c-MYC-(A) in the presence and absence of BTC lies between ∼44–42 Å. Interestingly, the FRET distribution of the K+-free c-MYC-(A) in the presence of BTC exhibits a narrow peak (Fig. 2b) with higher efficiency (εFRET ∼ 0.85) corresponding to a RDA value of ∼41 Å.
Fig. 2 Photon bursts of donor/acceptor (background corrected) (a), and FRET efficiency distributions (b) of 100 pM dual fluorescently labeled c-MYC-(A) G-quadruplex forming sequence in the presence and absence of K+ and BTC. For detailed information see Fig. S4 and S5, ESI.† |
The FRET histogram of K+-free h-TELO-(A) shows a wide distribution (mean εFRET ∼ 0.51) with a RDA of ∼54.6 Å (Fig. S4, ESI†). The K+-folded h-TELO-(A) shows a εFRET value of ∼0.8 with a RDA of ∼45 Å. Upon addition of BTC to either K+-free or the K+-folded h-TELO-(A), the εFRET distribution is shifted towards a higher value (∼0.9). The corresponding RDA values were determined to be ∼40.6 Å (BTC and K+) and ∼39 Å (only BTC). Considering the presence of extra nucleotides and the dye linkers present in the dual labeled sequences, the RDA values obtained for c-MYC-(A) and h-TELO-(A) in the presence of K+ and BTC (∼39–42 Å), are close to the size of a folded G-quadruplex [diameter ∼ 25 Å] conformation.
In unfolded structures, the distances between the donor and acceptor fluorophores (RDA ∼ 54–58 Å) are large and, consequently, the εFRET values are small. The low FRET peak observed for the K+-free c-MYC-(A) sequence at εFRET ∼0.4 may be assigned to the unfolded (single stranded) structure. The FRET state ∼0.6 indicates the existence of an intermediate state, presumably an ensemble of partially folded structures. The high FRET peak (∼0.8–0.9) indicates the formation of G-quadruplex structures.
BTC-induced quadruplex formation is further corroborated by carrying out sm-FRET experiments using a dual labeled mutant c-MYC sequence (Fig. S6 and S7, ESI†). Upon addition of BTC, the FRET histogram of c-MYC-mut doesn't show any significant change in FRET efficiency. These results suggest that G-rich sequences primarily remain as unfolded and partially folded structures in the absence of BTC and K+. BTC apparently induces and shifts the equilibrium to a folded G-quadruplex conformation as evidenced by the shift of FRET histograms towards higher values.
(1) |
The diffusion coefficients (Dt) were calculated from the diffusion time (τD) using equation S8, ESI† (Tables 1 and S4, ESI†). The Dt of K+-free c-MYC-(A) is 221 μm2 s−1 and increases to 296 μm2 s−1 in the K+-folded conformation. Similarly, the Dt of the K+-free h-TELO-(A) increases in the K+-folded conformations (Tables 1 and S4, ESI†). Addition of BTC to the K+-folded c-MYC-(A) and h-TELO-(A) results in further increase in the Dt values to 314 μm2 s−1 and 318 μm2 s−1, respectively. It is intriguing to observe that the binding of BTC to the K+-free c-MYC-(A) and h-TELO-(A) increases the Dt values to 318 μm2 s−1 and that is ∼43% and ∼32% higher compared to the K+-free c-MYC-(A) and h-TELO-(A) sequences. According to the Stokes–Einstein eqn (S10, ESI†), the diffusion coefficient is inversely related to the hydrodynamic radii (Rh) for a freely diffusing molecule. Therefore, the observed increase in Dt of the K+-free c-MYC-(A) and h-TELO-(A) sequences can be attributed to a ∼43% and ∼32% decrease in Rh upon the binding of BTC. The decrease in Rh estimated theoretically from the ratio of frictional coefficients (f/f0) indicates that the binding of BTC to the K+-free c-MYC-(A) and h-TELO-(A) leads to a ∼50% reduction in the Rh values. In addition, the diffusion parameters calculated using eqn S19, ESI† are quite similar to those obtained from the simplified eqn (1) (Table S5, ESI†). To further validate our findings, we carried out parallel diffusion experiments with single (FAM) labeled c-MYC-(B) and h-TELO-(B) G-quadruplex forming sequences (Fig. S9, ESI†). The τD values for single labeled c-MYC-(B) and h-TELO-(B) in the presence and absence of BTC and K+ ions are similar to those of the dual labeled c-MYC-(A) and h-TELO-(A) sequences. Hence the covalent attachment of the dyes does not significantly impact the diffusion parameters (Fig. S9 and Table S6, ESI†).
System | D t (μm2 s−1) | τ Cobs (μs) | k + (ms−1) | k − (ms−1) |
---|---|---|---|---|
c-MYC-(A) | 221 | 49 | 10.6 ± 1.1 | 9.8 ± 1 |
c-MYC-(A) + K+ | 296 | 9 | 37.0 ± 4 | 74.0 ± 7 |
c-MYC-(A) + K+ + BTC | 318 | 5 | 97.4 ± 10 | 102.6 ± 11 |
c-MYC-(A) + BTC | 318 | 6 | 88.0 ± 9 | 78.6 ± 8 |
h-TELO-(A) | 242 | 29 | 10.8 ± 1.1 | 23.7 ± 2.5 |
h-TELO-(A) + K+ | 283 | 9 | 41.7 ± 4.3 | 69.4 ± 7 |
h-TELO-(A) + K+ + BTC | 314 | 7 | 71.4 ± 7.4 | 71.5 ± 7 |
h-TELO-(A) + BTC | 318 | 7 | 66.8 ± 7 | 76.1 ± 8 |
D t, τCobs = ±10% |
Next, the kinetic parameters for the end-to-end contact formation (k+) and dissociation (k−) of the K+-free and K+-folded quadruplexes are calculated (Fig. 3b and d, Tables 1 and S4, ESI†). The end-to-end contact formation of K+-free DNA corresponds to the initial intra-chain contact formation during the G-quadruplex formation (Fig. 3e).42 In contrast, the k+/k− values of the folded quadruplexes indicate the motion of flanking sequences and the dye linkers (Fig. 3f).39 The K+-free c-MYC-(A) and h-TELO-(A) exhibit nearly identical k+ values of ∼10 ms−1. The k+ values are increased ∼3–4 fold for the K+-folded c-MYC-(A) and h-TELO-(A) quadruplexes. However, the K+-free c-MYC-(A) and h-TELO-(A) sequences show k− values of 9.8 ms−1 and 23.7 ms−1, respectively. The K+-folded c-MYC-(A) and h-TELO-(A) quadruplexes exhibit a 7-fold and a 3-fold increase in k− values, respectively. Upon binding to the BTC, c-MYC-(A) displays an ∼8–10 fold increase and the h-TELO-(A) displays a ∼6.0 and a ∼3.0 fold increase in the k+ and k− values, respectively. Interestingly, the k+/k− values obtained for single labeled c-MYC-(B) and h-TELO-(B) are also comparable to those of c-MYC-(A) and h-TELO-(A) (Fig. S9, Table S6, ESI†).
We presumed that ligand-induced changes in diffusion parameters are associated with the conformation of the DNA. The K+-free DNA sequences have a large surface area, which causes increased hydration on the surface of the DNA. This creates a constraint in flexibility in K+-free DNA and consequently, the rates of intra-chain contact formation and dissociation (k+/k−) have lower values (Fig. 3e). However, the BTC-mediated folding of K+-free DNA sequences into compact globular quadruplexes causes dehydration due to lower availability of the hydration sites resulting in high k+ and k− values.
(2) |
(3) |
The c-MYC-(A) exhibits a kFRET value of ∼0.8 ns−1 that increases to ∼1.3 ns−1 and ∼2.0 ns−1 in the presence of K+ and BTC, respectively. A similar increase in the kFRET value is observed for h-TELO-(A) after the addition of K+ and BTC, respectively. For further illustration, the RDA of c-MYC-(A) and h-TELO-(A) are estimated from the average FRET efficiency (εavg) using eqn S21, ESI† (Tables S7 and S8, ESI†). The calculated RDA is ∼49 Å for the c-MYC-(A), ∼44 Å for K+-folded c-MYC-(A), ∼42 Å for the BTC bound K+-folded c-MYC-(A) and ∼39 Å for the BTC folded c-MYC-(A) conformations. For h-TELO-(A), the RDA value decreases from ∼55 Å in the K+-free-state to ∼44 Å in K+-folded state and it reduces further to 41 Å in the BTC-bound conformation. These distances are in close agreement with the sm-FRET data and demonstrate that the binding of BTC decreases the distance between the donor and acceptor labels, as compared to the K+-free and K+-folded c-MYC-(A) and h-TELO-(A) conformations. Together, these results indicate that BTC can template the formation of compact quadruplexes from the K+-free G-rich sequences that facilitate faster donor decay via donor-to-acceptor energy transfer.
Fig. 4 NMR titration of c-MYC-(D) and h-TELO-(D) in the presence and absence of BTC and K+. (a) Scheme of the major conformation of c-MYC-(D), as determined using NMR from Ambrus et al.43 [PDB code: 1XAV], with sequence and numbering. Imino region of the 1D 1H NMR spectrum of c-MYC-(D) in the presence of an increasing amount of BTC, with (b) or without (c) 100 mM KCl. (d) Scheme of the proposed structure of h-TELO-(D), with sequence and numbering. Imino region of the 1D 1H NMR spectrum of h-TELO-(D) in the presence of an increasing amount of BTC, with (e) or without (f) 100 mM KCl. Experimental conditions: 100 μM DNA, 25 mM Tris–HCl buffer (pH 7.4), 298 K, 600 MHz. Guanine residues in anti and syn conformations are represented in cyan and magenta, respectively. |
Addition of BTC to c-MYC-(D) shows that the ligand selectively interacts and stabilizes a minor conformation, both in the presence and absence of K+ (Fig. 4c and d). The signals of the imino protons of the major conformation of c-MYC-(D) are not perturbed by BTC. However, a set of signals, already detectable in the absence of BTC (minor conformation), becomes more intense upon the addition of ligand, which suggests that the ligand can interact with the c-MYC-(D) quadruplex via conformational selection. As suggested by the CD data (Fig. S14, ESI†), we propose that the minor conformation stabilized by BTC preserves the all-parallel scaffold of the major conformation (Fig. 4a), but presents different arrangements of the capping structures.
In K+ containing buffer, the 1H 1D imino pattern revealed that h-TELO-(D) exists predominantly as a 2-tetrad basket quadruplex, probably with looser capping structures compared to the structure reported by Lim et al.17 (Fig. 4d and e). Other minor conformations are also observed in the 1H NMR spectra of the h-TELO-(D) quadruplex. Addition of BTC to the h-TELO pre-folded in the presence of K+ induces weak chemical shift perturbations (highlighted with grey shadows, Fig. 4e) of the imino signals from the minor conformation and, possibly, from the guanine residues involved in the capping structure (G9 and G21). In the absence of K+, the h-TELO-(D) is found to be partially folded, as suggested by the imino signals detectable in the Hoogsteen hydrogen bond region (Fig. 4f).
Interestingly, upon addition of BTC, at a ratio [Ligand]:[DNA] = 2, a broad signal at 13.5 ppm, indicating the formation of a Watson/Crick (WC) base pair, and a sharp signal at 12.5 ppm, possibly arising from the carbazole NH group or from the newly formed Hoogsteen imino bonds between guanines, were detected. These findings suggest that in the absence of K+, BTC is able to promote the restructuring of the capping structures via formation of intra- or inter-loop WC-A:T interactions.
In contrast with that observed by sm-FRET (Fig. S4, ESI†), the aromatic region of the 1H NMR spectra of h-TELO-(D) (Fig. S14c and d, ESI†) reveals that in the absence of K+, BTC is able to promote G-quadruplex folding only partially. Taken together, the NMR and CD data support the single-molecule results, confirming that the ligand is able to stabilize the c-MYC-(D) and h-TELO-(D) G-quadruplexes. However, the details of the interaction process, particularly in the absence of K+, might be different to the one reported by single-molecule studies, due to the different experimental conditions. In particular, the DNA concentration used in the NMR studies is 106-fold higher than the concentration used in the sm-FRET studies.
Footnote |
† Electronic supplementary information (ESI) available: Experimental details, synthetic procedures, characterization data of compounds, 1H NMR and 13C NMR spectra, sm-FRET, shot noise, FCS-FRET data, lifetime data, CD spectra, NMR spectra. See DOI: 10.1039/c6sc00057f |
This journal is © The Royal Society of Chemistry 2016 |