Dong Cai1, ZhiHua Zhang2, Yu Chen3, YanYan Zhang4, YuQi Sun5, YiXia Gong6. 1. College of Public Basic Sciences, Jinzhou Medical University, Jinzhou 121001, China. caidong0804@163.com. 2. School of Chemical and Environmental Engineering, Liaoning University of Technology, Jinzhou 121001, China. bridge1026@163.com. 3. School of Life Science and Biopharmaceutics, Shenyang Pharmaceutical University, Shenyang 110016, China. gzweishengwu@126.com. 4. School of Chemical and Environmental Engineering, Liaoning University of Technology, Jinzhou 121001, China. zyy19981231@sina.com. 5. School of Chemical and Environmental Engineering, Liaoning University of Technology, Jinzhou 121001, China. cpusyq@jzmu.edu.cn. 6. College of Public Basic Sciences, Jinzhou Medical University, Jinzhou 121001, China. gongyixia_2006@163.com.
Abstract
Novel 18β-glycyrrhetinic acid derivatives possessing a carbamate moiety and structurally similar ester derivatives were developed and evaluated for their efficacy as antitumor inhibitors. In the cellular assays, most of the N-substituted carbamate derivatives at the C3-position exhibited potent activities. The results of SAR investigation revealed that the introduction of the morpholine group at the C30-COOH led to a significant loss of the inhibitory potency. Among the ester derivatives, the ester group at C3-position also determined a noticeable reduction in the efficacy. Compound 3j exhibited the most prominent antiproliferative activity against six human cancer cells (A549, HT29, HepG2, MCF-7, PC-3, and Karpas299). Furthermore, compound 3j exerted a moderate inhibiting effect on the ALK. The results of molecular docking analyses suggested that it could bind well to the active site of the receptor ALK, which was consistent with the biological data. These results might inspire further structural optimization of 18β-glycyrrhetinic acid aiming at the development of potent antitumor agents. The structures 4d, 4g, 4h, 4j, and 4n were studied by X-ray crystallographic analyses.
Novel 18β-glycyrrhetinic acid derivatives possessing a carbamate moiety and structurally similar ester derivatives were developed and evaluated for their efficacy as antitumor inhibitors. In the cellular assays, most of the N-substituted carbamate derivatives at the C3-position exhibited potent activities. The results of SAR investigation revealed that the introduction of the morpholine group at the C30-COOH led to a significant loss of the inhibitory potency. Among the ester derivatives, the ester group at C3-position also determined a noticeable reduction in the efficacy. Compound 3j exhibited the most prominent antiproliferative activity against six humancancer cells (A549, HT29, HepG2, MCF-7, PC-3, and Karpas299). Furthermore, compound 3j exerted a moderate inhibiting effect on the ALK. The results of molecular docking analyses suggested that it could bind well to the active site of the receptor ALK, which was consistent with the biological data. These results might inspire further structural optimization of 18β-glycyrrhetinic acid aiming at the development of potent antitumor agents. The structures 4d, 4g, 4h, 4j, and 4n were studied by X-ray crystallographic analyses.
Glycyrrhetinic acid (GA) is a major bioactive triterpenoid metabolite of glycyrrhizin (GL) extracted from the roots of Glycyrrhiza uralensis Fisch (a particular species of licorice). GA exists as 18α and 18β-GA stereo-isomeric forms. Chemical structural differences between 18α-GA and 18β-GA lie in the spatial arrangement of hydrogen atom of C18-position. 18β-GA has been most extensively investigated and used not only for its abundance in root extract, but also the structural resemblance between 18β-GA and corticosteroids [1]. 18β-GA derivatives exhibit diverse pharmacological properties, such as anti-inflammation [2], antiulcer [3], antivirus [4], antitumor [5], antihepatotoxic [6], antibacterial [7], and antidiabetic activities [8].18β-GA derivatives exhibit prominent chemopreventive activities in various experimental cancer models [9,10,11,12]. The studies reported that 18β-GA derivatives have suggested protective effects against carcinogenic and tumorigenic factors by modulating the enzymatic antioxidant system and the attachment of carcinogenic factors to DNA or their receptors.Besides, the proapoptotic mechanisms of 18β-GA have been extensively studied over the past few decades. 18β-GA derivatives display anti-proliferative and pro-apoptotic effects against humanpituitary adenoma cells (GH3, MMQ) [13], breast cancer (MCF-7) [14], prostate cancer (DU-145) [15], ovarian cancer (SiHa, SK-OV-3, OVKAR-3) [16], lung cancer (A549, NCI-H460) [17], promyelotic leukemia (HL-60) [9], stomach cancer (KATO III) [18], hepatic cancer cells (HepG2, LX-2) [9,18], etc. The direct effects of 18β-GA derivatives can occur by suppressing tumor cells proliferation, with a noticeable accumulation of the tumor cells in the G1 phase, accompanied by a decrease in tumor cells in the S phase [18,19,20]. The antiproliferative activity transforms into cytotoxic effect when cell cycle arrest persists for long durations on several cancer lines [18]. There are also some 18β-GA derivatives that can exert anti-migratory and anti-invasive activities in humanbreast cancer cells (MDA-MB-231, MDA-MB-436) [21].18β-GA has been adopted as an attractive molecular scaffold to search for potential antitumor inhibitors. Current structural optimization of 18β-GA leading to antitumor agents primarily focused on modification of the C3-OH in ring A, 11-one in ring C, C30-COOH in ring-E and/or multi-fragment modified simultaneously (Figure 1). The results of SAR analyses revealed that the C3-OH is a critical structural feature. The modifications at the C3-OH, reducing the polarity of the entire molecule, resulted in the significant enhancement in the in vitro antiproliferative activity. Esterification of the C3-OH group induced an enhanced inhibition of chymotrypsin-like, trypsin-like, and caspase-like activities of the 20S proteasome [22,23]. Furthermore, the introduction of side chains containing substituted amino groups in the C3-OH position significantly affected the cytotoxic activities [24,25,26,27,28].
Figure 1
Structure of 18β-GA 1 and known derivatives A and B.
Carbamate derivatives (e.g., the steroid skeleton) have aroused scientific interest over the years for their antitumor activities [29,30,31,32,33]. This is because carbamate moiety can form extensive hydrophobic and hydrogen bonding interactions with binding sites. Bufalin-3-yl nitrogen containing carbamate derivative A exhibits robust antiproliferative activities. Oleanolic acid derivatives B partially act as dual inhibitors for both topoisomerase I and IIa [34]. According to the results, the carbamate moiety at C3-position had vital effect on the activity [35].To enhance antiproliferative activity of 18β-GA, a series of novel 18β-GA derivatives possessing a carbamate moiety was synthesized to delve into the effect of structural modifications at the positions of C3-OH and C30-COOH. Additional similar derivatives of esterification of the C3-OH were synthesized to explore the influence of introducing a substituted acetoxy moiety. The antiproliferative activities in vitro of the synthesized compounds were evaluated. Furthermore, docking simulation was also performed for exploring the binding mode of the active compound at the ALK active site.
2. Results and Discussion
2.1. Chemistry
The synthetic routes to compounds 2, 3a–3o, and 4a–4n are illustrated in Scheme 1. The 18β-GA 1 was activated by 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide hydrochlorid (EDCI), 1-hydroxybenzotriazole (HOBt), and triethylamine under reflux for 20 min, and then it underwent amidation reaction with morpholine to form amide 2. O-esterification of 1 or 2 with substituted isocyanates produced carbamate derivatives 3a–3o or 4a–4n, respectively.
When this amidation reaction was performed at ambient temperature, coupling of the commercially available 18β-GA with morpholine in the presence of EDCI, HOBt, and triethylamine initially formed the intermediate 2a, presumably because of the significant steric hindrance of the 1H-benzotriazol-1-yl group (Figure 2). Such intermediate 2a could be isolated from the reaction mixture in good yield. When the temperature rose to 80 °C, the intermediate 2a should not be isolated and then react with morpholine to give compound 2 by refluxing in CH3CN for 24 h. Results suggested that the reaction conditions (e.g., temperature, time, and solution) should be thoroughly controlled to make intermediate 2a completely transformed into amide 2.
Figure 2
(a) The proposed mechanism of amidation reaction. Energy minimization by MM2 was performed using ChemBio3D Ultra 14.0 (CambridgeSoft Corporation, 2014) software force field. (b) The picture of 3D structure was produced using PyMOL (Delano Scientific, USA).
Condensation of the carboxylic acids with substituted isocyanates can afford unstable carbamic carboxylic anhydrides, which were transformed into N-substituted amides after decarboxylation in the presence of bases (most often, 4-dimethylaminopyridine [36], N,N-diisopropylethylamine [37], trimethylamine [38]). This condensation reaction is a well-known method for a practical synthesis of N-substituted amides and peptide analogues. Nevertheless, compounds 3a–3o can be prepared in high yields from 18β-GA with substituted isocyanates in the absence of base catalyst. In this case, the competitive reactions of the C30-COOH and substituted isocyanates were not observed [38]. These methyl and cyclohexyl fragments in the α-position to C30-COOH provided steric hindrance that prevented steric hindrance that prevents the possible amidation reaction in this condensation. Obviously, an excess of substituted isocyanates should not be adopted in condensation for the possible amidation. Compounds 3a–3o were obtained at a molar ratio of 18β-GA/substituted isocyanates (1:1.2) in refluxing ethyl acetate.New derivatives of 18β-GA bearing a bulky ester moiety in position C3-OH were synthesized (Scheme 2). The hydroxyl group of 18β-GA underwent esterification with chloroacetic anhydride at 130 °C, and the resulting ester 5 was aminated with the secondary amines to obtain 6a–6d. Compounds 7a–7b can also be synthesized through the treatment of 18β-GA with a substituted acyl chlorides in the presence of bases (e.g., triethylamine). 1H-NMR and 13C-NMR spectra for all the prepared compounds can be seen in the Supplementary Materials.
Scheme 2
Reagents and conditions: (a) chloroacetic anhydride, 130 °C, 1 h; (b) substituted secondary amines, K2CO3/ I2, C2H5OH; (c) substituted acyl chlorides, ethyl acetate, NEt3, reflux, 24 h.
Due to the noticeable steric hindrance at the β-position of C3-OH (Figure 2), the esterification reactions of 18β-GA with substituted isocyanates required refluxing in ethyl acetate. For the identical reason, the reactions of 18β-GA with substituted benzyl chlorides or chloroacetic anhydride are easily are also affected by steric hindrance as well.On the whole, compound 5 bearing a weak electron withdrawing substituent (chloro group) at the α-position of ester group is stable, primarily existing in the form of keto tautomer. The 1H-NMR (400 MHz) spectrum in Chloroform-d of compound 5 (Figure 3a) displayed a quartet at δ 4.62 ppm attributed to the proton at C3-position and a doublet at 4.06 ppm attributed to CH2 protons at the α-position of ester group.
Figure 3
The 1H-NMR spectrum of (a): compound 5; (b): compound 7a.
In contrast, a potent conjugation between the phenyl ring and the enol moiety of compound 7a led to noticeable shift equilibrium toward the enol tautomer [39,40]. The 1H-NMR (400 MHz) spectrum in Chloroform-d of compound 7a (Figure 3b) is presented here to study enolate formation. A singlet at 3.58 ppm referred to the presence of the CH2 protons at position 1 of keto tautomer. The structure of the enol tautomer might form both the Z-configuration and E-configuration [41]. The chemical shifts of the proton 1′ on olefin of the enol tautomer were identified at 4.76 ppm (J = 7.1 Hz) and at 3.74 ppm (J = 4.1 Hz), respectively. Moreover, a noticeable singlet observed at δ 13.10 ppm was attributed to hydroxyl group of the enol tautomer. A small multiplet observed at 4.54–4.60 ppm (near proton at C3-position in the keto tautomer, of which the signal was observed at 4.46–4.52 ppm) could represent the proton at C3′-position of enol tautomer.
2.2. Subsection Crystal Structure Analysis of 4d, 4g, 4h, 4j, and 4n
The ORTEP of the compounds 4d, 4g, 4h, 4j, and 4n with thermal ellipsoids at 50% probability is shown in Figure 4. Table 1 represents crystal and experimental data of the molecules compound 4d, 4g, 4h, 4j, and 4n. Crystallographic data have been deposited in the Cambridge crystallographic data Center. “CCDC 1953934, 1953928, 1953941, 1913193, 1952936 contain the supplementary crystallographic data for this paper. These data can be obtained free of charge via http://www.ccdc.cam.ac.uk/conts/retrieving.html (or from the CCDC, 12 Union Road, Cambridge CB2 1EZ, UK; Fax: +44 1223 336033; E-mail: deposit@ccdc.cam.ac.uk)”.
Figure 4
ORTEP diagram of compounds 4d, 4g, 4h, 4j, and 4n at 50% probability.
Table 1
Details for the crystal structure determinations of compounds 4d, 4g, 4h, 4j, and 4n.
Entry
4d
4g
4h
4j
4n
CCDC number
1953934
1953928
1953941
1913193
1952936
Crystal color
colorless
colorless
colorless
colorless
colorless
Solution
CH3OH
CH3OH
CH3OH
CH3OH
CH3OH
Crystal system
Ortha
Orth
Orth
Orth
Orth
Space group
P212121
P212121
P212121
P212121
P212121
a/Å
14.5561(7)
14.5119(7)
14.5570(11)
14.1203(8)
13.7916(9)
b/Å
15.3392(8)
15.4564(8)
15.2530(11)
15.6010(9)
16.7118(11)
c/Å
17.0499(9)
17.0665(8)
17.0544(19)
18.3201(10)
17.5923(12)
α (°)
90
90
90
90
90
β (°)
90
90
90
90
90
γ (°)
90
90
90
90
90
Volume/Å3
3806.9(3)
3828.0(3)
3786.7(6)
4035.7(4)
4054.7(5)
Temperature/K
296
296
296
296
296
Z
4
4
4
4
4
Density (calculated)/g·cm−3
1.210
1.280
1.187
1.196
1.217
F000
1496
1568.0
1464.0
1560.0
1592.0
Independent reflections
9492
7350
6813
7157
7151
Goodness-of-fit on F2
1.053
1.029
1.017
1.023
1.035
Final R indices (I > 2σ (I))
R1 = 0.0478
R1 = 0.0478
R1 = 0.0505
R1 = 0.0615
R1 = 0.0509
wR2 = 0.1161
wR2 = 0.1097
wR2 = 0.1302
wR2 = 0.1622
wR2 = 0.1181
R indices (all data)
R1 = 0.0732
R1 = 0.0720
R1 = 0.0740
R1=0.1008
R1 = 0.0801
wR2 = 0.1287
wR2 = 0.1198
wR2 = 0.1459
wR2=0.1895
wR2 = 0.1343
Largest diff. peak and hole/ e.Å−3
0.351 and −0.436
0.801 and −0.857
0.408 and −0.280
0.267 and −0.180
0.395 and −0.251
a Orth, Orthorhombic.
The 3D superimposition of compounds 4d, 4g, 4h, 4j, and 4n (Figure 5) revealed that the geometries of the two structural scaffolds are virtually identical, and the main structural difference between these five structures is in orientation of the phenyl rings in the carbamate side chain.
Figure 5
3D superimposition crystal structure of compounds 4d (brown), 4g (blue), 4h (red), 4j (green), and 4n (yellow).
In contrast to the X-ray crystal structure of 18β-GA from acetone/H2O [42], the two structural scaffolds are perfectly superimposable, but the orientation of carbonyl groups at the 30-position are slightly different. The overlay diagrams of both conformers are depicted in Figure 6. These two pictures were produced using the ‘molecular overlay’ feature of Accelrys Discovery Studio, which gives an automated optimum overlay of a number of small molecules using equally weighted steric and electrostatic fields.
Figure 6
(a) Molecular structures and (b) overlay crystal structures of 18β-GA (Lines) compound 4d (Sticks).
2.3. In Vitro Cell Growth Inhibitoty Activity
The antiproliferative activities of all the synthesized compounds against HT-29 and A549 cells were evaluated by MTT assay. The growth inhibition of cancer cell line (%) from these tests are listed in Table 2.
Table 2
Antiproliferative activity of the target compounds (growth inhibition, %).
Entry
A549
HT29
2 μg/mL
20 μg/mL
2 μg/mL
20 μg/mL
3a
23.94
95.92
11.34
98.05
3b
13.96
77.35
14.12
87.72
3c
18.51
95.40
8.39
98.61
3d
9.93
55.43
13.73
56.25
3e
14.32
74.92
18.84
96.44
3f
11.36
80.21
15.58
89.22
3g
9.72
63.39
18.57
76.54
3h
34.44
56.72
11.78
58.59
3i
15.91
87.47
19.41
90.12
3j
17.74
88.31
24.79
99.56
3k
18.79
88.82
10.76
94.83
3l
9.74
53.73
6.06
47.44
3m
11.68
47.51
12.82
32.33
3n
20.08
97.07
18.99
97.65
3o
24.41
86.30
55.75
98.17
4a
0.00
0.00
8.88
18.34
4b
0.00
0.00
14.40
32.80
4c
0.00
0.00
7.41
2.76
4d
0.00
54.37
16.34
53.73
4e
0.00
0.00
7.22
13.25
4f
0.00
0.00
7.10
18.45
4g
0.00
1.02
14.34
27.12
4h
0.00
8.54
14.81
30.39
4i
0.00
4.46
13.70
27.32
4j
0.00
0.00
20.31
36.42
4k
4.47
0.00
13.04
14.91
4l
0.00
0.00
13.61
16.57
4m
0.00
8.11
15.69
23.22
4n
0.00
0.00
7.12
40.21
6a
0.00
0.00
0.00
0.00
6b
0.00
0.00
0.00
12.34
6c
0.00
0.00
0.00
0.00
6d
0.00
0.00
0.00
0.00
7a
11.35
27.45
33.27
57.45
7b
15.26
31.16
44.15
59.02
18β-GA
0.00
8.41
0.00
12.93
Crizotinib
67.61
97.29
68.80
99.79
In contrast to the results for the compounds 3a–3o, the C30-COOH was converted into a morpholine amide derivatives 4a–4n produced a marked loss of the inhibitory potency. The results listed in Table 1 demonstrated that the C30-COOH group in carbamate derivatives led to the potent antiproliferative activity of target compounds.Table 2 evidently shows that the compounds 3a–3o exhibited excellent antiproliferative activity against HT-29 and A549 cells, especially at the concentration of 20 μg/mL. Compared with Crizotinib, a powerful anticancer drug as positive control, compounds (3a, 3c, 3e, 3f, 3i, 3j, 3k, 3n, 3n, and 3o) exhibited cell growth inhibitory activity nearly the same as Crizotinib to HT-29 cells at the identical concentration of 20 μg/mL. It is noteworthy that these active compounds 3a–3l, and 3o either possess electron-withdrawing groups (e.g., F, Cl, Br, CF3O, and CF3) or weak electron-donating group (e.g., CH3) on their aromatic rings. In the meantime, compounds (3l, 3m), containing CH3O group, displayed a dramatic decrease in the activity, suggesting that the strong electron-donating group is not recommended.Besides, introducing halogen atom into the compounds can improve selectivity, intrinsic potency, and so on. In this paper, compounds containing halogen atom in their structure have shown better anticancer activity. Comparing the derivatives with various chloro-substitution positions on the phenyl ring, compounds (3a, 3b, 3c, 3e, 3f) containing 2,4-diCl, 4-Cl-3-CF3-, 3,5-diCl, 3-Cl, and 3-Cl-4-CH3 groups displayed better activity, whereas 4-Cl compounds 3d resulted in a slight decrease in potency. A similar phenomenon was observed for 4-bromo substituted compound 3g and 4-fluoro substituted compound 3h. It is therefore speculated that the presence of the chloro, bromo, and fluoro substituent in para-position of phenyl rings stimulated a detrimental effect on inhibitory activity of cancer cells.The carbamate moiety at the C3-position is considered to be critical to compounds 3a–3o. With the introduction of six-membered ring secondary amine groups to give compounds 6a–6d, the activity was enhanced dramatically as compared with that of compounds 3a–3o. In the meantime, the introduction of a substituted phenylacetoxy moiety resulted in a significant loss of cytotoxicity for the ester derivatives (7a, 7b) compared with the corresponding carbamate derivatives (3d, 3h).From the results obtained, the carboxyl group at the C30-position of 18β-GA is active essential group to improve the inhibitory potency. The introduction of the substituted phenyl carbamate moiety at the C3-position led to a consistent increase in the activity. The electronic effect of substituent group and the position of the substituent group on the phenyl ring significantly influenced anticancer activity. On the whole, the substituents on the phenyl ring, carbamate moiety, and carboxyl group are critical for inhibiting the growth of tumor cells.To deepen our research, eight compounds 3a–3c, 3e, 3j, 3k, 3n, and 3o were taken in assay for their antitumor potency indicated by IC50 values. Table 3 presents that compound 3j exhibited the greatest antiproliferative activity against six humancancer cells (A549, HT29, HepG2, MCF-7, PC-3, and Karpas299) with IC50 values of 2.81 μg/mL, 3.19 μg/mL, 5.55 μg/mL, 5.26 μg/mL, 5.96 μg/mL, and 5.59bμg/mL, respectively. The other seven compounds also exhibited a significant antiproliferative activity against six humancancer cells. Compared with 18β-GA, introduction of a carbamate moiety at the C3-position could significantly enhance inhibitory activity.
Table 3
Antiproliferative activity (IC50, μg/mL) of selected compound.
Entry
A549
HT29
HepG2
MCF-7
PC-3
Karpas299
3a
11.19
11.44
11.6
9.31
11.91
12.30
3b
11.91
11.97
9.25
8.92
12.07
18.57
3c
12.28
9.74
10.68
5.53
14.98
15.47
3e
11.34
10.91
11.95
8.51
9.75
7.31
3j
2.81
3.19
5.55
5.26
5.96
5.59
3k
2.93
6.47
11.25
7
10.33
9.16
3n
7.42
11.33
12.55
9.30
13.32
25.97
3o
3.24
7.30
4.75
7.86
10.69
11.44
18β-GA
>40
>40
>40
>40
>40
>40
Crizotinib
1.27
0.67
2.15
1.84
3.3
0.28
To effectively study whether compound 3j is concentration-dependent manner, the cell was treated using MTT methods, and five concentration gradients were selected. After 24, 48, or 72 h of interaction, the inhibition rate was determined. Compound 3j inhibited the proliferation of HepG2 cell in a significant concentration manner, whereas it did not occur in a time-dependent manner. The result was shown in Figure 7.
Figure 7
The relationship between different concentrations of compound 3j and cytotoxicity. Data are means ± SD of the inhibition (%) from three independent experiments.
Given the mentioned results, the most efficient compound 3j was selected for in depth study, and it was evaluated in different concentrations (0.064–40.0 μg/mL) towards HepG2 cell line and non-tumorigenic liver LO2 cell. A 48 h continuous drug exposure protocol was employed by the MTTcytotoxicity assay. The results verified that treatment with the upregulated dose of compound 3j had a significant inhibitory effect on HepG2 cell lines. In the meantime, compound 3j did not exhibit significant toxic action towards LO2 cells at relatively higher concentrations, suggesting that compound 3j might exhibit selective antiproliferative activity against humantumor cells. The result is presented in Figure 8.
Figure 8
Inhibitory effects of compound 3j on the proliferation of HepG2 and LO2 cells. Data are means ± SD of the inhibition (%) from three independent experiments.
2.4. Anaplastic Lymphoma Kinase (ALK) Activity
Based on the inhibitory activities in vitro, the compound 3j and 3k were taken for further evaluation of their anti-ALK potency. As shown in Table 4, the compound 3j displayed a moderate potency against ALK enzymatic activity with an IC50 value of 120.68 nM, comparable to that of the positive control Crizotinib (IC50 = 1.64 nM). Besides, compound 3k was less potent against the ALK with IC50 value of 278.21 nM. The mentioned results indicated that the inhibition of ALK is likely to be a mechanism for the antitumor effect of these carbamate derivatives, and this series of compounds further studies.
Table 4
Inhibitory effects of selected compounds on ALK activity.
Compound
3j
3k
Crizotinib
IC50 (nM)
120.68
278.21
1.64
2.5. Molecular Docking Studies
In order to gain the likely binding model of these carbamate derivatives with ALK, we performed docking analysis using Discovery Studio 3.5 software. The image files were generated by PyMOL. Using the crystal structure of the humanALK in complex with Crizotinib (PDB code: 2XP2), the overall binding interactions of the representative compound 3j was shown in Figure 9. Some differences were found in the binding mode between these two ligand. The pyridine-2-amine moiety in Crizotinib bound well to Met1199 and Glu1197 residues via two hydrogen bonds, respectively. Besides, the two phenyl rings of Crizotinib might form a π–π interaction with Leu1256 residue. When compound 3j was docked into the ALK active site, no interactions were observed with mentioned critical residues. As observed, the C30-COOH group was orientated toward the binding pocket, making a hydrogen bond with Lys1150 residue. These results comply with the moderate potency for compound 3j inhibition of ALK. The results of the docking study also implied that the C30-COOH group is conducive to the stability of binding conformation. Nevertheless, introducing an amine group (4a–4n) might break the hydrogen bond between the ligand and the receptor, resulting in compounds (4a–4n) with poor activities in cellular assay.
Figure 9
(a) Binding mode of Crizotinib (Stick) and the compound 3j (Line); (b) prediction of contacting residues in compound-protein interfaces. Hydrogen bonds are highlighted as yellow dashes.
3. Materials and Methods
3.1. Materials and Apparatus
Generally, all commercial reagents and solvents were used without additional purification. Analytical thin layer chromatography (TLC) was performed on Merck silica gel 60 F254, precoated on aluminum plates. Spot visualization was done by using use of UV light (254 nm, Shanghai Baoshan Gucun Electro-optical Instrument Factory, Shanghai, China). Melting points were recorded on a WRS-1B digital melting point apparatus (Shanghai Shenguang Instrument Co., Ltd., Shanghai, China) and were uncorrected. 1H- and 13C-NMR spectra were recorded on Agilent 400/54Premium Shielded NMR Magnet System (Agilent, Santa Clara, CA, USA) and chemical shifts were quoted in ppm, referenced to tetramethylsilane (TMS). Peak multiplicities are reported as follow: s (singlet), d (doublet), t (triplet), q (quartet), m (multiplet). The high-resolution mass spectrum (HRMS) were recorded on Agilent 6200 Series TOF and 6500 Series Q-TOF LC/MS System B.05.01. (B5125) in positive ion modes (Agilent).
Single crystals of compound 3j was recrystallized by slow evaporation from methanol. Data collection was performed at 296 on a Bruker SMART APEX II CCD diffractometer (Bruker BioSpin, Rheinstetten, Germeny) using Mo Ka (λ = 1.54178 Å) radiation. The crystal structure was solved by direct methods using SHELXS-97 and final refinement, based on F2, was carried out by a full matrix least squares with SHELXL-97. Refinement was performed anisotropically for all nonhydrogen atoms. In general, the hydrogen atoms were assigned to idealized positions and allowed to ride on the parent atom.
3.4. Primary Anticancer Assay
All the cells were seeded in 96-well plates and incubated with 5% CO2 at 37 °C for 24 h. Next, the compounds and the reference were dissolved into the culture medium. The final concentration of DMSO in the medium was less than 0.5%. After the cells were treated with compounds for 48 h, the supernatant was removed and 5 mg/mL of a fresh prepared solution of MTT was added to each well and incubated with the cells at 37 °C for another 4 h. The medium was removed, and 100 μL of testing solution was added to each well. After an overnight incubation with 5% CO2 at 37 °C, the absorbance was measured at 490/630 nm by the microplate reader. The IC50 was calculated using GraphPad Prism version 6.0 software (San Diego, CA, USA) from the non-linear curve.
3.5. Kinase Activity Determination
The effects of the compounds on the activities of the tyrosine kinase was determined using enzyme-linked immunosorbent assay (ELISA). Briefly, 20 μg/mL substrate [poly (Glu,Tyr)4:1], 50 μL aliquot of 10 μmol/L ATP solution diluted in kinase buffer (50 mmol/L HEPES pH 7.5, 50 mmol/L MgCl2, 0.5 mmol/L MnCl2, 0.2 mmol/L Na3VO4, and 1 mmol/L DTT) was added to each well; 1 μL of various concentrations of target compound and positive control drug diluted in 1% DMSO (v/v) were then added to each well. 1% DMSO was used as the negative control. After 5–10 min preincubation, the kinase reaction was initiated by the addition of purified ALK proteins diluted in 49 μL of kinase buffer. After 30 min preincubation, the anti-phosphotyrosine monoclonalantibody (100 μL; 1:500, diluted in 5 mg/mL BSA T-PBS) was then added. After a 30 min incubation at 37 °C, and 100 μL horseradish peroxidase conjugated goat anti-mouse IgG (1:2000, diluted in 5 mg/mL BSA T-PBS) was added. The plate was then incubated at 37 °C for 30 min. A 100 μL aliquot of a solution containing 0.03% H2O2 and 2 mg/mL o-phenylenediamine in 0.1 mol/L citrate buffer was added. The reaction was terminated by the addition of 50 μL of 2 mol/L H2SO4 as the color changed, and the plate was analyzed using a multi-well spectrophotometer at 490 nm.
3.6. Molecular Modeling
The humanALK in complex with Crizotinib (PDB code: 2XP2) was retrieved from the Protein Data Bank (http://www.rcsb.org). Molecular docking was performed using the CDocker protocol (the Discovery Studio 3.5 software package, Accelrys, Co. Ltd., San Diego, CA, USA). Protein preparation was carried out using the Prepare Protein protocol, and all crystallographic water was removed from the protein. The ligands preparation was carried out using the Prepare Ligand protocol. The docking parameters were set as default. The lowest binding energy was taken as the best-docked conformation of the representative compound for the protein. Molecular docking was validated by the docking of the co-crystallized inhibitor for enzyme, and root-mean-square deviation (RMSD) value for the backbone atoms between docked pose and crystallographic pose was below 1.5 Å.
4. Conclusions
18β-GA is considered an interesting scaffold for the development of potential antitumor inhibitors. Current structural optimization of 18β-GA primarily focused on alternation in position C3 or C30-position, exhibited remarkable chemopreventive activities in various experimental cancer models. In the present study, 35 derivatives of 18β-GA, altered in C3 and/or C30-position, were developed and evaluated for their efficacy as antitumor inhibitors. Among the mentioned derivatives, the carboxyl group at the C30-position of 18β-GA is beneficial to improve the inhibitory potency. Compound 3j exhibited the most excellent antiproliferative activity against six humancancer cells (A549, HT29, HepG2, MCF-7, PC-3, and Karpas299). Besides, compound 3j inhibited the proliferation of HepG2 cell in a significant concentration manner and exhibited selective antiproliferative activity against humantumor cells. Furthermore, compound 3j exerted moderate inhibitory effects against the ALK. The results of docking analysis suggested that the hydrogen bond interactions with kinase are conducive to the binding. These results may inspire further structural modifications of 18β-GA aimed at developing potent ALK inhibitors.
Authors: Laura Conesa-Milián; Eva Falomir; Juan Murga; Miguel Carda; Eef Meyen; Sandra Liekens; J Alberto Marco Journal: Eur J Med Chem Date: 2018-02-09 Impact factor: 6.514