Omaima M. Abdelhafez*a,
Hamed I. Alibc,
Kamelia M. Amind,
Mohamed M. Abdallae and
Eman Y. Ahmeda
aChemistry of Natural Products Dept, National Research Center, Dokki, Egypt. E-mail: dromaima45@gmail.com; Fax: +20 202 33370931; Tel: +20 202 37608284
bPharmaceutical Sciences Dept, Irma Lerma Rangel College of Pharmacy, Texas A&M Health Science Center, TX, USA
cPharmaceutical Chemistry Dept, Faculty of Pharmacy, Helwan University, Egypt
dPharmaceutical Chemistry Dept, Faculty of Pharmacy, Cairo University, Egypt
eResearch Unit, Mapco Pharmaceutical industries Balteim, Egypt
First published on 3rd March 2015
In continuation of our work concerning the relation between the anticancer and anti-vascular endothelial growth factor receptor (anti-VEGFR-2) activity of some synthesized compounds, we hereby designed and prepared three new series of furochromone and benzofuran derivatives (carbonitriles, sulfonyl hydrazides and imides). The prepared compounds were evaluated for their in vitro VEGFR-2 inhibitory activity, their cytotoxicity on fifteen human cancer cell lines and their in vivo antiprostate cancer activity. The highest anti-VEGFR-2 activity was demonstrated by 6-acetyl-4-methoxy-7-methyl-5H-furo[3,2-g]chromen-5-one (3), which exhibited the same IC50 value as the reference drug sorafenib (2.00 × 10−3 μM). On the other hand, most of the synthesized compounds showed potent cytotoxicity against most of the tested cell lines, in particular, the carbonitrile series (4a,b and 5a,d) which exhibited the best activity with IC50 values ranging from 3.56 × 10−13 to 4.89 × 10−7 μM. Moreover, the imide series (15–17) showed the most significant in vivo antiprostate cancer activity. An in silico GOLD molecular docking study has been done to explore the binding mode of interaction of the furochromone and benzofuran derivatives to VEGFR-2 kinase, and to reveal the correlation between IC50 (μM) of the enzymatic inhibition of VEGFR-2 kinase and the GoldScore fitness for further therapeutic application.
Some inhibitors of VEGFR-2 such as sunitinib13 and sorafenib,14 have been approved by the Food and Drug Administration for treating different kinds of tumors, also in literature chromone bioisosteres,15,16 benzofuran,17,18 carbonitrile,19,20 hydrazide,21,22 and imide derivatives,23,24 were identified as anti-VEGFR-2 and anticancer agents.
The previously mentioned studies encouraged us to design some furochromone and benzofuran derivatives in order to evaluate their biological activities and molecular docking.
Refluxing visnagin (1) in aqueous potassium hydroxide afforded the visnaginone derivative (2),25 which reacted with acetic anhydride in the presence of anhydrous sodium acetate under reflux to give the acetyl visnagin derivative (3).26 The reaction of compound (3) with the respective aromatic aldehyde namely, (benzaldehyde, 4-bromobenzaldehyde, 5-methyl furfural and 3,4,5-trimethoxybenzaldehyde),in the presence of malononitrile and/or ethylcyanoacetate and ammonium acetate afforded the iminocarbonitrile (4a–d) and the oxocarbonitrile (5a–d) derivatives respectively.
Moreover, khellinone (7) and acetyl khellin (9) were prepared using the same previously mentioned method of preparation of visnaginone (2) and acetyl visnagin (3) respectively. The Mannich bases (11a–c)27 were prepared by the reaction of visnaginone (2) with some secondary amines namely, (piperidine, N-methyl piperazine and morpholine) in ethanol in the presence of formaline. The different aryl sulfonyl hydrazides (8a,b, 10a,b and 12a–c) were synthesized by reacting the previously mentioned acetyl containing compounds (2, 3, 7, 9 and 11a–c) with p-toluene sulfonyl hydrazide in methanol and few drops of conc. HCL.
Furthermore, chlorosulphonation of visnagin (1) at 0 °C affored visnagin sulfonyl chloride (13),28 which upon stirring at room temperature with hydrazine hydrate 98% in 20 mL absolute ethanol and few drops of triethylamine gave the visnagin sulfonyl hydrazide (14),28 finally, the imide derivatives (15–17) were obtained by refluxing (14) with the appropriate anhydrides namely, (phthalic, succinic and maleic anhydrides) in glacial acetic acid.
Compound | Enzymatic inhibition IC50 μM |
---|---|
VEGFR-2 | |
3 | 2.00 × 10−3 |
4a | — |
4b | 2.70 × 10−3 |
4c | 2.90 × 10−3 |
4d | 8.40 × 10−3 |
5a | 1.01 × 10−2 |
5b | 2.20 × 10−3 |
5c | 6.30 × 10−3 |
5d | 1.22 × 10−2 |
8a | 2.62 × 10−2 |
8b | 5.70 × 10−3 |
10a | 3.20 × 10−3 |
10b | 3.20 × 10−3 |
11a | 1.20 × 10−2 |
11b | 5.20 × 10−3 |
11c | 1.63 × 10−2 |
12a | 4.65 × 10−2 |
12b | 1.48 × 10−2 |
12c | 1.79 × 10−2 |
14 | 3.49 × 10−2 |
15 | 2.62 × 10−2 |
16 | 4.22 × 10−2 |
17 | 2.17 × 10−2 |
Sorafenib | 2.00 × 10−3 |
Compound | IC50 μM | ||||
---|---|---|---|---|---|
KB | SK OV-3 | SF-268 | NCI H 460 | RKOP27 | |
3 | — | 5.66 × 10−7 | 2.56 × 10−9 | — | — |
4a | — | 4.35 × 10−3 | 4.56 × 10−7 | 2.67 × 10−3 | 4.78 × 10−1 |
4b | 3.67 × 10−7 | — | 3.45 × 10−11 | 2.39 × 10−10 | 4.67 × 10−11 |
4c | 3.33 × 10−6 | 8.36 × 10−3 | 9.70 × 10−4 | 8.80 × 10−3 | 9.32 × 10−3 |
4d | 2.36 × 10−4 | 4.40 × 10−4 | 3.30 × 10−5 | 6.40 × 10−3 | 8.80 × 10−5 |
5a | 2.55 × 10−3 | — | 2.33 × 10−11 | 4.62 × 10−11 | 3.67 × 10−13 |
5b | 6.30 × 10−5 | 3.42 × 10−3 | 8.86 × 10−4 | 8.56 × 10−6 | 8.00 × 10−4 |
5c | 7.40 × 10−3 | 6.34 × 10−3 | 5.30 × 10−3 | 8.70 × 10−4 | 8.97 × 10−5 |
5d | 7.54 × 10−7 | 2.67 × 10−7 | — | 1.45 × 10−10 | 3.56 × 10−13 |
8a | 6.40 × 10−4 | 4.60 × 10−5 | 8.77 × 10−3 | 2.27 × 10−3 | 7.27 × 10−3 |
8b | 2.10 × 10−4 | 7.23 × 10−5 | 5.80 × 10−4 | 8.56 × 10−4 | 7.70 × 10−6 |
10a | 6.00 × 10−4 | 3.00 × 10−4 | 4.40 × 10−5 | 6.60 × 10−5 | 7.80 × 10−5 |
10b | 4.70 × 10−4 | 4.00 × 10−3 | 5.50 × 10−3 | 4.40 × 10−3 | 3.30 × 10−5 |
11a | 4.30 × 10−4 | 5.50 × 10−5 | 5.50 × 10−3 | 3.30 × 10−3 | 4.80 × 10−3 |
11b | 7.74 × 10−3 | 5.50 × 10−5 | 6.50 × 10−3 | 5.50 × 10−3 | 9.40 × 10−3 |
11c | 6.80 × 10−3 | 7.70 × 10−3 | 5.70 × 10−3 | 9.47 × 10−3 | 5.50 × 10−3 |
12a | 7.74 × 10−3 | 4.40 × 10−4 | 8.70 × 10−6 | 6.60 × 10−5 | 9.40 × 10−3 |
12b | 2.27 × 10−3 | 6.60 × 10−4 | 4.40 × 10−5 | 4.79 × 10−3 | 9.00 × 10−3 |
12c | 7.98 × 10−5 | 3.33 × 10−2 | 5.20 × 10−3 | 9.00 × 10−4 | 9.59 × 10−4 |
14 | 6.98 × 10−3 | 6.21 × 10−3 | 4.60 × 10−3 | 4.20 × 10−3 | 1.25 × 10−5 |
15 | 3.00 × 10−5 | 7.00 × 10−3 | 5.00 × 10−4 | 6.00 × 10−4 | 8.00 × 10−3 |
16 | 4.50 × 10−4 | 6.50 × 10−5 | 4.30 × 10−5 | 7.60 × 10−3 | 8.90 × 10−3 |
17 | 7.00 × 10−3 | 5.00 × 10−3 | 4.00 × 10−3 | 3.00 × 10−3 | 5.00 × 10−5 |
Fluorouracil | 4.46 × 10−3 | — | — | — | — |
Doxorubicin | — | 4.16 × 10−3 | — | — | — |
Cytarabine | — | — | 7.68 × 10−3 | — | — |
Gemcitabine | — | — | — | 2.13 × 10−3 | — |
Capecitabine | — | — | — | — | 4.33 × 10−3 |
Compound | IC50 μM | ||||
---|---|---|---|---|---|
Leukemia | Melanoma | ||||
HL60 | U937 | K561 | G361 | SK-MEL-28 | |
3 | 3.67 × 10−8 | — | — | 3.76 × 10−8 | — |
4a | 4.67 × 10−8 | 4.39 × 10−7 | — | — | 8.55 × 10−7 |
4b | 4.67 × 10−9 | 4.56 × 10−7 | 4.56 × 10−9 | 5.33 × 10−8 | — |
4c | 8.00 × 10−4 | 9.86 × 10−3 | 2.70 × 10−4 | 6.20 × 10−4 | 6.20 × 10−4 |
4d | 4.40 × 10−3 | 6.60 × 10−3 | 7.46 × 10−3 | 634 × 10−3 | 7.50 × 10−5 |
5a | 3.67 × 10−9 | 5.87 × 10−7 | 5.32 × 10−7 | 7.43 × 10−8 | 5.67 × 10−7 |
5b | 4.80 × 10−3 | 8.00 × 10−3 | 8.60 × 10−3 | 7.00 × 10−3 | 6.00 × 10−5 |
5c | 4.80 × 10−5 | 4.80 × 10−3 | 4.80 × 10−4 | 4.80 × 10−5 | 4.80 × 10−3 |
5d | 4.55 × 10−7 | 5.44 × 10−7 | 7.88 × 10−8 | 3.34 × 10−7 | 4.89 × 10−7 |
8a | 8.70 × 10−6 | 5.40 × 10−3 | 4.90 × 10−4 | 6.60 × 10−5 | 7.70 × 10−3 |
8b | 4.55 × 10−7 | 5.44 × 10−7 | 7.88 × 10−8 | 3.34 × 10−7 | 4.90 × 10−7 |
10a | 4.00 × 10−4 | 6.00 × 10−6 | 9.60 × 10−5 | 7.70 × 10−3 | 8.80 × 10−3 |
10b | 3.00 × 10−5 | 9.00 × 10−3 | 5.60 × 10−3 | 8.89 × 10−3 | 9.90 × 10−3 |
11a | 2.48 × 10−4 | 9.60 × 10−5 | 8.63 × 10−3 | 7.50 × 10−3 | 3.00 × 10−3 |
11b | 9.70 × 10−4 | 4.00 × 10−4 | 6.60 × 10−5 | 9.90 × 10−5 | 8.00 × 10−5 |
11c | 9.60 × 10−4 | 9.70 × 10−4 | 5.50 × 10−5 | 7.70 × 10−3 | 2.00 × 10−5 |
12a | 6.30 × 10−4 | 5.00 × 10−5 | 5.50 × 10−4 | 8.80 × 10−4 | 5.87 × 10−7 |
12b | 9.00 × 10−4 | 6.00 × 10−3 | 6.77 × 10−3 | 8.90 × 10−5 | 3.00 × 10−5 |
12c | 9.86 × 10−3 | 4.40 × 10−3 | 8.50 × 10−3 | 5.50 × 10−3 | 6.60 × 10−5 |
14 | 8.00 × 10−3 | 7.50 × 10−4 | 4.30 × 10−3 | 7.35 × 10−3 | 3.40 × 10−5 |
15 | 2.60 × 10−5 | 5.00 × 10−3 | 7.64 × 10−4 | 8.00 × 10−5 | 8.50 × 10−3 |
16 | 7.53 × 10−4 | 5.96 × 10−5 | 6.00 × 10−6 | 8.00 × 10−4 | 8.00 × 10−4 |
17 | 7.00 × 10−3 | 7.70 × 10−3 | 7.70 × 10−3 | 8.60 × 10−3 | 5.70 × 10−5 |
Doxorubicin | 1.13 × 10−3 | 4.45 × 10−3 | 6.66 × 10−3 | — | — |
Aldesleukin | — | — | — | 6.66 × 10−3 | 3.45 × 10−3 |
Compound | IC50 μM (neuroblastoma) | |
---|---|---|
GOTO | NB-1 | |
3 | 4.89 × 10−7 | 5.32 × 10−7 |
4a | — | 5.78 × 10−7 |
4b | 6.55 × 10−9 | 3.44 × 10−7 |
4c | 8.60 × 10−6 | 7.90 × 10−5 |
4d | 3.00 × 10−4 | 5.68 × 10−4 |
5a | — | — |
5b | 8.60 × 10−6 | 9.70 × 10−5 |
5c | 8.00 × 10−4 | 3.64 × 10−4 |
5d | 4.66 × 10−9 | — |
8a | 9.00 × 10−4 | 3.00 × 10−5 |
8b | 8.60 × 10−4 | 7.20 × 10−5 |
10a | 7.50 × 10−4 | 7.20 × 10−4 |
10b | 6.41 × 10−4 | 7.20 × 10−3 |
11a | 7.60 × 10−4 | 5.00 × 10−5 |
11b | 8.60 × 10−3 | 7.25 × 10−4 |
11c | 6.80 × 10−3 | 7.00 × 10−3 |
12a | 8.60 × 10−3 | 7.25 × 10−4 |
12b | 8.60 × 10−3 | 7.20 × 10−3 |
12c | 8.00 × 10−5 | 5.30 × 10−3 |
14 | 8.50 × 10−3 | 3.87 × 10−3 |
15 | 4.63 × 10−5 | 8.00 × 10−4 |
16 | 7.40 × 10−4 | 3.80 × 10−5 |
17 | 9.60 × 10−3 | 8.63 × 10−3 |
Doxorubicin | 4.73 × 10−3 | 5.15 × 10−3 |
Compound | IC50 μM | ||
---|---|---|---|
PC-3 | HT1080 | HepG2 | |
3 | 6.40 × 10−3 | 5.43 × 10−7 | 3.37 × 10−7 |
4a | 4.78 × 10−10 | 8.90 × 10−7 | 5.64 × 10−7 |
4b | 8.60 × 10−3 | 5.65 × 10−7 | 2.88 × 10−7 |
4c | 8.00 × 10−6 | 4.64 × 10−3 | 8.00 × 10−3 |
4d | 6.40 × 10−4 | 8.08 × 10−5 | 8.00 × 10−3 |
5a | 9.70 × 10−3 | 7.68 × 10−7 | 5.55 × 10−7 |
5b | 8.70 × 10−5 | 7.00 × 10−3 | 5.70 × 10−3 |
5c | 6.88 × 10−3 | 8.70 × 10−3 | 3.60 × 10−3 |
5d | 8.80 × 10−6 | 2.34 × 10−7 | 2.99 × 10−7 |
8a | 8.00 × 10−4 | 9.70 × 10−3 | 7.50 × 10−3 |
8b | 4.74 × 10−4 | 9.70 × 10−4 | 8.60 × 10−4 |
10a | 5.70 × 10−4 | 9.76 × 10−5 | 9.70 × 10−5 |
10b | 5.00 × 10−4 | 9.70 × 10−3 | 9.00 × 10−4 |
11a | 7.25 × 10−5 | 8.00 × 10−3 | 1.46 × 10−3 |
11b | 6.50 × 10−3 | 3.46 × 10−4 | 6.35 × 10−4 |
11c | 8.65 × 10−3 | 3.57 × 10−3 | 3.60 × 10−3 |
12a | 4.60 × 10−3 | 4.97 × 10−3 | 6.34 × 10−3 |
12b | 7.00 × 10−3 | 8.60 × 10−3 | 6.30 × 10−3 |
12c | 9.70 × 10−5 | 7.90 × 10−3 | 1.34 × 10−3 |
14 | 1.25 × 10−5 | 6.00 × 10−3 | 7.60 × 10−3 |
15 | 8.00 × 10−3 | 7.80 × 10−4 | 7.00 × 10−5 |
16 | 8.90 × 10−3 | 8.00 × 10−5 | 7.00 × 10−3 |
17 | 5.00 × 10−5 | 4.60 × 10−3 | 9.86 × 10−3 |
Abiraterone acetate | 6.70 × 10−6 | — | — |
Imatinib | — | 13.24 × 10−5 | — |
Gemcitabine | — | — | 3.44 × 10−3 |
Screening the cytotoxicity of the tested compounds on cervical carcinoma (KB), leukemia (K561), neuroblastoma (NB-1), liver carcinoma (HepG2) cell lines, showed that compound (4b) has a very remarkable activity, with IC50 values ranging from (4.56 × 10−9 to 3.67 × 10−7 μM), compared to the reference drug used.
On ovarian carcinoma (SKOV-3), colonoadenocarcinoma (RKOP27), melanoma (SK-MEL-28), neuroblastoma (GOTO), fibrosarcoma (HT1080) cell lines, compound (5d) was significantly potent with IC50 values ranging from (3.56 × 10−13 to 4.89 × 10−7 μM), exhibiting more potency than the reference drugs.
Regarding the CNS cancer (SF-268), non-small lung cancer (NCI H460), leukemia (HL60) cell lines, it was found that compound (5a) (IC50 ranging from 2.33 × 10−11 μM to 3.67 × 10−9 μM) demonstrated a potent cytotoxicity among the other synthesized derivatives.
Estimation of the cytotoxicity of the prepared compounds on leukemia (U937), prostate cancer (PC-3) cell lines revealed that compound (4a) showed a very pronounced potency (IC50 = 4.39 × 10−7 μM) and (IC50 = 4.78 × 10−8 μM) respectively compared to the reference drug and the other furochromone and benzofuran derivatives.
Evaluation of the cytotoxicity on melanoma (G361) cell lines showed that compound (3) exhibited a significant IC50 (3.67 × 10−8 μM), better than the reference drug aldesleukin (IC50 = 6.66 × 10−3 μM).
The previous screening of the bioactivity of the tested compounds on the presented panel of cell lines indicated that, the majority of the synthesized compounds were significantly potent than the comparative reference drugs used. The carbonitrile series (4a,b and 5a,d) showed the best potency against all the cell lines compared to the reference drugs and to the rest of the newly synthesized furochromone and benzofuran derivatives.
Compound | ED50 μM |
---|---|
a Each ED50 value is the mean (three significant digits) ± SEM from five experiments. Level of statistical significance: P < 0.01 with respect to ED50 value of flutamide as determined by ANOVA/Dunnett's. | |
3 | 13.11 ± 0.004 |
4a | 21.49 ± 0.006 |
4b | 23.64 ± 0.005 |
4c | 10.64 ± 0.003 |
4d | 8.00 ± 0.004 |
5a | 26 ± 0.004 |
5b | 13.11 ± 0.004 |
5c | 7.21 ± 0.005 |
5d | 28.6 ± 0.004 |
8a | 5.66 ± 0.04 |
8b | 6.22 ± 0.04 |
10a | 6.85 ± 0.05 |
10b | 7.53 ± 0.06 |
11a | 3.40 ± 0.006 |
11b | 16.14 ± 0.006 |
11c | 19.53 ± 0.006 |
12a | 14.68 ± 0.005 |
12b | 17.76 ± 0.007 |
12c | 5.14 ± 0.005 |
14 | 3.73 ± 0.004 |
15 | 2.56 ± 0.005 |
16 | 2.82 ± 0.006 |
17 | 2.12 ± 0.007 |
Flutamide | 11.60 ± 0.09 |
In this article, we have implemented the GoldScore as a search algorithmic function for “drug-like” which was reported to give superior results than the Chemscore function as a scoring function for GOLD.31 This GoldScore fitness means the higher the values of fitness (Table 4) is the better binding interaction into the binding site.
GOLD fitness = Shb_ext + Svdw_ext + Shb_int + Svdw_int |
Comp. | GoldScore | External VDW | Hydrogen bonds between atoms of compounds and amino acids of VEGFR | RMSDa (Å) | |||
---|---|---|---|---|---|---|---|
Atoms of comp. | Amino acids | Distance (Å) | Angle (°) | ||||
a Root mean square deviation.b N-[4-({3-[2-(Methylamino)pyrimidin-4-yl]pyridin-2-yl}oxy)naphthalen-1-yl]-6-(trifluoromethyl)-1H-benzimidazol-2-amine. | |||||||
3 | 44.02 | 25.40 | 4-CH3O | HN of K868 | 1.97 | 166.6 | 6.96 |
6-C![]() |
HN of R1051 | 1.95 | 164.8 | ||||
4a | 72.00 | 52.43 | 5-C![]() |
HN of D1046 | 2.32 | 164.3 | 4.44 |
4b | 74.00 | 53.55 | 4′′-Br | HN of R1027 | 2.45 | 152.0 | 4.16 |
4c | 69.03 | 49.37 | 5-C![]() |
HN of D1046 | 2.07 | 164.8 | 4.73 |
4d | 79.85 | 55.78 | 5-C![]() |
HN of D1046 | 2.01 | 159.7 | 5.65 |
4′′-CH3O | HN of R1027 | 2.11 | 152.9 | ||||
5a | 71.80 | 52.78 | 5-C![]() |
HN of D1046 | 2.05 | 165.7 | 4.68 |
5b | 70.86 | 51.09 | 5-C![]() |
HN of D1046 | 2.05 | 165.4 | 4.63 |
5c | 70.66 | 53.09 | 5-C![]() |
HN of D1046 | 2.23 | 171.3 | 4.77 |
5d | 81.17 | 58.32 | 5-C![]() |
HN of D1046 | 2.48 | 162.0 | 5.48 |
8a | 61.23 | 35.49 | 6-O | HN of K868 | 1.85 | 163.8 | 3.34 |
8b | 63.93 | 37.98 | 4-O | HN of K868 | 1.48 | 167.4 | 4.70 |
6-O | HN of D1046 | 2.15 | 170.1 | ||||
10a | 68.08 | 45.62 | 4-O | 1HN of K868 | 2.25 | 174.2 | 4.55 |
6-O | 2HN of K868 | 2.09 | 166.1 | ||||
10b | 71.36 | 41.51 | 4-O | 1HN of K868 | 2.24 | 172.9 | 4.84 |
6-O | 2HN of K868 | 2.12 | 159.8 | ||||
11a | 57.89 | 39.98 | Furan O | HN of K868 | 2.37 | 171.5 | 3.86 |
6-O | HN of D1046 | 1.85 | 171.3 | ||||
11b | 59.91 | 36.94 | Furan O | HN of K868 | 2.21 | 171.9 | 4.00 |
6-O | HN of D1046 | 1.87 | 144.2 | ||||
11c | 58.27 | 38.14 | Furan O | HN of K868 | 2.27 | 174.0 | 3.42 |
6-O | HN of D1046 | 2.05 | 160.1 | ||||
12a | 74.71 | 54.22 | 4-CH3O | HN of K868 | 2.22 | 165.7 | 2.39 |
6-O | HN of D1046 | 2.16 | 151.4 | ||||
12b | 75.09 | 51.49 | 6-OH | O![]() |
1.94 | 141.6 | 6.27 |
5-CHNHNH | O![]() |
2.25 | 115.5 | ||||
12c | 78.28 | 48.43 | Furan O | HN of D1046 | 2.20 | 161.9 | 5.85 |
S![]() |
HN of K868 | 1.75 | 162.0 | ||||
14 | 57.17 | 27.14 | Furan O | HN of K868 | 1.64 | 154.2 | 2.40 |
5-C![]() |
HO of T916 | 1.92 | 139.7 | ||||
S![]() |
HN of K868 | 1.85 | 164.0 | ||||
S![]() |
O![]() |
1.48 | 127.6 | ||||
15 | 72.88 | 51.88 | Furan O | HN of K868 | 1.87 | 164.4 | 2.40 |
16 | 60.02 | 37.25 | Furan O | HN of K868 | 2.30 | 174.8 | 2.87 |
S![]() |
HN of D1046 | 2.09 | 152.5 | ||||
17 | 60.65 | 38.98 | 4-O | 1HN of K868 | 1.85 | 151.4 | 4.28 |
5-C![]() |
2HN of K868 | 1.54 | 150.2 | ||||
Pyran O | HN of D1046 | 2.15 | 154.9 | ||||
S![]() |
HN of D1046 | 2.48 | 114.2 | ||||
K111b | 80.74 | 60.92 | Pyrimidine-N1 | 1HN of K868 | 2.46 | 139.4 | 1.83 |
Naphthyloxy-O | 2HN of K868 | 2.34 | 114.7 | ||||
Naphthyl-NH | OH of T916 | 1.90 | 154.9 |
The fitness score is the negative sum of the component energy terms, such as protein–ligand external hydrogen bond energy (Shb_ext), protein–ligand external van der Waals energy (Svdw_ext), protein–ligand intramolecular hydrogen bond energy (Shb_int), and the ligand internal van der Waals energy (Svdw_int), so that larger fitness scores are better.
The implemented GoldScore algorithmic function was initially evaluated for its accuracy into VEGFR-2 kinase. As a function of RMSD (root mean square deviation) of the runs done by flexible docking using GOLD 5.1, the docking poses were nearest to the experimental binding mode of the co-crystallized ligands of VEGFR-2.32 Our docking results are matched with the reported requirements for successful scoring function,33 being of RMSD ≤ 2.0 Å from the experimental one.
In this study, the newly synthesized compounds were docked into the target crystal structure of human Protein tyrosine kinase of VEGFR-2 kinase domain (pdb code: 3EWH) in complex with a pyridyl-pyrimidine benzimidazole inhibitor (K111). Many compounds revealed good Gold fitness score (GoldScore) values namely, compounds (4a,b,d, 5a–d, 10b, 12a–c and 15) which revealed 72, 74, 79.85, 71.80, 70.86, 70.66, 81.17, 71.36, 74.71, 75.09, 78.28 and 72.88 respectively (Table 4). Almost all the docked compounds show the common hydrogen bond interactions (1–4) with D1046 (NH, CO), K868 (NH), E885 (C
O), and R1027 (NH). Fitness of the best docked compounds takes place within RMSD of 2.39–6.27 Å as cited in (Table 4). Such binding mode was identically similar to the native co-crystallized ligand (K111) which binds with similar binding orientation involving similar amino acids used by docked compounds. This reveals that the interactions between the K868 (NH) and the small molecules will be crucial to inhibit the VEGFR-2 kinase activity.
Compound (4a) possess a high potential fitness (GoldScore: 72) into the binding site of the 3D macromolecule (PDB: 3EWH). Its high affinity is presumably attributed to its hydrogen bond formed between its 5-CO group and D1046 (NH) amino acid and RMSD of 4.44 Å. In addition to the hydrophobic interaction of its 4′-phenyl pyridine moiety which occupies the hydrophobic pockets of the receptor site (I888, I889, I892, V898, V899, and I1044), as illustrated in Fig. 1. Compound (4b) has a better potential fitness (GoldScore: 74) into the binding site, forming one hydrogen bond between its terminal 4′′-Br and R1027 (NH) amino acid. In addition, both compounds (4a and 4b) revealed external vdw of 52.43 and 53.55, respectively.
As shown in (Fig. 2) and cited in (Table 4), compound (12a) (GoldScore: 74.71) has a promising antiproliferative activity through its higher kinase binding affinity. Where it binds into VEGFR2 kinase through two hydrogen bonds between its 4-CH3O and 6-O moieties and K868 (HN), and D1046 (HN) amino acids, respectively, where it revealed excellent RMSD of 2.39 Å. Furthermore its 7-piperidinyl moiety occupies the extended hydrophobic pocket similarly to the lipophilic trifluoromethyl-arene portion of the native cocrystallized (K111).
In addition to the role of the hydrogen bond and vdw interactions, the hydrophobic interaction plays a significant role to occupy the 3(4′-aryl pyridine) moiety of compounds (4a–d and 5a–d), and the 7-piperidinyl, N–Me piprazinyl, and morpholino moieties of compounds (11a–c and 12a–c) into the extended hydrophobic pocket similarly to the lipophilic trifluoromethyl-arene portion of the native co-crystallized (K111).
In the analysis of GOLD docking results, a high correlation (R2 = 0.712) between IC50 of (VEGFR-2) enzymatic inhibition and the GoldScore fitness for compounds (4b,c, 5a–c, 8a, 10a,b, 11c, 12a, 16 and 17) into (VEGFR-2) kinase is shown in (Fig. 3A).
Also well correlated result (R2 = 0.687) is shown between IC50 (μM) against prostate cancer cell line (PC-3) and GoldScore fitness of compounds (4a,d, 5d, 10a,b, 11b,c, 12a,c and 16) (Fig. 3B).
On other hand the correlations between the biological results IC50 (μM) against different tumor cell lines namely, (G361, HL60, K561, KB, RKOP27, SKOV-3, SK-MEL-28 and U937) were variably correlated and revealing correlation coefficients (R2) of 0.67, 0.74, 0.73, 0.764, 0.66, 0.658, 0.79 and 0.686, respectively.
On the other hand, the in vivo results of these compounds are correlated with this outcome, where compound (5b) revealed better in vivo antiprostate cancer activity (ED50 = 13.11 μM) than that of compounds (5a and 5d) (ED50 = 26.00 and 28.60 μM, respectively). This suggests that their prostate cancer antiproliferative activities are mediated through VEGFR-2 pathway.
Likewise within the in vivo prostate cancer experiment, compounds (15–17) have (ED50 of 3.73, 2.56, 2.82, and 2.12 μM, respectively). Similarly, they have equipotent inhibitory activities against VEGFR-2 of (3.49 x 10−2, 2.62 × 10−2, 4.22 × 10−2, and 2.17 × 10−2 μM, respectively). However, this correlation was deviated in the in vitro PC3 cell assay, where compound (14 and 17) have equipotent activities of IC50: 1.25 × 10−5 and 5.00 × 10−5 μM, respectively, revealing highly more potent cytotoxic effect than compounds (15 and 16) of IC50: 8.0 × 10−3 and 8.9 × 10−3 μM, respectively. This result may suggest the resistance of PC3 cell lines against compounds (15 and 16), and indicates that the in vitro data may not necessarily reflect the in vivo effect.
The reaction was initiated by adding VEGFR-2 kinase (5 μL), followed by addition of 100 μL of anti-phosphotyrosine antibody (PY99; 1:
500 dilution) then goat anti-mouse IgG horseradish peroxidase (100 μL; 1
:
2000 dilution) diluted in T-PBS containing 5 mg mL−1 BSA was added. After each addition there was an incubation period of 1 h at room temperature followed by washing.
Finally, 100 μL of (0.03% H2O2, 2 mg mL−1 o-phenylenediamine in citrate buffer 0.1 M, pH 5.5) was added and incubated at room temperature until color developed. The reaction was terminated by the addition of 100 μL of 2 M H2SO4. The inhibition rate (%) was calculated using the following equation:
Inhibition rate (%) = [1 − (A492/A492Control)] × 100%, where, A492 was measured using a multiwell spectrophotometer (VERSAmax™).
After being incubated for 24 h, cells were treated with different concentrations of each compound under test, where an aliquot of frozen concentrate was diluted to twice the desired final maximum test concentration and four 10 fold serial dilutions were made to provide a total of five drug concentrations plus control. The cells were incubated with the compounds for 48 h at 37 °C and in atmosphere of 5% CO2. Control cells were treated with vehicle alone.
After 48 h, cells were fixed, washed and stained with sulforhodamine B (SRB) (0.4% w/v) dissolved in 1% acetic acid for 10 min. Excess stain was washed with acetic acid and attached stain was solubilized with 10 mM trizma base. The absorbance was measured in an automated plate reader at wave length of 515 nm. The concentration required for 50% inhibition of cell viability (IC50) was calculated and compared with the reference drugs.
For pose selection and enrichment studies default values of speed settings and all other parameters were used. GOLD Score fitness algorithmic function was used with a default input and annealing parameters, which were managed within 4.0 and 2.5 Å for the van der Waals and H-bond interactions, respectively38,39
All the docked ligands were energetically minimized by using MOPAC with 100 iterations and minimum RMS gradient of 0.10.
In this study, 10 genetic algorithm (GA) docking runs were used with internal energy offset. For pose reproduction analysis, the radius of the binding pocket was set as the maximal atomic distance from the geometrical center of the ligand plus 3 Å. When the top three solutions attained RMSD values within 1.5 Å, docking was terminated.
The top ranked docking pose was retained for the 3D cumulative success rate analysis. Rescoring was conducted with the GOLD rescore option, in which poses would be optimized by the program. The top 10 ranked docking poses were selected based on their Gold fitness scores and their reported mode of interactions, which were analyzed using Accelrys Discovery studio to reveal the hydrogen bond interaction and binding mode within the binding domain.
With respect to ligand flexibility, flipping of all planar RNR1R2 ring, NH–R ring, carboxylic acids –COOH were included, and the torsion angle distribution and post process rotatable bonds were set as default. All other Genetic Algorithm parameters were set as default settings.
The crystal structure of the (VEGFR-2) kinase domain in complex with N-[4-({3-[2-(methylamino)pyrimidin-4-yl]pyridin-2-yl}oxy)naphthen-1-yl]-6-(trifluorometh-yl)-1H-benzimidazol-2-amine (K111) co-crystalized ligand (pdb code: 3EWH) was used as target protein structure, which was retrieved from the Protein Data Bank (PDB). PDB sum was accessed for identification of the key amino acids and the flexible residues (E885, T916, and D1046) in the binding site. All residues within 10 Å radius from the origin: 13.237, −2.016 and 11.23 of the co-crystallized ligand coordinates.
IC50 | Half maximal inhibitory concentration |
ED50 | Median effective dose |
TLC | Thin layer chromatography |
DMSO | Dimethyl sulfoxide |
MOPAC | Molecular orbital package |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4ra16228e |
This journal is © The Royal Society of Chemistry 2015 |