Richard J.
Mitchell
a,
Sarah M.
Kriger
d,
Alexander D.
Fenton
a,
Dmytro
Havrylyuk
a,
Ankit
Pandeya
a,
Yang
Sun
a,
Tami
Smith
b,
Jason E.
DeRouchey
a,
Jason M.
Unrine
b,
Viral
Oza
c,
Jessica S.
Blackburn
c,
Yinan
Wei
a,
David K.
Heidary
*d and
Edith C.
Glazer
*d
aDepartment of Chemistry, University of Kentucky, 505 Rose St., Lexington, KY 40506, USA. E-mail: ec.glazer@uky.edu; david.heidary@gmail.com
bDepartment of Plant and Soil Sciences, University of Kentucky, 1100 S. Limestone St, Lexington, KY 40546, USA
cDepartment of Molecular and Cellular Biochemistry, University of Kentucky, 741 S. Limestone St., Lexington, KY 40536, USA
dDepartment of Chemistry, North Carolina State University, 2620 Yarbrough DriveRaleigh, NC, 27695, USA. E-mail: eglazer@ncsu.edu
First published on 27th February 2023
Ruthenium complexes are often investigated as potential replacements for platinum-based chemotherapeutics in hopes of identifying systems with improved tolerability in vivo and reduced susceptibility to cellular resistance mechanisms. Inspired by phenanthriplatin, a non-traditional platinum agent that contains only one labile ligand, monofunctional ruthenium polypyridyl agents have been developed, but until now, few demonstrated promising anticancer activity. Here we introduce a potent new scaffold, based on [Ru(tpy)(dip)Cl]Cl (tpy = 2,2′:6′,2′′-terpyridine and dip = 4,7-diphenyl-1,10-phenanthroline) in pursuit of effective Ru(II)-based monofunctional agents. Notably, the extension of the terpyridine at the 4′ position with an aromatic ring resulted in a molecule that was cytotoxic in several cancer cell lines with sub-micromolar IC50 values, induced ribosome biogenesis stress, and exhibited minimal zebrafish embryo toxicity. This study demonstrates the successful design of a Ru(II) agent that mimics many of the biological effects and phenotypes seen with phenanthriplatin, despite numerous differences in both the ligands and metal center structure.
In contrast to cisplatin, carboplatin, and oxaliplatin, which can make two bonds to a biomolecule, pyriplatin and phenanthriplatin defied traditional structure–activity relationships (SAR) by exerting their anticancer properties as monoadduct generators.6–9 Recently, a variety of studies have demonstrated that phenanthriplatin acts via multiple biological mechanisms, including topoisomerase inhibition,10 damage to nuclear DNA,9 and, like oxaliplatin, ribosome biogenesis stress.3,4 As with oxaliplatin, the anticancer profile of phenanthriplatin is distinct from traditional platinum agents, suggesting potential application in tumor types that are refractory to the other drugs. It is also possible that monoadduct formers may evade resistance mechanisms, such as DNA repair, that are effective against biadduct formers. However, all platinum agents suffer from other mechanisms of resistance in cancer cells, such as thiol-mediated inactivation and decreased cellular accumulation.1,7 For this reason, other transition metals are intrinsically appealing if they can mimic the biologically efficacious molecular interactions of non-traditional Pt(II) agents, but overcome these limitations.
Polypyridyl ruthenium complexes with one labile ligand are potential replacements for platinum-based monoadduct-generators.11,12 This might seem counterintuitive, due to structural differences, as Ru(II) creates octahedral complexes and Pt(II) is square planar, and the platin agents do not contain polypyridyl ligands. However, both metals are relatively “soft” and exhibit aquation kinetics that typically occur over several hours, a similar time-scale to cell division processes.13 Ru(II) has a lower affinity for free thiols, making its complexes less susceptible to inactivation and efflux.13,14 Despite this fact, currently reported Ru(II) polypyridyl monofunctional agents have poor activity compared to phenanthriplatin.11,12 It is unclear if this is due to a failure to mimic appropriate mechanism of action (MoA), or other chemical or biological factors.
We questioned if it was possible to make a Ru(II) complex that recapitulated phenanthriplatin's activity, which is attributed to multiple different MoAs, but with an emphasis on ribosome biogenesis. This was anticipated to be a challenge, as recently, the DeRose group demonstrated very stringent structural requirements for Pt(II) compounds capable of disrupting ribosome biogenesis. To this date, there are only six molecules that have this effect: oxaliplatin, three oxaliplatin analogues, an oxaliplatin prodrug, and phenanthriplatin.4,15,16 However, we were encouraged by the fact that some Ru(II) complexes are able to localize to the nucleolus,17 as oxaliplatin and phenanthriplatin do, and other ruthenium complexes have been shown to alter translation.18,19 Moreover, specific Ru(II) complexes have also shown topoisomerase inhibition, and can either inhibit the catalytic site of the enzyme20 or generate cleavage complexes.21–23 These studies demonstrate individual Ru(II) compounds can exhibit the various biological activities of phenanthriplatin, and we hypothesized that this may be done in a single molecule if key structural features, such as polyaromatic lipophilic groups are incorporated into the system. Here we report the successful creation of an octahedral Ru(II) “mimic” of phenanthriplatin, which is a new addition to the select group of compounds that induce nucleolar stress.
Compound | IC50 (St. dev.) [μM] | Zebrafish viability (%) [20 μM] | ||||||
---|---|---|---|---|---|---|---|---|
HL-60 | DU145 | MIA PaCa2 | A549 | HEL299 | SI (HEL299:A549) | DU145 spheroids | ||
a Average values from triplicate measurements. b The percentage of viable embryos is representative of n = 6. c Reported values are following exposure to indigo (450 nm) light. d Values taken from prior work. | ||||||||
1 | 3.33 (0.04) | 1.50 (0.10) | 2.6 (0.2) | 5.73 (0.04) | 28.19 (1.85) | 4.9 | 27.33 (2.14) | 100 |
2 | 0.4 (0.1) | 0.17 (0.02) | 0.23 (0.06) | 0.57 (0.02) | 7.56 (0.37) | 13.3 | 4.91 (0.84) | 100 |
3 | 1.6 (0.2)d | 7.4 (1.0) | 26.5 (3.5) | 1.1 (0.3)d | 16.72 (4.27) | 15.2 | n.d. | n.d. |
Cisplatin | 3.1 (0.2)d | 0.62 (0.06) | 5.9 (0.4) | 3.4 (0.6) | 3.05 (0.44) | 0.9 | 5.28 (0.36) | 100 |
Phenanthriplatin | 1.8 (0.4) | 0.14 (0.01) | 0.15 (0.01) | 0.19 (0.01) | 2.63 (0.14) | 13.8 | 0.84 (0.06) | 86 |
In order to determine the selectivity for cancer cells compared to healthy cells, healthy lung fibroblasts (the HEL299 cell line) were treated with each compound and a selectivity index (SI) was determined for the ratio of IC50 values for HEL299:A549. The SI for cisplatin was 0.9, consistent with previous reports.26 Compound 1 had an SI of 4.9, while phenanthriplatin and 2 had an SI value of 13.8 and 15.2 respectively, displaying the capacity for this scaffold to selectively target cancerous cells over healthy cells. Zebrafish embryos were then used as a screening tool to assess tolerability for small molecules in vivo.35 We found that zebrafish embryos tolerated 20 μM treatment with compounds 1 and 2 (Table 1), demonstrating that these compounds exhibit selective toxic effects.
Different cancer types exhibit varying sensitivities to small molecules. This information can be utilized to understand mechanistic underpinnings of novel molecules.36 To estimate the selectivity of compounds for the various cancer cell lines, the difference in individual pIC50 values versus the average pIC50 was calculated (eqn (S1) and (S2), ESI†). From these data, the sensitivity of each cell line to each of the compounds was visualized, with negative and positive values corresponding to cell line resistance and sensitivity, respectively (Fig. S8, ESI†). The profiles for 1 and 2 were similar to phenanthriplatin, and distinct from cisplatin and 3.
This intriguing result motivated the submission of 1 and 2 to the National Cancer Institute (NCI) for the 60-cell line analysis at a single dose for percent growth inhibition (Fig. S9 and S10, ESI†), and five-dose response (Fig. S11 and S12, ESI†) to correlate the compounds to the Standard Agents GI50 library of small molecules with known mechanisms of action.36 The COMPARE analysis, which correlates specific cell line response profiles between a compound and the Standard Agents GI50 library,36 was used for 2. The resistance-sensitivity data for compounds with high correlations and distinct mechanisms of action were plotted in a heat map (Fig. 1). The COMPARE analysis identified chromomycin A3 as having the highest correlation to 2, with a Pearson correlation coefficient of 0.85 (Table S1, ESI†). Chromomycin A3 binds to DNA with GC-rich regions,37 with a preference for nucleolus organizer regions (NORs) of DNA,38 which have been implicated in nucleolar stress.39 Actinomycin D exhibited a correlation of 0.69; this compound was initially described as a DNA binding agent that interfered with elongating RNA strands,40,41 resulting in inhibition of transcription, but it is now established that actinomycin D is also a ribosome biogenesis stress inducing agent.42 There was also high correlations between 2 and microtubule-targeting molecules paclitaxel (0.71) and vinblastine (0.64). Very few coordination complexes have been reported to interact with microtubules, but hyperstabilization was shown with [Ru(dip)3]2+ in one study,43 and an organometallic osmium compound had a similar effect due to coordination to the sulfhydryl groups of microtubules.44
Notably, the correlation between 2 and DNA alkylating agents such as cisplatin and carboplatin was poor (0.12–(−0.16)), suggesting the primary mechanism of action is not reliant upon this kind of DNA damage (Table S1, ESI†). There was a positive correlation between the topoisomerase inhibitors etoposide and doxorubicin and 2 (0.55–0.39), although it was lower than that of transcription inhibitors. However, doxorubicin also induces hallmarks of nucleolar stress. Oxaliplatin and 1 were also plotted in the heatmap, but only for qualitative comparison, since they are not in the standard agents library. The five-dose data for 1 and oxaliplatin was collected independently for addition to Fig. 1. Oxaliplatin appeared to have a moderate correlation, however, non-small cell lung cancer cell lines exhibited a resistant profile for oxaliplatin but a sensitive profile for 2. Phenanthriplatin was not included in the heatmap because the 60-cell line data is not available on the public NCI DTP database. However, the data available from previous reports allowed for some comparisons.9 Non-small cell lung cancers and central nervous system (CNS) cancers were both sensitive to phenanthriplatin and 2.9 The response across ovarian cancers was comparable, suggesting that phenanthriplatin and 2 exhibit similar behaviors across different cell lines of the same cancer type. The response in renal and breast cancers, in contrast, differed between the two compounds. Renal cancers were sensitive to treatment with phenanthriplatin but breast cancer response was average.9 Renal cancers were resistant to 2 and breast cancers were generally sensitive. Since phenanthriplatin has been identified as having multiple MoAs, and 2 exhibits comparable sensitivity-resistance profiles to phenanthriplatin, we decided to evaluate 2 for each of the plausible mechanisms highlighted by the NCI analysis: DNA damage, topoisomerase inhibition, microtubule dynamics, and ribosome biogenesis stress.
Cellular uptake of the different metal complexes in E. coli was quantified by Inductively Coupled Plasma Mass Spectrometry (ICP-MS). An uptake of 1.24 × 10−15 and 1.67 × 10−15 moles of Ru per cell was found for 1 and 2, compared to 8.76 × 10−16 moles of Pt (Table S2 and Fig. S15, ESI†). Compound 3 was utilized as a Ru(II) positive control, and approximately 7.7% of the compound was observed in bacteria, similar to earlier reports.46 Next, cellular DNA was isolated and analyzed for metal content. The quantity of metal in the total genomic DNA was determined by ICP-MS, and the nucleotide to metal center ratio was calculated by dividing the number of nucleotides in the total genomic DNA present in each sample by the metal atoms. Phenanthriplatin produced the lowest nucleotide (nt) to metal center (mc) ratio, followed by cisplatin and 1, and then 2 (Table S3, ESI†). These results indicate that both compounds 1 and 2 irreversibly bind DNA. However, despite the fact that compound 2 and phenanthriplatin have the same MIC value in E. coli, the Ru(II) complex formed 10-fold fewer DNA adducts. This is consistent with an interpretation that an indirect mechanism of DNA damage, such as topoisomerase inhibition, or another mechanism, such as ribosome biogenesis stress, plays a greater role in the functional activity of 2. A ratio of 1 DNA-Ru(II) adduct to 106 atoms of Ru taken up into each E. coli cell also suggests a mechanism of cell death independent of DNA damage.
The cytological profile only provides supporting evidence of a MoA, as it reports solely on biological processes in bacteria. Compounds with MoAs that are specific to eukaryotic cells do not produce cytological profiles, as is the case with microtubule disrupting agents and members of the camptothecin family, which are Type I topoisomerase inhibitors that are ineffective in bacteria. The profile of the Ru(II) complex observed in E. coli does not eliminate the possibility of such mechanisms of action in mammalian cells, but it indicates that 2 has other MoA that are effective in bacteria that are also observed with phenanthriplatin and Topo II inhibitors.
Compounds 1 and 2 were also evaluated as DNA damaging agents via gel electrophoresis. The DNA damage EC50 for phenanthriplatin, 1, and 2 was 8.8 μM, 21.9 μM, and 12.2 μM, respectively (Table S4, Fig. S18 and S19, ESI†). DNA electrophoresis can indicate general DNA damaging capacity in an isolated system, but the shifts in DNA mobility are heavily dependent on distorting adducts being generated, and these studies are limited in their ability to differentiate between non-covalent and covalent adducts. To reconcile the results of the gel electrophoresis experiments with DNA, the reactivity of 1 and 2 for ligand exchange was examined to determine if a covalent or non-covalent MoA was more plausible. When 1 and 2 were dissolved in phosphate buffer, the labile chloride ligands could be exchanged with water. Compound 1 was converted to the aqua, while 2 exhibited some conversion, but some also precipitated out of solution (Fig. S20, ESI†). Both 1 and 2 did not precipitate in the presence of BSA or Opti-MEM, and conversion to the aqua was significantly slowed. The rates of ligand exchange reactions with DNA also differed. Compound 1 had a half-life (t1/2) of 2.5 hours, while 2 exhibited biphasic kinetics, with an average t1/2 of 38.6 hours in the presence of calf thymus (CT) DNA and BSA (Fig. S21, ESI;† compounds 1 and 2 were dosed at 100 μM, and a concentration of 0.66 μg μL−1 for each CT-DNA and BSA was used).
High performance liquid chromatography (HPLC) analysis was conducted to determine if 2 could covalently modify a DNA base over 16 h. First, 2 was evaluated in water alone. A new peak in the chromatogram appeared at 20.8 minutes, and increased in area over time, while the parent peak at 23.3 minutes decreased, indicating the conversion to the aqua complex (Fig. 3A). In the presence of deoxyguanosine (dG), the aqua peak increased initially, and then decreased as a new peak at 20.6 minutes emerged (Fig. 3C). These results indicate that 2 covalently modifies dG within 10 h, either directly or through an aqua intermediate. When phenanthriplatin was incubated at room temperature in water, the complex was stable, with only a small percentage exhibiting exchange of the chloride ligand (Fig. S22A, ESI†). This could be observed from the slight decrease in the peak at 9.5 minutes and a new peak observed at 8 minutes. However, in the presence of 5 mM dG, phenanthriplatin underwent full conversion to the dG-adduct in 16 h (Fig. S22B, ESI†). These data suggest phenanthriplatin does not require an aqua intermediate to react with dG, which is consistent with previous reports.50 When comparing the reaction kinetics of 2 and phenanthriplatin, the dG-conjugates were formed at similar times, but through different processes, as 2 converts to the aqua complex first, and full reaction with dG subsequently takes an additional 8 hours.
Although compound 2 can covalently modify DNA, it was unknown if this was a cause of cell death. Nucleotide excision repair (NER) is the DNA damage response mechanism utilized to repair distorting adducts, such as those seen with cisplatin. Xeroderma Pigmentosum Group A (XPA) is an initiator protein involved in global genome and transcription-coupled nucleotide excision repair mechanisms. When XPA is knocked out, distorting adducts are not repaired, resulting in cellular sensitization to DNA damaging agents.53,54 Thus, XPA knock out HeLa cells were treated with phenanthriplatin, cisplatin, 1, and 2. The XPA knockout was sensitized 4-fold versus the parental cell line to cisplatin (Fig. 3D). A 2-fold sensitization occurred for 1 (Fig. S24, ESI†), and no sensitization was observed for phenanthriplatin (Fig. S24, ESI†) or 2 (Fig. 3E), indicating that DNA adducts resolved by NER were not the cause of cell death for either compound. This result was consistent with previous reports for phenanthriplatin.9
Another commonly utilized reporter for DNA damage is phosphorylation of serine 139 on Histone 2A; phospho-H2AX (p-γH2AX) initiates DNA damage repair (DDR). Immunoblotting for p-γH2AX thus is a common approach to probe small molecules that covalently modify DNA and initiate DDR. To further investigate DNA damage, the lysates from A549 cells treated with cisplatin (5 μM), phenanthriplatin (0.5 μM), or 2 (0.5 μM) were probed for p-γH2AX at 4 and 24 h (Fig S25, ESI†). The 24 h treatment with phenanthriplatin yielded a significant p-γH2AX response. Compound 2-treated samples, however, did not exhibit p-γH2AX, indicating that another mechanism of action is involved.
In order to determine the topoisomerase inhibiting capacity of 2in vitro, the in vitro complex of enzyme (ICE) assay was conducted. If topoisomerase is unable to re-ligate the DNA, it stalls on the DNA as a cleavage complex (ccDNA). In the ICE assay, cells are treated acutely with high concentrations of the compound of interest, lysed, the DNA is extracted, and stalled topoisomerase is detected by immunoblotting. In contrast to phenanthriplatin, compound 2 did not generate ccDNA at 50 μM or 100 μM, indicating that 2 is not inhibiting topoisomerase function via cleavage complex generation.
To reconcile the results for the two topoisomerase assays, the capacity for 2 to directly inhibit topoisomerase IIα (TopIIα) was evaluated using a similar computational approach as described above, with Mitoxantrone docked into TopIIα as a control.60 Phenanthriplatin and 2 were both predicted to insert into the DNA-binding channel of TopIIα. A single binding mode was observed for both metal complexes, despite the entirety of the channel being included in the active site definition (Fig. S28 and S29, ESI†). However, compound 2 had fewer interactions with the binding channel than phenanthriplatin, suggesting that 2 has a lower affinity for the enzyme or is unable to compete with biological substrates. This data, in combination with the experimental results for phenanthriplatin and 2 indicates that the Ru(II) complex is unlikely to bind to and inhibit TopoIIα. It also highlights the limitations of an isolated system, as shown by the differences in the ICE assay and gyrase data.
Inspired by this, we used immunofluorescent microscopy to evaluate NPM localization. Concentrated NPM1 foci were observed in the localized phenotype of untreated cells (Fig. 5A and H), while the delocalized NPM1 phenotype was observed with treatment with actinomycin D (Fig. 5B), a positive control,42 consistent with prior reports.4,15,16 Analysis of the images revealed sharp and abundant spikes in fluorescent intensity across the diameter of the nucleus in the untreated population (Fig. 5J). By contrast, the translocated, diffuse NPM1 phenotype displayed overall reduced emission intensity and fewer spikes. Furthermore, the appearance of rounded nucleolar morphology, with NPM1 at the periphery in “cap-like” structures3,4 was reflected with a concave pattern for the pixel intensity across the nucleolus (Fig. 5I and J). Treatment with either phenanthriplatin or 2 induced the translocated NPM1 phenotype in 49% and 36% of the cell population respectively (Fig. 5D, F, and G). Moreover, compound 2 induced the translocated NPM1 distribution at the same concentration as phenanthriplatin (0.5 μM), very close to the IC50 values for cytotoxicity, supporting the implication that induction of ribosome biogenesis stress is correlated to cell death. Smaller changes in NPM1 distribution were observed in cells treated with cisplatin and 1 (Fig. 5C and E). The difference in phenotype between 1 and 2 suggest NPM1 translocation is not a general property of monofunctional Ru(II) compounds. The biological effect is likely a result of a specific combination of molecular features, including sterics and hydrophobicity, as is the case with monofunctional platinum compounds. Indeed, out of five monofunctional compounds, DeRose showed the NPM1 redistribution phenotype was unique to phenanthriplatin,4 highlighting the sensitivity of the cellular response to compound structure.
Using a semi-rational approach, we designed and characterized a potent Ru(II) compound that acts as a molecular mimic of oxaliplatin and phenanthriplatin. The ability of this Ru(II) complex to induce the hallmarks of ribosome biogenesis stress is noteworthy. To the best of our knowledge, there are only six inorganic compounds that have this effect, five of which are related to oxaliplatin, and phenanthriplatin. There are a few other Ru compounds that have been reported to interact with the nucleolus,19 but none have been reported to cause nucleolar stress. The incorporation of a phenyl-modified terpyridine and a dip ligand were key features in the design of the Ru(II) complex, and given the simplicity of this octahedral scaffold, there are many opportunities for additional synthetic modulation. Indeed, compound 2 opens the door to rational design of improved derivatives, and the use of phenotypic assessment streamlines planning of further biological evaluation. In sum, we believe that Ru(II) complex 2 provides a readily modifiable chemical scaffold for the development of improved ribosome biogenesis stress inducers.
What remains to be determined is if there is any particular significance to the ability of compound 2 to form monoadducts, or if the chemical structure, rather than the chemical reactivity at a single site, is the key to the biological effects observed. We also note the chemical or biological features of phenanthriplatin that prevented its advance into the clinic also unpublished, and are possibly unknown. This information is needed to determine whether Ru(II) analogues have a chance at a more successful ultimate outcome, and to allow for logically directed optimization of properties.
Footnote |
† Electronic supplementary information (ESI) available: Experimental approaches, additional figures and tables. See DOI: https://doi.org/10.1039/d2cb00247g |
This journal is © The Royal Society of Chemistry 2023 |