Next Article in Journal
A General Synthesis of C8-Arylpurine Phosphoramidites
Next Article in Special Issue
Cesium Carbonate-Catalyzed α-Phenylchalcogenation of Carbonyl Compounds with Diphenyl Dichalcogenide
Previous Article in Journal
Squalene Emulsions for Parenteral Vaccine and Drug Delivery
Previous Article in Special Issue
Theoretical Studies of [2,3]-Sigmatropic Rearrangements of Allylic Selenoxides and Selenimides
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 14
Issue 9
10.3390/molecules14093313
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Synthesis and Pharmacological Screening of Several Aroyl and Heteroaroyl Selenylacetic Acid Derivatives as Cytotoxic and Antiproliferative Agents

by
Carmen Sanmartín
1,*,
Daniel Plano
1,
Enrique Domínguez
1,
María Font
2,
Alfonso Calvo
3,
Celia Prior
3,
Ignacio Encío
4 and
Juan Antonio Palop
1,*
1
Synthesis Section, Department of Organic and Pharmaceutical Chemistry, University of Navarra, Irunlarrea, 1, E-31008 Pamplona, Spain
2
Molecular Modeling Section, Department of Organic and Pharmaceutical Chemistry, University of Navarra, Irunlarrea, 1, E-31008 Pamplona, Spain
3
Oncology Division, Center for Applied Medical Research, CIMA, University of Navarra, Pío XII, 53, E-31008 Pamplona, Spain
4
Department of Health Science, Public University of Navarra, Avda. Barañain, s/n, E-31008 Pamplona, Spain
*
Authors to whom correspondence should be addressed.
Molecules 2009, 14(9), 3313-3338; https://doi.org/10.3390/molecules14093313
Submission received: 1 July 2009 / Revised: 4 August 2009 / Accepted: 27 August 2009 / Published: 1 September 2009
(This article belongs to the Special Issue Selenium and Tellurium Chemistry)
Browse Figures
Versions Notes

Abstract

:
The synthesis and cytotoxic activity of a series of twenty six aroyl and heteroaroyl selenylacetic acid derivatives of general formula Ar-CO-Se-CH2-COOH or Heterar-CO-Se-CH2-COOH are reported. The synthesis was carried out by reaction of acyl chlorides with sodium hydrogen selenide, prepared in situ, and this led to the formation of sodium aroylselenides that subsequently reacted with α-bromoacetic acid to produce the corresponding selenylacetic acid derivatives. All of the compounds were tested against a prostate cancer cell line (PC-3) and some of the more active compounds were assessed against a panel of four human cancer cell lines (CCRF-CEM, HTB-54, HT-29, MCF-7) and one mammary gland-derived non-malignant cell line (184B5). Some of the compounds exhibited remarkable cytotoxic and antiproliferative activities against MCF-7 and PC-3 that were higher than those of the reference compounds doxorubicin and etoposide, respectively. For example, in MCF-7 when Ar = phenyl, 3,5-dimethoxyphenyl or benzyl the TGI values were 3.69, 4.18 and 6.19 μM. On the other hand, in PC-3 these compounds showed values of 6.8, 4.0 and 2.9 μM. Furthermore, benzoylselenylacetic acid did not provoke apoptosis nor did it perturb the cell cycle in MCF-7.

Introduction

Cancer is currently one of the leading causes of death worldwide and could become the most common cause in the future [1]. Our understanding of the biology of cancer has undoubtedly improved in the last decade. One characteristic of cancer cells is their highly proliferative nature. Consequently, inhibition of proliferative pathways is considered an effective strategy to fight cancer and a great deal of attention has recently been paid to the discovery and development of new, more selective anticancer agents. The results of numerous studies indicate that selenium is an essential trace element [2,3] with structural and enzymatic roles for viral suppression and against AIDS. This element is also implicated in delaying the aging process. Selenium deficiency has been linked to a number of disorders such as heart disease, diabetes and cancer. Selenium compounds have attracted renewed interest as chemopreventive agents for human cancer on the basis of the pioneering intervention study carried out by Clark and Combs [4]. Later epidemiological and clinical intervention studies support the protective role of selenium against the development of prostate cancer [5,6,7]. More recently, mounting evidence indicates that Se compounds inhibit neoplasm of lung, colon, mammary gland, sarcoma and glioblastoma [8,9]. In addition, selenium supplementation [10,11,12] could provide significant therapeutic potential for patients with gastrointestinal, liver and head and neck cancers.
Furthermore, literature reports have consistently shown that the dose and chemical forms of selenium [13,14] are determinant factors in anticancer activity, with organic Se compounds being more bioavailable than the inorganic forms.
Several chemical forms are under investigation for their potential as candidate drugs and these include methylseleno derivatives such as methylselenocysteine [15], p-xylylbis(methylselenide) [16] and methylseleninic acid [17], compounds with selenium in the heterocyclic ring as in selenazolidines [18], ebselen [19], D-501036 [20], BBSKE [21] or the presence of a selenide function in a symmetrical structure such as diphenyl-, dimethyl-, dibenzylselenide [22,23,24], bis(3-methyl-4-pyridyl)selenide [25] and di(3-indolyl)selenide [26].
On the other hand, the metabolism of selenium compounds is critical for activity and at least two selenium metabolite pools, hydrogen selenide and methylselenol [27,28], that induce distinct types ofantiangiogenesis, apoptosis, and cell cycle responses, have been implicated as active metabolites for the anticancer effects.
Based on the interesting and widely observed antitumoral activity of selenium compounds and our own experience in this field [29,30], we report here the synthesis, biological evaluation and molecular modelling study of a group of selenylacetic acid derivatives 1–26, (Figure 1) as potential agents in cancer therapy.
Molecules 14 03313 g001 550
Figure 1. General formula of compounds.
Figure 1. General formula of compounds.
Molecules 14 03313 g001
Among the strategies adopted for the design of these structures, particular consideration was given to the following approaches:
  • (a) The use of the selenide function due to the facile scission of Se from the organic moiety in these types of compounds. In addition, this proposed preliminary hypothesis concerning the action mechanism of these derivatives, related to the possible scission of Se from the organic moiety, allow us to determine the bond order. According to the Molecular Orbital Theory, the bond order (b.o.) is equivalent to the number of electrons in the antibonding molecular orbital minus the number of electrons in the bonding molecular orbital divided by two. This parameter can be taken as a quantitative descriptor for the bond strength and can be related to the aforementioned Se scission.
  • (b) The appropriate choice of the substituent on the phenyl ring was made according to a synthetic accessibility and with the aim of assessing the following aspects:
    (b.1)
    Influence of the substituents on the Se charge and bond stability, as well as on the acidic character of the derivatives, expressed as pKa.
    (b.2)
    Modulation of the electronic distribution over the aromatic moiety with a variety of electron-donating and electron-withdrawing substituents, placed at different position with respect to the lateral chain.
    (b.3)
    Molecular volume, conformational behaviour and hydrophobic character (expressed as AlogP), considering the presence of the keto moiety, the aromatic/heteroaromatic ring, the methylene bridge and the carboxylic moiety as the structural basic pattern.

Results and Discussion

Chemistry

The synthesis of the title compounds 1–26 was carried out as depicted in Scheme 1. In the first step sodium hydrogen selenide was produced by reaction of powdered grey selenium with metallic hydrides in an appropriate medium. In the second step the aroyl or heteroaroyl chloride was reacted with sodium hydrogen selenide to produce the intermediate sodium aroyl selenides. The last step involved the formation of selenylacetic acids 1–26 by treatment of this intermediate with α-bromoacetic acid. This synthesis was based on a published procedure with modifications [31] following our own protocol.
Molecules 14 03313 g009 550
Scheme 1. Synthesis of compounds 1–26.
Scheme 1. Synthesis of compounds 1–26.
Molecules 14 03313 g009
The compounds were obtained in yields from 3 to 78%. The analytical data are shown in Table 1. The purity levels of compounds 1–26 were assessed by TLC and elemental analyses and their structures were identified from spectroscopic data (IR, 1H NMR). The IR spectra of these compounds revealed intense absorption bands in the ranges 1,720–1,700 cm–1 and 1,690–1,665 cm–1, assigned to carboxyl and carbonyl groups, respectively. The frequencies of these absorptions varied according to the substitution pattern on the arylidene ring. As expected, a decrease in the carbonyl absorption frequencies was observed for compounds having electron-donating groups such as methylenedioxy, trimethoxy or tert-butyl attached to the benzylidene ring in comparison to electron-withdrawing groups (chloro, bromo).

Biological Evaluation

Cytotoxic activity in PC-3

