Jakub Żubertowskia,
Magdalena Rapp*a,
Jan Dolicherb,
Błażej Rubiś
b and
Anna Dembska
*a
aFaculty of Chemistry, Adam Mickiewicz University, Uniwersytetu Poznańskiego 8, Poznań 61-614, Poland. E-mail: magdrapp@amu.edu.pl; aniojka@amu.edu.pl
bDepartment of Clinical Chemistry and Molecular Diagnostics, University of Medical Sciences, Rokietnicka 3, 60-806 Poznań, Poland
First published on 17th February 2025
An important class of fluorescent dyes used in studying interactions and visualization of vital biomolecules are compounds with a skeleton origin 4,4-difluoro-4-bora-3a,4a-diaza-s-indacene, known as BODIPY. The objects of the presented study are the simple, polar and hydrophobic 3,7-dimethyl-substituted-core BODIPY dyes with the unmodified or modified phenyl aromatic ring at the meso position. Their optical properties as well as binding interactions with different DNA forms (i-motif, parallel G4, antiparallel G4, hybrid G4, dsDNA and ssDNA) were investigated by biophysical methods. The BODIPY derivatives interact more preferably with tetraplexes than other DNA forms. Especially, ligand 1 and 3 exibit tendency to destabilize parallel c-MYC G-quadruplex. The experiments with peroxidase-mimicking DNAzymes manifest that the main interaction between these BODIPY ligands and parallel G-quadruplex occurs via end-stacking mode. Moreover, their biological activity was evaluated by MTT assay.
Moreover, a conjugates of the fluorophore BODIPY and the macrocyclic heptaozaxole54 or Phen-DC355
were used to visualize G-quadruplexes (G4) in cell-free and cell-based assays systems. Lately, a few reports have been published studying the binding ability of BODIPY derivatives towards G4 DNA.56–59
G-quadruplex are DNA secondary structures formed in guanine-rich sequences60 and have been found to play an important role in regulating different biological processes such as replication, transcription, translation, and genomic stability. Their biological significance is a consequence of fact that guanine-rich sequences forming G-quadruplexes are present in different regions in the human genome, such as telomeres and the promoter region of different genes, including oncogene promoters.61–65 Therefore, the rational design of small molecules able to recognize, bind with high affinity and stabilize these non-canonical DNA become one of the key approaches in the development of new therapeutic and diagnostic strategies in cancer diseases.66 So far, a large number of G4-targeting ligands have been designed and synthesized, most of which exhibit common structural elements, such as a polycyclic aromatic core responsible for π-stacking interactions with the external G4 quartet and the presence of positively charged substituents that can interact electrostatically with the grooves of G4 DNA.67,68 Moreover, taking into account structural diversity of G-quadruplexes, an effort had been made to design and develop both stabilizing and destabilizing ligands which acquire specificity and selectivity without compromising affinity.69,70 For example, studies by Li et al. revealed high selectivity of BODIPY derivatives towards parallel c-MYC G4.59
Encouraged by these findings we decided to examine the potential of a few small BODIPY derivatives 1–6 to serve as a G4 ligands. The chosen BODIPY derivatives possess phenyl ring at meso position and 3,7-dimethyl-substituted or 3,7-dimethyl-2,8-disulfo-substituted-core. Our molecular target was the parallel G-quadruplex formed by the 5′-TGAGGGTGGGTAGGGT GGGTAA-3′ DNA sequence derived from the nuclease hypersensitive region of the c-MYC promoter (c-MYC G4).61 This choice was dictated not only by the fact that the c-MYC G4 is the most extensively studied and the best-characterized of all proto-oncogenes, due to its involvement in cellular proliferation and cell growth,71,72 but also by high selectivity of BODIPY derivatives towards c-MYC G4.59 All the ligands were tested against a parallel c-MYC G4 by means of UV-Vis, fluorescence and circular dichroism spectroscopy (CD). The CD technique is an essential tool providing information about the structure formation, stability and topology of G-quadruplexes, it is also commonly used to study the interaction between ligands and G-quadruplex structures.73–76 For example, the presence of induced exciton bands can be used as a method to distinguish between different binding modes G4 with ligand.73 The steady-state fluorescence measurements and DNAzyme assay were applied to explore binding properties between tested BODIPY dyes and c-MYC G4.
The BODIPY derivatives with best binding affinity were also tested against different DNA forms (i-motif, parallel G4, antiparallel G4, hybrid G4, dsDNA, and ssDNA) as well as their biological activity was evaluated.
Fluorine substituents diminish the π electron density of phenyl, enhancing the change in the quadrupole moment in aromatic ring facilitating participation in stacking79,80 and increasing the lipophilicity of compound. The ability of fluorine to participate in the formation of hydrogen bonds is crucial to biological activity of fluorochemicals.81 Therefore, we assumed that penta-fluorination of aromatic ring in 2 will increase the expected π–π interactions with G4 tetrad. The planar azide moieties in BODIPY dyes 341 and 5
41 offer the opportunity to interact with G-quadruplex structures acting as supramolecular template directing bioorthogonal macrocyclization.82 Moreover, the compounds 4 and 5 contain polar sulfonic groups, and were supposed to serve as a comparison to nonpolar compounds 1 and 2. Finally, the derivative 6
41 with aliphatic moiety containing positively charged (cationic) quaternary ammonium groups should presumably interact with the phosphate backbone to anchor the molecule in the DNA and to enhance their solubility. The role of “tails” substituents in better and more efficiently tuned selectivity with respect to G-tetrad stacking-only focused species was evidenced by many research.83 The boron-dipyrromethene (BODIPY) derivatives 1–6 were prepared accordingly to previously reported methods.23,41,77 Moreover, the single crystal of compound 2 was obtained and its structure confirms the presence of a pentafluorinated phenyl ring twisted in the relationship to BODIPY core containing boron centre (see: ESI†).
The compounds 1–3 are hardly soluble in aqueous solution, thus these compounds were dissolved in DMSO directly before spectroscopic measurements. Stock solutions of BODIPY derivatives 4–6 were prepared in high-purity water. The values of molar absorption coefficients for tested ligand (in Tris buffer) have been determined – both in presence and absence of potassium ions and are collected in Table S1.†
![]() | ||
Fig. 2 Fluorescence spectra of BODIPY derivatives 1–6 (1 μM). All spectra were measured in Tris–HCl buffer (10 mM, pH 7.2) containing 100 mM KCl. |
![]() | ||
Fig. 3 Normalized absorption spectra of 1 (A) and 3 (B) upon gradual concentration increase. All spectra were measured in Tris–HCl buffer (10 mM, pH 7.2) containing 100 mM KCl. |
From red-shifted absorption maximum of 1–3 and 6 the type of aggregates can be determined as J-aggregates – end-to-end stacking.87,88 Results obtained shows no evidence for aggregation of 4 and 5, as expected due to the repulsion between negatively charged sulfonic groups present in their structure (Fig. 1).
To verify the aggregation induced emission (AIE) properties90 of tested compounds we have examined their fluorescence in DMSO/water mixtures of different ratios (fw) (Fig. S5†). The significant increase in the fluorescence intensity of 1 and 2 was observed upon enhancing fw from 0 to 60% and 80%, respectively. With a further increase of fw to 99%, the AIE effect diminished which can be attributed to their precipitation.90 In contrast, the fluorescence intensity of 3 and 6 decreased gradually with increase of the water fraction. The latter changes in the fluorescence intensity suggested a significant aggregation quenched emission (ACQ) effect of 3 and 6 (Fig. S5C and F†).
At this point formation of aggregates can be compared to structural modification of fluorophores, introduced during molecules design, which hampers the rotation between donor and acceptor units and disables the charge transfer.84,86 Therefore, when 1–3 and 6 are introduced to the THF solutions with increasing water fractions the TICT is promoted by the solvent polarity, seen as initial decrease of emission intensity, until the point at which polarity of the mixture induces fluorophores aggregation. Disabled rotation reduces TICT process and aggregates revitalize their emission, until the water fraction reaches almost 100% and fluorescence intensity drops, which could be seen as the result of the formation of greater, non-emissive aggregates. Interestingly, despite being perfectly water-soluble, compound 6 displays properties similar to those of insoluble 1–3 rather than to 4 and 5. Introduction of quaternary ammonium groups guarantees water-solubility, however their location in long aliphatic chain causes 6 to be amphiphilic, which explains its partial similarity to behavior of 1–3 when introduced to mixtures of varying polarities.
Temperature is another variable known to impact emissive properties of systems with coexisting TICT and aggregation processes. It has been observed that cooling can cause the profound increase in emission intensity of many AIE systems by suppression of intramolecular rotation.91,92 In contrary, aggregation quenched emission (ACQ) systems should exhibit opposite behavior upon decrease of temperature. We observed the drop of fluorescence intensity when the temperature was increased from 0 to 30 °C for aqueous solutions of 1 and 4–6, while for 2 and 3 the opposite appeared and emission intensity increased (Fig. S6 and S7†). There are limited reports of BODIPY derivatives exhibiting AIE phenomenon which boost their emission intensity with the temperature increase and most of those concern measurements taken in nonpolar, organic solvents – like in series of benzodithiophene BODIPY derivatives synthesized by Sengupta et al.93 Searching for the cause of distinct and surprising behavior of 2, one should consider that introduction of fluorine atoms to the phenyl ring of 2 diminishes its electron density and could impact both TICT and AIE processes of this molecule. As fluorine is an EWG (electron-withdrawing group), phenyl part of a BODIPY becomes a weaker donor for charge transfer. Tendency for aggregates formation or their self-assembly's patterns could also be affected as electrons withdrawal would disfavor π–π interactions of aromatic rings. On the other hand, the presence of –CH2N3 moiety in 3 introduced another rotating junction – it has been noticed that additional sites for intramolecular rotations can result in less emissive fluorophores.84 On the other hand, rotation at the junction with azide moiety could also disturb overall molecule's planarity and as such alter the aggregation pattern or its dynamics. Such subtle changes are known to greatly impact both arrangement and spectral properties of molecular rotors.94,95
During the preparation of abovementioned experiments, we have also noticed that intensities of both absorption and fluorescence decreases in time if the sample solutions are vigorously mixed, but only for derivatives able to form aggregates – 1–3 and 6. Absorption spectra of all derivatives were measured for one hour, in 10 minutes intervals, with sample solutions being mixed on vortex mixer for 10 second after each measurement. Analogical measurement was taken for non-mixed samples. Obtained results (Fig. S8 and S9†) clearly indicates that mixing of fluorophore's solution decrease absorption intensity but only for aggregating ones, while there is no effect on absorption of non-aggregating compounds 4 and 5. Additionally, when used quartz cuvettes were emptied and field with DMSO as a wash, the absorption spectra characteristic for given BODIPY derivative still could be obtained. Once again, such observation was made only for aggregative compounds 1–3 and 6, while for 4 and 5 there were no traces of absorption after DMSO wash (Fig. S8 and S9†). We assume that nanosized aggregates are formed on the walls of a cuvette even thou they are macroscopically clear and transparent and the addition of DMSO causes them to come apart, which results in observed sharp spectra. Fluctuations of absorbance and fluorescence intensities while mixing are commonly observed among molecular rotors, as mixing can promote their rotation.96 Additionally, while vigorous mixing, the solution in the cuvette gets aerated. As the concentration of gasses such introduced to the solution increase, it decreases the density of a water-based solution.97 Less-dens media are also known to enhance molecular rotation and as such – favor TICT.
The variety of conditions and variables can influence on optical properties of BODIPY dyes. Consequently, data analysis could be challenging for systems in which multiple variables coexist.
Interestingly, we have also noticed that CD signal of c-MYC G4 at 264 nm did not decrease gradually after addition of BODIPY ligands. Therefore, we decided to investigate interaction between BODIPY ligands 1–6 and G4 by adding small amounts (0.25 equiv., 12 times) of BODIPY ligands to 5 μM solution of folded c-MYC G4. As shown in Fig. S11 and Table S2,† we can observe the fluctuations in intensity of the positive band peak at 265 nm as well as the negative band at range from 505–520 nm for ligands 1–5. In case of ligand 6 we rather noticed gradual decrease of the positive CD band. The UV-Vis spectra measured parallel with CD spectra also indicate the similar tendency for absorbance at 260 nm upon addition of ligands 1–6 (Fig. S11†), which also imply interactions between tested BODIPY ligands and folded c-MYC G4. The observed fluctuations and final decrease in CD signal indicated that BODIPY molecules 1–6 can be classified as G-quadruplex destabilizing compounds. However, caution should be exercised when making such conclusions based only on CD titrations used to study DNA/small molecule interactions because of possible ICD contributions and other possible artefacts (e.g., aggregation) that cannot be easily unravelled.98
What is important, none of the studied BODIPY-derived ligands did affect the position of the bands of the c-MYC G4, which indicate that the G4 retains the parallel topology. Observed fluctuations in the intensity of the CD signals after addition of BODIPY ligands suggests end-stacking interactions with G4 tetrad rather than process of disturbing the perfect arrangement of tetrads in the G-quadruplex by compounds 1–6.99,100 Moreover, in case of compounds 1 and 3 the induced negative signals (ICD) were clearly observed in the long-wavelength region where BODIPY-based ligands exhibit absorption band (around 500 nm) (Fig. 4). This observation rather excludes the groove binding to G-quadruplex structures, the binding mode with characteristic positive ICD signals in the long-wavelength absorption region of the ligand.76,101–106 The presence of ICD suggests close proximity between the BODIPY ring of ligand 1 or 3 and G-tetrad(s) and strong π–π interactions between them.107 Using CD, we also checked if ligands 1–6 were able to induce G-quadruplex structure formation. We found that the addition up to 7 equivalents of tested ligands did not result in the formation of the G-quadruplex structure (Fig. S12†). Based on our experiments we supposed that the level of fluctuations in CD spectra observed both for G4 and single-stranded DNA (Fig. S13†) after binding with BODIPY ligands comes from each ligand ability to form aggregates due to binding to DNA and complexity equilibria between free ligands and DNA/BODIPY complexes. As it has been proposed for porphyrins, there exist a few possibilities of ligand/G4 interactions, which are in equilibrium.108 Thus, different species are present in solution, among them there are free ligand monomers, monomers bound to DNA and ligand aggregates formed on the surface of DNA molecules. We supposed that this scenario is possible in case of studied BODIPY ligands and it is responsible for fluctuations observed in their CD and UV spectra.
To verify if the main way of interaction between BODIPY derivatives and G4 is through end-stacking interaction, we decided to use a G4/hemin peroxidase inhibition assay.57,109 G4-hemin complex behaves as DNAzyme and exhibits peroxidase activity, i.e. it is able to oxidize ABTS in the presence of H2O2.110 The hemin molecule binds to the DNA oligonucleotide through an end-stacking interaction between the pyrrole rings of porphyrin and DNA nucleobases. This interaction mode requires the easy access of hemin to the G-quartet; thus, the G4s forming parallel topologies are known to produce DNAzymes exhibiting higher peroxidase activity.111 One should remember that the examined c-MYC G4 structure has a sidewise-type of loops, which facilitates ligand/G-tetrad π–π stacking interaction. We assumed that our ligands would compete with hemin to bind on the G-tetrad surface, that should cause the decrease in the catalytic oxidation ability of DNAzyme. Thus, for our purpose we used two parallel G-quadruplexes: c-MYC G4 as well as Cat G4. The latter one possesses the sequence 5′-TGG GTA GGG CGG GTT GGG AAA-3′and it also adopts the parallel topology.112 Among tested ligands only ligands 1 and 3 inhibit the activity of both complexes: c-MYC G4/hemin as well as Cat G4/hemin (compare Fig. 5, S14 and S15†). The fact that such results were obtained for those two ligands is in a good agreement with previously discussed CD spectra, in which we observed ICD bands only for these ligands. However, the clear inhibition activity of c-MYC G4/hemin or Cat G4/hemin complexes was observed only in presence of BODIPY derivative 1 as indicated by their initial reaction rates 5.0 × 10−4 [1 s−1] and 10.0 × 10−4 [1 s−1] compared with values 8.0 × 10−4 [1 s−1] and 17.0 × 10−4 [1 s−1] for c-MYC G4/hemin or Cat G4/hemin, respectively. These results support the conclusion that binding process of tested BODIPY-based compounds on the parallel G-quadruplexes DNA can occur via end-stacking interaction with external G-tetrads of G-quadruplexes. According to literature61 the 5′-end G-tetrad is more accessible, providing a larger hydrophobic surface, where the ligand could more easily stack and interact stronger compared to the 3′-end G-tetrad. Taking into account the results of DNAzyme assay, we can assume that our BODIPY ligands are rather stacked on G-tetrad at 3′end as hemin has been demonstrated to bind selectively to the 3′-terminal G-tetrad of parallel G4s.113
We may also speculate the BODIPY molecules may induce the changes in conformation of the flanking bases at 3′ termini as was demonstrated for quindoline/c-MYC114 or 9CI/c-MYC G-quadruplex interactions.115 Moreover, such structural arrangement could also be explanation of the ICD band observed in case of ligand 1 and 3. We realize that such hypothesis should be further verified by molecular dynamic simulations or NMR studies.
To assess the affinity of BODIPY ligands 1–6 to G-quadruplex structure, we performed fluorescence titration experiment by adding incremental c-MYC G4 to solution of BODIPY and vice versa (reverse titration). Addition of c-MYC G4 to 1 μmolar solution of ligands causes quenching of their fluorescence except of BODIPY 4 and 5 (Fig. 6A). In case of 2 and 6 a dramatic decrease was observed in fluorescence intensity (around 90%) with no shifts for the maximum peaks. However, the reference experiments reveled that emission of ligands also decreases upon addition of the same amount (in ml) of KCl/Tris buffer (Fig. 6B). This effect we attributed to the aggregation of these ligands, which we demonstrated earlier (see part 2.2). As previously, the fluorescence emission of ligands 4 and 5 (with sulfonic group in BODIPY core) was stable and almost unaffected by addition of salt. Taking into account that observed changes in fluorescence of ligands 1–3 and 6 are not only caused by G4 but also by salt (or only mixing) it became clear that these data could not be used to obtain the reliable values Kb of binding constants. Thus, we decided to perform reverse fluorescence titrations to estimate the binding affinities of the tested BODIPY compounds. Such approach was used by other researchers studying ligands/G4 complexes,116,117 although it seems to be controversial due to fact that the fluorescence molecule is added to G4 solution and the inner filter effect (IFE) can contribute to level of final emission.
Anyway, we believe that obtained results at least can help to distinguish which of tested ligands binds with highest affinity to G-quadruplex. In our case, the ligands were added to 2 μM G4 solution in stepwise mode to reach the final concentration 10 or 14 μM (Fig. S16†). Interestingly, for emission of ligand 2 we observed the increase of weak fluorescence band located near 505 nm, which we attributed to aggregates formation. In the case of hydrophilic ligands 4–6, no plateau was observed for titration plots. Next, such obtained, reverse fluorescence titration data were fitted to the non-competitive DNA-binding model.118 The calculated binding constants, Kb values may be ranked as follows: 7.80 ± 0.03 × 106 for 3, 5.40 ± 2.10 × 106 for 1, 2.51 ± 0.80 × 105 for 6, 1.92 ± 0.10 × 105 for 5, 1.48 ± 0.03 × 105 for 4 and 1.30 ± 0.03 × 105 for 2. The obtained results indicated that BODIPY 1 and 3 bind with high affinity to G4 c-MYC, that is in good agreement with observed DNAzyme inhibition, as only ligands 1 and 3 bind to parallel G4 with affinity sufficient to compete with hemin molecule.
To investigate the binding stoichiometries of tested BODIPY derivatives with the G-quadruplex we applied a fluorescence-based Job plot. In this method, the total molar concentration of BODIPY ligand and c-MYC G4 were kept constant, but their mole fractions were varied.119 The valuable plots were obtained only for ligand 1 and water-soluble ligands 4–6. As shown in Fig. S17,† inflection points in the graphs are observed at the DNA fraction of 0.5 for all complexes. These results indicated that BODIPY ligands bind to c-MYC G4 DNA by 1:
1 binding stoichiometry. One should remember, that Job plot stoichiometry represents only the major binding event and the higher binding stoichiometries are difficult to determine precisely in case of ligands able to bind by multiple modes such as end-stacking, electrostatic interactions, groove binding.120 We cannot exclude such scenario, especially for ligands 2 and 3 due to their tendency to aggregation, and thus possibility of self-association on the DNA backbone, analogously as proposed by Gandini et al. for porphyrins.108 Nevertheless, the main interaction of our BODIPY ligands with G4 seems to be end-stacking process with 1
:
1 ratio.
Thus, we decided to check their ability to discrimination between parallel G-quadruplex topology c-MYC G4 and other structural DNA form. For our studies we used hybrid-type G-quadruplex 22HT/K+, antiparallel G-quadruplex 22HT/Na+, double-stranded DNA ds26, single-stranded DNA and i-motif (iM) c-MYC (Table S3†).60 Therefore, we acquired UV melting curves for different DNA tetraplexes with and w/o BODIPY ligands 1 and 3 at 295 nm (Table 1 and Fig. S18†). The most pronounced effect can be seen for ligand 3/c-MYC G4 complex with Tm 70.0 °C, which is 6.4 °C lower than for free c-MYC G4. Interestingly, changes in Tm indicate that ligand 1 and 3 exhibit destabilizing effect for G4 with different topologies, but in contrary exhibits stabilizing effect for iM. It is worth to mention that the rest among tested ligands seems to stabilize c-MYC G4 (Fig. S19†).
System | Tm ± SD [°C] | ΔTm [°C] | ΔH [kcal mol−1] | ΔS [cal mol−1 K−1] | ΔG (310 K) [kcal mol−1] | |
---|---|---|---|---|---|---|
c-MYC G4 (parallel) | DNA | 76.4 ± 0.5 | −52 | −149 | −5.8 | |
+Ligand 1 | 76.2 ± 0.5 | −0.2 | −38 | −108 | −4.5 | |
+Ligand 3 | 70.0 ± 0.3 | −6.4 | −28 | −83 | −2.3 | |
22HT/Na+ (antiparallel) | DNA | 57.5 ± 0.5 | −31 | −93 | −2.2 | |
+Ligand 1 | 56.1 ± 0.4 | −1.4 | −37 | −111 | −2.6 | |
Ligand 3 | 54.0 ± 0.2 | −3.5 | −32 | −99 | −6.3 | |
22HT/K+ (hybrid) | DNA | 63.2 ± 0.5 | −35 | −104 | −2.8 | |
+Ligand 1 | 61.3 ± 0.3 | −1.9 | −32 | −96 | −2.2 | |
+Ligand 3 | 62.5 ± 0.1 | −0.7 | −30 | −91 | −1.8 | |
iM c-MYC (i-motif) | DNA | 53.9 ± 0.5 | −60 | −184 | −3.0 | |
+Ligand 1 | 54.3 ± 0.2 | + 0.4 | −57 | −174 | −3.1 | |
+Ligand 3 | 55.2 ± 0.1 | + 1.3 | −56 | −170 | −3.3 |
A thermodynamic analysis of the denaturation profiles (shown in Fig. S18†) was performed according to two-state model, where the folded G-quadruplex is in equilibrium with the unfolded strand. The enthalpy and entropy were calculated from classical relation as follows: ΔG = −RTln
(K) = ΔH − TΔS. The equilibrium constant K can be written as K = f/(1 − f), where for intramolecular equilibrium f is the fraction of folded oligonucleotide. In all cases G-quadruplex or i-motif formation is enthalpically favorable and entropically unfavorable as observed previously.121 The thermodynamic data collected in Table 1 have shown that ligand/G4 binding events are the exotermic reactions. For all oligonucleotides the ΔG of tetraplex formation was negative in the appropriate buffer at physiological temperature (37 °C), indicating that the predominant species was the folded form in this temperature (Table 1). The interaction with ligands is associated with a net free energy change (ΔG). For example, ligand 3 binds to c-MYC G4 causing ΔG change from −5.8 to −2.3 kcal mol−1, whereas to 22HT/Na+ causing ΔG change from −2.2 to −6.3 kcal mol−1. Thus, ligand 3 pronounces destabilizing effect on c-MYC G4 and opposite on 22HT/Na+.
The CD spectra of pre-annealed solutions of G-quadruplexes of different topologies as well as i-motif and single-stranded and double-stranded DNA were recorded prior addition of BODIPY ligands in stepwise mode as previously (Fig. 7, S12A, C and S20†). In general, none of ligands (1 and 3) affect the positions of the CD signals, suggesting that the DNA topology is unchanged. Moreover, no significant changes in the long-wave part of CD spectra were observed even after the addition of 7 equivalents of ligand 3 to antiparallel G4 (22HT/Na+) or hybrid G4 (22HT/K+). The latter results suggest that binding mode of 3 to these G4s can be different than in complexes with parallel c-MYC G4. However, the ligand 1 induced CD band around 510 nm in case of antiparallel G4 (22HT/Na+) or hybrid G4 (22HT/K+), analogously to parallel G4 (c-MYC G4) (Fig. S20A and B†). Summarizing, the performed CD titration experiments of the different DNA structure with 1 and 3 also showed that these BODIPY ligands interact more preferably with tetraplexes than ds26 or ssDNA. The interaction of BODIPY derivatives 1–6 with iM will be explored deeply in the future.
Evaluation of dose-response to tested ligands 1 and 3 measured at rising concentrations indicated their non-cytotoxic effect. Both tested ligands slightly affected cell viability in a dose-dependent, but not significant manner. It was observed that ligand 1 after 24 h treatment increased T-47D viability from 12% to 18% for 5 μM or 50 μM, while viability reduction from 5–25% was observed for MRC-5 cells. To compare, the slight cytotoxic activity against T-47D and for MRC-5 cells after 48 h (viability >94% and >90%, respectively) and 72 h (viability >85%) was detected. Only in the case of the concentration 50 μM of BODIPY 1, the highest effect on MRC-5 cell after 48 h of incubation (viability 57%) was observed. Interestingly, the opposite tendency occurred in case of both cell lines ligand 3 treatment. Thus, the higher viability (from 107% to 117% for 1 μM ligand) of non-tumorigenic MRC-5 cells was observed after 24 h of incubation compared to 102% (1 μM) for T-47D cells. The results for T-47D and MRC-5 cells after 48 h treatment with 0.1 μM ligand 3 were 85%, compared to viability 82% after 72 h of incubation (Fig. 8).
Summarizing, in the ligand 1 and 3 concentration range of 0.1–50 μM, the chosen BODIPY derivatives indicated non-cytotoxicity against non-tumorigenic lung fibroblast cells, nor breast cancer cells. The interactions of ligands 1 and 3 with G4 structures demonstrated during our studies in vitro seems to be not associated with their anticancer activity, as it was shown by in cellulo evaluations. On the other hand, the non-cytotoxic properties of the BODIPY dyes 1 and 3 would be the assets in their potential applications in neurological disorders and age-related diseases investigations, as reported for other small molecules that destabilize G-quadruplex.122
However, more in-depth studies in the context of neurodegenerative diseases should be performed. G4 and i-motif-binding compounds are known for their ability to bind promoters and affect the metabolism of mitochondria that are critical for age-related and neurodegenerative diseases.123 Thus, further studies should involve evaluation of cortical neurons in the context of ROS levels assessment followed by evaluation of autophagy, inflammation and gene expression analysis. Similarly, such studies should be accompanied by 3D cell culture models and animal studies evaluation. Moreover, the evaluated compounds 1 and 3 seem to be the good candidates for the bio-imaging as the BODIPY-based compounds have already been found as fluorescent probes in Alzheimer disease and other related processes.124
Ligands 1–3 were dissolved in DMSO, whereas ligands 4–6 were dissolved in H2O for obtaining stock solutions of ∼1.5 mM concentration and were stored at 4 °C. The ligands 2 and 3 were stocked no longer than 48 h.
Among tested BODIPY derivatives, ligands 1 and 3 indicated the highest affinity to parallel c-MYC G4 and are able to compete with hemin to bind on the G-tetrad surface. The main mode of those interactions occurs via end-stacking and binding stoichiometry 1:
1. The end-stacking binding mode with external G-tetrads was also confirmed by circular dichroism spectroscopy titration. The strong negative ICD in the long wavelength absorption region of achiral BODIPY-based ligands 1 and 3 around 500 nm indicates that π–π end-stacking interactions may be the dominant interaction mode between compounds 1 or 3 and c-MYC G4. The CD spectra titrations also revealed that none of tested ligands with BODIPY skeleton were able to induce formation of a stable G4. Moreover, CD studies upon interactions of 1 and 3 with different DNA forms indicated that these ligands interact more preferably with tetraplexes (both G4s and i-motifs) than ssDNA and dsDNA. However, the ligand 1 induced CD band around 510 nm in case of antiparallel G4 (22HT/Na+) or hybrid G4 (22HT/K+), which suggest that binding mode of 1 to these G4s can be similar to complexes with parallel c-MYC G4. A decrease in CD signal intensity (the position of signals has not changed) recorded upon the addition of up to five equivalents of ligands, suggesting the presence of ligand-dependent perturbations in the studied G4 structure. However, this may be the result of the stacking interaction and not the destabilization of G4. To verify this hypothesis we performed a thermodynamic analysis of the melting profiles of different topology G4s with and w/o BODIPY derivatives. The latter results confirmed that the ligand 1 and 3 interacted the strongest with c-MYC G4 as changes in enthalpy and entropy are the most pronounced and a net free energy changes confirm the destabilizing effect of these ligands on c-MYC G4.
Summarizing, we showed that the simple BODIPY dyes with the unmodified or modified phenyl aromatic ring at meso position may serve as novel G4-recognizing ligands. In case of two ligands, we observed the distinct tendency toward disrupting the G4 structure and for this reason these ligands are good candidates for bio-imaging or to study neurodegenerative as well as age-related diseases.
Footnote |
† Electronic supplementary information (ESI) available: Experimental details, compounds characterization (NMR spectra, X-ray structure), spectroscopy measurements (CD, UV, FL), DNAzyme and MTT assays. CCDC 2241316. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4ra08464k |
This journal is © The Royal Society of Chemistry 2025 |