Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Berry Extracts and Their Bioactive Compounds Mitigate LPS and DNFB-Mediated Dendritic Cell Activation and Induction of Antigen Specific T-Cell Effector Responses
Next Article in Special Issue
Antioxidant Use after Diagnosis of Head and Neck Squamous Cell Carcinoma (HNSCC): A Systematic Review of Application during Radiotherapy and in Second Primary Cancer Prevention
Previous Article in Journal
Analysis of the Flavonoidome Reveals the Different Health-Promoting Flavonoid Characteristics in Fruit
Previous Article in Special Issue
The Contribution of Oxidative Stress to NF1-Altered Tumors
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First Generation of Antioxidant Precursors for Bioisosteric Se-NSAIDs: Design, Synthesis, and In Vitro and In Vivo Anticancer Evaluation

1
Department of Pharmaceutical Technology and Chemistry, University of Navarra, Irunlarrea 1, 31008 Pamplona, Spain
2
Instituto de Investigación Sanitaria de Navarra (IdiSNA), Irunlarrea 3, 31008 Pamplona, Spain
3
Department of Pharmacology, Penn State Cancer Institute, CH72, Penn State College of Medicine, 500 University Drive, Hershey, PA 17033, USA
4
Department of Health Sciences, Public University of Navarra, Avda. Barañain s/n, 31008 Pamplona, Spain
*
Authors to whom correspondence should be addressed.
Antioxidants 2023, 12(9), 1666; https://doi.org/10.3390/antiox12091666
Submission received: 18 July 2023 / Revised: 15 August 2023 / Accepted: 16 August 2023 / Published: 24 August 2023
(This article belongs to the Special Issue Oxidative Stress in Tumor Genesis, Progression and Therapy)
Graphical abstract
">
Figure 1
<p>Chemical structures of some NSAID derivatives bearing Se [<a href="#B60-antioxidants-12-01666" class="html-bibr">60</a>,<a href="#B61-antioxidants-12-01666" class="html-bibr">61</a>,<a href="#B62-antioxidants-12-01666" class="html-bibr">62</a>,<a href="#B63-antioxidants-12-01666" class="html-bibr">63</a>,<a href="#B64-antioxidants-12-01666" class="html-bibr">64</a>,<a href="#B65-antioxidants-12-01666" class="html-bibr">65</a>].</p> ">
Figure 2
<p>(<b>A</b>) Schematic diagrams of the chemical design of Se-modified NSAIDs, including previous possible cleavages according to where the bond is broken (1 or 2). (<b>B</b>) General structure of the new Se-NSAID derivatives <b>1</b>–<b>5</b> synthesized in this work and the parent NSAIDs used.</p> ">
Figure 3
<p>(<b>A</b>) DPPH radical scavenging activity of the Se-NSAID derivatives at 0.06 mg/mL and measured at eight time points. (<b>B</b>) DPPH radical scavenging activity of the Se-NSAID derivatives at 0.03 mg/mL and measured at eight time points. Data are presented as the mean ± SD of three independent experiments.</p> ">
Figure 4
<p>Induction of apoptosis in HCT-116 cells. (<b>A</b>) Cells were treated with five concentrations of <b>5</b> for 24 h and examined on a Muse™ automated cell analyzer with the Annexin V &amp; Dead Cell apoptosis assay. (<b>B</b>) Analogous independent experiment performed with the Caspase 3/7 apoptosis assay. (<b>C</b>,<b>D</b>) Quantification of the cell population with the Annexin V &amp; Dead Cell and Caspase 3/7 assays, respectively. Data are presented as the mean ± SEM of at least three independent experiments. ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05 comparing control and different concentrations of compound <b>5</b>.</p> ">
Figure 5
<p>In vivo toxicity of compound <b>5</b> in mice. Female nude mice (<span class="html-italic">n</span> = 5 in each group) were treated either with the vehicle or compound <b>5</b> for 27 days. (<b>A</b>) Biochemical analysis of the blood serum withdrawn at the endpoint. The treatment with a dose of 7.5 mg/Kg of <b>5</b> given thrice weekly showed no significant changes in the parameters measured compared to the control group. (<b>B</b>) Body weight changes of mice during the treatment. Data are presented as the mean ± SEM.</p> ">
Figure 6
<p>In vivo efficacy of compound <b>5</b> in a colon cancer xenograft mouse model. Female nude mice were inoculated with 2.5 × 10<sup>6</sup> HCT-116 cells on each flank. After 15 days, the HCT-116 tumor-bearing mice (<span class="html-italic">n</span> = 5 in each group) were treated either with the vehicle or compound <b>5</b> for additional 23 days. (<b>A</b>) Growth curves of tumor volume. The treatment with a dose of 10 mg/Kg of <b>5</b> given thrice weekly achieved significant tumor growth inhibition compared to the control group. (<b>B</b>) Tumor weight measured at the endpoint. (<b>C</b>) Body weight changes of mice during the treatment. Data are presented as the mean ± SEM. **** <span class="html-italic">p</span> &lt; 0.0001 comparing control and treatment with compound <b>5</b>.</p> ">
Scheme 1
<p>Synthesis of the Se-NSAID derivatives. Reagents and conditions: (i) ClCOCOCl, CH<sub>2</sub>Cl<sub>2</sub>, 12 h, room temperature; (ii) ClCOCOCl, CH<sub>2</sub>Cl<sub>2</sub>/<span class="html-italic">N</span>,<span class="html-italic">N</span>-DMF, 12 h, room temperature; (iii) H<sub>2</sub>O, room temperature; (iv) THF/H<sub>2</sub>O, 30–60 min, room temperature.</p> ">
Versions Notes

Abstract

:
The introduction of selenium (Se) into organic scaffolds has been demonstrated to be a promising framework in the field of medicinal chemistry. A novel design of nonsteroidal anti-inflammatory drug (NSAID) derivatives based on a bioisosteric replacement via the incorporation of Se as diacyl diselenide is reported. The antioxidant activity was assessed using the DPPH radical scavenging assay. The new Se-NSAID derivatives bearing this unique combination showed antioxidant activity in a time- and dose-dependent manner, and also displayed different antiproliferative profiles in a panel of eight cancer cell lines as determined by the MTT assay. Ibuprofen derivative 5 was not only the most antioxidant agent, but also selectively induced toxicity in all the cancer cell lines tested (IC50 < 10 µM) while sparing nonmalignant cells, and induced apoptosis partially without enhancing the caspase 3/7 activity. Furthermore, NSAID derivative 5 significantly suppressed tumor growth in a subcutaneous colon cancer xenograft mouse model (10 mg/kg, TGI = 72%, and T/C = 38%) without exhibiting any apparent toxicity. To our knowledge, this work constitutes the first report on in vitro and in vivo anticancer activity of an unprecedented Se-NSAID hybrid derivative and its rational use for developing precursors for bioisosteric selenocompounds with appealing therapeutic applications.

Graphical Abstract">

Graphical Abstract

1. Introduction

Selenium (Se) is an essential trace element that has been recognized to play an important role in human health [1]. Se exerts its biological activities mainly in the form of selenocysteine, a proteogenic amino acid bearing the selenol functionality (SeH), which is part of selenoproteins [2]. The unique features of this selenol moiety enable selenoproteins to accomplish a wide variety of different physiological functions [2], including protection against oxidative stress, immune response, cellular differentiation [3], and maintenance of the redox homeostasis [4], among others. Selenoprotein-mediated mechanisms are also involved in the prevention, onset, and clinical outcome of various diseases such as cancer [2]. Moreover, a positive relationship between Se status and a favorable prognosis of these diseases has been generally observed [2].
In this context, a large body of evidence supports the role of Se in preventing cancer incidence and blocking tumor metastasis [5]. Recently, a potential benefit of the dietary intake of certain microelements, including Se, on colorectal cancer was also reported [6]. Furthermore, epidemiological studies have shown that low Se intake and low plasma Se levels could lead to an increased incidence of several cancers [7,8,9]. Moreover, Se supplementation alongside traditional therapies has been shown to increase the efficacy of anticancer drugs, limit associated side outcomes, and improve the general conditions of the patients [10,11,12].
Accordingly, over recent decades Se-containing compounds have attracted growing interest as anticancer agents due to their high efficacy and selectivity against cancer cells [4]. In this context, several Se-containing compounds have been proven to exert potent anticancer activity in several xenograft models [13,14,15], and thus they have been considered appealing agents in the field of drug discovery for cancer therapy [1,4,16,17]. The proposed mechanisms involved in the anticancer activity of Se-containing compounds include involvement in angiogenesis, induction of cell cycle arrest, and the initiation of cell death triggered by diverse processes such as apoptosis, necrosis, autophagy, ferroptosis, necroptosis, or entosis, among others [16].
The efficacy of selenocompounds as anticancer agents is correlated with their metabolic route, dose, and the chemical form of Se [5,18]. Hence, the main Se entities that have been widely explored to date are selenides [19,20,21,22], selenocyanates [22,23,24,25,26], isoselenocyanates [27,28,29], selenazo heterocyclic compounds [13,30,31], selenoesters [32,33,34], selenoureas [35,36,37,38,39], and diselenides [19,22,23,24,25,39,40]. Additionally, several compounds bearing these functionalities have also been reported to have a potent radical scavenging capacity, allowing the design of molecules with an interesting dual profile as antioxidant and antiproliferative agents [19,25,37,38].
Furthermore, among the scaffolds that have attracted great attention due to their therapeutic properties, nonsteroidal anti-inflammatory drugs (NSAIDs) have recently emerged as promising frameworks to be modified through the incorporation of different Se functionalities [41]. NSAIDs are prominent agents with applications in cancer treatment and prevention [42] aside from their traditional analgesic, antipyretic, and anti-inflammatory prescriptions. Numerous meta-analyses and systematic reviews to date have supported their chemopreventive and chemotherapeutic potential [43,44,45,46], especially in colon cancer [47,48,49,50,51]. Hence, the use of aspirin (ASA) and ibuprofen (Ibup) for reducing the high risk of colorectal cancer with DNA mismatch repair (MMR) genes is reported [52]. In another study, a synergistic effect of sulindac and bexarotene was observed in familial adenomatous polyposis (FAP) patients [53]. Likewise, another recent cohort study suggested an association between NSAID exposure and reduced cancer recurrence without increasing the risk of anastomotic leakage [54]. The key mechanisms for the therapeutic activity of NSAIDs are the inhibition of cyclooxygenases (COX) and suppression of the production of prostaglandins (PGs) [55]. Interestingly, human colon cancers have elevated PGE2 levels, and overexpression of COX has been found in 45% of colon adenomas and 85% of colon carcinomas [56]. Additionally, several studies have shown that the anticancer activity of NSAIDs could also be exerted through COX-independent mechanisms [57], such as the regulation of the Wnt pathway [58] or the DLL1 and NOTCH1 expression [59].
In this context, the development of new NSAID analogs based on the combination of NSAIDs and Se entities has attracted great attention (Figure 1). The ASA analogs bearing Se through a selenazolidine heterocycle (AS-10) [60] or a selenocyanate group (Se-Aspirin) [61] showed potent growth inhibition via different mechanisms in pancreatic and colorectal cancer cells. Other Se-NSAID derivatives with a longer carbon chain elongation or even with trifluoromethyl selenides attached through an amide [62] or ester [63] group also displayed potent antiproliferative activity via induction of apoptosis. Remarkably, a derivative of celecoxib synthesized through the modification of the position-3 of the pyrazole ring with methylene selenocyanate, namely selenocoxib-1, showed better efficacy against prostate cancer cells than celecoxib by inducing apoptosis and decreasing HIF-1α, p-AKT, and Bcl-2 levels [64]. Likewise, the glutathione conjugate (selenocoxib-1-GSH) showed potent tumor growth inhibition associated with the suppression of COX-2 and PI3K/AKT signaling pathways [65]. Furthermore, we recently reported a library of Se-NSAID derivatives functionalized with a selenoester group. Among them, Se-indomethacin was identified as a potent antiproliferative agent and inducer of apoptosis in breast cancer cells [66].
Therefore, considering the chemotherapeutic activity of NSAIDs, the literature supporting the structural modification of NSAID scaffolds to create more potent anticancer drugs, and the reported antitumor profile of organoselenium compounds as chemopreventive and chemotherapeutic agents in preclinical models, we hypothesized that the combination of NSAIDs and a diacyl diselenide group could lead to novel Se-NSAID analogs with a therapeutic application as precursors of active molecules. To our knowledge, this is the first study that reports the therapeutic potential of precursors for bioisosteric Se-NSAIDs based on the diacyl diselenide moiety, and specifically their anticancer effect both in vitro and in vivo, along with their preliminary evaluation as antioxidant agents.
Herein, we designed and synthesized five new Se-NSAID hybrid analogs as precursors for the corresponding selenols, which isosterically replace the hydroxide (-OH) of the carboxylic acid functionality of the NSAIDs by -SeH, and evaluated their antiproliferative activity in several cancer cell lines. The selectivity of the new compounds was determined in two nonmalignant cell lines. The following five NSAIDs were selected to be structurally modified: aspirin (ASA), indomethacin (Ind), naproxen (Nap), ketoprofen (Ket), and ibuprofen (Ibup). Compared with our previously reported methylselenoester NSAIDs [66], Ibup was added to this set of newly synthesized diacyldiselenides. The capacity of the new diacyl diselenides to scavenge free radicals was also studied in vitro at different time points and concentrations. The most active and selective compound 5 was further evaluated to assess its capability to induce apoptosis and inhibit tumor growth in a subcutaneous colon cancer xenograft mouse model in order to determine if the novel combination of Se and NSAID scaffolds through a diacyl diselenide group is an appealing new framework for therapeutic development.