Initially, the new compounds 1–26 were evaluated for their in vitro cytotoxic activity against a human prostate cancer cell line (PC-3, ATCC, Manassas, VA, USA) using the MTT assay [32]. We selected this cell line because there are numerous clinical trials that show the activity of selenium compounds in the reduction of prostate cancer [5,9,15,16]. The results are tabulated as IC50 values. All experiments were independently performed at least three times and the values were calculated after 72 hours exposure (drug concentrations of 2, 5, 7 and 10 μM). The results are shown in Table 2.
The maximum cytotoxic activity was exhibited by compounds 1, 2, 11, 14 and 21. By comparing the activities of compounds with those of methylseleninic acid and etoposide, both effective drugs against prostate cancer that have been employed as controls for cytotoxic assays, the results revealed that three of the new compounds −1, 11, 14 − were more active than methylseleninic acid and etoposide and that five were better than etoposide. It is noteworthy that the introduction of a new methoxy groups in the 4-position in 14 to give 15 results in suppression of the antitumour activity on the PC-3 cell line. These findings, along with the observed antitumour activity of compound 14, suggest that the activity of the compounds depends on the location of the methoxy substituents on the phenyl ring. With the aim of evaluating the preliminary structure-activity relationships to gain an understanding of the activity shown by such compounds as antitumorals, we carried out a molecular modelling study on the aforementioned derivatives.
Table
Table 1. Physical constants for compounds 1–26.
Table 1. Physical constants for compounds 1–26.
Molecules 14 03313 i001
Ref.RYield (%)M. p. (ºC)Recryst. SolventCHNAnal (%) Calcd/Found
1phenyl7583aTolueneC9H8O3SeC, 44.44/44.47; H, 3.29/3.24; N, 0.00/0.00.
24-cyanophenyl78147aTolueneC10H7NO3SeC, 44.77/44.85; H, 2.61/2.67; N, 5.22/4.98.
34-(trifluoromethyl)phenyl62108aTolueneC10H7F3O3SeC, 38.58/38.73; H, 2.25/2.15; N, 0.00/0.00.
44-chlorophenyl19135-136b,cTolueneC9H7ClO3SeC, 38.95/38.71; H, 2.54/2.43; N, 0.00/0.00.
54-methylphenyl5392aTolueneC10H1003SeC, 44.69/44.95; H, 3.89/3.89; N, 0.00/0.00.
64-( tert-butyl)phenyl1999-103b,cChloroformC13H16O3SeC, 52.18/52.04; H, 5.39/5.31; N, 0.00/0.05.
74-methoxyphenyl16104-107bChloroform/Carbon tetrachlorideC10H10O4SeC, 43.97/43.49; H, 3.69/3.52; N, 0.00/0.00.
82-chlorophenyl44123-125bCarbon tetrachlorideC9H7ClO3SeC, 38.95/38.97; H, 2.54/2.52; N, 0.00/0.06.
92-bromophenyl32124-128bCarbon tetrachlorideC9H7BrO3SeC, 33.57/33.14; H, 2.19/2.02; N, 0.00/0.20.
102-iodophenyl3105-108bCarbon tetrachlorideC9H7IO3SeC, 29.29/29.14; H, 1.91/1.92; N, 0.00/0.05.
11benzyl4874aTolueneC9H10O3SeC, 46.69/46.87; H, 3.89/3.89; N, 0.00/0.00.
122-phenylethyl1165-69bHexaneC11H12O3S ¼ H2OC, 47.93/48.06; H, 4.57/4.47; N, 0.00/0.18.
133,5-dichlorophenyl17108-109bCarbon tetrachlorideC9H6Cl2O3Se·½ H2OC, 33.67/33.56; H, 2.20/1.84; N, 0.00/0.06.
143,5-dimethoxyphenyl69117aTolueneC11H12O5SeC, 43.56/43.70; H, 3.96/3.88; N, 0.00/0.00.
153,4,5-trimethoxyphenyl14107-110bEthanolC12H14O6SeC, 43.26/43.32; H, 4.24/3.97; N, 0.00/0.11.
163,4-methylendioxyphenyl47106-113bCarbon tetrachlorideC10H8O5SeC, 41.83/41.56; H, 2.81/2.67; N, 0.00/0.02.
17naphthyl37130aTolueneC13H10O3SeC, 53.24/52.99; H, 3.41/3.25; N, 0.00/0.00.
18diphenylmethyl16127-130bCarbon tetrachlorideC16H14O3SeC, 56.90/56.88; H, 4.33/4.11; N, 0.00/0.14.
194-pyridyl11119-121b,cEther / hexaneC8H7NO3SeC, 39.36/39.39; H, 2.89/2.78; N, 5.74/5.67.
203-pyridyl15147-150b,cMethanolC8H7NO3SeC, 39.36/39.57; H, 2.89/2.76; N,5.74/5.63.
213-(2-chloropyridyl) 157aTolueneC8H6ClNO3SeC, 34.47/34.68; H, 2.15/2.17; N, 5.01/5.26.
223-(2-propylthiopyridyl)36109-111bCarbon tetrachlorideC11H13NO3SSeC, 41.51/41.28; H, 4.12/3.90; N, 4.40/4.24.
232-thienyl3382-84bCarbon tetrachlorideC7H6O3SSeC, 33.75/33.45; H, 2.43/2.54; N, 0.00/0.07.
24pyrazinyl10138-140bIsopropanolC7H6N2O3SeC, 34.30/34.18; H, 2.47/2.49; N, 1.43/1.37.
252-quinolyl6131-132b,cTolueneC12H9NO3SeC, 49.00/49.33; H, 3.08/3.13; N, 4.76/4.75.
263-quinolyl9187-189b,cChloroformC12H9NO3SeC, 49.00/48.95; H, 3.08/3.23; N, 4.76/4.66.
a Determined by differential scanning calorimetry; b Determined by thermomicroscopy; c Fusion with degradation.
Table
Table 2. Cytotoxic activities of compounds 1–26 against the PC-3 cell line.
Table 2. Cytotoxic activities of compounds 1–26 against the PC-3 cell line.
Molecules 14 03313 i002
CompoundRPC-3 cell line
IC50 (μM)
1phenyl6.8
24-cyanophenyl10.0
34-(trifluoromethyl)phenylNEa
44-chlorophenylNE
54-methylphenylNE
64-tert-butylphenylNE
74-methoxyphenylNE
82-chlorophenylNE
92-bromophenylNE
102-iodophenylNE
11benzyl2.9
122-phenylethylNE
133,5-dichlorophenylNE
143,5-dimethoxyphenyl4.0
153,4,5-trimethoxyphenylNE
163,4-methylenedioxyphenylNE
17naphthylNE
18diphenylmethylNE
194-pyridylNE
203-pyridylNE
213-(2-chloropyridyl)10.0
223-(2-propylthiopyridyl)NE
232-thienylNE
24pyrazinylNE
252-quinolylNE
263-quinolylNE
MSAb 8.38 [33]
Etoposide 13.6 ± 2.2 [34]
a NE = no effect is observed; b methylseleninic acid.

Molecular modelling

