Mun Juinn
Chow
*a,
Mohammad
Alfiean
b,
Giorgia
Pastorin
cd,
Christian
Gaiddon
ef and
Wee Han
Ang
*ad
aDepartment of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore. E-mail: a0037979@u.nus.edu
bSchool of Chemistry and Life Sciences, Nanyang Polytechnic, Singapore
cDepartment of Pharmacy, National University of Singapore, 3 Science Drive 3, 117543 Singapore
dNUS Graduate School for Integrative Sciences and Engineering, Singapore. E-mail: chmawh@nus.edu.sg; Tel: +65 6516 5131
eUniversité de Strasbourg, Strasbourg, France
fU1113 INSERM, 3 Avenue Molière, Strasbourg 67200, France
First published on 28th February 2017
Multidrug resistance is a major impediment to chemotherapy and limits the efficacies of conventional anticancer drugs. A strategy to bypass multidrug resistance is to develop new drug candidates capable of inducing apoptosis-independent programmed cell death. However, cellular pathways implicated in alternative programmed cell death are not well-elucidated and multifactorial, making a target-based discovery approach a challenge. Here, we show that a coordination-directed three-component assembly and phenotypic screening strategy could be employed as a viable alternative for the identification of apoptosis-independent organoruthenium anticancer agents. Through an on-plate synthesis and screening of 195 organoruthenium complexes against apoptosis-sensitive and -resistant cancers, we identified two apoptosis-independent hits. Subsequent validation of the two hits showed a lack of induction of apoptotic biomarkers, a caspase-independent activity and an equal efficacy in both apoptosis-sensitive and -resistant colorectal cancers. This validated their apoptosis-independent modes-of-action, paving the way as potential candidates for the treatment of highly-refractory cancer phenotypes.
In recent years, several other types of programmed cell deaths have been identified involving cellular pathways and morphological changes distinct to that of apoptosis.12 Concomitantly, apoptosis-independent anticancer drug candidates have been explored as treatment options for apoptosis-resistant cancers.13–16 Several metallodrugs have shown potential for such applications. Cu(I) and Cu(II) complexes have been known to induce a form of programmed necrosis known as paraptosis, characterized by mitochondria swelling and cytoplasmic vacuolation.17–20 Two Re(V) oxo complexes were recently reported to induce necroptosis, a distinct form of programmed necrosis mediated by RIP1/RIP3 kinase complex.21 Several Cu, Fe, Au and Pt complexes have also exhibited the ability to induce autophagic cell death.22–25 A number of Ru(II) complexes induced non-apoptotic cell death, although the exact cell death types have yet to be characterized.26–29 For example, RAS-1H and RAS-1T could bypass apoptosis-resistance in an in vitro system (Fig. 1).28
Fig. 1 A phenotypic screening approach to identify apoptosis-independent Ru-arene Schiff-base complexes. (Left panel) Previous strategy with generic cytotoxicity screening approach; serendipitous discovery of apoptosis-independent compounds only after extensive mode-of-action studies.28,34,35 (Right panel) Phenotypic screening approach used for the current study; apoptosis-independent ‘hits’ could be identified right from the screening phase. Few steps required after to confirm validity of ‘hits’. |
In general, discovery of such compounds is sparse compared to the vast majority of apoptotic agents and few have been designed specifically for that goal. This is because, unlike the apoptosis pathway, which is extensively studied and contain many known druggable targets, cellular pathways implicated in other programmed cell deaths are less well-defined with fewer druggable targets. Thus, a target-based approach focused on the discovery of apoptosis-independent drugs is not entirely feasible. Furthermore, the efficiency of target-based drug discovery has been called into question, owing to the high attrition rates of candidates discovered through such an approach and the multifactorial nature of MDR.30–32 The low efficacy of these drug candidates has been attributed to false-positives in target validation and a failure to take into account redundancy and compensatory crosstalk in complex cellular signaling pathways.33
An alternative approach to target-based drug discovery would be the use of a phenotypic screening strategy for fast identification of hits. However, a generic phenotypic screening approach using only cytotoxicity as a standard-of-measure seldom yields drug candidates with novel mode-of-action.32 This means that identification of apoptosis-independent hits in this manner would rely on serendipitous discovery after extensive testing and mode-of-action studies (Fig. 1).28 Hence, a methodology combining the efficiency of the phenotypic screening strategy and the specificity of the target-based approach, coupled with the availability of a large compound library, would be required to accurately identify apoptosis-independent drug candidates in the early phases of the drug discovery process.
In the present study, we utilize a phenotypic screening strategy that would allow us to identify apoptosis-independent candidates right from the screening phase (Fig. 1). We coupled the use of (i) a previously established Coordination-Directed Three-Component Assembly (C3A) protocol as a convenient way to prepare a new library of 195 Ru-arene Schiff-base (RAS) complexes for screening,34 and (ii) the screening against two colorectal cancer cell lines with different sensitivity towards apoptosis i.e. apoptosis-sensitive HCT116 colorectal carcinoma and apoptosis-resistant TC7.36 We investigated the structural factors involved in self-assembly of RAS complexes and shed light on the structural-dependence of the C3A assembly process. Further analysis and comparison of the differential killing in both cell lines allowed for the identification of apoptosis-independent hits, which were validated by investigating their efficacy in resistant cell lines, caspase dependence and induction of apoptotic biomarkers. Although, similar phenotypic screening strategies exploiting differential cytotoxicity in paired cell lines for hit identification are known,37–39 this is the first one (to the best of our knowledge) developed specifically for the identification of apoptosis-independent anticancer drug candidates.
Briefly, stock solutions containing RAx (10 mM), PAy (40 mM) and ADz (40 mM) were prepared in DMSO and added in a predetermined order to a 96-well plate containing equal volume of ddH2O, such that each well would contain one final RAS complex (Fig. 2c). The reaction mixtures were incubated with shaking at rt for 36 h before the crude products were analyzed by ESI-MS and the assembly yields ascertained using RP-HPLC. In general, 166 out of 195 (85%) RAS complexes assembled via C3A had high assembly yield/purity >70% with a significant proportion achieving >90% yields (Fig. 3a; ESI Table S1†). Only 9 (5%) RAS complexes were unsuccessful with assembly yields <50%.
The success of C3A assembly was influenced by both steric and electronic factors. First, components that contributed to increased steric bulk around Ru resulted in lower assembly yields. For instance, RAS complexes synthesized using more sterically-encumbered RA5 had significantly lower yields compared to less sterically hindered RA3 or RA4 (Fig. 3b). Similarly, the bulkier PA2 and PA4 had lower assembly yields compared to less bulky PA1 (Fig. 3c). Moreover, structural isomers with functional moieties positioned in closer proximity to Ru directly increased steric bulk and resulted in lower assembly yield. For example, less RAS assembly was obtained from PA4 compared to isomeric PA5 despite other conditions being maintained (Fig. 3c and 4a). Similarly, meta-aniline derivatives (ADm), which had functional groups positioned nearer to the Ru, also resulted in lower assembly yields compared to their para-analogues (ADp) (Fig. 3d and 4b). Furthermore, steric bulk that hindered imine formation during ligand assembly, as is the case for PA2 (Fig. 3c and 4c), also affected reaction yields.
Beside the dominant steric factors, a secondary electronic effect also affected assembly yield. This was particularly evident in the case of an electron-donating R′ substituent on the AD component. R′ groups at the meta-positions resulted in steric hindrance that led to lower assembly while R′ groups at the para-positions had little influence on C3A. At the meta-position however, electron-donating R′ groups (AD6m, AD7m) consistently reduced assembly yields of RAS complexes (Fig. 3d) while electron-withdrawing R′ groups (AD2m, AD3m) improved yields marginally (Fig. 3e). The lower assembly levels by AD6m and AD7m were presumably due to competition of the substituents toward Ru-binding, which directly interfered with C3A. This was earlier observed for substituents comprising guanidine and thiol moieties which directly bound Ru. In the para-position (AD6p and AD7p), the substituents would be in conjugation with the aniline ring systems which reduced their nucleophilicity by mesomeric effect. In the meta-position however, the nucleophilicity of the substituents remained high and could compete for Ru-binding.
Taken together, C3A proved to be a highly reliable and efficient methodology for the synthesis of RAS complexes with highly predictable levels of RAS assembly. Complexes formed from the most favorable components gave the highest yield (e.g.437p, 98%) while those formed from the most steric encumbered components gave the lowest yield (e.g.547m, <10%).
First, 166 RAS complexes with assembly yield greater than 70% were selected for screening in HCT116. The RAS library was synthesized on-plate (96-well) without any further work-up or purification and exposed directly to cells at a fixed concentration of 25 μM. Cell viabilities were determined after 48 h (ESI Fig. S3 and Table S2†) and analyses were consistent with previous studies that showed a correlation between hydrophobicity and cytotoxicity (Fig. 5), despite the fact that the complexes were intrinsically hydrophilic as determined by logpow measurements.34,35 Variation in the position of the functional groups did not have any observable effect on cytotoxicities since RAS complexes prepared from PA4 and PA5 showed similar inhibition of cancer cell viabilities (Fig. 5c). The same was observed when comparing ADm and their ADp counterparts (Fig. 5d).
Next, efficacious RAS complexes (107 out of 166) were evaluated against TC7 under the same experimental conditions (ESI Fig. S4 and Table S3†). The differences in cell viabilities between TC7 and HCT116 treatment were determined and 22 initial hits with a difference <10% were chosen to ascertain dose-dependencies (Fig. 6a). Therefore, the estimated IC50 (eIC50) values of each initial hit were derived from their cell viability curves in the respective cell lines (ESI Fig. S5†) and resistance factors (RF) were calculated by taking the ratio of eIC50 in both cell lines (Fig. 6b). Two RAS complexes, 532m and 532p, with the lowest RF values were selected as lead candidates for apoptosis-independent compounds (Fig. 6c and d).
Preparative-scale syntheses of 532m and 532p were performed following similar reaction conditions to the earlier on-plate format (Fig. 2a). Crude products were purified using silica gel column chromatography and characterized by 1H NMR and ESI-MS. The 1H NMR spectra of 532m and 532p revealed resonances typical of RAS complexes (ESI Fig. S6 and S10†).34 The singlet resonance peaks corresponding to the imine proton at ca. 9.3–9.4 indicated chelation of the imino-pyridine ligands to Ru. ESI-MS spectra also showed the typical [M]+ molecular ion peaks and with isotopic distribution consistent with calculations (ESI Fig. S9 and S13†). The purity of the complexes was determined to be greater than 95% using elemental analyses, after accounting for residual solvent. UV-vis stability studies showed that the complexes were generally stable to ligand substitution in aqueous solution and that aquation of the Cl ligand could be suppressed in 0.9% w/v NaCl solution (ESI Fig. S14†).
The cytotoxic activities of 532m and 532p were determined in a panel of drug-sensitive (AGS, HCT116) and -resistant (HCT116 p53-/-, TC7) cell lines (Table 1, ESI Fig. S15†). Both compounds showed low micromolar IC50 values in AGS gastric cancer cells and were several times more cytotoxic than cisplatin, commonly used in gastric cancer treatment. In wild type HCT116 colorectal cancer cells, 532m and 532p had modest IC50 values and were less cytotoxic than oxaliplatin and 5-fluorouracil, which are clinical drugs used in the adjuvant treatment of colorectal cancers. However, 532m and 532p were more efficacious against p53-null HCT116 p53-/- and apoptosis-resistant TC7 cells. In particular, 532m was marginally more cytotoxic and 2 times more cytotoxic in HCT116 p53-/- cells than oxaliplatin and 5-fluorouracil, respectively. Against TC7 cells, the differences were even more evident, with 532m being 3 times and 68 times more active than oxaliplatin and 5-fluorouracil, respectively. The calculated RF values were also particularly revealing about the p53-independent and apoptosis-independent nature of these RAS complexes. For instance, when considering the p53-dependence variable in wild type HCT116 and HCT116 p53-/- cells, 532m and 532p exhibited RF of 0.9 and 0.8, respectively, compared to p53-dependent oxaliplatin and 5-fluorouracil with RF of 4.6 and 2.3. Similarly, taking into account the apoptosis-dependence when comparing wild type HCT116 and TC7 cells, 532m and 532p were equally efficacious against HCT116 and TC7 with RFs of 0.8, which was in good agreement with values obtained during initial screening. In contrast, apoptosis-inducing oxaliplatin and 5-fluorouracil were less potent in the apoptosis-resistant lineage, with RF values of 9.2 and 55, respectively (Fig. 7a).
Complex | IC50a [μM] | Resistance factorb | ||||
---|---|---|---|---|---|---|
AGS | HCT116 (wild type) | HCT116 p53-/- | TC7 | p53 dependence [p53-/-/wild type] | Apoptosis [TC7/HCT116] | |
a IC50 values is the concentration of Ru complexes required to inhibit 50% of cell growth with respect to control groups, measured by MTT assay after 48 h of incubation. Data obtained are based on the average of three independent experiments, and the reported errors are the corresponding standard deviations. The IC50 were corrected using actual [Ru] determined using ICP-OES. b Resistance factors were calculated by taking the ratio of IC50. | ||||||
532m | 4.13 ± 0.57 | 10.1 ± 0.3 | 9.20 ± 0.44 | 8.03 ± 0.78 | 0.9 | 0.8 |
532p | 6.75 ± 2.06 | 17.8 ± 1.5 | 14.5 ± 1.18 | 13.6 ± 1.3 | 0.8 | 0.8 |
Oxaliplatin | n.d. | 2.62 ± 0.63 | 12.1 ± 1.8 | 24.0 ± 12.4 | 4.6 | 9.2 |
5-Fluorouracil | n.d. | 9.94 ± 2.3 | 22.6 ± 5.5 | 548 ± 218 | 2.3 | 55 |
Cisplatin | 31.5 ± 3.5 | n.d. | n.d. | n.d. | n.d. | n.d. |
We validated the apoptosis-independent activity of 532m and 532p by examining several hallmarks of apoptosis. Anticancer drugs such as oxaliplatin and 5-fluorouracil that act via apoptosis are dependent on apoptosis-executor caspases and inhibition of these caspases would result in decreased efficacies. Treatment with these compounds would result in cleavage and activation of caspase, and increase expression of other related apoptotic biomarkers. Indeed, co-incubation with a pan caspase inhibitor, z-VAD-fmk, conferred a cytoprotective effect to HCT116 cells treated with oxaliplatin but had no effect on cells treated with 532m or 532p (Fig. 7b and ESI Fig. S16†). Western blot analysis of HCT116 cells treated with oxaliplatin and 5-fluorouracil revealed significant induction of p53 expression and cleaved caspase-3 at both 24 h and 48 h time-points, and downstream cleavage of PARP-1 at 48 h. In contrast, cells treated with 532m and 532p did not display such expression profile regardless of duration or concentration of treatment (Fig. 7c). Further investigation also showed a lack of induction of other apoptotic biomarkers such as cleaved caspase-7 or -9, or DNA-damage biomarker pH2AX by 532m or 532p (ESI Fig. S17†).
Taken together, the similar activities in both apoptosis-sensitive (HCT116) and apoptosis-resistant (TC7) cell lines, caspase-independent activity and the lack of induction of apoptotic biomarkers confirmed the apoptosis-independent activity of 532m and 532p, validating the accuracy and efficiency of the phenotypic screening strategy.
As a first line of investigation for the identification of apoptosis-independent candidates, determination of the exact mode of cell death would be less relevant than ascertaining that cell death induction was indeed non-apoptotic. In that regard, phenotypic screening designed to delineated non-apoptotic from apoptotic mechanisms would be more revealing than a generic efficacy-based approach. In order to fully realized the potential and suitability of apoptosis-independent candidates such as 532m and 532p for the treatment of resistant cancer progenies, more throughout scrutiny including target identification, mode-of-action and in vivo studies would have to be applied. Nevertheless, the combined use of a convenient compound-library generation methodology, and a specific phenotypic screening strategy presented in this study is a reliable first step to identifying effective treatment for apoptosis-resistant cancers and development of similar strategies may be a way to accelerate the discovery of such compounds.
Footnote |
† Electronic supplementary information (ESI) available: Experimental protocols, assembly yield/purity of RAS complexes from the on-plate synthesis and cell viability data from phenotypic screening assay, 1H NMR, ESI-MS and UV-vis spectra for 532m and 532p, and western blot data in PDF format. See DOI: 10.1039/c7sc00497d |
This journal is © The Royal Society of Chemistry 2017 |