2. Materials and Methods

2.1. Chemistry

2.1.1. General Remarks

Starting materials, reagents, and solvents were purchased from commercial suppliers and used as received without further purification. Reaction courses were monitored by thin-layer chromatography (TLC) performed on precoated silica gel 60 F254 aluminum sheets (Merck, Darmstadt, Germany), and the spots were visualized under UV light. The crude reaction products were purified by silica gel column chromatography using silica gel 60 Å (0.040–0.063 mm, Merck, Darmstadt, Germany) or by flash chromatography, using hexane/ethyl acetate as the elution solvent. 1H, 13C, and 77Se nuclear magnetic resonance (NMR) spectra were recorded on a Bruker Avance Neo 400 MHz operating at 400, 100, and 76 MHz, respectively, using CDCl3 as solvent and TMS as the internal standard. Chemical shifts are reported in δ values (ppm) and coupling constants (J) values are reported in hertz (Hz). Melting points (mp) were determined with a Mettler FP82 + FP80 apparatus (Greifensee, Switzerland). High-resolution mass spectrometry (HRMS) was performed on a Micromass Q-TOF mass spectrometer. The mass spectrum of compound 5 was obtained by the direct injection of a methanolic solution of 5 in an ExpressionS CMS Advion mass spectrometer with a MX-Class pump operating in a negative ionization mode. For purity testing, quantitative NMR (qNMR) spectroscopy was performed in concordance with a standardized protocol [67]. Dimethyl sulfone (TraceCERT certified reference compound, lot no. BCBQ8884V, purity: 99.96%) purchased from Sigma-Aldrich (Alcobendas, Madrid, Spain) was used as an internal calibrant. All the newly synthesized diacyl diselenides were >95% pure.

2.1.2. General Procedure for the Preparation of the Diacyl Diselenide Derivatives

The NSAID acyl chloride (2 mmol) was added to a water solution of sodium hydrogen selenide (2 mmol) formed in situ by the reaction of powdered selenium (2 mmol) with sodium borohydride (4 mmol) in water, following the procedure previously reported by Sanmartin et al. [68] with modifications. THF (5 mL) was added, and the mixture was stirred at room temperature for 30–60 min. The reaction mixture was then extracted with methylene chloride (3 × 20 mL). The organic layers were dried with sodium sulfate and concentrated in vacuo. The resulting product was purified using silica gel chromatography or flash chromatography with a hexane/ethyl acetate gradient. For NSAIDs Ind, Ket, and Ibup, the synthesis of the acyl chloride derivatives was necessary prior to the formation of the diacyl diselenide compounds. These derivatives were obtained by the reaction of the carboxylic acid reagent (1 mmol) with oxalyl chloride (3 mmol) in methylene chloride (25 mL). In the case of the reactions involving Ket or Ibup acids, N,N-dimethylformamide (0.1 mL) was also added. The mixture was stirred overnight at room temperature. The crude product was isolated by the evaporation of the solvent in vacuo and washed with methylene chloride (2 × 15 mL). The acyl chlorides synthesized were then used without further purification.
2-Acetoxybenzoic diselenoperoxyanhydride (1). The title compound was synthesized from O-acetylsalicyloyl chloride, selenium, and sodium borohydride according to the general procedure described above. A yellow solid was obtained. Yield: 43%; mp: 55–58 °C. 1H NMR (400 MHz, CDCl3) δ 2.40 (s, 6H, CH3, CH3′); 7.21 (dd, 2H, J = 8.2 and 0.9 Hz, H3, H3′); 7.39 (td, 2H, J = 7.7 and 1.1 Hz, H5, H5′); 7.63 (td, 2H, J = 8.1, 7.6 and 1.6 Hz, H4, H4′); 8.02 (dd, 2H, J = 7.8 and 1.6 Hz, H6, H6′). 13C NMR (100 MHz, CDCl3) δ 21.5 (CH3, CH3′), 124.3, 126.6, 129.8, 130.4, 134.6, 148.0 (Caryl), 169.1 (CH3-C=O), 184.5 (Se-C=O). 77Se NMR (76 MHz, CDCl3) δ 659. HRMS m/z: calcd for C18H14O6Se2 [M + Na]+, 508.9013; found, 508.9025. Purity (qNMR: 25 °C, CDCl3, ms = 1.8 mg, mIC = 4.0 mg): 97.1%.
2-(1-(4-Chlorobenzoyl)-5-methoxy-2-methyl-1H-indol-3-yl)acetic diselenoperoxyanhydride (2). The title compound was synthesized from 1-(4-chlorobenzoyl)-5-methoxy-2-methyl-3-indoleacetic acid, oxalyl chloride, selenium, and sodium borohydride according to the general procedure described above. A yellow solid was obtained. Yield: 13%; mp: 78–79 °C. 1H NMR (400 MHz, CDCl3) δ 2.44 (s, 6H, CH3, CH3′); 3.84 (s, 6H, OCH3, OCH3′); 4.00 (s, 4H, CH2, CH2′); 6.71 (dd, 2H, J = 9.0 and 2.4 Hz, H6A, H6′A); 6.91 (d, 2H, J = 2.2 Hz, H4A, H4′A); 6.93 (d, 2H, J = 9.1 Hz, H7A, H7′A); 7.48 (dd, 4H, J = 8.4 Hz, H3B + H5B, H3′B + H5′B); 7.68 (dd, 4H, J = 8.4 Hz, H2B + H6B, H2′B + H6′B). 13C NMR (100 MHz, CDCl3) δ 13.7 (CH3, CH3′), 41.2 (CH2, CH2′), 55.9 (OCH3, OCH3′), 100.9, 110.8, 112.7, 115.3, 129.4, 130.4, 131.0, 131.5, 133.6, 138.2, 139.8, 156.5 (Caryl), 168.4 (Ph-C=O), 193.4 (Se-C=O). 77Se NMR (76 MHz, CDCl3) δ 657. HRMS m/z: calcd for C38H30Cl2N2O6Se2 [M + H]+, 840.9884; found, 840.9842. Purity (qNMR: 25 °C, CDCl3, ms = 2.4 mg, mIC = 4.0 mg): 96.6%.
2-(6-Methoxynaphthalen-2-yl)propanoic diselenoperoxyanhydride (3). The title compound was synthesized from (S)-2-(6-methoxy-2-naphthyl)propionyl chloride, selenium, and sodium borohydride according to the general procedure described above. A pink solid was obtained. Yield: 20%; mp: 113–115 °C. 1H NMR (400 MHz, CDCl3) δ 1.63 (d, 6H, J = 7.1 Hz, CH3, CH3′); 3.92 (s, 6H, OCH3, OCH3′); 4.17 (q, 2H, J = 7.1 Hz, CH, CH′); 7.13 (d, 2H, J = 2.3 Hz, H5, H5′); 7.16 (dd, 2H, J = 8.9 and 2.5 Hz, H7, H7′); 7.37 (dd, 2H, J = 8.6 and 1.7 Hz, H8, H8′); 7.73 (m, 6H, H1 + H3 + H4, H1′ + H3′ + H4′). 13C NMR (100 MHz, CDCl3) δ 18.3 (CH3, CH3′), 55.5 (OCH3, OCH3′), 56.6 (CH, CH′), 105.8, 119.5, 126.9, 127.8, 128.2, 129.0, 129.6, 133.3, 134.5, 158.3 (Caryl), 197.2 (Se-C=O). 77Se NMR (76 MHz, CDCl3) δ 636. HRMS m/z: calcd for C28H26O4Se2 [M + Na]+, 609.0054; found, 609.0061. Purity (qNMR: 25 °C, CDCl3, ms = 1.8 mg, mIC = 4.0 mg): 96.5%.
2-(3-Benzoylphenyl)propanoic diselenoperoxyanhydride (4). The title compound was synthesized from 2-(3-benzoylphenyl)propionic acid, oxalyl chloride, selenium, and sodium borohydride according to the general procedure described above. A white solid was obtained. Yield: 36%; mp: 65–66 °C. 1H NMR (400 MHz, CDCl3) δ 1.60 (dd, 6H, J = 7.1 and 1.5 Hz, CH3, CH3′); 4.16 (m, 2H, CH, CH′); 7.50 (m, 6H, H3B+H4B+H5B, H3′B+H4′B+H5′B); 7.59 (m, 4H, H5A + H6A, H5′A + H6′A); 7.79 (m, 8H, H2A + H4A + H2B + H6B, H2′A + H4′A + H2′B + H6′B). 13C NMR (100 MHz, CDCl3) δ 18.3 (CH3, CH3′), 56.4 (CH, CH′), 128.5, 128.6, 129.2, 130.2, 130.3, 130.4, 130.5, 132.6, 132.7, 132.8, 137.4, 138.4, 138.5, 138.6 (Caryl), 195.9 (Se-C=O), 196.3 (Ph-C=O). 77Se NMR (76 MHz, CDCl3) δ 633.8 and 634.1. HRMS m/z: calcd for C32H26O4Se2 [M + Na]+, 657.0054; found, 657.0062. Purity (qNMR: 25 °C, CDCl3, ms = 1.7 mg, mIC = 4.0 mg): 97.5%.
2-(4-Isobutylphenyl)propanoic diselenoperoxyanhydride (5). The title compound was synthesized from 2-(4-isobutylphenyl)propanoic acid, oxalyl chloride, selenium, and sodium borohydride according to the general procedure described above. A yellow oil was obtained. Yield: 44%. 1H NMR (400 MHz, CDCl3) δ 0.89 (d, 12H, J = 6.6 Hz, 2CH3, 2CH3′); 1.46 (d, 6H, J = 7.0 Hz, CH3, CH3′); 1.84 (sept, 2H, J = 6.8 Hz, CH, CH′); 2.44 (d, 4H, J = 7.2 Hz, CH2, CH2′); 3.98 (q, 2H, J = 7.2 Hz, CH-Ph, CH′-Ph); 7.09 (m, 8H, Hz, H2 + H3 + H5 + H6, H2′ + H3′ + H5′ + H6′). 13C NMR (100 MHz, CDCl3) δ 18.0, 18.1 (CH3, CH3′), 22.5 (2CH3, 2CH3′), 30.3 (CH, CH′), 45.2 (CH2, CH2′), 59.0, 59.1 (CH-Ph, CH′-Ph), 128.30, 128.33, 129.7, 135.28, 135.32, 141.57, 141.58 (Caryl), 198.6, 198.7 (Se-C=O). 77Se NMR (76 MHz, CDCl3) δ 791 and 792. HRMS m/z: calcd for C26H34O2Se2 [M + Na]+, 561.0781; found, 561.0748. Purity (qNMR: 25 °C, CDCl3, ms = 3.0 mg, mIC = 4.0 mg): 98.7%.