From the molecular modelling point of view, three approaches were applied in an effort to gain an insight into the structural requirements for activity in these compounds (Figure 2).
Molecules 14 03313 g002 550
Figure 2. Overview of structural variations carried out for the studied derivatives: (a) aryl or heteroaryl monocyclic derivatives, (b) aryl or heteroaryl bicyclic derivatives.
Figure 2. Overview of structural variations carried out for the studied derivatives: (a) aryl or heteroaryl monocyclic derivatives, (b) aryl or heteroaryl bicyclic derivatives.
Molecules 14 03313 g002
Firstly, it was planned to obtain some descriptive parameter at a molecular level, such as the log of the partition coefficient expressed as ALogP and the molecular volume, particularly the volume of the aromatic moiety, that would allow discrimination between active and inactive compounds.
A second approach was carried out. This was based on the pharmacophoric hypothesis that an active compound can be characterised by a set of chemical functions with a certain spatial arrangement and molecular shape that could be related to the observed biological activity but would not be associated with the inactive compounds. Thus, the presence of certain chemical functions is not only necessary, but these functions must be exposed in a manner that can be recognised by the associated binding site — a property that requires specific conformational behaviour [35]. In agreement with this approach, a compound can be considered as a collection of multiple instances, one for each conformation. An active compound contains at least one active instance, while an inactive compound does not contain any active instance. On the assumption that all of the active compounds bind to the same site, we would expect their binding conformations to have similar shapes. For this second approach it is necessary to carry out a study of the conformational behaviour of the compounds and subsequently select and analyse the supposed bioactive conformations of the compounds.
The third approach considered here concerned the structural characteristics of designed molecules, especially the presence of the selenium atom, the proposed structural variations related to the length of the chain, which acts like a bridge between the keto group and the aromatic or heteroaromatic moiety, and the electronic characteristics of the groups located in these rings. These structural variations led us to carry out a study from the mechano-quantic point of view in order to obtain data on parameters such as the value of the partial charge on the selenium, the energy and distribution of the HOMO and LUMO orbitals as well as the electronic density (ED) and molecular electrostatic potential (MEP) distribution. The aim was to try to relate these values with the biological activity of our derivatives.
The calculations were performed on a Dell Precision 380 workstation, provided with the software package Discovery Studio v1.7 [36], and on a SGI Virtu VS100 workstation, provided with MOPAC2009 [37] and Mercury [38] software packages.
Once the models for the compounds had been constructed, the initial geometries were fully minimized to an energy gradient below 10–3 kcal mol–1Å–1. The minimum energy conformers were superimposed, with the ring moiety taken as adjusting atoms. The effectiveness of the superimposed models was evaluated in terms of the Root Mean Square (RMS) values obtained. The energy differences between the different conformations analysed for each trajectory were in the range 2–5 kcal.
Analysis of the resulting data (Table 3, Table 4) enabled the following general observations to be made in an attempt to correlate the data with the cytotoxic activity levels found for the compounds on the PC-3 cell line.
1. The permissible maximum volume for the aromatic fragment is approximately 111–112 Å3. Thus, the bicyclic compounds (naphthyl, diphenylmethyl and quinolyl derivatives) and the monocyclic derivatives with voluminous substituents are inactive.
2. In relation to the conformational behaviour, it can be observed that, in general, in the conformational trajectory of active compounds the number of extended conformations is greater than that in the inactive ones, with a better general superposition of the conformations, and the total volume of the conformational trajectory is also greater (rms data not shown for the sake of brevity). An example is shown in Figure 3.
3. With respect to the ALogP98 values determined on the preliminary geometries obtained after the first minimization (Table 3), although a direct correlation could not be established, the resulting data did shed some light on the biological data, since the active products have a value between 1.518 (compound 21, IC50 = 10.0 μM) and 1.829 (compound 11, the most active compound, IC50 = 2.9 μM). As expected, some structural modifications give rise to a marked reduction in the AlogP values, as is the case for compounds 19 and 20, derived from pyridine, with a value of 0.644, or for compound 24, derived from pyrazine, with a negative value for this parameter. On the other hand, it can be seen that the compounds with an AlogP greater than 1.90 are inactive in this assay.
Table
Table 3. Molecular descriptors obtained for the analyzed compounds (aryl or heteroaryl monocyclic derivatives) a.
Table 3. Molecular descriptors obtained for the analyzed compounds (aryl or heteroaryl monocyclic derivatives) a.
Molecules 14 03313 i003
Ref.XYZnR2R3R4R5R6AlogPVolbμ (D)cIC50(μM)e
1 CCC0HHHHH1.79568.583.4656.8
2 CCC0HHCNHH1.67382.162.86510.0
3 CCC0HHCF3HH2.73792.641.137NE
4 CCC0HHClHH2.45983.192.628NE
5CCC0HHCH3HH2.28182.764.430NE
6CCC0HHC(CH3)3HH3.195125.554.584NE
7 CCC0HHOCH3HH1.77889.515.337NE
8 CCC0ClHHHH2.45982.972.874NE
9 CCC0BrHHHH2.54391.024.032NE
10 CCC0IHHHH2.373100.493.259NE
11 CCC1HHHHH1.82982.733.4182.9
Ref. X Y Z n R2 R3 R4 R5 R6 AlogP Volb μ ( D)c IC50( μ M)d
12 CCC2HHHHH2.28693.294.243NE
13 CCC0HClHClH3.12398.152.441NE
14 CCC0HOCH3HOCH3H1.762111.054.5914.0
15 CCC0HOCH3OCH3OCH3H1.745131.495.090NE
16 CCC0H-O-CH2-O-HH1.56388.784.352NE
19 CCN0HH-HH0.64464.451.081NE
20 NCC0HHHHH0.64464.643.543NE
21 NCC0ClHHHH1.51878.984.09410.0
22 NCC0S(CH2)2CH3HHHH2.597126.664.964NE
23 2-thienyl1.52064.023.315NE
24 NNC0HHHH--0.07860.353.096NE
MSA ------------8.38
Etoposide------------13.6 ± 2.2
a General structure for the analysed compounds showing the bonds (a–e) selected for the conformational analysis. b Volume (average value obtained from the lowest energy conformations) of the cyclic fragment in Å 3. c Dipolar moment (in Debyes) calculated for the representative low-energy. d Cytotoxic activity in PC-3 cell line, NE= no effect.
Table
Table 4. Molecular descriptors obtained for the analyzed compounds (aryl or heteroaryl bicyclic derivatives) a.
Table 4. Molecular descriptors obtained for the analyzed compounds (aryl or heteroaryl bicyclic derivatives) a.
Molecules 14 03313 i004
Ref.RAlogPVolbμ (D)cIC50(μM)d
17naphthyl2.703106.653.759NE
18diphenylmethyl3.324144.343.632NE
252-quinolyl2.409102.414.451NE
263-quinolyl1.981102.222.685NE
MSA----8.38
Etoposide----13.6 ± 2.2
a General structure for the analysed compounds showing the bonds (a–d) selected for the conformational analysis. b Volume (average value obtained from the lowest energy conformations) of the cyclic fragment in Å 3. c Dipolar moment (in Debyes) calculated for the representative low-energy conformation. d Cytotoxic activity in PC-3 cell line, NE= no effect.
Molecules 14 03313 g003 550
Figure 3. Conformational models for some representative compounds: (a) active compounds, left 11; right 1; (b) inactive compounds: left 12; right 6. Carbon in green; hydrogen in grey; oxygen in red; selenium in violet.
Figure 3. Conformational models for some representative compounds: (a) active compounds, left 11; right 1; (b) inactive compounds: left 12; right 6. Carbon in green; hydrogen in grey; oxygen in red; selenium in violet.
Molecules 14 03313 g003
4. With respect to the length of the bridging chain between the aromatic moiety and the keto group, the introduction of a methylene fragment significantly improves the activity (IC50 = 6.8 μM for compound 1 and IC50 = 2.9 μM for compound 11; Table 3 and Figure 4c). However, the introduction of a second methylene leads to the disappearance of the activity. This change could be related to the conformational behaviour. For example, in compound 12 the preferred conformations are the folded ones (see Figure 4c).
Molecules 14 03313 g004 550
Figure 4. Influence of the alkyl chain and ring on the HOMO 0 (a), LUMO 0 (b) distribution, valuated as descriptive parameters, showed on a low-energy representative conformation (c). Carbon in cyan; hydrogen in white; oxygen in red; selenium in violet; sulphur in yellow.
Figure 4. Influence of the alkyl chain and ring on the HOMO 0 (a), LUMO 0 (b) distribution, valuated as descriptive parameters, showed on a low-energy representative conformation (c). Carbon in cyan; hydrogen in white; oxygen in red; selenium in violet; sulphur in yellow.
Molecules 14 03313 g004
5. As regards the dipolar moment, a correlation could not be established and great variability was observed for this parameter, even within the conformational trajectory of each individual compound.
6. With respect to the descriptor parameters obtained by quantum chemistry semi-empirical calculations, the structural variations that lead to a negative value for the partial charge on the selenium lead to the loss of activity. Concerning the distribution of the frontier molecular orbitals (see 4a, 4b, 5a and 5b), it was also found that in the most active compounds the HOMO orbital is preferentially located on the selenium, with a significant contribution from carbons of the aromatic nucleus. In general it can be observed that the presence of substituents in the 4-position brings about the displacement of the HOMO towards the ring, with the contribution of selenium decreasing along with the activity. As an example, the distribution of this orbital for compounds 1 (IC50 = 6.8 μM, Figure 4a) and 2 (IC50 = 10.0 μM, Figure 5a) can be compared with that of compound 7, which contains a methoxy group and is inactive. Nevertheless a correlation with the data of energy for orbital HOMO and LUMO, nor with the gap (LUMO-HOMO energy) cannot be established (Table 5, Table 6).
Molecules 14 03313 g005 550
Figure 5. Influence of the substituents placed in 4-position on the HOMO 0 (a), LUMO 0 (b) distribution, valuated as descriptive parameters, showed on a low-energy representative conformation (c). Carbon in cyan; hydrogen in white; oxygen in red; selenium in violet; nitrogen in blue; chlorine in green.
Figure 5. Influence of the substituents placed in 4-position on the HOMO 0 (a), LUMO 0 (b) distribution, valuated as descriptive parameters, showed on a low-energy representative conformation (c). Carbon in cyan; hydrogen in white; oxygen in red; selenium in violet; nitrogen in blue; chlorine in green.
Molecules 14 03313 g005
Table
Table 5. Mechano-quantic descriptive parameters (semiempirical: PM6) obtained for the analyzed compounds (aryl or heteroaryl monocyclic derivatives).
Table 5. Mechano-quantic descriptive parameters (semiempirical: PM6) obtained for the analyzed compounds (aryl or heteroaryl monocyclic derivatives).
Molecules 14 03313 i005
Ref.XYZnR2R3R4R5R6pKabBond orderHOMOcLUMOdEL-MQ_SeeIC50
abc(M)f
1 CCC0HHHHH2,1191,83150,90680,9618-9,210-1,5527,6580.05776.8
2 CCC0HHCNHH2,0571,84780,90500,9608-9,422-2,0227,4000.077810.0
3 CCC0HHCF3HH2,1481,84920,90480,9608-9,399-1,8947,505-0.0303NE
4 CCC0HHClHH2,1881,79140,96590,9624-9,382-1,6667,716-0.0394NE
5 CCC0HHCH3HH2,1841,82300,90370,9613-9,002-1,4877,5150.0491NE
6 CCC0HHC(CH3)3HH2,1821,82560,90270,9611-9,108-1,4037,7050.0485NE
7 CCC0HHOCH3HH2,1961,81590,89400,9592-9,039-1,2087,8310.0433NE
8 CCC0ClHHHH2,3601,85080,93200,9673-9,258-1,3847,8740.0292NE
9 CCC0BrHHHH2,3461,86110,91720,9524-9,295-1,2448,0510.0134NE
10 CCC0IHHHH2,3881,87210,90690,9626-8,864-1,5827,282-0.0182NE
11 CCC1HHHHH2,3701,84190,97420,9841-9,439-0,9368,5030.06812.9
12 CCC2HHHHH2,3781,83600,97990,9622-9,485-1,1388,3470.0580NE
13 CCC0HClHClH1,8581,84820,91230,9626-9,410-1,9617,449-0.0331NE
14 CCC0HOCH3HOCH3H2,7911,92050,91260,9673-9,275-1,0638,2120.06934.0
15 CCC0HOCH3OCH3OCH3H2,4311,80210,95280,9633-8,770-1,3297,4410.0652NE
16 CCC0H-OCH2O-HH2,0831,84130,90680,9620-9,223-1,6307,5930.0720NE
19 CCN0HH-HH1,9851,78910,99180,9766-9,546-1,7177,8290.0825NE
20 NCC0HHHHH2,2301,77750,98400,9767-9,483-1,6637,8200.0712NE
21 NCC0ClHHHH2,1451,82790,95400,9686-9,368-1,7127,6560.012910.0
22 NCC0S(CH2)2CH3HHHH2,4861,87220,89690,9628-8,871-1,5857,286-0.0093NE
23 2-thienyl 2,1811,78720,94170,9598-9,204-1,6967,508-0.0438NE
24 NNC0HHHH-1,8131,78291,01080,9746-9,406-2,0307,3760.1126NE
Table
Table 6. Mechano-quantic descriptive parameters (semiempirical: PM6) obtained for the analyzed compounds (aryl or heteroaryl bicyclic derivatives) a.
Table 6. Mechano-quantic descriptive parameters (semiempirical: PM6) obtained for the analyzed compounds (aryl or heteroaryl bicyclic derivatives) a.
Molecules 14 03313 i006
Ref.RpKabBond orderHOMOcLUMOdΔEL-MQ_SeeIC50
abc(μM)f
17Naphtyl2,6551,84320,89210,9578-9,076-1,6477,4290.0644NE
18Diphenylmethyl2,5921,88290,94920,9689-9,456-1,1108,3460.0598NE
252-quinolyl2,3061,78841,00120,9606-9,214-1,9187,2960.0918NE
263-quinolyl2,3751,84530,89720,9589-9,222-1,8817,3410.0619NE
Notes for Table 5, Table 6: a General structure for the analyzed compounds showing the bonds (a-b) selected for the bond order determination. b pKa value calculated for the low-energy representative conformation c HOMO 0 orbital energy, in eV, calculated for the low-energy representative conformation d LUMO 0 orbital energy, in eV, calculated for the low-energy representative conformation e Charge on Selenium atom, value obtained for the low-energy representative conformation. fCytotoxic activity in PC-3 cell line, NE= no effect.
7. With respect to the bond order (Table 5, Table 6) a correlation with the activity cannot be established, which contradicts our starting hypothesis related to the stability of the C-Se bond.
8. In relation to the acidic character evaluated by the pKa values (Table 5, Table 6), the values of pKa for the most active compounds oscillate between 2.05 and 2.80, although the proposed structural modifications apparently do not affect in a significant manner this parameter.
In the present state of our research on these selenium derivatives, we do not try to propose a pharmacophore applying only a conformacional analysis. We have carried out a conformacional analysis in order to obtain data about the conformacional behavior of these molecules, with a high degree of conformational freedom, and to sample what are the preferred conformations, between that we hoped to find the bioactive one.

