DNA Interactions and Biological Activity of 2,9-Disubstituted 1,10-Phenanthroline Thiosemicarbazone-Based Ligands and a 4-Phenylthiazole Derivative

1. Introduction

Since the discovery of the double structure of DNA [1,2], numerous efforts have been made to understand the different biological functions and structures that nucleic acids can adopt. Currently, it is widely recognized that DNA is a structurally dynamic molecule, with the capacity to adopt a variety of alternative secondary structures apart from the double helix, such as DNA hairpins, Holliday junctions, cruciforms, triplexes or G-quadruplexes (G4) [3], just to mention a few. Among them, DNA and RNA quadruplexes have been gaining prominence in the last years for their biological relevant roles, including their establishment as therapeutic targets for cancer treatment [4], but also for many other applications in biomedicine and nanotechnology [5,6,7]. Regarding cancer, it is a leading cause of death worldwide, accounting for almost 10 million of deaths in 2020, or nearly one in six deaths, with an estimated 19.3 million new cancer cases that year [8].

The studies conducted by Bang and Gellert revealed that guanylic acid tends to aggregate in solution, resulting in jelly-like structures [9,10]. Subsequent studies showed that four guanylic acid molecules can be arranged to form a planar tetrad [11,12]. In this tetrad, each guanine can function as a donor and acceptor of two pairs of hydrogen bonds, called Hoogsteen H-bonds. In addition, these tetrads tend to stack on top of each other through π-π stacking interactions and are stabilized by positive physiological cations (mainly K⁺ and Na⁺), thus giving rise to three-dimensional non-canonical secondary structures designated as G-quadruplexes (G4s) [13]. G4s can be unimolecular or intermolecular and can adopt a wide variety of topologies arising from several combinations of strand direction and loop sequence composition [14]. G4s have been associated to many different cellular processes and are present in a vast number of relevant genomic regions: in telomeres, in replication origins, the 5'-UTR regions, and in several oncogenic promoters. Among their biological applications, the consideration of telomeric quadruplexes as potential anticancer targets deserves special attention. Telomeres are structures found at the ends of chromosomes, made up of regions of non-coding DNA with highly guanine-rich repetitive sequences. Their main function is to conserve genetic information, preserving genes' integrity and functionality over time. Specifically, telomeric DNA is composed of a repeat of nucleotides rich in guanines (TTAGGG for eukaryotic cells) and has a length of between 9,000 and 15,000 nucleotides [15].

The activity of ribonucleoprotein telomerase is intimately ligated to the type of secondary DNA and RNA structures formed in the telomeric region. It is responsible for maintaining telomere length and DNA synthesis through reverse transcription of an RNA template contained as part of its structure. In most somatic cells there is a progressive shortening of their telomeres, while in embryonic, stem and epithelial cells, telomerase activity is upregulated. Interestingly, this enzyme is overexpressed in approximately 80-90% of human tumor cells. Therefore, inhibition of this enzyme has been postulated as a therapeutic strategy against cancer [16]. Several telomerase inhibitors have already entered clinical or preclinical trials [17,18].

Stabilization of telomeric G4 structures through small molecules can prevent telomerase function, triggering apoptosis in cancer cells. In general, these G4 ligands fulfil a series of structural characteristics, possessing a) an aromatic π-system with enough surface to favor end stacking interactions with the quadruplex, b) an electronic deficiency near the center of the G-quadruplex, with a stabilizing effect similar to that produced by K⁺ or Na⁺ ions and c) one or more flexible side-chains positively charged at physiological pH, providing additional interactions with the grooves and/or loops [19,20,21,22]. However, it should be noted that there are some examples in the literature, such as the case of the molecule funtumin which, without fulfilling the mentioned characteristics, can interact with G4s [23]. Considering these generally accepted structural features, a wide variety of molecules have been designed and studied as G4 ligands, mainly those belonging to three families: fused polyaromatic derivatives, macrocyclic compounds and non-fused aromatic systems. The most representative examples of compounds that have shown relevant antitumor properties [24,25,26] include 2,6-diamidoanthraquinone, BisA, BRACO-19 (trisubstituted acridine), PIPER (perylene diimide derivative), Telomestatin (macrocyclic compound), Pyridostatin, TmPyP4 (cationic porphyrin derivative), Quarfloxacin, CX-5461, 360A, PhenDC3, benzofuran derivatives, etc. From these, it is worth mentioning that Quarfloxacin and CX-5461 are currently in Phase I/II clinical trials [27,28].

We have recently reported a new G4 ligand derived from the 1,10-phenanthroline heterocycle (Figure 1, A), compound PhenQE8 [29], which has proven cytotoxic in several types of tumor cells (PC-3, HeLa and MCF-7), with IC₅₀ values in the low micromolar concentration range. At the same time, it exhibited lower cytotoxicity in the normal cells evaluated, HFF-1 and RWPE-1, demonstrating a good selectivity potential. As the interaction assays demonstrated its ability to selectively interact with telomeric DNA G-quadruplex compared to DNA duplex, our interest is now focused on the study of analogues of PhenQE8. We pursue in this work the synthesis and biological study of 1,10-phenanthroline derivatives substituted with hydrazinecarbothiamide (Thiosemicarbazone derivatives, TSC) or hydrazine 4-phenylthiazole chains.

It is important to highlight that thiosemicarbazones (TSCs) have been and continue to be used as versatile scaffolds for ligand synthesis in medicinal chemistry (Figure 1, B). They have a wide range of biological activities including antibacterial, antifungal, antituberculosis, anticonvulsant, antioxidant and anticancer. In addition, they constitute the precursor structures for the synthesis of several types of heterocycles, such as 1,2,4-triazole-3-thiones, 4-thiazolidinones, oxadiazoles, and 1,3,4-thiadiazoles, all of which are of enriched biological interest. Their mechanisms of action have been extensively studied but remain unclear. Some hypotheses that are worth mentioning include, among others, the inhibition of cellular iron uptake from Tf, iron mobilization from cells, the inhibition of ribonucleotide reductase activity, the up-regulation of the metastasis suppressor protein, and the formation of redox active metal complexes that give rise to reactive oxygen species (ROS) [31,32]. Thiosemicarbazones have also been used in the form of bidentate ligands or Schiff bases with various metals, e.g. Pt(II) and Au(III), as antitumor agents, including G-quadruplex DNA or thioredoxin reductase as therapeutic targets [33].