2.2. Biology

2.2.1. DPPH Radical Scavenging Assay

The antioxidant activity of the novel diacyl diselenide-based compounds was determined by the 2,2-diphenyl-1-picrylhydrazyl (DPPH) colorimetric assay following a previously published procedure [22]. Briefly, the compounds were dissolved at a concentration of 1 mg/mL in absolute methanol, and a dilution was prepared. The stability of the DPPH radical through the time of analysis was checked by absorbance measurements of a methanolic solution (100 µM) of DPPH protected from light. Then, 100 µL of the DPPH solution was added to 100 µL of the compounds’ solution to a final concentration of 0.03 mg/mL or 0.06 mg/mL, and the decolorization of the purple radical to the yellowish reduced form was followed by recording the absorbance at 517 nm. Ascorbic acid (Asc) and Trolox were used as positive controls. Both of them, Asc and Trolox, were purchased from Sigma-Aldrich. Determinations were recorded at eight different time intervals from 0 to 120 min on a BioTeck PowerWave XS spectrophotometer, and the data were collected using KCJunior version 1.41. software. All the measurements were carried out in triplicate. Results are expressed as the percentage of the radical scavenged, calculated using the following formula: % DPPH radical scavenging = [(AcontrolAsample)/Acontrol] × 100, wherein Acontrol refers to the absorbance of the negative control and Asample refers to the absorbance of the tested compounds.

2.2.2. Cell Culture Conditions

The cell lines were purchased from the American Type Culture Collection (ATCC). The cancer cell lines DU-145, PC-3, T-47D, MDA-MB-231, H1299, and A549, and the non-tumorigenic cells (184B5 and BEAS-2B), were maintained in RPMI 1640 medium (Gibco, Alcobendas, Madrid, Spain) supplemented with 10% fetal bovine serum (FBS; Gibco) and 1% antibiotics (10.00 units/mL penicillin and 10.00 mg/mL streptomycin; Gibco). Colon cancer cells HT-29 and HCT-116 were cultured in McCoy’s 5A medium (Gibco) supplemented with 10% fetal bovine serum (FBS; Gibco) and 1% antibiotics (10.00 units/mL penicillin and 10.00 mg/mL streptomycin; Gibco). Cells were maintained in tissue culture flasks at 37 °C and 5% CO2. The culture medium was replaced every three days.

2.2.3. Cell Viability Assay

The antiproliferative activity and selectivity in vitro of the Se-NSAID derivatives were determined by the (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) (MTT) assay [37]. Compounds were dissolved in DMSO at a concentration of 0.01 M and serial dilutions were prepared. The Se-NSAID derivatives were tested at seven concentrations ranging between 1 and 100 µM in the colon (HT-29 and HCT-116), prostate (DU-145 and PC-3), breast (MDA-MB-231 and T-47D), and lung (H1299 and A549) cancer cell lines. The non-tumorigenic breast (184B5) and lung (BEAS-2B) cell lines were also tested to assess the selectivity of the compounds following the same procedure. A total of 1 × 104 cells were seeded per well in 96-well plates and incubated at 5% CO2 and 37 °C for 24 h. Then, cells were treated with either DMSO or increasing concentrations of the corresponding Se-NSAID derivatives for 48 h. A volume of 20 µL of MTT (5 mg/mL) was added to each well 2.5 h before the termination point. At the end of the experiment, the medium was removed, and the resultant formazan crystals formed by mitochondrial reduction of MTT were dissolved in 50 µL of DMSO. Absorbance was measured at 550 nm. IC50, GI50, TGI, and LC50 values were calculated using OriginPro version 8.5.1. software by nonlinear curve fitting. Selectivity indexes (SI) were calculated as the ratio of the IC50 values determined for the nonmalignant and the tumoral cells in the breast (IC50 (184B5)/IC50 (MDA-MB-231) and IC50 (184B5)/IC50 (T-47D)) and lung (IC50 (BEAS-2B)/IC50 (H1299) and IC50 (BEAS-2B)/IC50 (A549)) cell lines. Data were obtained from at least three independent experiments performed in triplicates.

2.2.4. Apoptosis Assays

The capacity of the Se-NSAID derivatives to induce apoptosis was analyzed using both a Caspase 3/7 assay kit and Annexin V & Dead Cell assay kit (EMD Millipore, Darmstadt, Germany). HCT-116 cells were seeded in 6-well plates at a density of 5 × 105 cells per well. Cells were treated either with DMSO (control) or with different concentrations of compound 5 and incubated for 24 h. At the end of the treatment, cells were collected and treated with the respective dyes according to the manufacturer’s protocol. For the Caspase 3/7 assay, cells were stained with 5 µL of Muse™ Caspase 3/7 working solution and incubated at 37 °C for 30 min. After incubation, 150 mL of Muse™ Caspase 7-AAD working solution was added and the samples were incubated for another 5 min in the dark at room temperature prior to analysis. For Annexin V assay, cells were stained with 100 µL pf Muse™ Annexin V & Dead Cell Reagent and incubated for 20 min at room temperature in the dark. Samples for both assays were analyzed on a Muse™ Cell Analyzer (Merck Millipore, Darmstadt, Germany) using Muse™ version 1.4. software.

2.2.5. In Vivo Xenograft Studies

Female mice of 4–5 weeks of age (The Jackson Laboratory, Bar Harbor, ME, USA) were kept on the daily cycle of 12 h light/dark and received food and water ad libitum throughout the studies. All animal procedures were performed following the National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals and protocols were approved by the Animal Care and Use Committee of the Penn State College of Medicine.

2.2.6. Toxicity Study

Mice were treated either with intraperitoneal doses of 7.5 mg/Kg of compound 5 or with DMSO (vehicle) thrice weekly and kept under observation for 27 days. Body weight, behavior, and food and water intake were monitored every alternate day. At the end of the study, animals were sacrificed by euthanizing the animals by CO2 asphyxiation. The blood serum was withdrawn and subjected to the biochemical analysis of creatinine (mg/dL), AST(SGOT) (U/L), ALT(SGPT) (U/L), Alk. Phos (U/L), and BUN (mg/dL).

2.2.7. Assay of Tumor Growth Inhibition on HCT-116 Colon Cancer Xenograft Model

All mice were housed in quarantine for one week prior to the inoculation of the cancer cells. Female nude mice (4–5 weeks old) were subcutaneously injected (100 µL of PBS/flank) with 2.5 × 106 HCT-116 colon cancer cells in each flank. Two weeks after inoculation, mice were sorted randomly into treatment and control groups (five mice/group) such that the mean tumor volumes were similar between groups. Mice from the treatment group were then injected intraperitoneally with 10 mg/Kg of compound 5 every 2 days over a 5-week treatment period. Mice that were injected with only DMSO (vehicle) were used as negative controls. The body weight of animals and the tumor volume were monitored thrice weekly. Tumor volumes were calculated using the formula (L × W2)/2, wherein L and W are the length and width of the tumor mass measured with a vernier caliper, respectively. At the end of the experiment, mice were sacrificed by euthanizing the animals by CO2 asphyxiation, and the tumors were excised and weighed. Inspection of major morphological changes in organs was also carried out. The antitumor effect was assessed by the tumor growth inhibition (TGI) and the relative increment ratio (T/C) parameters in the aspects of tumor volume and tumor weight, respectively. The TGI was calculated using the following formula: TGI (%) = [1 − (TT0)/(CC0)] × 100, wherein T and C are the final tumor volumes measured at the end of the treatment of compound-treated mice and vehicle-treated mice, respectively, whereas T0 and C0 are the tumor volumes measured at the beginning of the treatment of compound-treated mice and vehicle-treated mice, respectively. The T/C value was calculated as WT/WC × 100, wherein WT is the mean of tumor weight of the compound-treated mice, and WC is the mean of tumor weight of the vehicle-treated mice.

2.2.8. Statistical Analysis

Data were expressed as the mean ± SD (standard deviation) unless otherwise specified, and experiments were performed at least thrice in triplicates. Non-linear curve regression analysis calculated by OriginPro version 8.5.1. software was used to assess the IC50, GI50, TGI, and LD50 values. For animal experiments, all animal studies were conducted using five animals per group (n = 5) and two tumors per animal. The two-way analysis of variance (ANOVA) was used to calculate the statistical significance of differences comparing compound 5 and control. Data were analyzed using GraphPad Prism version 8.0.1., and the statistically significant values (p-value) for ANOVA analysis were taken as **** p < 0.0001, ** p < 0.01, and * p < 0.05.

3. Results and Discussion

3.1. Structural Design

In recent years, the role of NSAIDs as single agents or in combination with current therapies in the prevention and treatment of cancer has attracted great interest, and their anticancer effect has been widely reported [41]. The modifications made to date in the literature reporting Se-NSAID analogs (Figure 1) considered the inclusion of the Se atom through well-known functional groups susceptible to being later metabolized to form active fragments. However, the part of the structure that retained the Se atom was never the NSAID itself, given the structural characteristics of these functional groups [60,61,62,63,64]. In this context, only selenocoxib-1 and selenocoxib-1-GSH [65] retained the Se moiety, but the NSAIDs were already modified and thus the parent celecoxib was not directly involved in the biological activity of the Se moiety once the compounds were metabolized. Consequently, there is currently no study reported to date regarding the biological activity of therapeutic Se-NSAID hybrids in which the entire NSAID could potentially be maintained without any structural modifications other than isosteric Se replacement of the oxygen atom of the NSAID-carboxylic acid group (Figure 2A). In this regard, given that the direct attempt to obtain the selenoic acid derivatives of NSAIDs could lead to unstable molecules, a bioisosteric replacement approach is proposed instead. Thus, we suggest that the modification of the NSAID scaffold via the novel inclusion of a diacyl diselenide functional group, which has recently been reported as carbonic anhydrase inhibitors [69], could lead to precursors for bioisosteric selenoic NSAID acids (Figure 2A). The rationale of our design resides in the liability and possible hydrolysis of the unique Se–Se bond present in the diacyl diselenide group. We hypothesized that this Se–Se bond would be liable to be broken in the metabolic process of the drug. Thus, the parent NSAID would retain the Se without any further modifications. Furthermore, the NSAID would also maintain the Se moiety in the form of a selenoic acid, in which -OH has been substituted with -SeH. SeH is a precursor of hydrogen selenide (H2Se), one of the main metabolites of Se [70] as it is the common intermediate of the reductive pathway and the catabolism of selenoamino acids [2] and with a prominent role in the biological activity of Se compounds [71].
Hence, five new selenocompounds were synthesized comprising diacyl diselenides decorated with two moieties of the corresponding NSAID forming a symmetric molecule, according to the general structure shown in Figure 2B. The NSAID scaffolds employed in the synthesis of the compounds covered a variety of chemical structures with reported chemotherapeutic activity, such as salicylate (ASA), acetic (indomethacin, Ind), and arylpropionic (naproxen, Nap, ketoprofen, Ket, and Ibup) acid derivatives. To our knowledge, this is the first time that a series of diacyl diselenide-based compounds has been reported to possess potent anticancer activity.

3.2. Chemistry