Cytotoxic activity in CCRF-CEM, HTB-54, HT-29, MCF-7 and 184B5

In order to investigate the effect of the active compounds in more detail, we examined the activity on proliferation in other cancer cell lines. The most active compounds in PC-3 were tested, along with the reference compound doxorubicin, for cytotoxic and antiproliferative activities. Tests were carried out in the Department of Health Science, Public University of Navarra against a panel of four human tumour cell lines: lung (HTB-54), colon (HT-29), leukaemia (CCRF-CEM) and breast adenocarcinoma (MCF-7). Cytotoxicity assays were performed based on the reactivity of MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide], as described by NCI [39]. The results are expressed as GI50 values, i.e., the concentration that reduces by 50% the growth of treated cells with respect to untreated controls, TGI, the dose that completely inhibits cell growth, and LD50, the concentration that kills 50% of the cells. The cytotoxic effect of each substance was tested at five different doses between 0.01 and 100 μM, or at lower levels when the GI50 was less than 10 nM. Mean GI50, TGI, and LD50 values are summarized in Table 7. Doxorubicin was used as a control. As guidance with regard to selectivity, all of the compounds were further examined for toxicity in a mammary gland cell culture derived from non-malignant cells (184B5). Drug concentrations ranged from 0.01 to 100 μM.
Table
Table 7. Cytotoxic activities (average GI50a, TGIb and LD50c values) for compounds against tumour cell lines.
Table 7. Cytotoxic activities (average GI50a, TGIb and LD50c values) for compounds against tumour cell lines.
Comp.C.Pa (μM)Cell lines
CCRF-CEMHTB-54HT-29MCF-7184B5
1bGI50>1000.587.390.091.97
cTGI>10042.4755.983.6922.50
dLD50>100>100>10058.0176.73
2GI50>1001.908.952.642.29
TGI>1009.9151.956.5420.14
LD50>100>100>10023.1278.91
11GI50>10011.1416.190.0061.05
TGI>100>10072.194.189.64
LD50>100>100>10053.4972.34
14GI50>10016.825.991.570.0009
TGI>100>10047.426.198.21
LD50>100>100>10052.6562.55
21GI50>10011.797.530.0030.05
TGI>100>10060.948.271.89
LD50>100>100>10089.217.20
Doxorub.GI500.033<0.01nde0.88nd
TGI0.0711.25nd>100nd
LD500.293.45nd>100nd
a Cytotoxic parameters; b Dose that inhibits 50% of cell growth; c Dose that completely inhibits cell growth; d Dose that kills 50% of cells; e n.d.: not determined.
The data show that the compounds under investigation influenced tumour cell growth differently depending on the cell line. The human breast adenocarcinoma cell line (MCF-7) was most sensitive to the antiproliferative effects of the investigated compounds. In particular, 1 (GI50 = 0.09 μM), 11 (GI50 = 0.006 μM), 14 (GI50 = 1.57 μM ) and 21 (GI50 = 0.003 μM) were strongly antiproliferative. Compounds 1, 11 and 21 were 10, 147 and 293 times more active, respectively, than standard doxorubicin (GI50 = 0.88 μM). Besides, if we compare the TGI values for these compounds (TGI = 3.69, 4.18, 6.19 and 8.27 μM, respectively), all of them are lower than that of doxorubicin. In addition, 11 and 21 exhibited a better antitumoral profile than paclitaxel (GI50 = 0.010 ± 0.04 and TGI = 20 ± 1.9 μM) and vinorelbine (GI50 = 5 ± 1.3 and TGI = 100 ± 10.2 μM), both drugs that are widely used in clinic. The least sensitive cell line was CCRF-CEM leukemia with cytotoxic parameters >100 μM for all the compounds. HT-29 and HTB-54 showed moderate levels of susceptibility to the compounds with the exceptions of 1 (GI50 = 7.39 and 0.58 μM) and 2 (GI50 = 8.95 and 1.90 μM).
As an example, curves with the original data from which the GI50, TGI and LD50 values for compounds 1 and 2 were calculated are shown in Figure 6. Fortunately, the LD50 in 184B5 is higher than in MCF-7 and the selectivity indexes are 1.3 and 3.4, respectively (184B5/MCF-7 LD50 ratio).
Molecules 14 03313 g006 550
Figure 6. Cytotoxic effects of 1 and 2 on CCRF-CEM, HTB-54, HT-29, MCF-7 and 184B5 cells. Cells were incubated in the presence of each compound at the indicated concentration for 72 hours. Cytotoxicity was then determined by a colorimetric microassay based on the use of MTT. Data are expressed as percentage growth ± SEM from at least 3 independent experiments performed in quadruplicate.
Figure 6. Cytotoxic effects of 1 and 2 on CCRF-CEM, HTB-54, HT-29, MCF-7 and 184B5 cells. Cells were incubated in the presence of each compound at the indicated concentration for 72 hours. Cytotoxicity was then determined by a colorimetric microassay based on the use of MTT. Data are expressed as percentage growth ± SEM from at least 3 independent experiments performed in quadruplicate.
Molecules 14 03313 g006
It is evident that compound 1 displayed the most prominent antiproliferative activity against all treated cell lines compared to the untreated control. This compound is therefore a representative candidate for these compounds for preliminary studies in order to explore a possible mechanism of action.

Apoptosis