Based on the previous considerations, we report herein the preparation and characterization of three 2,9-disubstituted 1,10-phenanthroline compounds that contain hydrazinecarbothiamide in their side chains (thiosemicarbazones, TSCs), one as a neutral ligand and two as its corresponding PF₆¯ and BF₄¯ ionic salts (compounds 1-3, Figure 2). The neutral ligand 1 is a known compound [34], and it was used as a precursor to obtain a heterocyclic derivative (compound 4), in which the 1,10-phenanthroline is substituted with two hydrazine 4-phenylthiazole motifs. Quadruplex and dsDNA interactions have been preliminarily studied, especially for the neutral derivative 1, by DNA FRET-based melting assays, equilibrium dialysis, circular dichroism (CD) and viscosity assays. In addition, cytotoxic activity against several tumor cells (PC-3, DU145, HeLa, MCF-7 and HT29), and two non-tumor cell lines (HFF-1 and RWPE-1) has been tested. Cell cycle analysis and Annexin V/PI assays have been carried out to preliminary assess effects on cell cycle phases and whether the novel compounds can produce cell death by apoptosis. Their biological activity is compared to other known 1,10-phenanthroline G4 ligands endowed with antitumor properties, such as PhenDC3 [35] or PhenQE8.

2. Materials and Methods

2.1. General methods

All reagents and solvents were acquired from commercial sources in high purity (Merck, BLD Pharm, Carlo Erba Reagents, PanReac AppliChem-ITW Reagents, Fisher, Scharlab, Acros Organics, Thermo Scientific Chemicals) and used as received. 1,10-phenanthroline-2,9-carbaldehyde was synthesized following a reported procedure [29]. Analytical TLC was performed on precoated aluminum silica gel sheets (Merck 60 F₂₅₄ 0.25 mm). Melting points were measured on an Electrothermal Digital IA9100 equipment and are reported uncorrected. IR spectra from powdered solids were recorded on an IR Bruker Alpha II spectrophotometer, under ambient conditions between 4000-400 cm^-1 on a diamond ATR unit. ¹H NMR and ¹³C NMR spectra were registered on a Bruker 400 Ultrashield in DMSO-d₆ at 298 K. Chemical shifts (δ) are reported considering DMSO-d₆ (δ_H 2.50 ppm), DMSO-d₆ (δ_C 39.5 ppm). Coupling constants values (J) are reported in Hz. Signal multiplicity is designated with the notations s (singlet), bs (broad singlet), d (doublet), t (triplet), at (apparent triplet) or hept (heptuplet). Electron spray ionization mass spectra (ESI-MS) were recorded on a MAXISII instrument (SiDI, Universidad Autónoma de Madrid, UAM). Combustion elemental analysis was conducted on a LECO CHNSO-932 apparatus.

Short oligonucleotides sequences were acquired in high purity from IDT® (Integrated DNA Technologies) and molecular biology grade reagents from Sigma-Aldrich (Merck). DNA stock solutions were prepared using BPC-grade water. In all DNA interaction studies, Milli-Q water was employed for the preparation of the working buffers. Stock solutions of compounds 1-4 at ~ 1 mM concentration were made freshly immediately before each assay in water and a maximum 10% DMSO, and then gradually diluted to desired final concentrations in the corresponding buffer or cell culture medium.

UV-visible experiments were conducted using a Perkin-Elmer Lambda 35 spectrophotometer. CD spectra were recorded employing a Jasco J-715 spectropolarimeter. FRET DNA melting assays were performed in microplate-format on an QuantStudio 6 PRO 0.2 mL, 96-well package (Applied Biosystems). Flow times for viscosity determinations were recorded at 25 °C using a Visco System AVS 470, with a microUbbelohde (K= 0.01018) capillary viscometer. MTT viability assays were performed by specialized personnel at the Cell Culture Unit (UAH) following standard reported protocols and using an ELISA plate reader (ELX 800 Biotech Instruments). Cell cycle and apoptosis assays were performed by flow cytometry on a MACSQuant® Analyzer 10 Flow Cytometer (Miltenyi Biotech).

2.2. Synthesis and characterization of chemical compounds

2. ,2’-(1,10-phenanthroline-2,9-diyl)-bis(methylene)-bis(hydrazine-1-carbothioamide) (1)

1,10-phenanthroline-2,9-carbaldehyde (98 mg, 0.42 mmol) was dissolved in dry MeOH (8 mL) with stirring, and then thiosemicarbazide (76 mg, 0.84 mmol) and 3 drops of 96% H₂SO₄ were added. The reaction mixture was then heated at reflux for 1 h, cooled at room temperature and filtered. The resulting solid was washed once with cold water, then with MeOH and finally three times with Et₂O, and dried under vacuum. The solid was recrystallized using a mixture of 9:1 AcCN/MeOH, providing compound 1 (128 mg, 81%) as an orange solid. M.p.: > 180 °C (decompn.). IR: v (cm^-1): 3250, 3143, 2944, 1592, 1513, 1491, 1200, 1115, 1043. ¹H-NMR (400 MHz, DMSO-d₆): δ (ppm) 11.99 (s, 2H, NH), 8.77 (d, J = 8.6 Hz, 2H, H-4,7), 8.57 (d, J = 8.6 Hz, 2H, H-3,8), 8.51 (bs, 2H, NH), 8.48 (s, 2H, N=CH), 8.42 (bs, 2H, NH), 8.05 (s, 2H, H-5,6). ¹³C-NMR (101 MHz, DMSO-d₆): δ (ppm) 179.1 (-C=S), 153.7 (C-2,9), 143.1 (C-10a,10b), 141.4 (CH=N), 138.4 (CH, C-4,7), 129.4 (C-4a,6a), 127.6 (CH, C-5,6), 121.1 (CH, C-3,8). HRMS (ESI, m/z) for C₁₆H₁₅N₈S₂ [M + H]⁺: calculated: 383.0861; found: 383.0857. HRMS (ESI, m/z) for C₁₆H₁₄N₈S₂Na [M + Na]⁺: calculated: 405.0680; found: 405.0673. Anal. Calcd. (%) for C₁₆H₁₄N₈S₂·6H₂O: C, 39.17; H, 5.34; N, 22.84. Found: C, 38.77; H, 4.82; N, 17.47. 1,10-phenanthroline derivatives are known to crystallize as hydrates and/or solvates.

2. ,2’-(1,10-phenanthroline-2,9-diyl)-bis(methylene)-bis(hydrazine-1-carbothioamide) hexafluorophosphate salt (2)