Some methods for obtaining diacyl diselenides have been reported previously. These methods include the reaction of acyl chlorides either with bis(trimethylsilyl)selenide (Me3SiSeSiH2Me3) [72]; with reduced Se from the reaction with samarium diiodide (SmI2) [73]; or with LiAlHSeH prepared from elemental Se and LiAlH4 and followed by quenching with the addition of I2 and KI [74,75,76]. The oxidation of potassium selenocarboxylates with xenon difluoride (XeF2) [77] has also been described, among others. However, these methods often include some inconveniences, such as the number of steps, the difficulty of purification, and the limited availability or expensive cost of some of the starting materials. Herein, we propose a new synthetic method for obtaining diacyl diselenide compounds in a single step by means of affordable reagents and using a non-contaminant solvent such as water. Thus, the syntheses of the novel Se derivatives of ASA, Ind, Nap, Ket, and Ibup were carried out by the reaction of the corresponding NSAID acyl chloride with an aqueous solution of sodium hydrogen selenide (NaHSe) formed in situ. In the case of NSAIDs Ind, Ket, and Ibup, the acyl chloride was not commercially available and had to be previously obtained from carboxylic acid. The synthesis route of these derivatives (15) is described in Scheme 1. The compounds were purified by silica gel column chromatography or flash chromatography with a gradient of hexane and ethyl acetate as eluent. All the diacyl diselenide NSAID derivatives were obtained with feasible yields (20–50%) in a one-pot synthesis.
The structures of all the compounds were confirmed by HRMS and by 1H, 13C, and 77Se NMR, as described in the previous section. The purity of all the compounds was determined by qNMR spectroscopy (see Supplementary Material). Concerning 1H NMR, a shift in the H peak in the groups close to the diacyl diselenide moiety could be observed due to the modification in the chemical environment of the carbonyl group. Thus, when the carbonyl of the diacyl diselenide moiety is attached to a methylene group as in the case of compound 2, or to a methine group for the arylpropionic acid derivatives 3, 4, and 5, a shift downfield of the corresponding H signal was found compared to the position of the same peak in the parent NSAID spectra. The methylene protons of compound 2 appeared at 4.00 ppm (Figure S6), whereas the same methylene attached to the carboxylic acid group of Ind appeared at 3.69 ppm. Likewise, the peak belonging to the single H of the methine group linked to the carbonyl appeared in the range of 3.98–4.17 ppm for compounds 3, 4, and 5, while these signals appeared to be shifted upfield at 3.69–3.88 ppm in the 1H NMR spectra of the parent NSAIDs Nap, Ket, and Ibup, respectively, although the J values in all the cases were maintained at 7.2 Hz. Additionally, the two characteristic peaks of the ‘para’ substitution in the benzene ring of Ibup (7.09 and 7.21 ppm) merged into one signal at 7.09 ppm for compound 5 (Figure S15) marked as a multiplet. Inspection of 13C NMR spectra also revealed the displacement of the position of the group linked to the carbonyl of the diacyl diselenide moiety. As in the case of 1H NMR spectra, the peaks belonging to the methylene or the methine groups appeared to shift downfield when compared to the position of the signals in the NSAID spectra. For compound 2, the methylene group appeared at 41.2 ppm (Figure S7), and this same peak appeared at 30.1 ppm in the case of Ind. The chemical environment of the methine group in Nap, Ket, and Ibup is highly similar, and so the signals in the corresponding 13C NMR spectra appeared almost in the same position, around 45.3 ppm. As expected, the incorporation of the diacyl diselenide moiety also modified the position of the methine in the same manner, these signals being shifted downfield in the range of 56.4–59.1 ppm. When comparing the new Se-NSAID derivatives with their corresponding parent NSAIDs, an evident shift downfield in the position of the carbonyl group attached to the double Se bond is observed. In this regard, the carboxylic acid groups of the parent NSAIDs appeared in the range of 170.3–181.1 ppm, ASA being the NSAID with the peak appearing as the most upfield. Accordingly, the carbonyl of the Se-NSAID derivatives showed an interval of 184.5–198.7 ppm, in which compound 1 (Figure S4) was the only derivative with the signal appearing below 190 ppm. This could probably be due to the aromatic ring attached to the carbonyl of 1 in comparison with the rest of the NSAID derivatives included in this study. Regarding 77Se NMR, overall, the Se peak for the diacyl diselenides appeared as one sharp peak in the range of 634–792 ppm. Interestingly, in the case of derivatives 4 (Figure S14) and 5 (Figure S17), two peaks located in close proximity to one another could be defined separately, corresponding to the two Se atoms of the double bond. These chemical shifts observed for the Se atom are in concordance with other ranges reported in the literature concerning diacyl diselenides [75]. Thus, compounds 3 and 4 showed a similar chemical shift for the Se peak in the range of 634–636 ppm, whereas the Se signals of derivatives 1 and 2 appeared slightly shifted downfield to 659 and 657 ppm, respectively, despite not having the same chemical environment as in the case of 3 and 4. Likewise, the Se peaks of compound 5 appeared downfield at 791 and 792 ppm, although the chemical environment of the Se was expected to be similar to that of the other Se-NSAID derivatives of arylpropionic acids 3 and 4.

3.3. Evaluation of the Radical Scavenging Activity of the Se-NSAID Derivatives

We first evaluated the radical scavenging capacity of the diacyl diselenides as a preliminary study to assess their potential antioxidant activity. Thus, the DPPH assay was performed for all the compounds at two different concentrations (0.06 and 0.03 mg/mL) and measured at eight different time points. Asc and Trolox were used as positive controls. The eight different time points and the two different concentrations were selected with the aim of achieving a significant DPPH radical scavenging activity for the diacyldiselenides. The same conditions, for comparison purposes, were used for the positive controls (Asc and Trolox), even though the maximal scavenging activity was demonstrated by both positive controls for each time point and concentration. The results obtained are summarized in Figure 3 and presented as the percentage of DPPH radical scavenging activity.
As shown in Figure 3, at the lowest concentration tested (Figure 3A), the compounds exhibited in general low activity, with the sole exception of compound 5 having an antioxidant capacity of 20%. Nevertheless, the Se-NSAID derivatives displayed an increase in their antioxidant activity over time, this effect is clearly observed at the highest concentration tested (Figure 3B). Thus, compounds 14 showed values of around 20–30% by the end of the experiment, surpassing the radical scavenging activity of other reported selenocompounds with therapeutic properties [22,66,78]. Remarkably, compound 5 was found to display the most potent radical scavenging activity among all the Se-NSAID derivatives, achieving values even comparable to those of the positive controls by the end of the assay (Figure 3B). Interestingly, no antioxidant activity could be ascribed to several selenocompounds that included the NSAID scaffolds [66], so the effect observed for these derivatives could be highly related to the presence of the diacyl diselenide group in their structures.
To the best of our knowledge, the results outlined in Figure 3 represent the first report of the antioxidant capacity of this type of molecule, therefore revealing the diacyl diselenide moiety as a new scaffold worth further exploration for developing compounds with antioxidant properties.

3.4. Evaluation of the Antiproliferative Activities of the Se-NSAID Derivatives

As a first approach to evaluate the possible anticancer properties of the new diacyl diselenides, the effect of the Se-NSAID derivatives on the cell viability was determined in vitro across a panel of cancer cell lines using the MTT assay as previously described [37]. The antiproliferative activity of compounds 15 was thereby assessed against two cell lines each of four different types of cancer (colon, prostate, breast, and lung) at seven concentrations ranging between 1 and 100 µM for 48 h. The calculated IC50 values are outlined in Table 1. The GI50, TGI, and LD50 values were also determined and are included in Table S1. The IC50 values obtained for each parent NSAID (ASA, Ind, Nap, Ket, and Ibup) are included for reference.
As reported in Table 1, compounds 15 displayed distinct antiproliferative profiles. Since the difference between compounds lies in the NSAID scaffold attached to the diacyl diselenide group at both sides of the molecule, it can be implied that the NSAID is a determinant for the antiproliferative activity. In this context, ASA derivative 1 showed potent antiproliferative activity in one of the colon cancer cell lines (HT-29) and in a breast cancer cell line (T-47D). Additionally, 1 showed IC50 values below 10 µM in both the prostate cancer cells (DU-145 and PC-3) and the lung cancer cells (H1299 and A549). In contrast, 1 presents relatively less antiproliferative effect towards triple-negative breast cancer cells (MDA-MB-231) and particularly to HCT-116 colon cancer cells, with an IC50 value of almost 40 µM. Compound 4 including a Ket scaffold displayed a similar profile as compound 1, as it showed higher antiproliferative activity towards certain cell lines. In this case, 4 was active in at least one of the cell lines for the different types of cancer tested, with low IC50 values in HCT-116, PC-3, MDA-MB-231, and A549 cells. Interestingly, Ind derivative 2 showed no activity against all the eight cell lines evaluated, this being the only compound in which the inclusion of Se did not improve the antiproliferative capacity of the parent NSAID. The GI50, TGI, and LD50 values calculated also demonstrated poor activity of compound 2 in comparison with the rest of the series (Table S1).
Compound 3, featuring a Nap scaffold, showed potent activity in both the prostate (DU-145 and PC-3) and breast (MDA-MB-231 and T-47D) cancer cell lines, with an average IC50 value of around 5 µM in these cells. On the contrary, this compound displayed selective activity in HCT-116 and H1299 cells, with no activity at all in the other cancer cells of the same type (HT-29 and A549, respectively). Ibup derivative 5 was the only compound among the Se-NSAID derivatives which displayed potent antiproliferative activity in all the cancer cells tested, with IC50 values below 10 µM in every colon, prostate, breast, and lung cancer cell line assessed. The GI50, TGI, and LD50 values obtained for this compound showed the same trend observed for the IC50 values (Table S1). Interestingly, only compounds 3 and 5 showed GI50 values below 10 µM across all the cancer cell lines tested, but compound 5 was the only Se-NSAID derivative with GI50 values even below 5 µM in at least six out of the eight cell lines.
Since unwanted side effects are usually associated with a lack of selectivity in the antiproliferative activity towards healthy cells, we then decided to test the Se-NSAID derivatives in mammary gland (184B5) and bronchial epithelium (BEAS-2B) nonmalignant cell lines. Thus, the selectivity indexes (SI) were determined as the ratio of the IC50 values obtained for the nonmalignant cells and the homolog breast and lung cancer cells. The IC50 values obtained in these nonmalignant cell lines and the calculated SI values are summarized in Table 2, together with the results already shown for MDA-MB-231, T-47D, H1299, and A549 cells for clarity. The GI50, TGI, and LD50 values were also calculated in the nonmalignant cells and are included in Table S1.
As shown in Table 2, compounds 1 and 5 showed the highest IC50 values among the Se-NSAID derivatives in the range of 14–19 µM in the nonmalignant cell lines. Hence, 1 showed SI values > 2 in the lung cells, but it lacked selectivity in breast cells when comparing the triple-negative breast cancer cell line (MDA-MB-231) with the nonmalignant cells (184B5). On the contrary, compound 5 displayed significant selectivity towards all the cell lines tested, with fair SI values of almost five- and four-fold in breast and lung cancer cells, respectively. Moreover, 5 was not only the best compound in terms of IC50 values in breast cancer cells, but it also inhibited cell growth without leading to cell death in breast nonmalignant cells (Table S1). On the other hand, similar to the antiproliferative activity tested in cancer cells, compound 2 was completely inactive in the two nonmalignant cell lines, and no IC50, GI50, TGI, and LD50 values could be determined (Table 2 and Table S1). Compound 3 was fairly selective in breast cells with SI values around two in both cases, but it showed poor selectivity in lung cancer cells, especially since this compound was not active in A549 cells. Compound 4 displayed IC50 values around 5 µM in the two nonmalignant cells, and thus was the less selective derivative among the compounds, mainly in lung cancer cells (SI values < 0.5). Therefore, compound 5 was identified as the most effective and selective Se-NSAID derivative overall.
Although the number of synthesized compounds is relatively small to make solid structure-activity relationships, some conclusions can be drawn from an overview analysis of the IC50 and SI values outlined in Table 1 and Table 2: (1) the Se atom incorporated in the novel form of a diacyl diselenide group into NSAID derivatives is a valid approach to obtain potent antiproliferative compounds; (2) the NSAID is also determinant for the biological activity, as the presence of an Ind scaffold at both sides of the diacyl diselenide group resulted in an inactive compound 2 with no activity in the panel of cancer cell lines tested; (3) NSAID Ibup yielded the best NSAID derivative in terms of both the activity and selectivity, as it stood out with a good selectivity towards cancer cells while evincing the most potent antiproliferative activity in all the cancer cell lines tested.
Considering the selectivity and anticancer activity displayed by compound 5, this Ibup derivative was selected for further biological studies. This compound showed outstanding antiproliferative activity with IC50 values below 10 µM in all the cancer cell lines evaluated (Table 1). Additionally, 5 proved to be the most selective compound among all the Se-NSAID derivatives (Table 2). Interestingly, compound 5 also displayed the most potent antioxidant activity in the DPPH radical scavenging assay (Figure 3), thus suggesting a dual behavior for this analog.