Mounting evidence indicates that apoptosis is a critical mechanism for cancer prevention by Se compounds [40,41]. For this reason, we investigated whether apoptosis was involved in cell growth inhibition in the MCF-7 cells by 1.
The apoptotic status of the cells after 48 hours of treatment with 25 μM of the corresponding compound was determined using the Apo-Direct kit (BD Pharmingen) [42] based on the TUNEL technique. Camptothecin was used as a positive control. The results obtained are shown in Figure 7. As can be seen, for compound 1 the induction of cell death was independent of the apoptotic process.
Effects on cell cycle progression
Cell cycle arrest is one of the targets of many anticancer drugs, including doxorubicin and camptothecin. In an effort to ascertain whether compound 1 could affect cell cycle progression in MCF-7 cells, these cells were treated with 25 μM of the corresponding compound for 48 hours and cell cycle progression was determined by flow cytometry analysis [42]. As shown in Figure 8, DNA flow cytometric analysis indicated that treatment of the cells with compound 1 did not induce any specific phase arrest of the cell cycle.
Molecules 14 03313 g007 550
Figure 7. MCF-7 cells were incubated either with 25 µM of the indicated compound or vehicle (control cells) for 48 hours. The results are presented as the mean ± SEM of three independent experiments (duplicate wells). * p < 0.01 with respect to the control.
Figure 7. MCF-7 cells were incubated either with 25 µM of the indicated compound or vehicle (control cells) for 48 hours. The results are presented as the mean ± SEM of three independent experiments (duplicate wells). * p < 0.01 with respect to the control.
Molecules 14 03313 g007
Molecules 14 03313 g008 550
Figure 8. Cell cycle phase distribution of MCF-7 cells after 48 hours of treatment with 25 μM of the indicated compound or vehicle (control). Results are expressed as percentages of total cell counts. Each bar represents mean ± SEM of three independent experiments (duplicate wells).
Figure 8. Cell cycle phase distribution of MCF-7 cells after 48 hours of treatment with 25 μM of the indicated compound or vehicle (control). Results are expressed as percentages of total cell counts. Each bar represents mean ± SEM of three independent experiments (duplicate wells).
Molecules 14 03313 g008

Conclusions

In summary, we have synthesized a series of 26 aroyl and heteroaroyl selenylacetic acids and evaluated the cytotoxicity against a prostate cancer cell line (PC-3). Five compounds (1, 2, 11, 14 and 21) showed interesting activity and three of them (1, 11 and 14) were more powerful than etoposide and methylseleninic acid. Exploration of the SAR study suggested that the permissible maximum volume for the aromatic fragment is approximately 111–112 Å. A study of the ALogP data indicated that the compounds with values much greater than 1.90 were inactive. In terms of the length of the bridging chain between the aromatic moiety and the keto group, the presence of a methylene fragment improves the activity. We also tested some selected compounds against a panel of four human tumoral cell lines (CCRF-CEM, HTB-54, HT-29, MCF-7) and one mammary gland-derived non-malignant cell line (184B5). Compounds 1, 11 and 21 exhibited excellent growth inhibition activity against MCF-7, with values in the nanomolar range, and moderate values in HTB-54 and HT-29. However, these compounds did not show inhibitory potency against CCRF-CEM. Comparable activity was observed on two of the cell lines studied, one a hormone-independent prostate (PC-3) and the other a hormone-dependent MCF-7 breast cancer cell line. This could suggest that the mode of action is to a large extent unrelated to the expression of the androgen receptor (PC-3 negative) and oestrogen receptor (MCF-7 positive). Compound 1, which exhibited the best profile, was examined as an apoptotic agent and cell cycle modulator in MCF-7. The preliminary results revealed that the mechanism of action is independent of apoptosis and additional studies are required to determine the mechanism by which this compound exerts its activity. This work is currently in progress in our laboratory. The results reported here open up new perspectives for future investigations into the synthesis of new compounds of this family and these are potentially useful for the modulation of pharmacological potential with enhanced and selective anticancer activity.

Experimental

General

Melting points were determined by differential scanning calorimetry with a Perkin-Elmer DSC Diamond and by thermomicroscopy with a Mettler FP82+FP80 apparatus (Greifense, Switzerland) and have not been corrected. The 1H NMR spectra (Table 8) were recorded on a Bruker 400 UltrashieldTM (Rheinstetten, Germany), using TMS as the internal standard (Table 8). The IR spectra (Table 8) were performed on a Thermo Nicolet FT-IR Nexus in KBr pellets (Table 8). Elemental microanalyses were carried out on vacuum-dried samples (Table 1) using an Elemental Analyser (LECO, CHN-900 Elemental Analyser). Silica gel 60 (0.040–0.063 mm) 1.09385.2500 (Merck KGaA, 64271 Darmstadt, Germany) was used for Column Chromatography and Alugram® SIL G/UV254 (Layer: 0.2 mm) (Macherey-Nagel GmbH & Co. KG. Postfach 101352. D-52313 Düren, Germany) was used for Thin Layer Chromatography. Chemicals were purchased from E. Merck (Darmstadt, Germany), Scharlau (F.E.R.O.S.A., Barcelona, Spain), Panreac Química S.A. (Montcada i Reixac, Barcelona, Spain), Sigma-Aldrich Química, S.A. (Alcobendas, Madrid, Spain), Acros Organics (Janssen Pharmaceuticalaan 3a, 2440 Geel, België), Maybridge (Cambridge, CB5 8BZ, UK) and Lancaster (Bischheim-Strasbourg, France).
Table
Table 8. Spectroscopic data for compounds 1–26.
Table 8. Spectroscopic data for compounds 1–26.
Ref.IR (KBr, υ in cm-1)1H NMR (400 MHz, δ/ppm, J in Hz)
11702, 1678DMSO-d6, 3.81 (s, 2H, Se-CH2-COOH), 7.59 (m, 2H, H3 + H5), 7.75 (tt, 1H, H4, J4-5=J4-3=7.9, J4-6=J4-2=1.2), 7.90 (dd, 2H, H2 + H6, J6-5=J2-3=8.4), 12.80 (br s, 1H, COOH).
22232, 1681CDCl3, 3.91 (s, 2H, Se-CH2-COOH), 7.81 (d, 2H, H3 + H5, J3-2=J5-6=8.7), 8.00 (d, 2H, H2 + H6).
31694, 1712CDCl3 3.91 (s, 2H, Se-CH2-COOH), 7.78 (dd, 2H, H3+ H5, J3-2=J5-6=8.1, J3-CF3=J5-CF3=0.5), 8.02 (dd, 2H, H2 + H6, J2-CF3=J6-CF3=0.6).
41697, 1686DMSO-d6, 3.81 (s, 2H, Se-CH2-COOH), 7.66 (d, 2H, H3+H5; J3-2= J5-6=8.5), 7.91 (d, 2H, H2+H6), 12.89 (br s, 1H, COOH).
51691CDCl3, 2.43 (s, 3H, CH3), 3.85 (s, 2H, Se-CH2-COOH), 7.29 (d, 2H, H3+ H5, J3-2=J5-6= 8.1), 7.81 (d, 2H, H2 + H6).
61715, 1666CDCl3, 1.36 (s, 9H, C-(CH3)3), 3.85 (s, 2H, Se-CH2-COOH), 7.51 (d, 2H, H3+H5, J3-2= J5-6=8.7), 7.86 (d, 2H, H2+H6).
71707, 1692CDCl3, 3.83 (s, 2H, Se-CH2-COOH), 3.91 (s, 3H, CH3O), 6.97 (d, 2H, H3+H5, J3-2= 9.0), 7.90 (d, 2H, H2+H6).
81707, 1688CDCl3, 3.89 (s, 2H, Se-CH2-COOH), 7.37-7.43 (m, 1H, H5), 7.49-7.50 (m, 2H, H4+H3), 7.76-7.79 (m, 1H, H6).
91703, 1684CDCl3, 3.89 (s, 2H, Se-CH2-COOH), 7.40 (tt, 1H, H5, J5-6=1.8, J5-4=7.5 J5-3=8.9), 7.44 (tt, 1H, H4, J4-3=1.4 Hz, J4-6=7.6), 7.70 (dd, 1H, H3), 7.72 (dd, 1H, H6).
101701, 1677CDCl3, 3.90 (s, 2H, Se-CH2-COOH), 7.23 (dt, 1H, H5, J5-3=1.5, J5-6=7.9, J5-4=7.7), 7.48 (dt, 1H, H4,J4-6=0.4, J4-3=7.7), 7.69 (dd, 1H, H3); 8.00 (dd, 1H, H6).
111694, 1685CDCl3, 3.62 (s, 2H, Se-CH2-COOH), 3.90 (s, 2H, Ar-CH2-COSe), 7.32 (m, 2H, H3 + H5), 7.38 (m, 3H, H2 + H4 + H6).
121708, 1686CDCl3, 3.03 (s, 4H, Ar-CH2-CH2-COSe + Ar-CH2-CH2-COSe), 3.69 (s, 2H, Se-CH2-COOH), 7.20-7.26 (m, 3H, H3+H4+H5), 7.30-7.34 (m, 2H, H2+H6).
131699, 1667CDCl3, 3.84 (s, 2H, Se-CH2-COOH), 7.57 (s, 1H, H4), 7.70 (s, 2H, H2+H6), 8.44 (br s, 1H, COOH).
141714, 1696CDCl3, 3.85 - 3.86 (s + s, 8H, 2 OCH3 + Se-CH2-COOH, JCH2-Se=72.2), 6.72 (dt, 1H, H4, J4-2=J4-6=0.8 Hz, J4-OCH3=2.3), 7.04 (dd, 2H, H2 + H6, J2-OCH3=J6-OCH3=2.2).
151703, 1671CDCl3, 3.86 (s, 2H, Se-CH2-COOH), 3.94 (s, 6H, 3,5-diCH3O), 3.95 (s, 3H, 4-CH3O), 7.16 (s, 2H, H2).
161701, 1675DMSO-d6, 3.76 (s, 2H, Se-CH2-COOH), 6.18 (s, 2H, O-CH2-O), 7.09 (dd, 1H, H5, J5-2= 0.9 Hz, J5-6=7.8), 7.34 (dd, 1H, H2, J2-6=2.9), 7.55 (dd, 1H, H6), 12.79 (br s, 1H, COOH).
171719, 1671CDCl3, 3.92 (s, 2H, Se-CH2-COOH), 4.40 (br s, 1H, COOH), 7.63 (m, 2H, H4+ H5), 7.91 (m, 3H, H3 + H6 + H7), 8.00 (d, 1H, H8, J8-7=8.1), 8.48 (s, 1H, H2).
181720, 1693CDCl3, 3.68 (s, 2H, Se-CH2-COOH), 5.23 (s, 1H, (Ph)2-CH-COSe), 7.31-7.40 (m, 10 H, 2H2, 2H6, 2H3, 2H5, 2 H4).
191717, 1659 DMSO–d6, 3,87 (s, 2H, Se-CH2-COOH), 7,78 (td, 2H, H3+H5), 8,85 (td, 2H, H2+H6).
201713, 1673DMSO-d6, 3.87 (s, 2H, Se-CH2-COOH), 7.61-7.65 (m, 1H, H5), 8.26-8.30 (m, 1H, H4), 8.87-8.90 (m, 1H, H6), 9.02-9.03 (m, 1H, H2), 12.84 (br s, 1H, COOH).
211724, 1690DMSO-d6, 3.87 (s, 2H, Se-CH2-COOH), 7.63 (ddd, 1H, H5, J5-4=7.7, J5-6=4.8, J5-Cl= 1.1), 8.24 (ddd, 1H, H4, J4-Cl=1.1), 8.63 (ddd, 1H, H6, J6-Cl=1.1), 12.90 (br s, 1H, COOH).
221700, 1668DMSO-d6, 0.97 (t, 3H, CH3), 1.58-1.70 (m, 2H, S-CH2-CH2-CH3), 3.11 (dt, 2H, S-CH2-CH2-CH3), 3.81 (s, 2H, Se-CH2-COOH), 7.32 (dd, 1H, H4, J4-5=4.8, J4-6=7.8), 7.97 (dt, 1H, H5, J5-6=1.7), 8.29 (dd, 1H, H6), 12.89 (br s, 1H, COOH).
231712, 1643DMSO-d6, 3.80 (s, 2H, Se-CH2-COOH), 7.30 (dd, 1H, H4, J4-5=3.9 Hz,J4-3=4.9 Hz), 8.02 (dd, 1H, H5, J5-3=1.1 Hz), 8.16 (dd, 1H, H3), 12.82 (bs, 1H, COOH).
241703, 1674DMSO-d6, 3.74 (s, 2H, Se-CH2-COOH), 8.87 (dd, 1H, H6, J6-3=1.2 Hz, J6-5=2.4 Hz), 9.05 (dd, 1H, H5, J5-3=0.2 Hz), 9.06 (dd, 1H, H3), 12.73 (br s, 1H, COOH).
251714, 1690DMSO-d6, 3.73 (s, 2H, Se-CH2-COOH), 7.82 (t, 1H, H7, J7-6=7.5, J7-8=8.0), 7.94 (t, 1H, H6, J6-5=8.3), 7.97 (d, 1H, H5), 8.15 (d, 1H, H8), 8.19 (d, 1H, H3, J3-4=8.4), 8.66 (d, 1H, H4), 12.68 (br s, 1H, COOH).
261709, 1675DMSO-d6, 3,91 (s, 2H, Se-CH2-COOH), 7,77 (t, 1H, H6, J6-7=7,6 ), 7,98 (t, 1H, H7), 8,13 (d, 1H, H5, J5-8=8,3), 8,31 (d, 1H, H8), 9,07 (d, 1H, H4, J4-2=1,8), 9,23 (d, 1H, H2),12,75 (bs, 1H, COOH).