1,10-phenanthroline-2,9-carbaldehyde (50 mg, 0.21 mmol) was dissolved in dry MeOH (4 mL) with stirring, and then thiosemicarbazide (39 mg, 0.42 mmol) and 3 drops of 96% H₂SO₄ were added. The reaction mixture was heated at reflux for 1 h, and then KPF₆ (0.42 mmol, 78 mg) was added. The mixture was then stirred at reflux for another 30 min and 30 min at room temperature, after which the solid obtained was filtered and washed first with cool MeOH once, and then three times with Et₂O. The product was recrystallized using a mixture of 9:1 AcCN/MeOH, and filtered again to provide compound 2 (51 mg, 46%) as an orange solid. M.p.: > 240 °C (decompn.). IR: v (cm^-1): 3346, 3249, 3154, 2944, 1592, 1517, 1200, 1053, 845, 751. ¹H-NMR (400 MHz, DMSO-d₆): δ (ppm) 12.02 (s, 2H, NH), 8.79 (d, J = 8.6 Hz, 2H, H-4,7), 8.60 (d, J = 8.6 Hz, 2H, H-3,8), 8.52 (bs, 2H, NH), 8.49 (s, 2H, N=CH), 8.42 (bs, 2H, NH), 8.07 (s, 2H, H-5,6). ¹³C-NMR (101 MHz, DMSO-d₆): δ (ppm) 179.0 (-C=S), 154.0 (C-2,9), 144.5 (C-10a, 10b), 142.4 (CH=N), 137.6 (CH, C-4,7), 129.3 (C-4,6a), 127.5 (CH, C-5,6), 120.7 (CH, C-3,8). ¹⁹F-NMR (376 MHz, DMSO-d₆): δ (ppm) -70.14 (d, J = 710.4 Hz). ³¹P-NMR (162 MHz, DMSO-d₆): δ (ppm) -144.20 (hept, J = 710.4 Hz). HRMS (ESI, m/z) for C₁₆H₁₄N₈S₂Na [M + Na]⁺: calculated: 405.0680; found: 405.0666. HRMS (ESI, m/z) for PF₆^¯ [M]^¯: calculated: 144.9641; found: 144.9646. Anal. Calcd. (%) for C₁₆H₁₅N₈S₂PF₆: C, 36.37; H, 2.86; N, 21.20. Found: C, 38.11; H, 3.98; N, 19.77.

2. ,2’-(1,10-phenanthroline-2,9-diyl)-bis(methylene)-bis(hydrazine-1-carbothioamide) tetrafluoroborate salt (3)

An analogous procedure than the one used to obtain compound 2 was used for the preparation of compound 3, using the same starting materials (1,10-phenanthroline-2,9-carbaldehyde (50 mg, 0.21 mmol) in dry MeOH (4 mL), thiosemicarbazide (39 mg, 0.42 mmol) and 3 drops of 96% H₂SO₄) heating at reflux for 1 h and then adding NaBF₄ (46 mg, 0.42 mmol). Compound 3 was obtained as an orange solid (51 mg, 52%). M.p.: > 240 °C (decompn.). IR: v (cm^-1): 3349, 3249, 3151, 2953, 1592, 1521, 1446, 1205, 1050, 861. ¹H-NMR (400 MHz, DMSO-d₆): δ (ppm) 12.01 (s, 2H, NH), 8.78 (d, J = 8.6 Hz, 2H, H-4,7), 8.59 (d, J = 8.6 Hz, 2H, H-3,8), 8.52 (bs, 2H, NH), 8.49 (s, 2H, N=CH), 8.43 (bs, 2H, NH), 8.08 (s, 2H, H-5,6). ¹³C-NMR (101 MHz, DMSO-d₆): δ (ppm) 179.0 (-C=S), 153.9 (C-2,9), 144.2 (C-10b, 10a), 142.2 (CH=N), 137.7 (CH, C-4,7), 129.3 (C-4, 6a), 127.5 (CH, C-5, 6), 120.7 (CH, C-3,8). ¹⁹F-NMR (376 MHz, DMSO-d₆): δ (ppm) -148.22, -148.28. HRMS (ESI, m/z) for C₁₆H₁₄N₈S₂Na [M + Na]⁺: calculated: 405.0680; found: 405.0669. HRMS (ESI, m/z) for BF₄^¯ [M]^¯: calculated: 87.0029; found: 87.0070. Anal. Calcd. (%) for C₁₆H₁₅N₈S₂BF₄ ·3H₂O: C, 36.65; H, 4.04; N, 21.37. Found: C, 38.09; H, 3.80; N, 19.81.

2. ,9-bis((E)-(2-(4-phenylthiazol-2-yl)hydrazineylidene)methyl)-1,10-phenanthroline (4)

To a solution of 2,2’-(1,10-phenanthroline-2,9-diyl)-bis(methylene)-bis(hydrazine-1-carbothioamide) (30 mg, 0.16 mmol) in EtOH (8 mL) 2-bromoacetophenone (31 mg, 0.16 mmol) was added, and the resulting mixture was stirred under reflux for 24 h. Afterwards, the crude was cooled at room temperature, and centrifuged to 1500 rpm for 5 min. The supernatant was extracted, and the resulting precipitate was filtered, washed with Et₂O and dried under vacuum, obtaining compound 4 (34 mg, 75%) as an orange solid. M.p.: > 250 °C (decompn.). IR: v (cm^-1): 3245, 2839, 1540, 1439, 1267, 1147, 1037. ¹H-NMR (400 MHz, DMSO-d₆): δ (ppm) 8.73 (d, J = 8.0 Hz, 2H, H-4,7-phenanthroline), 8.60 (s, 2H, H-5,6-phenanthroline), 8.38 (d, J = 8.0 Hz, 2H, H-3,8-phenanthroline), 8.15 (s, 2H, H-thiazole), 7.90 (d, J = 7,6 Hz, 4H, H-2,6-phenyl), 7.50 (s, 2H, HC=N), 7.44 (at, J = 7,8 Hz, 4H, H-3,5-phenyl), 7.34 (at, J = 7.3 Hz, 2H, H-4-phenyl). ¹³C-NMR (101 MHz, DMSO-d₆): δ (ppm) 167.8 (S-C=N), 152.8 (C=N, thiazole), 151.2 (C-2,9-phenanthroline), 140.0 (C-4a,6a-phenanthroline), 138.5 (CH=N), 134.7 (C-4,7-phenanthroline), 129.2 (C-3,5-phenyl), 128.8 (C-4-phenyl), 128.3 (C-5,6-phenanthroline), 128.3 (C-10a,10b-phenanthroline), 127.5 (C-3,8-phenanthroline), 126.0 (C-2,6-phenyl), 105.8 (C-S, thiazole). HRMS (ESI, m/z) for C₃₂H₂₃N₈S₂ [M + H]⁺: calculated: 583.1487; found: 583.1472. HRMS (ESI, m/z) for C₃₂H₂₂N₈S₂Na [M + Na]⁺: calculated: 605.1306; found: 605.1293.

2.3. DNA interaction studies and biological activity

2.3.1. Fret-based DNA melting assays

FRET DNA melting assays were conducted following standard reported procedures [29]. Briefly, the employed oligonucleotides F21T (5'- FAM – (GGG TTA)₃GGG – TAMRA - 3') and F10T (5'- FAM – TAT AGC TA TA/Sp18/TA GCT ATA – TAMRA - 3') were dissolved in buffer at 0.25 μM concentration and annealed by heating at 90 °C for 5 min, and then immediately cooling to 0 °C, in the case of F21T, or slowly to rt (F10T) during a period of three hours. The solutions were allowed to stand at 4 °C at least overnight.

The tested compounds 1-4 were dissolved in water with up to a maximum 10% DMSO (50 μM stock solution). Each stock solution was then gradually diluted to reach the final concentrations in the assay from 0 to 10 μM and a fixed 0.2 μM concentration of target sequences F21T or F10T. The buffer used for G4 sequence F21T was 10 mM KCl, 90 mM LiCl and 10 mM Li-cacodylate, pH 7.3. The buffer employed with F10T was 100 mM LiCl and 10 mM Li-cacodylate.

The melting process program consisted of an incubation step at 24 °C for 5 min and then a temperature ramp (1°C/min rate) with fluorescence recordings at every degree up to 95 °C. Subsequently, the microplate with the samples (50 μL) was incubated for 5 min at 96 °C and then gradually cooled by using a temperature ramp of -1 °C/min, with fluorescence measurements at every degree down to 25 °C. The melting and refolding (or renaturing) curves were registered based on the change in emission of the FAM (6-carboxyfluorescein) fluorophore, excited at 492 nm and emitting at 516 nm. Melting temperature values (T_m) were determined from replicate experiments using normalized curves and considering mid-transition T_1/2 temperatures.

2.3.2. Tel22 CD titration spectra

CD experiments were performed as previously reported [29]. Unlabeled Tel22 telomeric quadruplex sequence (5'-A(GGGTTA)₃GGG-3') at 4 μM concentration was folded in the same conditions than those employed for FRET-based melting experiments, using the buffering system than for G4 F21T. Tel22 and ligand 1 solutions at DNA/ligand ratios ranging from 1:1 to 1:5 were freshly prepared before the assay and allowed to stand for 2 hours prior measurements. The temperature used in the titrations was kept at 25 °C during the experiments and CD spectra were registered using a 0.5 cm-path cell, 1 nm band width, and 0.5 nm increments.

2.3.3. Non-competition and competition equilibrium dialysis experiments

Equilibrium dialysis experiments were conducted with unlabeled quadruplex DNA and/or with dsDNA sequences as previously reported [29]. G4 sequence Tel22 is 5’-A(GGGTTA)₃GGG- 3’ and ds17 is 5’-GGG TTA CTA CGA ACT GG- 3’/5’-CCA GTT CGT AGT AAC CC-3’. The buffering system employed contained 10 mM potassium phosphate buffer and 100 mM of KCl (pH = 7.2). DNA solution concentrations were determined at 90 °C by UV-visible spectrophotometry using the λ_max values and extinction coefficients provided by IDT®. Extinction coefficients of TSC ligand 1 were calculated in water and in buffer in the absence and in the presence of 1% (w/v) triton X-100. Competition and non-competition dialysis experiments were carried out using the protocol described by Chaires [36], with minor modifications. 0.2 mL volume dialysis bags (Biotech Regenerated Cellulose (RC) membrane, part number 133192, Spectrum Laboratories, Inc.) were used, with G4 Tel22 or duplex ds17 DNA at a concentration of 75 μM, in monomeric units. The dialysis bags were placed in a beaker containing 200 mL of a solution of compound 1 (2 μM ) in buffer. The beaker was sealed with Parafilm and the solution was allowed to equilibrate at room temperature (22 °C) for 24 h with continuous stirring. Then, the contents of the dialysis bags were transferred into microcentrifuge tubes and triton X-100 was added (1% (w/v) final concentration). After 2 h equilibration at room temperature, the total concentration of the ligand (C_t) was determined by UV-visible spectroscopy. The concentration of free compound 1 (C_f) was also determined spectrophotometrically from the dialysate solution. The concentration of DNA-bound compound (C_b) was then calculated by the equation C_b = C_t - C_f and K_app were determined.

2.3.4. CT viscosity assay

The DNA viscosity titration was conducted as previously reported [29]. For that purpose, Calf thymus CT DNA was used. DNA and compound 1 solutions were freshly prepared in 10 mM NaH₂PO₄/Na₂HPO₄ buffer, pH = 7.2. CT DNA solutions (~0.35 mM, nt) were equilibrated for 20 min at 25 °C in a thermostated bath and 20 flow times were then recorded. Several aliquots (ca. 25–60 μL each) of a solution of ligand 1 (~1.8 mM) were sequentially added. After a 20-min equilibration at 25 °C, flow times were registered. Finally, viscosities (μ) were calculated at each titration point, from the averaged flow times recorded and the viscometer constant. The viscosity graphs were obtained by representing (μ/μ₀)^1/3 versus r, the molar ratio of bound ligand to DNA. The parameter μ₀ represents the DNA solution viscosity in the absence of compound.

2.3.5. Cell culture conditions and MTT viability assay

A panel of five tumor cell lines and two normal cell lines was used to preliminary test the cytotoxicity of compounds 1-4, and the experiments were conducted as previously reported [29]. All human cell lines were acquired from American Type Culture Collection (ATCC): tumor PC-3 and DU145 (prostate), HeLa (cervix), MCF7 (breast) and HT29 (colorectal) and healthy HFF-1 (fibroblast) and RWPE-1 (prostate). The cell catalogue reference and culture media used in the experiments were CRL-1435^TM, RPMI (PC-3), HTB-81™, RPMI (DU145), CCL-2^TM, DMEM (HeLa), HTB-22^TM, DMEM (MCF-7), HTB-38^TM, McCoy's 5a Medium Modified, Catalogue No. 30-2007 (HT29), SCRC-1041^TM, DMEM (HFF-1) and CRL-11609^TM, Keratinocyte-SFM (RWPE-1). Culture media were purchased from Merck: RPMI-1640, DMEM + 10% FBS (Fetal Bovine Serum) + 10% antibiotic (penicillin/streptomycin/amphotericin B). Media and supplements for RWPE-1 were from ThermoFisher (17005075). MCF7 cells were additionally supplemented with insulin (0.001 mg/mL in 500 mL media). Cells were kept at 37°C (5% CO₂) and culture media was renewed three times per week. When the cells reached 80-100% confluence, they were trypsinized, centrifuged and the concentration was adjusted to each experiment. Cells were seeded at a density of ~10,000 cells/well into 24-well plates. The cells were treated with ligands 1-4 or the cytotoxic agents PhenQE8, PhenDC3 or cisplatin at different concentrations (10 nM-100 μM). Compound solutions were freshly prepared using the appropriate culture medium. Triplicate experiments were run. For MTT viability assays, after 72 h incubation, 50 μL of MTT (5 mg/mL) was added and an additional incubation of 4 h was allowed. The medium in each well was then replaced by DMSO and the absorbance measured at 570 nm using an ELISA plate reader. Absorbances were normalized using negative controls (untreated cells, in triplicate, with a maximum 0.9% DMSO content). Dose-response curved were obtained by plotting the percentage of viable cells versus logarithmic compound concentration. IC₅₀ values were determined using Kaleidagragh^TM (3.52) software by adjusting the data to equation 1/(1+10^(m2*(log(m1)-x))); m1=.000003; m2=1

2.3.6. Cell cycle and Annexin-V apoptosis assays

For cell cycle experiments, exponentially growing HeLa cells in DMEM were seeded in 6-well plates, at a density of 2.5×10⁵/well. The cells were treated with compounds 1, 2 or reference compounds with known antitumor activity (PhenDC3, 360A or PhenQE8), under several conditions: 72 h at ½ IC₅₀ concentrations, 24 or 48 h at a 20 μM. Untreated cells were used as controls. After washing with PBS, cells were detached with 0.25% trypsin/0.2% EDTA and centrifuged at 1500 rpm for 5 min at rt. The resulting pellets were mixed with cold 70% ethanol (0.5 mL) and kept at 4°C for 30 min. Ethanol was removed by centrifugation and the pellets were washed with 2 mL PBS + 2% BSA and centrifuged. After discarding the supernatants, a propidium iodide solution (PI/RNAse, FxCycle^TM, Invitrogen^TM ref. F10797, 0.5 mL) was added, mixed thoroughly and incubated for 15 min at room temperature prior to analysis by flow cytometry. 10,000 events were registered for each sample. The obtained results were analyzed with the MacsQuantify 2.13.1 program and represented as histograms.

For the Annexin V apoptosis experiments [37], approximately half of the cells used in combined experiments with cell cycle analysis were harvested by trypsinization, PBS washed and stained with the Annexin-V FITC & PI kit (Invitrogen^TM ref. V 13242) following a standard protocol [38]. 10,000 events were recorded, and the results obtained by flow cytometry were represented as dot plots.

3. Results

3.1. Synthesis and characterization of 1,10-phenanthroline TSCs and derived 4-phenylthiazole analogue.

The synthesis of 2,2'-((1,10-phenanthroline-2,9-diyl)bis(methanelylidene))-bis(hydrazine-1-carbothioamide) (1), and its corresponding PF₆^¯ and BF₄^¯ionic salts, compounds 2 and 3, respectively, was conducted via a straight-forward and fast procedure from precursor dialdehyde, 1,10-phenanthroline-2,9-dicarboxaldehyde. The dialdehyde was first prepared following established protocols, as previously reported [29]. Subsequently, through a condensation reaction with thiosemicarbazide (2 equiv., Scheme 1) in an acidic methanolic solution at reflux for 1h, in the absence of inorganic salts, compound 1 was obtained satisfactorily [34], in neutral form, with 81% yield.

In addition to the neutral compound 1, two derived salts were prepared, the PF₆^¯ (2) and the BF₄^¯(3) ionic salts. This was accomplished to study whether there is an effect on DNA interactions and/or biological activity, dependent on the existence and the nature of a counterion. As a matter of fact, it has been previously shown that altering the nature of the counterion in positively charged compounds can have significant effects, for instance, in cytotoxic activity [39,40,41]. Thus, the synthesis of TSC derivatives 2 and 3 was approached, and an analogous experimental procedure was followed by slightly varying the reaction conditions (Scheme 2). The process included an additional reaction step. After 1 h at reflux, a KPF₆ or NaBF₄ salt was added in stoichiometric excess (2 equiv.) to the reaction mixture, at then cooled to room temperature, isolating the desired ionic derivative. Reactions yields were considerably lower for ionic TSCs, compared to the neutral compound: 46% for 2 (PF₆^¯ salt) and 52%, for 3 (BF₄^¯ salt), respectively.

Finally, as thiosemicarbazones can be used as synthetic precursors to render heterocycles of biological relevance, compound 1 was used as the starting material for the synthesis of a thiazole derivative substituted with a benzene ring, compound 4 (2,9-bis((E)-(2-(4-phenylthiazol-2-yl)hydrazineylidene)methyl)-1,10-phenanthroline). In the absence of x-Ray diffraction data, the isomer E is inferred by analogy to previous results with analogue PhenQE8. From a mechanistic perspective, the compound was obtained in good yield (75%) by reaction of 1 with 2-bromoacetophenone in ethanol at reflux for 24 h (Scheme 2), through a Hantszch-type cyclodehydration reaction [42].

After purification by crystallization, the structural characterization of compounds 1-4 in solution was mainly conducted by ¹H and ¹³C NMR, using (CD₃)₂SO as solvent (in the case of compounds 2 and 3, additionally by ³¹P- and ¹⁹F-NMR), FTIR-spectroscopy, high-resolution MS spectrometry and combustion elemental analysis (Materials and methods and Figures S1–S18).

3.2. DNA interactions with telomeric G4

3.2.1. FRET-based DNA melting experiments

Initial assessment of telomeric G4 thermal stabilization by compounds 1-4 was carried out by FRET-based DNA melting experiments. The sequence F21T [43], 5’-labeled with the fluorophore FAM and 3'-labeled with TAMRA was utilized for that purpose, and the ds-DNA forming sequence F10T [44], labeled in a similar fashion, was used as a comparison for duplex stabilization. Both DNA sequences were titrated with compounds 1-4 at a concentration ranging from 1 μM to 10 μM.

At the reference concentration of 1 μM, compounds 1-4 produced negligible stabilization effects. However, as concentration of compound augmented, all derivatives produced a dose-dependent increase in the T_m of the telomeric G4 sequence. Compounds 1-3 showed a similar behavior, with notable stabilization effects (ΔT_m ~ 7 °C, 9 °C and 6 °C at 5 μM; ΔT_m ~ 12 °C, 15 °C and 10 °C at 10 μM, respectively, Figure 3, top panels A-C) whereas compound 4 hardly stabilized the quadruplex structure (ΔT_m ~ 3 °C at 10 μM; Figure 3, top panel D). All the reported ΔT_m values refer to the effects produced by compounds 1-4 on the telomeric G-quadruplex renaturing process, starting from its denatured state to reach the folded conformation.

The thermal stabilizing effects produced by the TSCs compounds were appreciable both in the DNA melting curves and in the opposite renaturing process, although the latter yielded better quality sigmoidal curves. In general, these results highlight that compounds 1-3 are good stabilizing agents of the telomeric G-quadruplex structure, while the derivative 4 produces only subtle stabilization of the structure. In addition, when parallel experiments were run with dsDNA sequence F10T, none of the compounds showed stabilization of duplex structure (ΔT_m ~ 0 °C) up to the highest concentration tested, 10 μM (Figure 3, bottom panels A-D). These results support, for compounds 1-3, the existence of a certain binding preference and higher stabilization effects towards the quadruplex structure than towards dsDNA. However, the establishment of binding selectivity would require additional experiments performed under competitive conditions.

3.2.2. Circular dichroism

Having detected that compounds 1-3, the TSCs compounds, can stabilize the telomeric quadruplex sequence, we were interested in assessing whether the core heterocyclic structure of these derivatives can induce structural changes on the quadruplex structure. Compound 1 was selected for further DNA-ligand interaction studies, starting with circular dichroism (CD) experiments employing the well-known telomeric sequence Tel22 (5’-A(GGGTTA)₃GGG- 3’). Although CD spectra do not provide precise information of G4 folded structures and their binding complexes, they can be utilized as a qualitative estimation of induced changes in the characteristic G4 folding topologies under a wide range of experimental conditions [45,46,47].

The effects of compound 1 on telomeric quadruplex structure were analyzed by recording the Tel22 CD spectra, folded under K⁺ ionic conditions, using various molar [DNA]/[ligand 1] ratios, from 1:1 to 1:5. The main modifications of Tel22 CD signature are depicted in Figure 4. The intensity of the 295 nm band increased in a dose-dependent fashion, up to the maximum ligand/DNA ratio evaluated, molar ratio 1:5, whereas the intensity of the 270 nm band diminished. Overall, the spectral changes mentioned, including the appearance at ~ 260 nm of a negative signal, indicate that compound 1 stabilizes the telomeric G4 and, most importantly, the TSC derivative produces conformational changes in the G4 structure, increasing its antiparallel character. This observation is in good agreement with previously results obtained with compound PhenQE8 [29].

However, two main differences were detected in the CD experiments, when compared to compound PhenQE8. On the one hand, unlike PhenQE8, spectral changes associated to stabilization of the structure with compound 1 are observed at DNA/ligand molar ratios higher than 2:1. On the other hand, an apparent ICD band (induced circular dichroism) was observed at 350 nm, at 1:5 molar ratio DNA/ligand. This suggests that, besides the expected end-stacking binding mode of compound 1, it is feasible that, when added in excess, the compound may bind the quadruplex structure in the grooves or by partial intercalation.

3.2.3. Non-competition and competition equilibrium dialysis

We were then interested in establishing the apparent binding association constant of compound 1 and telomeric sequence Tel22 [48]. To this end, equilibrium dialysis experiments under non-competitive and competive conditions were carried out. For comparison purposes and to estimate quadruplex/dsDNA selectivity, ds17 sequence [49], was also studied, in both conditions. For the competitive assays, the experimental protocol developed by Chaires [36] was followed, with slight modifications (Materials and methods).

The results obtained with compound 1 with telomeric G4 Tel22 and duplex DNA ds17 are presented in Table 1. Assays were performed at room temperature with ~ 2 μM solutions of 1, equilibrated with 75 μM (in monomeric units, tetrads or base pairs, respectively) of the tested DNA sequence for 24 h, in a buffer containing 100 mM KCl. Afterwads, the surfactanct triton X-100 was added to the solutions contained in the dialysis bags, and recording of UV-visible spectra was used to determine the concentrations of free and DNA-bound ligand for each situation: Tel22, ds17 or competition experiment Tel22-ds17. The concentration of DNA-bound ligand was determined based on at least two different experiments, performed in duplicate or tiplicate. The concentration data was then used to determine the apparent association constants of compound 1 with Tel22, with ds17 or with both structures simultaneously, in the competitive version, using the equation K_app = C_b/(C_f)(S_total-C_b)[32], where C_b is the concentration of ligand bound, C_f is the free ligand concentration and S_total= 75 μM, total DNA concentration, in monomeric units.

The determined K_app values are shown in Table 1. Compound 1 binds the telomeric sequence with an association constant in the order of 10⁵ M^-1, and dsDNA ds17 with diminished affinity, in the order of 10³ M^-1. In both competitive and non-competitive conditions, the compound demonstrated higher binding affinity for the quadruplex sequence than for dsDNA sequence ds17, displaying a certain degree of selectivity. The determined apparent association constants under competitive conditions are in good agreement with those determined under non-competitive conditions. Concerning binding selectivity (determined by the quotient K_app (Tel22)/K_app (ds17) in the competitive assay), compound 1 exhibited preferential binding to the telomeric G4 over dsDNA (ca. 24-fold selectivity). This selectivity is considerably lower than the one previously reported for compound PhenQE8 (90-fold) [29]. Overall, the structural molecular difference that implies the substitution of two nitrogen atoms by two sulfur atoms, with changes in compound net charge, is reflected by a lower binding affinity to both G-quadruplex and dsDNA secondary structures.

3.2.4. CT DNA viscosity titration

As the previous interaction assays revealed that the TSCs derivatives can recognize both quadruplex and dsDNA structures, we set out to study whether the binding mode towards double-stranded DNA occurs either by DNA intercalation, groove and/or electrostatic external binding or both. Viscosity titrations assays are known to provide useful information, mainly in reversible dsDNA-ligand non-covalent binding modes [50] and, especially, they are normally used to differentiate between classical DNA intercalants from groove or external binders. The theory of Cohen and Eisenberg [51], established that the existence (or not) of changes of the viscosity of a DNA solution upon binding to small-molecule ligand could be correlated to the ligand binding mode. Therefore, DNA-ligand viscosity titrations can be employed to produce graphs representing the cubed root of the relative DNA viscosity (η/η_o)^1/3 against the molar ratio of bound ligand to DNA (r). If a linear behavior is observed, the obtained slope value can be broadly associated to the DNA-ligand interaction mode. Groove and electrostatic external binders do not normally effect DNA viscosity and commonly display viscosity slope values ~0.0. On the other hand, classical monointercalating agents, such as ethidium bromide, which increase the viscosity of the DNA solution upon binding, give rise to slope values ~1.0 [50,51]. In practice, the slope values experimentally obtained for a great variety of known minor-groove binders oscilate from -0.3 to 0.2 [52].

Compound 1 was chosen as the representative compound from the series, and viscometric titrations were performed in a thermostated bath (25 ± 0.01 °C) with solutions of Calf Thymus (CT) DNA and gradual additions of small aliquots of 1. After an initial 20 min equilibration and further equilibrations following each subsequent addition, at least twenty flow times were registered in the absence (η_o) or in the presence (η) of ligand. Flow time results were converted to viscosity values, which were plotted as the relative viscosity values, (η/η_o)^1/3 relative to r (ratio [ligand]/[DNA (nt)]), as displayed in Figure 5.