3.5. Compound 5 Induced Cell Death by Partially Triggering Apoptosis in Colon Cancer Cells

The assays performed on cell viability evinced that compound 5 showed potent antiproliferative activity. Thus, we evaluated this derivative to further confirm if 5 induced cell death through apoptosis. Since sequential activation of caspases is essential for beginning the process of apoptosis [79], we tested the effect of compound 5 on caspase-3 and caspase-7 with the Caspase 3/7 assay. The Annexin V & Dead Cell assay was also performed to complement the results.
NSAIDs have been widely studied and their chemopreventive activity in colon cancer has been well documented in the literature [48,49]. Several studies also reported their therapeutic effect and their mechanism of action in this cancer [47,58,80,81], which encouraged us to select the colon cancer HCT-116 cell line as representative cells to focus our experiments on. As shown in Table 1 and Table S1, compound 5 showed potent antiproliferative activity in all the cancer cell lines tested, including HCT-116. Therefore, HCT-116 cells were treated either with DMSO (vehicle) or compound 5 for 24 h at five different concentrations ranging from 2.5 to 25 µM. The results were obtained following the manufacturer’s protocol and are shown in Figure 4. In addition, a quantitative comparison of the difference in the cell population induced by compound 5 is included in Figure 4C,D.
In general, cells treated only with DMSO were located as expected mainly in the lower left quadrant in both assays. The Annexin V & Dead Cell assay evinced that the treatment with compound 5 induced a shift from healthy cells towards an apoptotic state (Figure 4A,C). However, even at a concentration close to the IC50 value (10 µM), the cells remained mainly located in the lower left quadrant after 24 h and no discernible difference was observed compared to the control cells. As the concentration increased up to 15 µM, the live cell population decreased to ~63%, while the percentage of early apoptotic cells was consequently increased to ~23%. The population of cells in a late apoptosis state was also enhanced to ~13% (Figure 4A). Interestingly, at 25 µM, the percentage of viable cells treated with 5 was maintained at a value close to ~60%, while more than 30% of cells were detected to be apoptotic. In this context, a slight population of cells shifted to a necrotic or late apoptotic state, as there were ~10% or ~20% of cells, respectively, with 7-AAD positive (Figure 4A). Nevertheless, the presence of outer PS residues determined by the Annexin V & Dead Cell assay seems to be insufficient to explain the potent antiproliferative capacity observed in the cell viability assays (Table 1 and Table S1).
Considering the Caspase 3/7 assay (Figure 4B,D), the percentage of the live cell population was maintained even at high concentrations of compound 5 with almost no variations, although at 15 and 25 µM a significant change could be observed (Figure 4D). There was only a slight increase in the population of dead cells up to ~11% eventually at 25 µM, but the distribution of cells in an apoptotic state was not significantly different with respect to the control cells treated only with DMSO. Thus, these results imply that compound 5 might not exert its anticancer activity through significant involvement in the caspase 3/7 activity.
Therefore, the results shown in Figure 4, along with the anticancer activity displayed by compound 5 in HCT-116 cells (Table 1), could suggest that apoptosis is not the only mechanism by which this compound induces cell death in this type of cancer cells.

3.6. Compound 5 Effectively Inhibited Tumor Growth in a Colon Cancer Xenograft Model without Any Apparent Systemic Toxicity

The overall appealing profile of compound 5 regarding its antiproliferative activity in vitro prompted us to further assess its toxicity and efficacy in vivo.
For chronic toxicity assessment, the mice were given a dose of 5 (7.5 mg/Kg body weight) intraperitoneally (IP) thrice per week, and related changes in body weight, behavior, food, and water uptake were monitored for almost a month. The blood serum was also collected after the sacrifice and analyzed. The results are displayed in Figure 5. As shown in Figure 5A, the blood biochemistry profile measured at the end of the treatment with 5 revealed no significant changes in the compound-treated mice compared to the vehicle-treated animals. Moreover, the body weight of the mice was maintained with no appreciable variations during the treatment with 5 (Figure 5B) and we observed no changes in behavior and food/water uptake.
The low toxicity of compound 5 warranted the evaluation of its antitumor activity in a subcutaneous xenograft tumor mouse model bearing HCT-116 colon cancer cells. Nude mice were randomized and treated IP with compound 5 (10 mg/Kg body weight) for 23 days. The results of the in vivo efficacy assay are outlined in Figure 6. A higher dose for the efficacy study in vivo was chosen based on the favorable results at 7.5 mg/Kg (Figure 5). The TGI and the T/C parameters were calculated at the end of the treatment to assess the antitumor effects in terms of tumor volume and tumor weight, respectively.
As shown in Figure 6, compound 5 significantly inhibited tumor growth in the colon cancer xenograft model. The pronounced inhibitory effect was maintained throughout the duration of the treatment (Figure 6A), and on day 38 the TGI reached 72% compared to the vehicle-treated mice. Furthermore, the comparison of the weight of the tumors excised at the end of the treatment (Figure 6B) with a good T/C ratio of 38% also evinced the significant inhibition of tumor growth by compound 5. Besides, a slight change in the body weight of the mice was observed in the 5-treated group compared to the vehicle-treated mice during the treatment (Figure 6C), but no evident toxicity signs as morphological abnormalities in major organs were detected. Taken together, our in vivo results demonstrated that compound 5 remarkably produced a significant inhibition of the tumor growth in a colon cancer xenograft model with no noticeable toxicity, suggesting that 5 is a promising candidate worthy of further preclinical studies. Nevertheless, the in vivo efficacy of compound 5 seems to present a limited off-target effect considering its different behavior in vitro and in vivo. This discrepancy could be attributed to several factors, such as instability of the selenocompound or inadequate ADME properties. Thus, different approaches can be explored to further develop this new class of anticancer agents, including structural modifications to improve solubility and biodistribution features, or employing a targeted delivery to both increase the stability and ADME properties.

4. Conclusions

In summary, we synthesized for the first time a small library of NSAID derivatives containing Se in the novel form of a diacyl diselenide group as precursors with the aim of retaining the Se moiety in the NSAID scaffolds. The work presented here evinces that this design based on the incorporation of Se as diacyl diselenide moieties into traditional NSAIDs, led to analogs with potent antioxidant and/or anticancer activity. To our knowledge, this constitutes the first report of the appealing in vitro and in vivo therapeutic profile of diacyl diselenide-based compounds. The results described herein demonstrate that our design approach to maintain the entire NSAIDs along with Se provides a feasible new frame to develop novel biologically active agents. Compound 5 was identified as the most potent antiproliferative agent among the Se-NSAID derivatives with IC50 values below 10 µM in a panel of eight cancer cell lines in vitro. Furthermore, 5 proved to be the most selective analog when compared to the nonmalignant cell lines. This compound also displayed the most potent antioxidant activity within the series, achieving values comparable to those of the positive controls as demonstrated in the DPPH radical scavenging assay. The capability of 5 to induce cell death through apoptosis was also evaluated in HCT-116 colon cells. Interestingly, no clear enhancement on the caspase 3/7 activity was observed, but 5 was found to partially induce apoptosis determined by an Annexin V assay. Notably, compound 5 also markedly inhibited the tumor growth in an HCT-116 colon cancer xenograft model at a dose of 10 mg/Kg administrated IP, with an impressive TGI value of 72% and a T/C value of 38%. Overall, these findings qualify compound 5 and the diacyl diselenide functionality as appealing candidates to explore the potential and utility of diacyl diselenide-based derivatives as leads to design precursors for bioisosteric selenocompounds as therapeutic agents, especially in colon cancer treatment.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/antiox12091666/s1, Table S1: GI50, TGI and LD50 values; Figures S1–S20: 1H, 13C and 77Se NMR spectra, qNMR spectra; Figure S21: Mass spectrum of compound 5.

Author Contributions

Conceptualization, D.P. and C.S.; investigation and data curation, S.R.-I.; animal studies, C.A. (Cesar Aliaga); formal analysis, S.R.-I., I.E., A.R. and D.P.; writing—original draft preparation, S.R.-I.; writing—review and editing, A.K.S., C.A. (Carlos Aydillo), N.M.-S., D.P. and C.S.; supervision, A.K.S., D.P. and C.S.; funding acquisition, A.K.S., D.P. and C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Plan de Investigación de la Universidad de Navarra, PIUNA (2018–19), and the Department of Pharmacology and Penn State Cancer Institute of the Penn State College of Medicine.

Institutional Review Board Statement

The animal study protocol was approved by the Animal Care and Use Committee of the Penn State College of Medicine (protocol code: PROTO202102111—Development of Small Molecules Therapeutics, and date of approval: 22 August 2022).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon reasonable request from the corresponding authors.

Acknowledgments

Sandra Ramos-Inza acknowledges the support of the FPU program from the Spanish Ministry of Universities for a Ph.D. fellowship (FPU18/04679) and a mobility grant (EST19/00898).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ASA, aspirin; Asc, ascorbic acid; ATCC, American Type Culture Collection; COX, cyclooxygenase; DPPH, 2,2-diphenyl-1-picrylhydrazyl; FAP, familial adenomatous polyposis; FBS, fetal bovine serum; HRMS, high-resolution mass spectrometry; Ibup, ibuprofen; Ind, indomethacin; IP, intraperitoneally; Ket, ketoprofen; MMR, DNA mismatch repair; mp, melting point; MTT, 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide; Nap, naproxen; NMR, nuclear magnetic resonance; NSAIDs, nonsteroidal anti-inflammatory drugs; PG, prostaglandin, Se, selenium; SI, selectivity index; T/C, relative increment ratio; TGI, tumor growth inhibition; TLC, thin-layer chromatography.