General procedure for the preparation of selenylacetic acids 1–6, 8–26

A solution of sodium borohydride (1.00 g, 26.3 mmol) in distilled water (12.5 mL) was added to a stirred suspension of grey selenium (1.00 g, 12.6 mmol) in distilled water (12.5 mL) at room tempera-ture. The reaction mixture was stirred until an almost colourless solution of NaHSe was formed. The acyl chloride (12.6 mmol) was added in small portions and the reaction mixture was magnetically stirred for 1 hour. A yellow solution was formed and the bromoacetic acid (1.75 g, 12.6 mmol) was added. Within 30 minutes a solid was formed. The product was filtered off and washed with distilled water (3 × 25 mL). The products were recrystallised from an appropriate solvent (Table 1). The following compounds were synthesized according to this procedure: (Benzoylselenyl)acetic acid (1). (4-cyanobenzoylselenyl)acetic acid (2). (4-trifluoromethylbenzoylselenyl)acetic acid (3). (4-chlorobenzoylselenyl)acetic acid (4). (4-methylbenzoylselenyl)acetic acid (5). (4-tert-butylbenzoyl-selenyl)acetic acid (6). (2-chlorobenzoylselenyl)acetic acid (8). (2-bromobenzoylselenyl)acetic acid (9). (2-iodobenzoylselenyl)acetic acid (10). Phenylacetoylselenylacetic acid (11). (3-phenyl-propanonylselenyl)acetic acid (12). (3,5-dichlorobenzoylselenyl)acetic acid (13). (3,5-dimethoxy-benzoylselenyl)acetic acid (14). (3,4,5-trimethoxybenzoylselenyl)acetic acid (15). (3,4-methylenedioxy-benzoylselenyl) acetic acid (16). (2-naphthoylselenyl)acetic acid (17). (diphenylacetylselenyl)acetic acid (18). (4-pyridoylselenyl)acetic acid (19). (3-pyridoylselenyl)acetic acid (20). (3-(2-chloro)-pyridoylselenyl)acetic acid (21). (3-(2-propylthio)pyridoylselenyl)acetic acid (22). (2-thienoylselenyl)-acetic acid (23). (pyrazinoylselenyl)acetic acid (24). (2-quinoloylselenyl)acetic acid (25). (3-quinoloylselenyl)acetic acid (26).

Procedure for the preparation of (4-methoxyphenylselenyl)acetic acid (7)

An equivalent of grey-powder selenium was suspended in tetrahydrofuran (20 mL) and two equivalents of lithium aluminium hydride were added slowly. The carbonyl chloride was added to the reaction mixture and the mixture was stirred for 1 h. The suspension was filtered and the solid was discarded. An equivalent of bromoacetic acid was added to the liquid filtrate. After 20 min the mixture was poured onto 500 g of ice. The resulting solid was filtered off, washed and recrystallised (Table 1).

Cytotoxic activity in PC-3

PC-3 cells were seeded in 96-well plates (Millipore, Eschborn, Germany) at a density of 5 × 103 cells per well. The samples were incubated at 37 ºC under 5% CO2 overnight prior to the addition of the compounds. Compounds were diluted in complete medium. After 3 days of incubation, 10 μL MTT solution (5 mg/mL in PBS) was added to the cells in each well and these were stored for an additional 4 h at 37 ºC. The absorbance of formazan at λ = 570 nm was measured on a Polarstar Galaxy plate reader (BMG LabTechnologies GmbH, Offenburg, Germany). The percentage of viable cells was calculated to obtain IC50-values.
PC-3 are human tumorigenic and metastatic prostate cancer cells and these were obtained from American Type Culture Collection (ATCC, Manassas, VA, USA) (passage 36). The cells were cultured under standard conditions (Dulbecco’s RPMI 1640 medium, with GlutamaxTM 1, Invitrogen supplemented with 10% fetal bovine serum, Fetalclone III, SH30109.03, HYCLONE and 1% Penicillin-Streptomycin, Invitrogen, Carlsbad, CA, USA).

Molecular modelling