Compound 1 barely affected the CT DNA viscosity, showing a minor and linear effect at the [ligand]/[DNA] tested ratios, with a mean slope value of ca. 0.13, slightly lower than the one determined for compound PhenQE8 in previous studies, ca. 0.24. Similarly to what was observed with PhenQE8 [29], and consisting with their similar structure, these results indicate that compound 1 does not bind CT dsDNA through intercalation and insertion of the 1,10-phenanthroline ring between base pairs, but it interacts with dsDNA as a groove or as an electrostatic external binder.

3.3. Biological activity. MTT, cell cycle and apoptosis assays

Finally, the cytotoxic activity and potential antitumor properties of compounds 1-4 in cultured cells were preliminary evaluated by the MTT viablility assay and, in the case of compounds 1-2, by additional cell cycle and apoptosis assays. A panel of five human tumor cells was selected for this purpose: PC-3 and DU145 (prostate), HeLa (cervix), MCF-7 (breast) and HT29 (colorectal). As comparison and controls for cytotoxicity in non-tumor cells, two additional cell lines were also tested: fibroblasts (HFF-1) and prostate (RWPE-1). The cytotoxic activity of previously reported structural analogue PhenQE8 derivative, the highly selective G-quadruplex ligand and structurally-related compound PhenDC3, and metallo-drug cisplatin was also studied, to preliminary assess the biological activity of compounds 1-4.

The results are shown in Figure 6 and Table 2. With the exception of compound 4, which proved to be non-cytotoxic in all cell lines tested, the rest of compounds, after 72 h, exhibited moderate cytotoxicity at micromolar concentrations. Interestingly, the neutral derivative 1 did not showed cytotoxic activity in HeLa and in HT29 up to 100 μM concentration and, overall, its determined cytotoxicity in cell lines evaluated, with the exception of PC-3, was lower than the one observed for compounds 2 and 3, representing the same core structure, but in the form of ionic salts. As a matter of fact, compounds 2 and 3 were found to be highly cytotoxic in normal RWPE-1 cells, especially compound 3. Compound 1, on the contrary, was the only derivative in the series with a reduced cytotoxicity in a non-tumor cell line, RWPE-1. Overall, compounds 1-3 are less cytotoxic by the MTT assay than PhenQE8 and cisplatin in PC-3, DU145 and HeLa cells, but 2 and 3 showed similar or even improved activity in MCF7 and HT29 cells. When compared to PhenDC3, derivatives 1-3 were more cytotoxic in PC-3 cells, and 2 and 3 in HeLa cells (Figure 6, Table 2). Finally, regarding the cytotoxicity in normal cells, derivatives 1-3 revealed lower cytotoxicity than PhenDC3 and cisplatin in the normal cells tested, especially in HFF-1 and, in particular, compound 1. However, none of the derivatives demonstrated the selectivity previously reported for compound PhenQE8 [29].

The neutral derivative 1 and a representative example of an ionic salt, compound 2, were selected for further preliminary studies of their potential antitumor properties. Thus, cell cycle experiments using flow cytometry were performed first. Tumor HeLa and PC3 cells were treated with compounds 1, 2 or selected antitumor controls (PhenQE8, 360A and/or PhenDC3) at ½ IC₅₀ equivalent concentrations for 72 h, and changes in the cell cycle distribution (SubG0, G0/G1, S or G2/M populations) were quantitated. In HeLa cells, compounds 1 and 2 notably increased the number of cells in SubG0, which suggests DNA fragmentation and apoptosis, while they induced a decrease in the number of cells in the G0/G1 phase (Figure S19). Interestingly, very minor alterations in cell cycle phase distribution were found for the tested compounds and controls in PC-3 cells, under the experimental conditions previously used in the MTT viability assays.

Cell cycle analysis experiments were also carried out in combination with the Annexin V-PI FITC apoptosis assay, in HeLa cells. This time, cells were treated with compounds 1, 2 or antitumor controls for 48 h at a fixed, common concentration of 20 μM, and cell cycle effects and apoptosis induction were analyzed in parallel from the same batch of treated cells, using flow cytometry. As shown in Figure 7 and Figure 8, compounds 1 and 2 were found to possess pro-apoptotic properties in HeLa cells at 20 μM concentration.

Cell cycle experiments corroborated that compounds 1 and 2 produce an important increase in the number of cells in SubG0 (Figure 7, panels D and E), while the Annexin V-PI assay was used to established whether apoptosis was involved in the cell death produced by the TSC derivatives. Dot plots also allowed to distinguish percentage variations in the cell populations in early apoptosis from late-stage apoptosis, in comparison to untreated cells employed as controls (Figure 8). Compounds 1 and 2 were found to produce significant apoptosis after 48 h, with a notable increase (~ 20-25%) in the number of early apoptotic and late apoptotic cells (Figure 8, panels D and E). As a matter of fact, the effect produced by compound 1 and 2, especially the latter, is equivalent to the activity observed with antitumor compound PhenDC3 (Figure 8, panel C). Interestingly, although compounds 1 a 2 were found to be less cytotoxic in MTT assays than analogue PhenQE8, their preliminary biological action by cell cycle and apoptosis assays in HeLa cells (20 μM, 48 h) is slightly superior (Figure 7 and Figure 8, panels B, D and E). This highlights that the potential antitumor properties of the core structure of this family of TSCs derivatives are maintained when the nitrogen atom of the guanidine group is replaced by a sulfur in a thiourea moiety.

4. Discussion

We have carried out the synthesis of four 1,10-phenanthroline derivatives disubstituted at 2,9-positions, as potential DNA ligands endowed with antitumor properties in vitro. Compounds 1-3 are thiosemicarbazone-based compounds (TSCs), either in neutral form or as ionic salts (PF₆¯ and BF₄¯, respectively) and compound 4 is a 4-phenylthiazole analogue, synthetically accessible from derivative 1. To the best of our knowledge, we describe for the first time the DNA binding properties of known compound 1 and these three analogues, providing experimental results that support their potential applications as G-quadruplex DNA ligands and as cytotoxic agents with antitumor properties. Compared to our previous work on the closest structurally related analogue, the PhenQE8 G-Quadruplex ligand, the 1,10-phenanthroline heterocyclic core is substituted with two thiosemicarbazone units (TSCs), replacing two nitrogen atoms of the parent ligand for two sulfur atoms in compounds 1-3. The synthesized derivatives are poly-nitrogenated compounds that lack additional aromatic or heteroaromatic rings, apart from the heterocyclic 1,10-phenantrholine system, rings that are normally present in other well-known G4 ligands and make an important contribution to maximize end stacking interactions with the quadruplex structure. Likewise PhenQE8 [29], X-ray diffraction studies [34] indicate that these molecules can adopt extended conformations, and this may be the reason why they are found to be optimal for quadruplex recognition, besides the fact that they possess only one heteroaromatic ring in their structure.