References

  1. Chuai, H.; Zhang, S.Q.; Bai, H.; Li, J.; Wang, Y.; Sun, J.; Wen, E.; Zhang, J.; Xin, M. Small molecule selenium-containing compounds: Recent development and therapeutic applications. Eur. J. Med. Chem. 2021, 223, 113621. [Google Scholar] [PubMed]
  2. Barchielli, G.; Capperucci, A.; Tanini, D. The role of selenium in pathologies: An updated review. Antioxidants 2022, 11, 251. [Google Scholar] [CrossRef]
  3. Steinbrenner, H.; Speckmann, B.; Klotz, L.O. Selenoproteins: Antioxidant selenoenzymes and beyond. Arch. Biochem. Biophys. 2016, 595, 113–119. [Google Scholar] [CrossRef] [PubMed]
  4. Gandin, V.; Khalkar, P.; Braude, J.; Fernandes, A.P. Organic selenium compounds as potential chemotherapeutic agents for improved cancer treatment. Free Radic. Biol. Med. 2018, 127, 80–97. [Google Scholar]
  5. Kim, S.J.; Choi, M.C.; Park, J.M.; Chung, A.S. Antitumor effects of selenium. Int. J. Mol. Sci. 2021, 22, 11844. [Google Scholar] [CrossRef] [PubMed]
  6. Vernia, F.; Longo, S.; Stefanelli, G.; Viscido, A.; Latella, G. Dietary factors modulating colorectal carcinogenesis. Nutrients 2021, 13, 143. [Google Scholar] [CrossRef] [PubMed]
  7. Cai, X.; Wang, C.; Yu, W.; Fan, W.; Wang, S.; Shen, N.; Wu, P.; Li, X.; Wang, F. Selenium exposure and cancer risk: An updated meta-analysis and meta-regression. Sci. Rep. 2016, 6, 19213. [Google Scholar] [CrossRef]
  8. Vinceti, M.; Filippini, T.; Cilloni, S.; Bargellini, A.; Vergoni, A.V.; Tsatsakis, A.; Ferrante, M. Health risk assessment of environmental selenium: Emerging evidence and challenges (Review). Mol. Med. Rep. 2017, 15, 3323–3335. [Google Scholar] [CrossRef]
  9. Kuršvietienė, L.; Mongirdienė, A.; Bernatonienė, J.; Šulinskienė, J.; Stanevičienė, I. Selenium anticancer properties and impact on cellular redox status. Antioxidants 2020, 9, 80. [Google Scholar] [CrossRef]
  10. Wrobel, J.K.; Power, R.; Toborek, M. Biological activity of selenium: Revisited. IUBMB Life 2016, 68, 97–105. [Google Scholar] [CrossRef]
  11. Song, M.; Kumaran, M.N.; Gounder, M.; Gibbon, D.G.; Nieves-Neira, W.; Vaidya, A.; Hellmann, M.; Kane, M.P.; Buckley, B.; Shih, W.; et al. Phase I trial of selenium plus chemotherapy in gynecologic cancers. Gynecol. Oncol. 2018, 150, 478–486. [Google Scholar] [CrossRef] [PubMed]
  12. Stolwijk, J.M.; Garje, R.; Sieren, J.C.; Buettner, G.R.; Zakharia, Y. Understanding the redox biology of selenium in the search of targeted cancer therapies. Antioxidants 2020, 9, 420. [Google Scholar] [CrossRef]
  13. Ruberte, A.C.; Sanmartin, C.; Aydillo, C.; Sharma, A.K.; Plano, D. Development and therapeutic potential of selenazo compounds. J. Med. Chem. 2020, 63, 1473–1489. [Google Scholar] [CrossRef] [PubMed]
  14. Lin, T.; Ding, Z.; Li, N.; Xu, J.; Luo, G.; Liu, J.; Shen, J. Seleno-cyclodextrin sensitises human breast cancer cells to TRAIL-induced apoptosis through DR5 induction and NF-κB suppression. Eur. J. Cancer 2011, 47, 1890–1907. [Google Scholar] [CrossRef]
  15. Zeng, H.; Cheng, W.H.; Johnson, L.K. Methylselenol, a selenium metabolite, modulates p53 pathway and inhibits the growth of colon cancer xenografts in Balb/c mice. J. Nutr. Biochem. 2013, 24, 776–780. [Google Scholar] [CrossRef] [PubMed]
  16. Radomska, D.; Czarnomysy, R.; Radomski, D.; Bielawski, K. Selenium compounds as novel potential anticancer agents. Int. J. Mol. Sci. 2021, 22, 1009. [Google Scholar] [CrossRef]
  17. Ali, W.; Benedetti, R.; Handzlik, J.; Zwergel, C.; Battistelli, C. The innovative potential of selenium-containing agents for fighting cancer and viral infections. Drug Discov. Today 2021, 26, 256–263. [Google Scholar] [CrossRef]
  18. Genchi, G.; Lauria, G.; Catalano, A.; Sinicropic, M.S.; Carocci, A. Biological activity of selenium and its impact on human health. Int. J. Mol. Sci. 2023, 24, 2633. [Google Scholar] [CrossRef]
  19. Álvarez-Pérez, M.; Ali, W.; Marć, M.A.; Handzlik, J.; Domínguez-Álvarez, E. Selenides and diselenides: A review of their anticancer and chemopreventive activity. Molecules 2018, 23, 628. [Google Scholar] [CrossRef]
  20. Pang, Y.; Lin, H.; Ou, C.; Cao, Y.; An, B.; Yan, J.; Li, X. Design, synthesis, and biological evaluation of novel benzodiazepine derivatives as anticancer agents through inhibition of tubulin polymerization in vitro and in vivo. Eur. J. Med. Chem. 2019, 182, 111670. [Google Scholar] [CrossRef]
  21. Ali, W.; Spengler, G.; Kincses, A.; Nové, M.; Battistelli, C.; Latacz, G.; Starek, M.; Dąbrowska, M.; Honkisz-Orzechowska, E.; Romanelli, A.; et al. Discovery of phenylselenoether-hydantoin hybrids as ABCB1 efflux pump modulating agents with cytotoxic and antiproliferative actions in resistant T-lymphoma. Eur. J. Med. Chem. 2020, 200, 112435. [Google Scholar] [CrossRef]
  22. Ramos-Inza, S.; Henriquez-Figuereo, A.; Moreno, E.; Berzosa, M.; Encío, I.; Plano, D.; Sanmartín, C. Unveiling a new selenocyanate as a multitarget candidate with anticancer, antileishmanial and antibacterial potential. Molecules 2022, 27, 7477. [Google Scholar] [CrossRef]
  23. Huang, Y.M.; Cheng, Y.; Peng, Z.N.; Pang, L.P.; Li, J.Y.; Xiao, J.A.; Zhang, Y.F.; Cui, J.G. Synthesis and antitumor activity of some cholesterol-based selenocyanate compounds. Steroids 2023, 194, 109217. [Google Scholar] [CrossRef]
  24. Nie, Y.; Zhong, M.; Li, S.; Li, X.; Zhang, Y.; Zhang, Y.; He, X. Synthesis and potential anticancer activity of some novel selenocyanates and diselenides. Chem. Biodivers. 2020, 17, e1900603. [Google Scholar] [CrossRef]
  25. Roldán-Peña, J.M.; Puerta, A.; Dinić, J.; Jovanović Stojanov, S.; González-Bakker, A.; Hicke, F.J.; Mishra, A.; Piyasaengthong, A.; Maya, I.; Walton, J.W.; et al. Biotinylated selenocyanates: Potent and selective cytostatic agents. Bioorg. Chem. 2023, 133, 106410. [Google Scholar] [PubMed]
  26. Huang, Y.; Wei, M.; Peng, Z.; Cheng, Y.; Zhang, Y.; Li, J.; Xiao, J.; Gan, C.; Cui, J. Synthesis of estrone selenocyanate compounds, anti-tumor activity evaluation and structure-activity relationship analysis. Bioorg. Med. Chem. 2022, 76, 117086. [Google Scholar] [CrossRef]
  27. Crampsie, M.A.; Jones, N.; Das, A.; Aliaga, C.; Desai, D.; Lazarus, P.; Amin, S.; Sharma, A.K. Phenylbutyl isoselenocyanate modulates phase I and II enzymes and inhibits 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone-induced DNA adducts in mice. Cancer Prev. Res. 2011, 4, 1884–1894. [Google Scholar] [CrossRef] [PubMed]
  28. Wu, W.; Karelia, D.; Pramanik, K.; Amin, S.G.; Sharma, A.K.; Jiang, C.; Lu, J.X. Phenylbutyl isoselenocyanate induces reactive oxygen species to inhibit androgen receptor and to initiate p53-mediated apoptosis in LNCaP prostate cancer cells. Mol. Carcinog. 2018, 57, 1055–1066. [Google Scholar] [CrossRef]
  29. Milczarek, M.; Cierpiał, T.; Kiełbasiński, P.; Małecka-Giełdowska, M.; Świtalska, M.; Wietrzyk, J.; Mazur, M.; Wiktorska, K. An organofluorine isoselenocyanate analogue of sulforaphane affects antimetabolite 5-fluorouracil’s anticancer activity: A perspective for new combinatory therapy in triple-negative breast cancer. Molecules 2023, 28, 5808. [Google Scholar] [CrossRef]
  30. Marković, S.B.; Maciejewska, N.; Olszewski, M.; Višnjevac, A.; Puerta, A.; Padrón, J.M.; Novaković, I.; Kojić, S.; Fernandes, H.S.; Sousa, S.F.; et al. Study of the anticancer potential of Cd complexes of selenazoyl-hydrazones and their sulfur isosters. Eur. J. Med. Chem. 2022, 238, 114449. [Google Scholar] [CrossRef] [PubMed]
  31. Chen, Z.; Li, D.; Xu, N.; Fang, J.; Yu, Y.; Hou, W.; Ruan, H.; Zhu, P.; Ma, R.; Lu, S.; et al. Novel 1,3,4-selenadiazole-containing kidney-type glutaminase inhibitors showed improved cellular uptake and antitumor activity. J. Med. Chem. 2019, 62, 589–603. [Google Scholar] [CrossRef]
  32. Radomska, D.; Czarnomysy, R.; Szymanowska, A.; Radomski, D.; Domínguez-Álvarez, E.; Bielawska, A.; Bielawski, K. Novel selenoesters as a potential tool in triple-negative breast cancer treatment. Cancers 2022, 14, 4304. [Google Scholar] [CrossRef]
  33. Khalkar, P.; Díaz-Argelich, N.; Antonio Palop, J.; Sanmartín, C.; Fernandes, A.P. Novel methylselenoesters induce programed cell death via entosis in pancreatic cancer cells. Int. J. Mol. Sci. 2018, 19, 2849. [Google Scholar] [CrossRef] [PubMed]
  34. Szemerédi, N.; Dobiasová, S.; Salardón-Jiménez, N.; Kincses, A.; Nové, M.; Habibullah, G.; Sevilla-Hernández, C.; Benito-Lama, M.; Alonso-Martínez, F.J.; Viktorová, J.; et al. Cyano- and ketone-containing selenoesters as multi-target compounds against resistant cancers. Cancers 2021, 13, 4563. [Google Scholar] [CrossRef]
  35. Moussa, Z.; Kaddoura, R.; Saadeh, H.A.; Abutaha, N.; Ahmed, S.A. Highly bioactive novel aryl-, benzyl-, and piperazine-selenoureas: Synthesis, structural characterization and in vitro biological evaluation. Heliyon 2022, 8, e10709. [Google Scholar] [CrossRef]
  36. Barbosa, F.A.R.; Siminski, T.; Canto, R.F.S.; Almeida, G.M.; Mota, N.; Ourique, F.; Pedrosa, R.C.; Braga, A.L. Novel pyrimidinic selenourea induces DNA damage, cell cycle arrest, and apoptosis in human breast carcinoma. Eur. J. Med. Chem. 2018, 155, 503–515. [Google Scholar] [CrossRef]
  37. Ruberte, A.C.; Ramos-Inza, S.; Aydillo, C.; Talavera, I.; Encío, I.; Plano, D.; Sanmartín, C. Novel N,N’-disubstituted acylselenoureas as potential antioxidant and cytotoxic agents. Antioxidants 2020, 9, 55. [Google Scholar] [CrossRef]
  38. De Franco, M.; Saab, M.; Porchia, M.; Marzano, C.; Nolan, S.P.; Nahra, F.; Van Hecke, K.; Gandin, V. Unveiling the potential of innovative gold(I) and silver(I) selenourea complexes as anticancer agents targeting TrxR and cellular redox homeostasis. Chemistry 2022, 28, e202201898. [Google Scholar] [CrossRef]
  39. Zhou, J.N.; Zhang, B.; Wang, H.Y.; Wang, D.X.; Zhang, M.M.; Zhang, M.; Wang, X.K.; Fan, S.Y.; Xu, Y.C.; Zeng, Q.; et al. A functional screening identifies a new organic selenium compound targeting cancer stem cells: Role of c-Myc transcription activity inhibition in liver cancer. Adv. Sci. 2022, 9, e2201166. [Google Scholar] [CrossRef] [PubMed]
  40. Krasowska, D.; Iraci, N.; Santi, C.; Drabowicz, J.; Cieslak, M.; Kaźmierczak-Barańska, J.; Palomba, M.; Królewska-Golińska, K.; Magiera, J.; Sancineto, L. Diselenides and benzisoselenazolones as antiproliferative agents and glutathione-S-transferase inhibitors. Molecules 2019, 24, 2914. [Google Scholar] [CrossRef] [PubMed]
  41. Ramos-Inza, S.; Ruberte, A.C.; Sanmartín, C.; Sharma, A.K.; Plano, D. NSAIDs: Old acquaintance in the pipeline for cancer treatment and prevention─Structural modulation, mechanisms of action, and bright future. J. Med. Chem. 2021, 64, 16380–16421. [Google Scholar] [CrossRef]
  42. Zappavigna, S.; Cossu, A.M.; Grimaldi, A.; Bocchetti, M.; Ferraro, G.A.; Nicoletti, G.F.; Filosa, R.; Caraglia, M. Anti-inflammatory drugs as anticancer agents. Int. J. Mol. Sci. 2020, 21, 2605. [Google Scholar] [CrossRef]
  43. Sauer, C.M.; Myran, D.T.; Costentin, C.E.; Zwisler, G.; Safder, T.; Papatheodorou, S.; Mucci, L.A. Effect of long term aspirin use on the incidence of prostate cancer: A systematic review and meta-analysis. Crit. Rev. Oncol. Hematol. 2018, 132, 66–75. [Google Scholar] [CrossRef]
  44. Simon, T.G.; Duberg, A.S.; Aleman, S.; Chung, R.T.; Chan, A.T.; Ludvigsson, J.F. Association of aspirin with hepatocellular carcinoma and liver-related mortality. N. Engl. J. Med. 2020, 382, 1018–1028. [Google Scholar] [CrossRef]
  45. Ma, S.; Guo, C.; Sun, C.; Han, T.; Zhang, H.; Qu, G.; Jiang, Y.; Zhou, Q.; Sun, Y. Aspirin use and risk of breast cancer: A meta-analysis of observational studies from 1989 to 2019. Clin. Breast Cancer 2021, 21, 552–565. [Google Scholar] [CrossRef] [PubMed]
  46. Liu, J.; Zheng, F.; Yang, M.; Wu, X.; Liu, A. Effect of aspirin use on survival benefits of breast cancer patients: A meta-analysis. Medicine 2021, 100, e26870. [Google Scholar] [CrossRef]
  47. Mohammed, A.; Yarla, N.S.; Madka, V.; Rao, C.V. Clinically relevant anti-inflammatory agents for chemoprevention of colorectal cancer: New perspectives. Int. J. Mol. Sci. 2018, 19, 2332. [Google Scholar] [CrossRef]
  48. Kuo, C.N.; Pan, J.J.; Huang, Y.W.; Tsai, H.J.; Chang, W.C. Association between nonsteroidal anti-inflammatory drugs and colorectal cancer: A population-based case-control study. Cancer Epidemiol. Biomarkers Prev. 2018, 27, 737–745. [Google Scholar] [CrossRef]
  49. Tomić, T.; Domínguez-López, S.; Barrios-Rodríguez, R. Non-aspirin non-steroidal anti-inflammatory drugs in prevention of colorectal cancer in people aged 40 or older: A systematic review and meta-analysis. Cancer Epidemiol. 2019, 58, 52–62. [Google Scholar] [CrossRef] [PubMed]
  50. Figueiredo, J.C.; Jacobs, E.J.; Newton, C.C.; Guinter, M.A.; Cance, W.G.; Campbell, P.T. Associations of aspirin and non-aspirin non-steroidal anti-inflammatory drugs with colorectal cancer mortality after diagnosis. J. Natl. Cancer Inst. 2021, 113, 833–840. [Google Scholar] [CrossRef] [PubMed]
  51. Zhang, Y.; Chan, A.T.; Meyerhardt, J.A.; Giovannucci, E.L. Timing of aspirin use in colorectal cancer chemoprevention: A prospective cohort study. J. Natl. Cancer Inst. 2021, 113, 841–851. [Google Scholar] [CrossRef]
  52. Ait Ouakrim, D.; Dashti, S.G.; Chau, R.; Buchanan, D.D.; Clendenning, M.; Rosty, C.; Winship, I.M.; Young, J.P.; Giles, G.G.; Leggett, B.; et al. Aspirin, ibuprofen, and the risk of colorectal cancer in Lynch Syndrome. J. Natl. Cancer Inst. 2015, 107, djv170. [Google Scholar] [CrossRef]
  53. Bowen, C.M.; Walter, L.; Borras, E.; Wu, W.; Ozcan, Z.; Chang, K.; Bommi, P.V.; Taggart, M.W.; Thirumurthi, S.; Lynch, P.M.; et al. Combination of sulindac and bexarotene for prevention of intestinal carcinogenesis in familial adenomatous polyposis. Cancer Prev. Res. 2021, 14, 851–862. [Google Scholar] [CrossRef]
  54. Grahn, O.; Lundin, M.; Chapman, S.J.; Rutegård, J.; Matthiessen, P.; Rutegård, M. Postoperative nonsteroidal anti-inflammatory drugs in relation to recurrence, survival and anastomotic leakage after surgery for colorectal cancer. Colorectal Dis. 2022, 24, 933–942. [Google Scholar] [CrossRef]
  55. Knights, K.M.; Mangoni, A.A.; Miners, J.O. Defining the COX inhibitor selectivity of NSAIDs: Implications for understanding toxicity. Expert Rev. Clin. Pharmacol. 2010, 3, 769–776. [Google Scholar] [CrossRef]
  56. Kolawole, O.R.; Kashfi, K. NSAIDs and cancer resolution: New paradigms beyond cyclooxygenase. Int. J. Mol. Sci. 2022, 23, 1432. [Google Scholar] [CrossRef]
  57. Gurpinar, E.; Grizzle, W.E.; Piazza, G.A. COX-independent mechanisms of cancer chemoprevention by anti-inflammatory drugs. Front. Oncol. 2013, 3, 181. [Google Scholar] [CrossRef]
  58. Feng, Y.; Tao, L.; Wang, G.; Li, Z.; Yang, M.; He, W.; Zhong, X.; Zhang, Y.; Yang, J.; Cheung, S.; et al. Aspirin inhibits prostaglandins to prevents colon tumor formation via down-regulating Wnt production. Eur. J. Pharmacol. 2021, 906, 174173. [Google Scholar] [CrossRef]
  59. Farhangian, M.; Azarafrouz, F.; Fallahi, H.; Akrami, H. Investigating the effect of ibuprofen on DLL1 and NOTCH1 expression in gastric cancer stem cells derived from MKN-45 cell line. Middle East J. Cancer 2019, 10, 292–298. [Google Scholar]
  60. Karelia, D.N.; Kim, S.; Pandey, M.K.; Plano, D.; Amin, S.; Lu, J.; Sharma, A.K. Novel seleno-aspirinyl compound AS-10 induces apoptosis, G1 arrest of pancreatic ductal adenocarcinoma cells, inhibits their NF-κB signaling, and synergizes with gemcitabine cytotoxicity. Int. J. Mol. Sci. 2021, 22, 4966. [Google Scholar] [CrossRef]
  61. Plano, D.; Karelia, D.N.; Pandey, M.K.; Spallholz, J.E.; Amin, S.; Sharma, A.K. Design, synthesis, and biological evaluation of novel selenium (Se-NSAID) molecules as anticancer agents. J. Med. Chem. 2016, 59, 1946–1959. [Google Scholar] [CrossRef]
  62. He, X.; Zhong, M.; Li, S.; Li, X.; Li, Y.; Li, Z.; Gao, Y.; Ding, F.; Wen, D.; Lei, Y.; et al. Synthesis and biological evaluation of organoselenium (NSAIDs-SeCN and SeCF3) derivatives as potential anticancer agents. Eur. J. Med. Chem. 2020, 208, 112864. [Google Scholar] [CrossRef]
  63. He, X.; Nie, Y.; Zhong, M.; Li, S.; Li, X.; Guo, Y.; Liu, Z.; Gao, Y.; Ding, F.; Wen, D.; et al. New organoselenides (NSAIDs-Se derivatives) as potential anticancer agents: Synthesis, biological evaluation and in silico calculations. Eur. J. Med. Chem. 2021, 218, 113384. [Google Scholar] [CrossRef]
  64. Desai, D.; Sinha, I.; Null, K.; Wolter, W.; Suckow, M.A.; King, T.; Amin, S.; Sinha, R. Synthesis and antitumor properties of selenocoxib-1 against rat prostate adenocarcinoma cells. Int. J. Cancer 2010, 127, 230–238. [Google Scholar] [CrossRef]
  65. Gowda, R.; Madhunapantula, S.V.; Desai, D.; Amin, S.; Robertson, G.P. Simultaneous targeting of COX-2 and AKT using selenocoxib-1-GSH to inhibit melanoma. Mol. Cancer Ther. 2013, 12, 3–15. [Google Scholar] [CrossRef] [PubMed]
  66. Ramos-Inza, S.; Encío, I.; Raza, A.; Sharma, A.K.; Sanmartín, C.; Plano, D. Design, synthesis and anticancer evaluation of novel Se-NSAID hybrid molecules: Identification of a Se-indomethacin analog as a potential therapeutic for breast cancer. Eur. J. Med. Chem. 2022, 244, 114839. [Google Scholar] [CrossRef]
  67. Pauli, G.F.; Chen, S.N.; Simmler, C.; Lankin, D.C.; Godecke, T.; Jaki, B.U.; Friesen, J.B.; McAlpine, J.B.; Napoitano, J.G. Importance of purity evaluation and the potential of quantitative 1H NMR as a purity assay. J. Med. Chem. 2014, 57, 9220–9231. [Google Scholar] [CrossRef]
  68. Sanmartín, C.; Plano, D.; Domínguez, E.; Font, M.; Calvo, A.; Prior, C.; Encío, I.; Palop, J.A. Synthesis and pharmacological screening of several aroyl and heteroaroyl selenylacetic acid derivatives as cytotoxic and antiproliferative agents. Molecules 2009, 14, 3313–3338. [Google Scholar] [CrossRef]
  69. Tanini, D.; Capperucci, A.; Locuoco, M.; Ferraroni, M.; Costantino, G.; Angeli, A.; Supuran, C.T. Benzoselenoates: A novel class of carbonic anhydrase inhibitors. Bioorg. Chem. 2022, 122, 105751. [Google Scholar] [CrossRef]
  70. Weekley, C.M.; Harris, H.H. Which form is that? The importance of selenium speciation and metabolism in the prevention and treatment of disease. Chem. Soc. Rev. 2013, 42, 8870–8894. [Google Scholar] [CrossRef] [PubMed]
  71. Lü, J.; Zhang, J.; Jiang, C.; Deng, Y.; Özten, N.; Bosland, M.C. Cancer chemoprevention research with selenium in the post-SELECT era: Promises and challenges. Nutr. Cancer 2016, 68, 1–17. [Google Scholar] [CrossRef]
  72. Capperucci, A.; Degl’Innocenti, A.; Tiberi, C. Organoselenosilane-mediated selective mild access to selenolesters, selenoanhydrides and diacyl diselenides. Synlett 2011, 2011, 2248–2252. [Google Scholar] [CrossRef]
  73. Jia, X.; Zhang, Y.; Zhou, X. Samarium diiodide-induced reduction of elemental selenium leading to a selenolate anion species. A selective synthesis of diacyl diselenides. Synth. Commun. 1993, 23, 1403–1408. [Google Scholar] [CrossRef]
  74. Ishihara, H.; Koketsu, M.; Fukuta, Y.; Nada, F. Reaction of lithium aluminum hydride with elemental selenium: Its application as a selenating reagent into organic molecules. J. Am. Chem. Soc. 2001, 123, 8408–8409. [Google Scholar] [CrossRef] [PubMed]
  75. Koketsu, M.; Nada, F.; Hiramatsu, S.; Ishihara, H. Reactions of acyl chlorides with LiAlHSeH. Preparation of diacyl selenides, diacyl diselenides, selenocarboxylates and cyclic selenoanhydrides. J. Chem. Soc. Perkin Trans. 2002, 1, 737–740. [Google Scholar] [CrossRef]
  76. Kumar, S.; Tripathi, S.K.; Singh, H.B.; Wolmershäuser, G. Synthesis, reactivity, electrochemical and crystallographic studies of diferrocenoyl diselenide and ferrocenoyl selenides. J. Organomet. Chem. 2004, 689, 3046–3055. [Google Scholar] [CrossRef]
  77. Niyomura, O.; Tani, K.; Kato, S. A facile synthesis of potassium selenocarboxylates and their oxidation with XeF2 to diacyl diselenides: An X-ray structural analysis of di(4-methoxybenzoyl) diselenide. Heteroat. Chem. 1999, 10, 373–379. [Google Scholar] [CrossRef]
  78. Astrain-Redin, N.; Talavera, I.; Moreno, E.; Ramírez, M.J.; Martínez-Sáez, N.; Encío, I.; Sharma, A.K.; Sanmartín, C.; Plano, D. Seleno-analogs of scaffolds resembling natural products a novel warhead toward dual compounds. Antioxidants 2023, 12, 139. [Google Scholar] [CrossRef]
  79. Kesavardhana, S.; Malireddi, R.K.S.; Kanneganti, T.D. Caspases in cell death, inflammation, and pyroptosis. Annu. Rev. Immunol. 2020, 38, 567–595. [Google Scholar] [CrossRef]
  80. Arisan, E.D.; Ergül, Z.; Bozdağ, G.; Rencüzoğulları, Ö.; Çoker-Gürkan, A.; Obakan-Yerlikaya, P.; Coşkun, D.; Palavan-Ünsal, N. Diclofenac induced apoptosis via altering PI3K/Akt/MAPK signaling axis in HCT 116 more efficiently compared to SW480 colon cancer cells. Mol. Biol. Rep. 2018, 45, 2175–2184. [Google Scholar] [CrossRef]
  81. Din, F.V.; Valanciute, A.; Houde, V.P.; Zibrova, D.; Green, K.A.; Sakamoto, K.; Alessi, D.R.; Dunlop, M.G. Aspirin inhibits mTOR signaling, activates AMP-activated protein kinase, and induces autophagy in colorectal cancer cells. Gastroenterology 2012, 142, 1504–1515. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Chemical structures of some NSAID derivatives bearing Se [60,61,62,63,64,65].