The initial computational work was performed on a Dell Precision 380 workstation, provided with the software package Discovery Studio v1.7. The three-dimensional models of the studied compounds were constructed, in the vacuum phase, using atoms and structural fragments from the Viewer module (Discovery Studio) and using the Dreiding force field [43]. Once the models were constructed, a preliminary conformational analysis was carried out. The applied protocol (Diverse Conformational Generation integrated in the Pharmacophore protocol Discovery Studio) can be summed up as follows: (a) Initial construction of the model and first minimization by application of the Dreiding minimize protocol (Steepest descent algorithm with a convergence criterion of 10–6). The AlogP98 [44,45] descriptor (an implementation of the atom-based ALogP method) was calculated for each compound. (b) Application of the BEST routine for conformation generation (First: Conjugate-gradient minimization in torsion space; second: conjugate-gradient minimization in Cartesian space; third: Quasi_Newton minimization in Cartesian space). (c) Elimination of those conformations whose relative energy is greater than 10 Kcal/mol at a global minimum. (d) Analysis of conformational trajectory and selection of representative lowest energy conformations. Root mean square (rms) deviations of the structures were monitored. The energy differences between the different conformations analysed for each trajectory were in the range 2–5 Kcal.
For each of the compounds, ten lowest energy conformations were selected and a new minimization cycle was applied. The volumes of the ring moiety and the whole molecule were calculated for each of the new representative low-energy conformations selected.
The mechano-quantic analysis of the conformations obtained in the previous step was carried out with the package Mopac2009, PM6 [46] (or PM3 [47] for halogen and/or sulfur-containing derivatives) semi-empirical approaches, with the geometry optimized using an eigenvector following algorithm. The energy and distribution of the HOMO and LUMO orbitals, the selenium charge, the electronic density (ED) and molecular electrostatic potential (MEP) distribution and the Dipolar moment were obtained for each of the conformations obtained. The data corresponding to the representative lowest energy conformation for each compound was selected and used in the establishment of the preliminary structure-activity relationships.

Cytotoxic activity in CCRF-CEM, HTB-54, HT-29, MCF-7 and 184B5

The cytotoxic effects of each substance were tested at five different doses between 0.01 and 100 μM. Each substance was initially dissolved in DMSO at a concentration of 0.1 M, and serial dilutions were prepared using culture medium. The plates with cells from the different lines, to which medium containing the substance under test were added, were incubated for 72 h at 37 °C in a humidified atmosphere containing 5% CO2. European Collection of Cell Cultures (ECACC) or American Type Culture Collection (ATCC) provided human tumour cell lines. Four cell lines were used: one human lymphocytic leukemia (CCRF-CEM) and three human solid tumours, one colon carcinoma (HT-29), one lung carcinoma (HTB-54) and one breast adenocarcinoma (MCF-7). CCRF-CEM, HT-29 and HTB-54 cells were grown in RPMI 1640 medium (Life Technologies, Carlsbad, CA, USA) supplemented with 10% foetal calf serum, 2 mM L-glutamine, 100 units/mL penicillin, 100 μg/mL streptomycin and 10 mM HEPES buffer (pH = 7.4). MCF-7 cells were grown in EMEM medium (Clonetics) supplemented with 10% foetal calf serum, 2 mM L-glutamine, 100 units/mL penicillin and 100 μg/mL streptomycin. 184B5 cells were grown Hams F-12/DMEM (50:50) supplemented as described by Li et al. [48]. Cytotoxicity was then determined by the MTT method. Results are expressed as GI50 values, the concentration that reduces by 50% the growth of treated cells with respect to untreated controls, TGI, the dose that completely inhibits cell growth, and LD50, the concentration that kills 50% of the cells. Data were obtained from at least three independent experiments performed in quadruplicate.

Apoptosis and cell cycle

For breast adenocarcinoma MCF-7 cells, the apoptotic status and cell cycle analysis of the cells were determined using the Apo-Direct kit (BD Pharmingen), based on the TUNEL technique, under the conditions described by the manufacturer.

Acknowledgements

The authors wish to express their gratitude to the University of Navarra Research Plan (Plan de Investigación de la Universidad de Navarra, PIUNA) and CAN Foundation for financial support for this project.
  • Sample Availability: Samples of the compounds 1–26 are available from the authors.