In this work, the human telomeric sequence Tel22 was chosen as the target G4 structure, as it constitutes one of the most representative examples of biologically relevant DNA quadruplex. In particular, we have focused on the two hybrid structures, which are present in solution when the quadruplex is folded in K⁺ ionic conditions. The study of the in vitro interactions of compounds 1-4 with the telomeric quadruplex has been approached by using several techniques, such as FRET-based DNA melting assays, CD, and competitive and non-competitive equilibrium dialysis. Furthermore, the binding mode of compound 1 with CT DNA has been investigated through viscometric titrations.

These DNA interaction studies have shown that compounds 1-3 can thermally stabilize the telomeric quadruplex DNA to a greater extent than compound 4, the 4-phenythiazole derivative, whereas they produce negligible stabilization effects on double stranded DNA. The degree of stabilization of the telomeric structure was very similar in the three compounds, although 2 exhibited slightly higher ΔTm values. Compound 1, representing the core TSC structure of the family, was selected and studied in further detail. Circular dichroism spectra revealed that compound 1 stabilizes the hybrid structure of the telomeric G4 DNA, at least up to a ratio 1:5 Tel22/ligand and, most importantly, it induces a conformational change similar to that produced by ligand PhenQE8, increasing the antiparallel character of the structure. However, unlike PhenQE8, an apparent ICD band was observed at 350 nm at the 1:5 Tel22/ligand ratio, where the compound absorbs, suggesting than the derivative, besides end-stacking interactions with the telomeric quadruplex, may have other additional binding sites as a groove binder, or even as a partial intercalating agent.

We have determined by equilibrium dialysis experiments that compound 1 interacts quadruplex and dsDNA with different binding affinities, with an apparent association constant to the telomeric quadruplex ~10⁵ M^-1, while its binding affinity to duplex DNA is significantly lower (~5·10³ M^-1). Under competitive conditions, the selectivity of 1 to preferentially bind the telomeric quadruplex Tel22 over ds17 dsDNA was estimated to be ~ 24-fold. This represents a notable decrease in quadruplex versus dsDNA selectivity, compared to ligand PhenQE8. Overall, replacement of two nitrogen atoms in PhenQE8 by two sulfur atoms seems to produce a profound effect on DNA binding affinity, although the derivatives still retain certain degree of selectivity for the quadruplex structure. In addition, viscosity titration assays with CT DNA demonstrated that compound 1 does not behave as a DNA intercalating agent, but it recognizes the double helix by groove or electrostatic external binding, similarly to PhenQE8. Future molecular modelling studies, including docking and molecular dynamics simulations, will be used to shed light on the binding mechanism of these derivatives with the telomeric and other G4 DNA sequences.

The in vitro antitumor properties of compounds 1-4 have been preliminary assessed in cultured cells by using the MTT viability assay, cell cycle analysis and the annexin V-PI apoptosis assay. The TSCs compounds (1-3) have shown modest cytotoxic activity in the cell lines selected for this study (PC-3, DU-145, HeLa, MCF-7 and HT-29) with micromolar IC₅₀ values, whereas compound 4 did not prove active in any of the cell lines tested. Overall, compounds 1-3 exhibited lower cytotoxicities than the G4 ligand PhenQE8 and cisplatin metalodrug in all cell lines tested, but similar or higher than PhenDC3 in the case of PC3 and HeLa cells. The decreased cytotoxicity with respect to PhenQE8 is more pronounced for the neutral derivative which, in contrast, was the only compound which display a modest degree of selectivity in prostate cells (prostate tumor PC3 versus prostate non-tumor RWPE-1, selectivity index (SI) ~ 5). In general, within the TSC compounds, those in the form of ionic salts were found to be more cytotoxic to tumor and normal cells than the neutral derivative.

Finally, compounds 1 and 2 have been shown to alter cell cycle in HeLa cells after 48 h of treatment, with a considerable increase of the number of cells in SubG0 phase, which suggests that the compounds can produce cell death by apoptosis. This have been further confirmed by the Annexin-PI apoptosis assay, which also revealed a significant increase of number of cells in early and late apoptosis stages, more pronounced in the case of compound 2. The observed effect is similar of that produced by antitumor compound PhenDC3, and slightly superior to the activity observed for analogue PhenQE8. Further investigation on the mechanism of action of this family of compounds will be needed to precisely determine the molecular targets of PhenQE8 and TSC analogues that account for their biological activity.

5. Conclusions

In this work, a series of four 2,9-disubstituted 1,10-phenanthroline derivatives has been synthesized, three of them representing novel compounds. The compounds represent thiosemicarbazones (TSCs), the neutral compound 1 and two ionic derivatives (PF₆¯ and BF₄¯ salts, 2 and 3, respectively), and a fourth is a 4-phenylthiazole derivative (4). Their binding interactions with telomeric quadruplex and with dsDNA, especially for derivative 1, have been assessed by several techniques, from which we conclude that the TSCs derivatives retain the ability to bind the telomeric quadruplex with more binding affinity than with dsDNA, as other known 1,10-phenanthroline-based 2,9-disubstituted quadruplex ligands, although they have demonstrated only a moderate selectivity. Compounds 1-3 have proven cytotoxic against various tumor cell lines (PC-3, DU145, HeLa, MCF-7 and HT29), and also in the two normal cell lines tested (HFF-1 and RWPE-1), with only compound 1 showing minor selectivity in prostate cells (SI =5). On the other hand, modification of the structure with incorporation of a thioazole ring in the side chain of these 2,9-substituted heterocyclic systems clearly diminishes its thermal stabilization of the quadruplex DNA structure and, most importantly, it abolishes cytotoxic activity below 100 μM concentration. Although the cytotoxicity of derivatives 1-3 in tumor cells is moderate compared to other well-known G4 ligands, the compounds are able to produce cell death by apoptosis, an effect that compares well to, for example, PhenDC3 or PhenQE8 biological activities.

The poly-nitrogenated nature of the reported compounds and the extended conformations that these molecules are known to adopt may account for their capacity to recognize and bind G4 DNA structures with superior affinities than their double-stranded counterparts, despite the lack of additional aromatic rings besides the 1,10-phenanthroline system. However, this work brings to the fore that replacement of two nitrogen atoms in the related PhenQE8 analogue by two sulfurs in the TSCs compounds results in a decrease of overall binding affinity for DNA secondary structures and a diminished cytotoxicity and selectivity, although it preserves the ability to induce apoptosis in tumor cell lines. Further investigations will be required to assess the scope and applications of the reported derivatives as potential antitumor agents.