Figure 1. Chemical structures of some NSAID derivatives bearing Se [60,61,62,63,64,65].
Antioxidants 12 01666 g001
Figure 2. (A) Schematic diagrams of the chemical design of Se-modified NSAIDs, including previous possible cleavages according to where the bond is broken (1 or 2). (B) General structure of the new Se-NSAID derivatives 15 synthesized in this work and the parent NSAIDs used.
Figure 2. (A) Schematic diagrams of the chemical design of Se-modified NSAIDs, including previous possible cleavages according to where the bond is broken (1 or 2). (B) General structure of the new Se-NSAID derivatives 15 synthesized in this work and the parent NSAIDs used.
Antioxidants 12 01666 g002
Scheme 1. Synthesis of the Se-NSAID derivatives. Reagents and conditions: (i) ClCOCOCl, CH2Cl2, 12 h, room temperature; (ii) ClCOCOCl, CH2Cl2/N,N-DMF, 12 h, room temperature; (iii) H2O, room temperature; (iv) THF/H2O, 30–60 min, room temperature.
Scheme 1. Synthesis of the Se-NSAID derivatives. Reagents and conditions: (i) ClCOCOCl, CH2Cl2, 12 h, room temperature; (ii) ClCOCOCl, CH2Cl2/N,N-DMF, 12 h, room temperature; (iii) H2O, room temperature; (iv) THF/H2O, 30–60 min, room temperature.
Antioxidants 12 01666 sch001
Figure 3. (A) DPPH radical scavenging activity of the Se-NSAID derivatives at 0.06 mg/mL and measured at eight time points. (B) DPPH radical scavenging activity of the Se-NSAID derivatives at 0.03 mg/mL and measured at eight time points. Data are presented as the mean ± SD of three independent experiments.
Figure 3. (A) DPPH radical scavenging activity of the Se-NSAID derivatives at 0.06 mg/mL and measured at eight time points. (B) DPPH radical scavenging activity of the Se-NSAID derivatives at 0.03 mg/mL and measured at eight time points. Data are presented as the mean ± SD of three independent experiments.
Antioxidants 12 01666 g003
Figure 4. Induction of apoptosis in HCT-116 cells. (A) Cells were treated with five concentrations of 5 for 24 h and examined on a Muse™ automated cell analyzer with the Annexin V & Dead Cell apoptosis assay. (B) Analogous independent experiment performed with the Caspase 3/7 apoptosis assay. (C,D) Quantification of the cell population with the Annexin V & Dead Cell and Caspase 3/7 assays, respectively. Data are presented as the mean ± SEM of at least three independent experiments. ** p < 0.01, * p < 0.05 comparing control and different concentrations of compound 5.
Figure 4. Induction of apoptosis in HCT-116 cells. (A) Cells were treated with five concentrations of 5 for 24 h and examined on a Muse™ automated cell analyzer with the Annexin V & Dead Cell apoptosis assay. (B) Analogous independent experiment performed with the Caspase 3/7 apoptosis assay. (C,D) Quantification of the cell population with the Annexin V & Dead Cell and Caspase 3/7 assays, respectively. Data are presented as the mean ± SEM of at least three independent experiments. ** p < 0.01, * p < 0.05 comparing control and different concentrations of compound 5.
Antioxidants 12 01666 g004
Figure 5. In vivo toxicity of compound 5 in mice. Female nude mice (n = 5 in each group) were treated either with the vehicle or compound 5 for 27 days. (A) Biochemical analysis of the blood serum withdrawn at the endpoint. The treatment with a dose of 7.5 mg/Kg of 5 given thrice weekly showed no significant changes in the parameters measured compared to the control group. (B) Body weight changes of mice during the treatment. Data are presented as the mean ± SEM.
Figure 5. In vivo toxicity of compound 5 in mice. Female nude mice (n = 5 in each group) were treated either with the vehicle or compound 5 for 27 days. (A) Biochemical analysis of the blood serum withdrawn at the endpoint. The treatment with a dose of 7.5 mg/Kg of 5 given thrice weekly showed no significant changes in the parameters measured compared to the control group. (B) Body weight changes of mice during the treatment. Data are presented as the mean ± SEM.
Antioxidants 12 01666 g005
Figure 6. In vivo efficacy of compound 5 in a colon cancer xenograft mouse model. Female nude mice were inoculated with 2.5 × 106 HCT-116 cells on each flank. After 15 days, the HCT-116 tumor-bearing mice (n = 5 in each group) were treated either with the vehicle or compound 5 for additional 23 days. (A) Growth curves of tumor volume. The treatment with a dose of 10 mg/Kg of 5 given thrice weekly achieved significant tumor growth inhibition compared to the control group. (B) Tumor weight measured at the endpoint. (C) Body weight changes of mice during the treatment. Data are presented as the mean ± SEM. **** p < 0.0001 comparing control and treatment with compound 5.
Figure 6. In vivo efficacy of compound 5 in a colon cancer xenograft mouse model. Female nude mice were inoculated with 2.5 × 106 HCT-116 cells on each flank. After 15 days, the HCT-116 tumor-bearing mice (n = 5 in each group) were treated either with the vehicle or compound 5 for additional 23 days. (A) Growth curves of tumor volume. The treatment with a dose of 10 mg/Kg of 5 given thrice weekly achieved significant tumor growth inhibition compared to the control group. (B) Tumor weight measured at the endpoint. (C) Body weight changes of mice during the treatment. Data are presented as the mean ± SEM. **** p < 0.0001 comparing control and treatment with compound 5.
Antioxidants 12 01666 g006
Table 1. IC50 values (in µM) 1 for the Se-NSAID analogs and parent NSAIDs in colon (HT-29, HCT-116), prostate (DU-145, PC-3), breast (MDA-MB-231, T-47D), and lung (H1299, A549) cancer cell lines.
Table 1. IC50 values (in µM) 1 for the Se-NSAID analogs and parent NSAIDs in colon (HT-29, HCT-116), prostate (DU-145, PC-3), breast (MDA-MB-231, T-47D), and lung (H1299, A549) cancer cell lines.
Comp.Colon Cell LinesProstate Cell LinesBreast Cell LinesLung Cell Lines
HT-29HCT-116DU-145PC-3MDA-MB-231T-47DH1299A549
16.6 ± 2.239.2 ± 7.35.0 ± 0.13.3 ± 0.618.3 ± 6.95.9 ± 1.48.6 ± 0.87.0 ± 1.3
2>50>50>50>50>50>50>50>50
3>509.1 ± 5.96.1 ± 0.65.2 ± 0.25.7 ± 0.85.6 ± 0.18.1 ± 2.5>50
439.7 ± 3.35.9 ± 1.127.9 ± 3.25.8 ± 0.64.3 ± 1.226.0 ± 2.031.6 ± 1.711.7 ± 0.9
58.5 ± 2.79.0 ± 3.07.0 ± 1.37.4 ± 0.58.3 ± 2.53.8 ± 0.29.3 ± 0.64.5 ± 0.6
ASA 2>50>50>50>50>50>50>50>50
Ind>50>50>50>50>50>50>50>50
Nap>50>50>50>50>50>50>50>50
Ket>50>50>50>50>50>50>50>50
Ibup>50>50>50>50>50>50>50>50
1 IC50 values are presented as the mean ± SD of at least three independent experiments determined by the MTT assay. 2 Some IC50 values for parent NSAIDs were obtained from Plano et al. [61], He et al. [62,63], and Ramos-Inza et al. [66].
Table 2. IC50 values (in µM) 1 for the Se-NSAID analogs in breast (MDA-MB-231, T-47D) and lung (H1299, A549) cancer cell lines along with the IC50 values (in µM) obtained in nonmalignant cells (184B5 and BEAS-2B) and the respective SI calculated.
Table 2. IC50 values (in µM) 1 for the Se-NSAID analogs in breast (MDA-MB-231, T-47D) and lung (H1299, A549) cancer cell lines along with the IC50 values (in µM) obtained in nonmalignant cells (184B5 and BEAS-2B) and the respective SI calculated.
Comp.Breast Cell LinesLung Cell Lines
184B5MDA-MB-231SI 2T-47DSI 3BEAS-2BH1299SI 4A549SI 5
114.4 ± 1.518.3 ± 6.90.85.9 ± 1.42.418.5 ± 3.08.6 ± 0.82.27.0 ± 1.32.6
2>100>100->100->100>100->100-
311.9 ± 1.05.7 ± 0.82.15.6 ± 0.12.111.2 ± 1.08.1 ± 2.51.4>100>0.1
45.9 ± 0.84.3 ± 1.21.426.0 ± 2.00.25.2 ± 0.431.6 ± 1.70.211.7 ± 0.90.4
518.2 ± 0.88.3 ± 2.52.23.8 ± 0.24.815.9 ± 2.49.3 ± 0.61.74.5 ± 0.63.5
1 IC50 values are presented as the mean ± SD of at least three independent experiments determined by the MTT assay. 2 SI calculated in breast cells as IC50 (184B5)/IC50 (MDA-MB-231). 3 SI calculated in breast cells as IC50 (184B5)/IC50 (T-47D). 4 SI calculated in lung cells as IC50 (BEAS-2B)/IC50 (H1299). 5 SI calculated in lung cells as IC50 (BEAS-2B)/IC50 (A549).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ramos-Inza, S.; Aliaga, C.; Encío, I.; Raza, A.; Sharma, A.K.; Aydillo, C.; Martínez-Sáez, N.; Sanmartín, C.; Plano, D. First Generation of Antioxidant Precursors for Bioisosteric Se-NSAIDs: Design, Synthesis, and In Vitro and In Vivo Anticancer Evaluation. Antioxidants 2023, 12, 1666. https://doi.org/10.3390/antiox12091666