References

  1. Varmus, H. The new era in cancer research. Science 2006, 312, 1162–1165. [Google Scholar] [CrossRef]
  2. Jariwalla, R.J.; Gangapurkar, B.; Nakamura, D. Differential sensitivity of various human tumour-derived cell types to apoptosis by organic derivatives of selenium. Br. J. Nutr. 2009, 101, 182–189. [Google Scholar] [CrossRef]
  3. Latreche, L.; Chavatte, L. Selenium incorporation into selenoproteins, implications in human health. In 10th International Symposium on Metal Ions in Biology and Medicine, Bastia, France; 2008; pp. 731–737. [Google Scholar]
  4. Clark, L.C.; Combs, G.F. Selenium-compounds and the prevention of cancer-research needs and public health implications. J. Nutr. 1986, 116, 170–173. [Google Scholar]
  5. Sonn, G.A.; Aronson, W.; Litwin, M.S. Impact of diet on prostate cancer: A review. Prostate Cancer Prostatic Dis. 2005, 8, 304–310. [Google Scholar] [CrossRef]
  6. El-Bayoumy, K.; Sinha, R.; Pinto, J.T.; Rivlin, R.S. Cancer chemoprevention by garlic and garlic-containing sulfur and selenium compounds. J. Nutr. 2006, 136, 864–869. [Google Scholar]
  7. Gonzalez, C.A.; Salas-Salvado, J. The potential of nuts in the prevention of cancer. Br. J. Nutr. 2006, 96, S87–S94. [Google Scholar] [CrossRef]
  8. Sharma, A.K.; Sharma, A.; Desai, D.; Madhunapantula, S.V.; Huh, S.J.; Robertson, G.P.; Amin, S. Synthesis and anticancer activity comparison of phenylalkyl isoselenocyanates with corresponding naturally occurring and synthetic isothiocyanates. J. Med. Chem. 2008, 51, 7820–7828. [Google Scholar] [CrossRef]
  9. Sanmartín, C.; Plano, D.; Palop, J.A. Selenium compounds and apoptotic modulation: A new perspective in cancer therapy. Mini Rev. Med. Chem. 2008, 8, 1020–1031. [Google Scholar] [CrossRef]
  10. Bardia, A.; Tleyjeh, I.M.; Cerhan, J.R.; Sood, A.K.; Limburg, P.J.; Erwin, P.J. Montori, V.M. Efficacy of antioxidant supplementation in reducing primary cancer incidence and mortality: Systematic review and meta-analysis. Mayo Clin. Proc. 2008, 83, 23–34. [Google Scholar] [CrossRef]
  11. Bjelakovic, G.; Nikolova, D.; Simonetti, R.G.; Gluud, C. Systematic review: Primary and secondary prevention of gastrointestinal cancers with antioxidant supplements. Alimen. Pharmacol. Ther. 2008, 28, 689–703. [Google Scholar] [CrossRef]
  12. Connelly-Frost, A.; Poole, C.; Satia, J.A.; Kupper, L.L.; Millikan, R.C.; Sandler, R.S. Selenium, folate, and colon cancer. Nutr. Cancer 2009, 61, 165–178. [Google Scholar] [CrossRef]
  13. Smith, M.L.; Lancia, J.K.; Mercer, T.I.; Ip, C. Selenium compounds regulate p53 by common and distinctive mechanisms. Anticancer Res. 2004, 24, 1401–1408. [Google Scholar]
  14. Rikiishi, H. Apoptotic cellular events for selenium compounds involved in cancer prevention. J. Bioenerg. Biomembr. 2007, 39, 91–98. [Google Scholar] [CrossRef]
  15. Hurst, R.; Elliot, R.M.; Goldson, A.J.; Fainweather-Tait, S.J. Se-methylselenocysteine alters collagen gene and protein expression in human prostate cells. Cancer Lett. 2008, 269, 117–126. [Google Scholar] [CrossRef]
  16. Pinto, J.T.; Sinha, R.; Papp, K.; Facompre, N.D.; Desai, D.; El-Bayoumy, K. Differential effects of naturally occurring and synthetic organoselenium compounds on biomarkers in androgen responsive and androgen independent human prostate carcinoma cells. Int. J. Cancer 2007, 120, 1410–1417. [Google Scholar]
  17. Singh, U.; Null, K.; Sinha, R. In vitro growth inhibition of mouse mammary epithelial tumor cells by methylseleninic acid: Involvement of protein kinases. Mol. Nutr. Food Res. 2008, 52, 1281–1288. [Google Scholar] [CrossRef]
  18. Poerschke, R.L.; Franklin, M.R.; Moos, P.J. Modulation of redox status in human lung cell lines by organoselenocompounds: Selenazolidines, selenomethionine, and methylseleninic acid. Toxicol. Vitro 2008, 22, 1761–1767. [Google Scholar] [CrossRef]
  19. Talas, Z.S.; Ozdemir, I.; Yilmaz, I.; Gok, Y. Antioxidative effects of novel synthetic organoselenium compound in rat lung and kidney. Ecotoxicol. Environ. Saf. 2009, 72, 916–921. [Google Scholar] [CrossRef]
  20. Shiah, H.S.; Lee, W.S.; Juang, S.H.; Hong, P.C.; Lung, C.C.; Chang, C.J.; Chou, K.M.; Chang, F.Y. Mitochondria-mediated and p53-associated apoptosis induced in human cancer cells by a novel selenophene derivative, D-501036. Biochem. Pharmacol. 2007, 73, 610–619. [Google Scholar]
  21. Xing, F.; Li, S.L.; Ge, X.Y.; Wang, C.; Zeng, H.H.; Li, D.; Dong, L. The inhibitory effect of a novel organoselenium compound BBSKE on the tongue cancer Tca8113 in vitro and in vivo. Oral Oncol. 2008, 44, 963–969. [Google Scholar] [CrossRef]
  22. Barbosa, N.B.D.; Nogueira, C.W.; Guecheva, T.N.; Bellinaso, M.D.; Rocha, J.B.T. Diphenyl diselenide supplementation delays the development of N-nitroso-N-methylurea-induced mammary tumors. Arch. Toxicol. 2008, 82, 655–663. [Google Scholar] [CrossRef]
  23. Xiao, H.; Parkin, K.L. Induction of Phase II enzyme activity by various selenium compounds. Nutr. Cancer 2006, 55, 210–223. [Google Scholar] [CrossRef]
  24. Rosa, R.M.; Picada, J.D.; Saffi, J.; Henriques, J.A.P. Cytotoxic, genotoxic, and mutagenic effects of diphenyl diselenide in Chinese hamster lung fibroblast. Mut. Res. 2007, 628, 87–98. [Google Scholar] [CrossRef]
  25. Zhao, G.H.; Lin, H.K.; Zhu, S.R.; Sun, H.W.; Chen, Y.T. Dinuclear palladium(II) complexes containing two monofunctional [Pd(en)(pyridine)Cl]+ units bridged by Se or S. Synthesis, characterization, cytotoxicity and kinetic studies of DNA-binding. J. Inorg. Biochem. 1998, 70, 219–226. [Google Scholar] [CrossRef]
  26. Abele, E.; Popelis, J.; Shestakova, I.; Domracheva, I.; Arsenyan, P.; Lukevics, E. Synthesis and cytotoxicity of di(3-indolyl) selenide derivatives. Khimiya Geterotsiklicheskikh Soedinenii 2004, 6, 868–872. [Google Scholar]
  27. Zeng, H.W.; Briske-Anderson, M.; Idso, J.P.; Hunt, C.D. The selenium metabolite methylselenol inhibits the migration and invasion potential of HT1080 tumor cells. J. Nutr. 2006, 136, 1528–1532. [Google Scholar]
  28. Ohta, Y.; Suzuki, K.T. Methylation and demethylation of intermediates selenide and methylselenol in the metabolism of selenium. Toxicol. Appl. Pharmacol. 2008, 226, 169–177. [Google Scholar] [CrossRef]
  29. Plano, D.; Sanmartín, C.; Moreno, E.; Prior, C.; Calvo, A.; Palop, J.A. Novel potent organoselenium compounds as cytotoxic agents in prostate cancer cells. Bioorg. Med. Chem. Lett. 2007, 17, 6853–6859. [Google Scholar] [CrossRef]
  30. Plano, D.; Moreno, E.; Font, M.; Encío, I.; Palop, J.A.; Sanmartín, C. Synthesis and biological evaluation of selenadiazole derivatives as anticancer agents. Bioorg. Med. Chem. 2009. (Submitted for publication).. [Google Scholar]
  31. De Athayde, P.F.; De Souza, A.G.; De Morais, S.A.; Botelho, J.R.; Barbosa, J.M.; Miller, J.; Lira, B.F. Synthesis and characterization of three new organo-selenium compounds. A convenient synthesis of aroylselenoglycolic acids. Arkivoc 2004, 6, 22–26. [Google Scholar]
  32. Mosmann, T. Rapid colorimetric assay for cellular growth and survival: Application to proliferation and cytotoxicity assays. J. Inmunol. Methods 1983, 65, 55–63. [Google Scholar] [CrossRef]
  33. Dong, Y.; Zhang, H.; Hawthorn, L.; Ganther, H.E.; Ip, C. Delineation of the molecular basis for selenium-induced growth arrest in human prostate cancer cells by oligonucleotide array. Cancer Res. 2003, 63, 52–59. [Google Scholar]
  34. Van Brussel, J.P.; Oomen, M.A.; Vosselbeld, P.J.M.; Wiemer, E.A.C.; Sonneveld, P.; Mickisch, G.H.J. Identification of multidrug resistance-associated protein 1 and glutathione as multidrug resistance mechanisms in human prostate cancer cells: chemosensitization with leukotriene D4 antagonists and buthionine sulfoximine. BJU Int. 2004, 93, 1333–1338. [Google Scholar] [CrossRef]
  35. Font, M.; Ardaiz, E.; Cordeu, L.; Cubedo, E.; García-Foncillas, J.; Sanmartín, C.; Palop, J.A. Structural characteristics of novel symmetrical diaryl derivatives with nitrogenated functions. Requirements for cytotoxic activity. Bioorg. Med. Chem. 2006, 14, 1942–1948. [Google Scholar]
  36. Accerlys Software Inc. Discovery Studio Modeling Environment Release 1.7.; Accerlys Software Inc.: San Diego, USA, 2007.
  37. MOPAC MOPAC2009, Stewart, J.P. Stewart Computational Chemistry, Colorado Springs, CO, USA. Available online: http://OpenMOPAC.net,2008.
  38. Macrae, C.F.; Edgington, P.R.; McCabe, P.; Pidcock, E.; Shields, G.P.; Taylor, R.; Towler, M.; Van de Streek, J. Mercury: Visualization and analysis of crystal structures. J. Appl. Crystallogr. 2006, 39, 453–457. [Google Scholar] [CrossRef]
  39. Available online: http://dtp.nci.nih.gov/branches/btb/ivclsp.html accessed on 27 June 2009.
  40. Philchenkov, A.; Zavelevich, M.; Khranovskaya, N.; Surai, P. Comparative analysis of apoptosis induction by selenium compounds in human lymphoblastic leukemia MT-4 cells. Exp. Oncol. 2007, 29, 257–261. [Google Scholar]
  41. Naithani, R. Organoselenium compounds in cancer chemoprevention. Mini Rev. Med. Chem. 2008, 8, 657–668. [Google Scholar] [CrossRef]
  42. Fimognari, C.; Nüsse, M.; Cesari, R.; Iori, R.; Cantelli-Forti, G.; Hrelia, P. Growth inhibition, cell-cycle arrest and apoptosis in human T-cell leukemia by the isothiocyanate sulforaphane. Carcinogenesis 2002, 23, 581–586. [Google Scholar] [CrossRef]
  43. Mayo, S.L.; Olafson, B.D.; Goddard, W.A. Dreiding: A generic force field for molecular simulations. J. Phys. Chem. 1990, 94, 8897–8909. [Google Scholar] [CrossRef]
  44. Viswanadhan, V.N.; Ghose, A.K.; Revankar, G.R.; Robins, R.K. Atomic physicochemical parameters for three-dimensional structure-directed quantitative structure–activity relationships. Additional parameters for hydrophobic and dispersive interactions and their application for an automated superposition of certain naturally-occurring nucleoside antibiotics. J. Chem. Inf. Comput. Sci. 1989, 29, 163–172. [Google Scholar] [CrossRef]
  45. Ghose, A.K.; Viswanadhan, V.N.; Wendoloski, J.J. Prediction of hydrophobic (lipophilic) properties of small organic molecules using fragmental methods: An analysis of ALOGP and CLOGP methods. J. Phys. Chem. A 1998, 102, 3762–3772. [Google Scholar]
  46. Stewart, J.J.P. Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements. J. Mol. Model. 2007, 13, 1173–1213. [Google Scholar] [CrossRef]
  47. Stewart, J.J.P. Optimization of parameters for semi-empirical methods I-Method. J. Comp. Chem. 1989, 10, 209–220. [Google Scholar] [CrossRef]
  48. Li, Y.Z.; Pan, J.; Li, J.L.; Lee, J.H.; Tunkey, C.; Saraf, K.; Garbe, J.C.; Whitley, M.Z.; Jelinsky, S.A.; Stampfer, M.R.; Haney, S.A. Transcriptional changes associated with breast cancer occur as normal human mammary epithelial cells overcome senescence barriers and become immortalized. Mol. Cancer 2007, 6, Article 7. [Google Scholar]

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/molecules14093313

AMA Style

Sanmartín C, Plano D, Domínguez E, Font M, Calvo A, Prior C, Encío I, Palop JA. Synthesis and Pharmacological Screening of Several Aroyl and Heteroaroyl Selenylacetic Acid Derivatives as Cytotoxic and Antiproliferative Agents. Molecules. 2009; 14(9):3313-3338. https://doi.org/10.3390/molecules14093313

Chicago/Turabian Style

Sanmartín, Carmen, Daniel Plano, Enrique Domínguez, María Font, Alfonso Calvo, Celia Prior, Ignacio Encío, and Juan Antonio Palop. 2009. "Synthesis and Pharmacological Screening of Several Aroyl and Heteroaroyl Selenylacetic Acid Derivatives as Cytotoxic and Antiproliferative Agents" Molecules 14, no. 9: 3313-3338. https://doi.org/10.3390/molecules14093313

APA Style

Sanmartín, C., Plano, D., Domínguez, E., Font, M., Calvo, A., Prior, C., Encío, I., & Palop, J. A. (2009). Synthesis and Pharmacological Screening of Several Aroyl and Heteroaroyl Selenylacetic Acid Derivatives as Cytotoxic and Antiproliferative Agents. Molecules, 14(9), 3313-3338. https://doi.org/10.3390/molecules14093313

Article Metrics

Back to TopTop