AMA Style

Ramos-Inza S, Aliaga C, Encío I, Raza A, Sharma AK, Aydillo C, Martínez-Sáez N, Sanmartín C, Plano D. First Generation of Antioxidant Precursors for Bioisosteric Se-NSAIDs: Design, Synthesis, and In Vitro and In Vivo Anticancer Evaluation. Antioxidants. 2023; 12(9):1666. https://doi.org/10.3390/antiox12091666

Chicago/Turabian Style

Ramos-Inza, Sandra, Cesar Aliaga, Ignacio Encío, Asif Raza, Arun K. Sharma, Carlos Aydillo, Nuria Martínez-Sáez, Carmen Sanmartín, and Daniel Plano. 2023. "First Generation of Antioxidant Precursors for Bioisosteric Se-NSAIDs: Design, Synthesis, and In Vitro and In Vivo Anticancer Evaluation" Antioxidants 12, no. 9: 1666. https://doi.org/10.3390/antiox12091666

APA Style

Ramos-Inza, S., Aliaga, C., Encío, I., Raza, A., Sharma, A. K., Aydillo, C., Martínez-Sáez, N., Sanmartín, C., & Plano, D. (2023). First Generation of Antioxidant Precursors for Bioisosteric Se-NSAIDs: Design, Synthesis, and In Vitro and In Vivo Anticancer Evaluation. Antioxidants, 12(9), 1666. https://doi.org/10.3390/antiox12091666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop