Mona M Abdel-Atty1, Nahla A Farag1, Rabah A T Serya2, Khaled A M Abouzid2,3, Samar Mowafy1,4. 1. Pharmaceutical Chemistry Department, Faculty of Pharmacy, Misr International University, Cairo, Egypt. 2. Pharmaceutical Chemistry Department, Faculty of Pharmacy, Ain Shams University, Cairo, Egypt. 3. Organic and Medicinal Chemistry Department, Faculty of Pharmacy, University of Sadat City, Sadat City, Egypt. 4. Department of Chemistry, University of Washington, Seattle, WA, USA.
Abstract
A series of thieno[2,3-d]pyrimidine-based hydroxamic acid hybrids was designed and synthesised as multitarget anti-cancer agents, through incorporating the pharmacophore of EGFR, VEGFR2 into the inhibitory functionality of HDAC6. Three compounds (12c, 15b and 20b) were promising hits, whereas (12c) exhibited potent VEGFR2 inhibition (IC50=185 nM), potent EGFR inhibition (IC50=1.14 µM), and mild HDAC6 inhibition (23% inhibition). Moreover, compound (15c) was the most potent dual inhibitor among all the synthesised compounds, as it exhibited potent EGFR and VEGFR2 inhibition (IC50=19 nM) and (IC50=5.58 µM), respectively. While compounds (20d) and (7c) displayed nanomolar selective kinase inhibition with EGFR IC50= 68 nM and VEGFR2 IC50= 191 nM, respectively. All of the synthesised compounds were screened in vitro for their cytotoxic effect on 60 human NCI tumour cell lines. Additionally, molecular docking studies and ADMET studies were carried out to gain further insight into their binding mode and predict the pharmacokinetic properties of all the synthesised inhibitors.
A series of thieno[2,3-d]pyrimidine-based hydroxamic acid hybrids was designed and synthesised as multitarget anti-cancer agents, through incorporating the pharmacophore of EGFR, VEGFR2 into the inhibitory functionality of HDAC6. Three compounds (12c, 15b and 20b) were promising hits, whereas (12c) exhibited potent VEGFR2 inhibition (IC50=185 nM), potent EGFR inhibition (IC50=1.14 µM), and mild HDAC6 inhibition (23% inhibition). Moreover, compound (15c) was the most potent dual inhibitor among all the synthesised compounds, as it exhibited potent EGFR and VEGFR2 inhibition (IC50=19 nM) and (IC50=5.58 µM), respectively. While compounds (20d) and (7c) displayed nanomolar selective kinase inhibition with EGFR IC50= 68 nM and VEGFR2 IC50= 191 nM, respectively. All of the synthesised compounds were screened in vitro for their cytotoxic effect on 60 humanNCI tumour cell lines. Additionally, molecular docking studies and ADMET studies were carried out to gain further insight into their binding mode and predict the pharmacokinetic properties of all the synthesised inhibitors.
The development of multitarget drug therapy has become an important strategy for cancer treatment. Recently, the modulation of receptor tyrosine kinase (RTK) pathways by the inhibition of histone deacetylases (HDACs) has become a promising approach for cancer therapy. HDACs are a family of numerous epigenetic enzymes that are important therapeutic targets for cancer. HDACs have been classified into four distinct classes, in which class I (HDAC1–3 and 8), class II (HDAC4–7, 9 and 10), and class IV (HDAC11). HDAC6 is considered a therapeutically important target for cancer treatment due to its interaction with proteins involved in cell growth, migration, protein degradation, and apoptosis. It is reported that HDAC inhibitors with large or rigid hydrophobic surface recognition moiety (SRM) and bulky aromatic or short linkers, are more efficient to achieve HDAC6 selectivity, which could achieve a closer approach of the hydroxamate group to Zn2+. Numerous HDAC6 inhibitors have been developed having large SRM connected to the zinc-binding group (ZBG) with diverse linkers (Figure 1).
Figure 1.
The structures of histone deacetylase 6 (HDAC6) inhibitors, colour codes show essential pharmacophore for HDAC inhibitors composed of zinc-binding group (blue), linker (red), and SRM (purple).
The structures of histone deacetylase 6 (HDAC6) inhibitors, colour codes show essential pharmacophore for HDAC inhibitors composed of zinc-binding group (blue), linker (red), and SRM (purple).Despite the efficiency of HDAC inhibition as monotherapy in haematologic malignancies, it did not show significant competency against solid tumours. As a result, multitarget therapy is developed as a challenging approach for increasing the efficacy of HDAC inhibitors. The combination of HDAC and TKIs is considered a promising approach to enhance therapeutic effect and repress resistance. For example, a synergistic effect is observed by a combination of gefitinib (EGFR inhibitor) with vorinostat (HDAC inhibitor) by enhancing tumour growth inhibition and apoptosis compared with monotherapy of selective EGFR inhibitors,. Another example, combination therapy of pazopanib (VEGFR2 inhibitor) and abexinostat (HDAC inhibitor) showed that HDAC inhibition could improve response and holdback resistance to pazopanib in patients with solid tumour malignancies and renal cell carcinoma. Thus, it is suggested that a single molecule that simultaneously inhibits HDAC and RTK activities targeting multiple biological molecules and multiple signal pathways can not only enhance drug efficacy but also can have additive or even synergistic antitumor effects and can overcome the problems of cancer resistance and relapse. Additionally, it can overcome unfavourable pharmacokinetic properties, drug–drug interactions, poor patient compliance, and high drug cost involved in multicomponent drug cocktails,. CUDC-101 is a first-in-class multi-targeted hybrid that was designed to simultaneously inhibit HDACs, EGFR, and HER2. It was found that it induces significantly synergistic tumour regression in breast, NSCLC, liver, colon, head, neck, and pancreatic cancers. Series of hybrids of 4-anilinoquinaziline based hydroxamic acid derivatives were synthesised and discovered as dual VEGFR-2/HDAC inhibitors and EGFR/HDAC inhibitors (Figure 2),.
Figure 2.
The structures of several reported dual inhibitors HDAC/EGFR or VEGFR2.
The structures of several reported dual inhibitors HDAC/EGFR or VEGFR2.Moreover, the development of dual-inhibitors targeting EGFR and VEGFR2 has recently become a promising strategy for the treatment of diverse types of cancer to overcome problematic drug resistance and low response to selective tyrosine kinase inhibitors (TKIs) and to produce a synergistic effect. Inhibition of both EGFR and VEGFR2 produces a more significant effect on angiogenesis than selective VEGFR2 inhibitors, where EGFR-induced the production of angiogenic growth factors, which indirectly affects angiogenesis,. Several EGFR/VEGFR-2 dual inhibitors have been designed and synthesised (Figure 3),.
Figure 3.
Examples of dual inhibitors acting on EGFR and VEGFR2.
Examples of dual inhibitors acting on EGFR and VEGFR2.In this study, six novel series of thieno[2,3-d]pyrimidine-based hydroxamic acid derivatives were rationally designed. The design is based on the strategy of multitarget drug therapy and the clinical advantage of combining HDAC, VEGFR2, or/and EGFR inhibition in one single drug. Investigating the structural requirements of EGFR, VEGFR2, and HDAC inhibitors led us to design several series of novel leads of multitarget or dual-target inhibitors by combining the pharmacophore of EGFR, VEGFR2 inhibitors with that of HDAC inhibitors.
Rationale of the design
Design strategy is based on the study of several reported dual inhibitors such as CUDC101 and compound B, which were reported as dual EGFR/HDAC and VEGFR2/HDAC inhibitors respectively; (bearing 4-anilinoquinaziline fragment) (Figure 2),. Additionally, novel dual EGFR/VEGFR2 inhibitors such as compound D and compound E; (bearing thieno[2,3-d] pyrimidine fragment;) were recently reported (Figure 3).Based on these reported dual inhibitors and by investigating the structural requirements for EGFR, VEGFR2, and HDAC inhibition, we designed six series of hybrid molecules (7, 11, 12, 15, 19 and 20) with potent multi acting EGFR/VEGFR/HDAC inhibition in an attempt to synergise inhibitory activity, decrease resistance and relapse caused by single-target therapy and also avoid multicomponent drug cocktails therapy.The design aimed at replacing the quinazoline fragment in CUDC101 and compound B with thieno[2,3-d]pyrimidine fragment as in compound D and compound E, representing SRM of HDAC inhibitors and also showing suitable interactions in the hinge region of ATP-binding site of EGFR and VEGFR2 (Figure 4). This bioisostere replacement is also supported by recent reports of novel HDAC inhibitors bearing thieno[2,3-d]pyrimidine fragments representing the SRM that displayed excellent inhibitory activities on HDACs, such as compound A as shown in (Figure 1). Addition of aniline, 5-aminobenzimidazole, (4-phenylthioureido)phenylamine) and (4-phenylureido)phenylamine) fragments to 4-thieno[2,3-d]pyrimidine ring will extend the structure to the hydrophobic pocket of EGFR and VEGFR2 (Figure 4). Next, a hydroxamic acid functionality was directly attached to the thieno[2,3-d]pyrimidine fragment (SRM) or connected with a short linker through a connecting amide unit as the previously reported HDAC6 inhibitors (Figure 1),, in which it served as a bridge to properly display HDAC6 inhibitory pharmacophore features and the solvent-accessible area of EGFR- and VEGFR2-binding sites (Figure 4).
Results and discussion
Chemistry
The route adopted for the synthesis of key intermediates (3) and (5a-d) were prepared according to literature as illustrated in (Schemes 1 and 2), respectively. The synthetic route of the first three series of the target novel hydroxamic acid derivatives with different urea and thiourea moieties (7a-d, 11a-c, 12a-c) is depicted in (Scheme 3). Moreover, the synthesis of another three series of the novel target hydroxamic acid derivatives with different aniline moieties (15a-d, 19 b-d, 20 b-d) is depicted in (Scheme 4).
Synthesis of. Ethyl (4-chloro-5-methylthieno[2,3-d]pyrimidine)-6-carboxylate. Reagents and conditions: (a) Morpholine, S, EtOH, 70 °C, 4 h, 0 °C, 24 h; (b) Formamide, acetic acid, 40 h., 150 °C; (c) POCl3, reflux, 3 h.Synthesis of 1–(4-aminophenyl)-3-phenylurea and 1–(4-aminophenyl)-3-phenylthiourea derivatives. Reagents and conditions: (a) DCM, rt, 24 h. (b) H2, Pd/C, MeOH, rt, 4 h.Synthesis of phenyl ureido phenyl amino thieno [2,3-d] pyrimidine hydroxamic acid derivatives. Reagents and conditions: (a) EtOH, TEA, reflux, 18–48 h; (b) NH2OH.HCl, MeOH,/Na, 0 °C; (c) Li(OH)2, THF, reflux, 7 h; (d) EDCI/HOBt, NMM, DMF, 0 °C, rt, 24 h; (e) NH2OH.HCl, EtOH/Na, 0 °C.Synthesis of 5-amino benzimidazole thieno [2,3-d] pyrimidine hydroxamic acid derivative and 4-anilino thieno [2,3-d] pyrimidine hydroxamic acid derivatives. Reagents and conditions: (a) H2, Pd/C, MeOH, rt, 6 h; (b) EtOH, TEA, reflux, 18–48 h; (c) NH2OH.HCl, MeOH,/Na, 0 °C; (d) Li(OH)2, THF, reflux, 7 h; (e) EDCI/HOBt, NMM, DMF, 0 °C, rt, 24 h; (f) NH2OH.HCl, EtOH/Na, 0 °C.In (Scheme 1), intermediate (1) was synthesised by applying the Gewald reaction, in which ethyl cyanoacetate, ethyl acetoacetate, and sulphur were dissolved in absolute ethanol in the presence of morpholine. The cyclisation of (1) into (2) was achieved through refluxing (1) in a mixture of formamide and acetic acid for 40 h. The structure was confirmed by its reported melting point. Compound (3) was synthesised through chlorination of (2) by refluxing it with phosphorus oxychloride for 3 h, (Scheme 1).
Scheme 1.
Synthesis of. Ethyl (4-chloro-5-methylthieno[2,3-d]pyrimidine)-6-carboxylate. Reagents and conditions: (a) Morpholine, S, EtOH, 70 °C, 4 h, 0 °C, 24 h; (b) Formamide, acetic acid, 40 h., 150 °C; (c) POCl3, reflux, 3 h.
The synthesis of the intermediates (4a–d) was carried out by reacting the appropriate isocyanate with p-nitroaniline in dry dichloromethane for 48 h at room temperature to yield compounds (4a–d), which are then reduced to their corresponding amino derivatives (5a-d) using 10% Pd-C in methanol, (Scheme 2).
Scheme 2.
Synthesis of 1–(4-aminophenyl)-3-phenylurea and 1–(4-aminophenyl)-3-phenylthiourea derivatives. Reagents and conditions: (a) DCM, rt, 24 h. (b) H2, Pd/C, MeOH, rt, 4 h.
In (Scheme 3), the synthesised key intermediates (5a–d) were then coupled with the previously prepared chloro derivative (3) to yield (6a-d). The first novel series of the target hydroxamic acid derivatives (7a-d) was obtained through the reaction of the ethyl ester group in compounds (6a-d) with freshly prepared hydroxylamine in absolute ethanol to give (7a-d),, (Scheme 3). Compounds (7a-d) were confirmed by many spectral data, in which 1HNMR signals were consistent with protons of the targeted compounds (7a-d).
1HNMR signals revealed equally integrated signals between δ 8.23–9.82 ppm representing D2O exchangeable protons of the NH linker, the urea protons, and the NH of hydroxamic acid moiety in all target compounds (7a-d). Noticeably, 1HNMR signals revealed the absence of the triplet and the quartette peak around δ 1.31 − 1.36 ppm and δ 4.30–4.38 ppm, which confirms the absence of the ester group and the entrance of hydroxamic acid moiety.Moreover, the synthesis of another two novel series of hydroxamic acid derivatives with urea and thiourea moieties (11a-c and 12a-c) was depicted in (Scheme 3), through multistep reactions. The synthesis starts with the hydrolysis of the carboxylate ethyl ester of the previously prepared compounds (6a-c) (Scheme 3), through refluxing it with Li(OH)2 in THF and water to prepare the intermediate carboxylic acid derivatives (8a-c). Furthermore, the synthesis of novel two series of intermediate ethyl ester derivatives (9a-c, 10a-c) (Scheme 3) was achieved through coupling of the previously prepared carboxylic acid derivatives (8a-c) with the appropriate amino acid esterhydrochloride salt, utilising EDCI with HOBt in dry DMF. This method offered a high yield of novel intermediate compounds (9a-c and 10a-c) confirmed by many spectral data. 1HNMR signals confirmed the prepared structures of the novel esters, in which the aliphatic protons were shown as quartette signals for synthesised novel esters, where it appeared around δ 4.05–4.25 ppm representing the two protons of OCH2-CH3 with J = 8 Hz. Another quartette signal around δ 3.26–3.78 ppm appeared in all compounds referring to two protons of CH2-NH with J = 8 Hz. Triplet signal appeared around δ 2.36–2.71 ppm for all compounds, representing the two protons of CH2-CO with J = 8 Hz. Moreover, another triplet signal appeared for all compounds around δ 1.17–1.34 ppm, representing three protons of OCH2-CH3 with J = 8 Hz. An additional multiplet signal appeared around δ 1.76–2.04 ppm, representing two protons of CH2-CH2-CH2 in all prepared ethyl ester derivatives of series (10a-c) (Scheme 3).Finally, the ethyl ester group in the previously prepared novel compounds (9a-c, 10a-c) reacted with freshly prepared hydroxylamine in absolute ethanol to obtain target novel hydroxamic acid derivatives (11a-c, 12a-c) (Scheme 3). The two target novel series (11a-c, 12a-c) were confirmed by many spectral data, where 1HNMR spectrum revealed the disappearance of the quartette signal of the two protons of OCH2-CH3 around δ 4.05–4.25 and the triplet signal of the three protons OCH2-CH3 around δ 1.17–1.34 ppm, in which it confirms the absence of the ester group and the entrance of the hydroxamic acid moiety.In (Scheme 4), the purchased compound (5-nitrobenzimidazole) was reduced to afford the corresponding amino derivative (13a) using 10% Pd-C in methanol. Prepared compound (13a) and purchased aniline derivatives (13b-d) were also coupled with the previously prepared chloro derivative (3) to give the corresponding thieno[2,3-d]pyrimidine derivatives bearing aniline moieties (14a-d) (Scheme 4). Finally, the route adopted for the synthesis of the three target novel series of hydroxamic acid derivatives bearing aniline moiety (15a-d, 19 b-d and 20 b-d) (Scheme 4) is the same route depicted in (Scheme 3) with the same reagents and conditions. The first novel series of the target hydroxamic acid derivatives (15a-d) was obtained through the reaction of the ethyl ester group in compounds (14a-d) with freshly prepared hydroxylamine in absolute ethanol to give (15a-d),, (Scheme 4). Then, the hydrolysis of the carboxylate ethyl ester of the previously prepared compounds (14b-d) (Scheme 4), through refluxing it with Li(OH)2 in THF and water to prepare the intermediate carboxylic acid derivatives (16b-d). The previously prepared carboxylic acid derivatives (16b-d) are coupled with the appropriate amino acid esterhydrochloride salt, utilising EDCI with HOBt in dry DMF, to obtain a high yield of novel intermediate compounds (17b-d and 18b-d), (Scheme 4). Finally, the ethyl ester group in the previously prepared novel compounds (17b-d, 18b-d) reacted with freshly prepared hydroxylamine in absolute ethanol to obtain target novel hydroxamic acid derivatives (19b-d, 20b-d) confirmed by many spectral data, (Scheme 4).
Biological evaluation
The target compounds were evaluated for their in vitro inhibitory activities against EGFR, VEGFR-2, and HDAC6. The EGFR and VEGFR-2 inhibitory assays were performed at Thermo Fischer Scientific, Madison, WI, USA (SelectScreenServices@thermofisher.com). The HDAC inhibitory assays were performed at BPS Bioscience, San Diego, CA, USA (www.bpsbioscience.com). Initially, all of the synthesised target compounds were evaluated for their percent inhibitory activity against EGFR, VEGFR-2, and HDACs at 10 µM single dose concentration, in which their mean percent inhibition is summarised in (Figure 5, Tables 1 and 2). Representative compounds were evaluated for their percent inhibitory activity against EGFR and VEGFR-2 at 6 different concentrations (0.1 nM, 1 nM, 10 nM, 100 nM, 1 µM, 10 µM) to subsequently calculate their IC50 values (Table 3). Additionally, all of the synthesised target compounds were submitted to the National Cancer Institute “NCI” Bethesda, Maryland, USA (www.dtp.nci.nih.gov), to evaluate theirs in vitro antiproliferative activities against NCI 60- cancer cell line panel.
Figure 5.
Mean % inhibition of the investigated compounds on EGFR, VEGFR-2 and HDAC6 at 10 µM concentration.
Table 1.
Mean % inhibition of the investigated compounds on EGFR, VEGFR-2 and HDAC6 at 10 µM concentration.
Code
R1
X
Y
EGFR
VEGFR-2
HDAC6
7a
H
O
OH
48%
89%
5%
7b
OCH3
O
OH
39%
78%
6%
7c
CF3
O
OH
11%
100%
8%
7d
H
S
OH
78%
45%
6%
11a
H
O
12%
47%
7%
11b
OCH3
O
5%
75%
13%
11c
CF3
O
11%
91%
8%
12a
H
O
11%
48%
7%
12b
OCH3
O
5%
83%
7%
12c
CF3
O
66%
100%
23%
Red-coloured data represent potent inhibition 70–100%, green-coloured data represent moderate inhibition 30–70%, blue-coloured data represent mild inhibition 20–30%, and the black-coloured data represent very weak inhibition <20%.
Table 3.
The IC50 value of representative synthesised compounds against EGFR and VEGFR-2.
Compound
IC50 in enzyme assays
Code
R1
R2
X
Y
EGFR
VEGFR-2
7c
CF3
–
O
OH
ND
191 nM
12c
CF3
–
O
1.14 µM
185 nM
15c
H
Cl
–
OH
19 nM
5.58 µM
20b
H
H
–
155 nM
ND
20d
Cl
F
–
68 nM
ND
Gefitnib [27]
1.08 µM
ND
Vandetanib [17]
0.5 µM
40 nM
ND: not determined.
Mean % inhibition of the investigated compounds on EGFR, VEGFR-2 and HDAC6 at 10 µM concentration.Mean % inhibition of the investigated compounds on EGFR, VEGFR-2 and HDAC6 at 10 µM concentration.Red-coloured data represent potent inhibition 70–100%, green-coloured data represent moderate inhibition 30–70%, blue-coloured data represent mild inhibition 20–30%, and the black-coloured data represent very weak inhibition <20%.The IC50 value of representative synthesised compounds against EGFR and VEGFR-2.ND: not determined.
In vitro EGFR inhibitory assay
By closely investigating the structure activity relationship, thieno[2,3-d]pyrimidine derivatives with aniline moiety showed an excellent inhibitory profile compared to derivatives with urea and benzimidazole moieties (Tables 1 and 2). Moreover, thieno[2,3-d]pyrimidine derivatives with aniline moiety, mono-substitutions, or di-substitutions at 3- or 4- positions on its terminal phenyl ring exhibited excellent inhibitory profile, as previously reported.In general, EGFR inhibition increases by increasing the number of carbons in the linker chain between the hydroxamic acid moiety and the amide group to three carbons or by directly attaching the hydroxamic acid with the thienopyrimidine fragment. So that compounds (20b and 20c) showed 92 and 96% inhibition, compared to inhibition of 76 and 78% for the two carbons linker of (19b and 19c), respectively (Table 2).
Table 2.
Mean % inhibition of the investigated compounds on EGFR, VEGFR-2, and HDAC6 at 10 µM concentration.
Code
R1
R2
Y
EGFR
VEGFR-2
HDAC6
15a
OH
45%
36%
8%
15b
H
H
OH
98%
32%
30%
15c
H
Cl
OH
97%
87%
7%
15d
Cl
F
OH
97%
54%
4%
19b
H
H
76%
13%
10%
19c
H
Cl
78%
21%
9%
19d
Cl
F
96%
26%
0%
20b
H
H
92%
25%
56%
20c
H
Cl
91%
35%
7%
20d
Cl
F
96%
20%
9%
Red-coloured data represent potent inhibition 70–100%, green-coloured data represent moderate inhibition 30–70%, blue-coloured data represent mild inhibition 20–30%, and the black-coloured data represent very weak inhibition <20%.
Mean % inhibition of the investigated compounds on EGFR, VEGFR-2, and HDAC6 at 10 µM concentration.Red-coloured data represent potent inhibition 70–100%, green-coloured data represent moderate inhibition 30–70%, blue-coloured data represent mild inhibition 20–30%, and the black-coloured data represent very weak inhibition <20%.Remarkably, compound (12c) bearing urea moiety with 3-CF3 at its terminal phenyl ring showed an exceptional significant inhibitory activity towards EGFR kinase activity (66%), compared to other thieno[2,3-d]pyrimidine derivatives bearing urea moiety. Also, compound (7d) bearing thiourea demonstrated noticeable higher inhibitory activity (78% inhibition) than compound (7a) bearing urea moiety (48% inhibition).Promising candidates, (12c, 15c, 20b and 20d) were chosen as representatives for evaluation of potential enzyme inhibitory activity (IC50) (Table 3). Compound (15c) demonstrated the most potent EGFR inhibition with IC50 value of 19 nM since its hydroxyl amine group is directly attached to the thieno[2,3-d]pyrimidine ring and also it has 4-Cl substitution on its terminal phenyl ring. Compound (20d) exhibited more potent inhibition than (20b) with IC50 value of 68 nM compared to 155 nM. Since both (20d and 20b) have a linker of 3 carbons, thus the (4-F, 3-Cl) substitution at the terminal phenyl ring of compound (20d) is suggested to be responsible for the increase of its inhibitory activity compared to (20b), with no substitution at its terminal phenyl ring.Noticeably, compound (12c) bearing urea moiety with 3-CF3 at its terminal phenyl ring exhibited IC50 value of 1.14 µM, in which it showed the highest EGFR inhibition among all other synthesised derivatives bearing urea moiety.
In vitro VEGFR-2 inhibitory assay
On the contrary of EGFR inhibitory assay, most of the thieno[2,3-d]pyrimidine derivatives bearing urea moiety (7a-c, 11b,c, 12b,c) have demonstrated potent inhibition ranging from (75-100%) for VEGFR-2 kinase activity at an initial screening of 10 μM concentration. Compounds with 3-CF3 substitution at its terminal phenyl ring (7c, 11c and 12c) demonstrated the most potent VEGFR-2 inhibition by 100, 91 and 100%, respectively (Table 1). Moreover, by observing the effect of the hydrocarbon linker on the VEGFR-2 inhibitory activity, it was found that the most potent inhibitory activities were observed with derivatives bearing urea moieties having hydroxylamine moiety directly attached to the thieno[2,3-d]pyrimidine ring, such as (7a, 7b and 7c) exhibited 89, 78 and 100% inhibition, respectively. Noticeably, compound (7a) bearing a urea moiety showed higher VEGFR-2 inhibition (89% inhibition) than compound (7d) bearing thiourea moiety (45% inhibition) (Table 1).As for the aniline series, compound (17c) with no linear hydrocarbon linker between the hydroxamic acid moiety and the thieno[2,3-d]pyrimidine fragment and 4-chloro substitution on its terminal phenyl ring showed exceptional potent inhibition (87% inhibition) (Table 2).Promising candidates, (7c, 12c, and 15c) were chosen as representatives for the evaluation of potential enzyme inhibitory activity (IC50) (Table 3). All of the evaluated target compounds exhibited potent VEGFR-2 inhibitory activity with IC50 values ranging from (185 nM to 5.58 µM). Compounds (7c and 12c) bearing urea moiety exhibited potent VEGFR-2 inhibition with IC50 values of 191 and 185 nM, respectively. These results suggest that -CF3 substitution at the terminal phenyl ring of compounds (7c and 12c) and the urea moiety is responsible for their high inhibitory activity, regardless of having 3 carbons in the linker chain (12c) or without having a linker (7c). Compound (15c) with 4-chloro aniline moiety exhibited significant inhibition to VEGFR-2 with a higher IC50 value of 5.58 µM than other derivatives bearing urea moiety.
In vitro HDAC inhibitory assay
Initially, all synthesised compounds were screened against HDAC6 at 10 μM concentration. HDAC6 enzyme inhibition assay was selected since our newly synthesised hydroxamic acid compounds have a large capping group suitable for fitting in the large hydrophobic channel of HDAC6, which is unique to the HDAC6 receptor. The newly designed compounds were synthesized either having short hydrocarbon linkers between the zinc-binding group (ZBG) and the connecting amide unit or having no linear hydrocarbon linker between the hydroxamic acid moiety and the thieno[2,3-d]pyrimidine fragment, as the previously reported potent HDAC6 inhibitors,.Screening of the synthesised compounds against HDAC6 at 10 μM concentration is shown in (Tables 1 and 2). For derivatives bearing urea moiety, compound (12c) with three carbons in the linker region showed considerable HDAC6 inhibition by 23% (Table 1). For derivatives bearing aniline moiety, compound (20b) exhibited the most significant HDAC6 inhibition by 56% (Table 2). Compound (20b); bearing unsubstituted aniline moiety; is the smallest capping group among all synthesized compounds. The small capping group may allow it to fit in the HDAC6 hydrophobic channel. Additionally, compound (20b) had 3 carbons in the linker region compared to (15b) with no linear hydrocarbon linker and (19b) with 2 carbons in the linker region showing 30% and 10% inhibition, respectively (Table 2). Thus, having a small capping group of unsubstituted aniline and a longer carbon chain in the linker region between the hydroxamic acid moiety and the connecting amide is recommended for optimum HDAC6 inhibition.In conclusion, potential novel multitarget hybrid lead inhibitors acting on EGFR, VEGFR-2 and HDAC6 are considered, such as compounds (12c, 15 b and 20 b). Compound (12c) exhibited potent VEGFR-2 inhibition (100%, IC50 = 185 nM), significant EGFR inhibition (66%, IC50 =1.14 µM) and mild HDAC6 inhibition (23%). Compound (15 b) exhibited potent EGFR inhibition (98%), moderate VEGFR-2 (32%) and moderate HDAC6 inhibition (30%). Compound (20 b) exhibited potent EGFR inhibition (92%, IC50 = 155 nM), moderate HDAC6 inhibition (56%) and mild VEGFR-2 inhibition (25%), (Figure 5).Moreover, many of the synthesised novel compounds can be considered as promising leads of EGFR and VEGFR-2 dual inhibitors, (Table 2 and Figure 5). Where, compound (15c) is considered to be the most potent dual inhibitor among all the synthesised compounds, since it exhibits potent EGFR and VEGFR-2 inhibition (97%, IC50 = 19 nM) and (87%, IC50 =5.58 µM), respectively. Compound (15d) exhibited potent EGFR and moderate VEGFR-2 inhibition (97%) and (54%), respectively. Two selective potent inhibitors were also considered, such as (7c and 20d). Compound (7c) is considered as a potent novel VEGFR-2 inhibitor (100% inhibition, IC50 = 191 nM), while compound (20d) is considered as a potent novel EGFR inhibitor (96% inhibition, IC50 = 68 nM).
In vitro antiproliferative activity
All of the synthesised structures were selected and submitted to NCI based on computer modelling techniques and the degree of diversity of structures for the evaluation of their antiproliferative activity. The 20 synthesised compounds were selected with their respective NCI codes [NCS D-820157- NCS D-820176] for the inhibition percent assay at an initial single dose of 10 μM on the full NCI 60 cell panel, where the mean graph of the percent growth for all compounds on the full NCI cell panel is evaluated, compared to the untreated control cells (Supplementary Material). In vitro cell line inhibition % of all synthesised compounds is represented in (Table 1 of Supplementary Material).First, in series (7a-d) bearing urea moiety with no linear hydrocarbon linker between the hydroxamic acid moiety and the thieno[2,3-d]pyrimidine fragment, compound (7c) with urea moiety and 3-CF3 at the terminal phenyl ring showed the most significant anti-proliferative activity among all of the synthesised compounds, where it exhibited inhibitory activity against the central nervous system (CNS) cancer [SNB-75 cancer cell line (46% inhibition) and SF-539 cancer cell line (38.6% inhibition)] and colon cancerKM12cancer cell line (33.7% inhibition). (7c) exhibited inhibition activity also against CNS cancer [U251cancer cell line (26.9% inhibition) and SF-268 cancer cell line (10.8% inhibition)], renal cancerUO-31 (24.7% inhibition), leukaemia MOLT-4 (11.8% inhibition), NSCLC [NCI-H522cancer cell line (11.7% inhibition) and NCI-H23 (11.1% inhibition)] and breast cancerMCF7cancer cell line (11.6% inhibition). While in series (15a-d) bearing aniline moiety with no linear hydrocarbon linker between hydroxamic acid moiety and thieno[2,3-d]pyrimidine fragment, compound (15c) with 4-Claniline exhibited inhibitory activity mostly against renal cancer UO-31 cancer cell line (28.8% inhibition), NSCLC [NCI-H522cancer cell line (25.4% inhibition) and EKVX cancer cell line (20.6% inhibition)] and CNS cancer SNB-75 cell line (20.4% inhibition).In series (12a-c) bearing urea moiety with three carbons in the linker chain between the hydroxamic acid moiety and the amide group, compound (12c) with urea moiety and CF3 at the terminal phenyl ring showed moderate inhibitory activity mostly against CNS cancer [SNB-75 cancer cell line (29.2% inhibition) and SF-539 cell line (28.5% inhibition)]. In series (20 b-d) bearing aniline moiety with three carbons in linker chain between the hydroxamic acid moiety and the amide group, compound (20d) exhibited inhibitory activity mostly against renal cancer UO-31 cancer cell line (23.7% inhibition), breast cancerMDA-MB-468 (16.2% inhibition) and NSCLCNCI-H226 (11.3% inhibition).In series (11a-c) and (19 b-d) with two carbons in the linker chain between the hydroxamic acid moiety and the amide group, most of the compounds showed lower inhibitory activities against most of the investigated cell lines.Recent studies revealed that VEGFR-2 is a major contributor to the growth of CNS tumours and its interruption may effectively suppress tumour growth in CNS cancer. It was also reported that VEGFR-2 is a promising platform for anti-angiogenesis cancer therapy for colon cancer cells. Also, EGFR was reported to be overexpressed in different cancer cell lines, including NSCLC, renal cancer, and brain cancer. Noticeably, it is observed that most of the compounds showed weak to moderate inhibition towards Renal UO-31 and CNS cancer SNB-75 cell lines (20–46% inhibition).From these findings, it is suggested that the potent in vitro inhibition of compound (7c) towards VEGFR-2 (100% inhibition, IC50 = 191 nM) contributed to the significant inhibitory activity against CNS [SNB-75 cancer cell line (46% inhibition) and SF-539 cancer cell line (38.6% inhibition)] and colonKM12cancer cell line (33.7% inhibition). It is also suggested that the potent in vitro inhibition of compound (15c) towards EGFR (97% inhibition, IC50 = 19 nM) and VEGFR2 (87% inhibition, IC50 = 5.58 µM) contributed to the significant inhibitory activity against renal, CNS and NSCLC cancer cell lines as discussed. Additionally, it is also observed that (15c) is the only compound that showed considerable inhibition towards NSCLCH522 (25.4% inhibition) and EKVX (20.6% inhibition) cancer cell lines, which is contributed to its potent EGFR inhibition with (IC50 = 19 nM).It is obvious that despite most of our newly synthesised compounds exhibited potent EGFR and/or VEGFR-2 in vitro inhibitory activities (90–100% inhibition) but unexpectedly some of the compounds showed moderate in vitro inhibitory activity against NCI cell lines at 10 µM concentration. It is suggested that maybe the increased polarity of our synthesised hydroxamic acids, might contribute to poor cell membrane permeability towards NCI cancer cell lines. Also, VEGFR-2 inhibitors showed less sensitivity towards NCI cancer cell lines than human vascular endothelial cells, which contributed to their moderate antiproliferative activity towards NCI cancer cell lines.
Molecular modelling studies
Docking study
A molecular docking study using Discovery Studio 4.5 was performed using the C-DOCKER algorithm. A docking study is performed to analyse binding affinities, binding modes, and orientation of the newly synthesised target compounds into the active binding site of three different receptors such as EGFR, VEGFR-2, and HDAC6. The best docking pose is selected based on its binding mode in comparison to that of the reference ligand. The crystal structure of EGFR, VEGFR-2, and HDAC6 co-crystallised with their reference ligands were downloaded from protein data bank (www.rcsb.org) with PDB codes: 4HJO, 3VHE, and 5G0H, respectively. Re-docking the reference compounds into the X-ray crystal structure of the active site of its receptor was performed for validation.
Docking study on EGFR
The first docking study was done on EGFR (PDB code 4HJO) co-crystallised with the reference ligand (erlotinib). The root mean square difference (RMSD) between the re-docked and the co-crystallised conformers of erlotinib is 0.5 Å, which signifies the validity of the C-DOCKER algorithm. Besides, re-docked erlotinib retrieved the reported binding mode into the active binding site of EGFR (-CDOCKER energy = 27.18 Kcal/mol) as depicted in (Figure 6). Where nitrogen of pyrimidine ring forms a hydrogen bond with NH of MET769, pyrimidine ring forms hydrophobic interactions with ALA719, MET769 and LEU820, phenyl ring forms hydrophobic interaction with LEU694 and LEU820 and the terminal phenyl ring forms hydrophobic interaction with VAL702, ALA719 and LYS721 (Figure 6). Docking of all synthesised compounds was done on EGFR (PDB code 4HJO), where the binding mode of representative thieno[2,3-d]pyrimidine derivatives were investigated to justify the potent or significant EGFR inhibition (high % inhibition or low IC50) exhibited by some of the synthesised compounds (Table 2 of Supplementary Material). In general, docking results revealed that thieno[2,3-d]pyrimidine derivatives with aniline moiety (15c, 15d, 19d, 20 b-d) showed higher binding affinities than erlotinib, fulfilled the reported binding mode and showed extra interactions with EGFR active site. Although thieno[2,3-d]pyrimidine derivatives with urea moiety showed higher binding affinities than erlotinib but it did not fulfil all the key interactions with the active site of EGFR. The binding mode and energies of three representatives docked compounds are presented in (Figure 7). Compound (15c) showed the most potent EGFR inhibitory activity (97%, IC50 = 19 nM) and revealed better binding affinity than erlotinib (-CDOCKER energy = 28.34 Kcal/mol). The binding mode of (15c) was consistent with erlotinib as shown in (Table 2 of Supplementary Material). Additionally, the methyl group of (15c) showed extra hydrophobic interactions with VAL702 and LEU820 and chlorine showed extra hydrophobic interactions with adjacent hydrophobic pockets including LEU764 and LEU834. Compound (20d) exhibited potent inhibitory activity towards EGFR (97%, IC50 = 68 nM) and revealed better binding affinity than erlotinib (-CDOCKER energy = 42.52 Kcal/mol). The binding mode was consistent with erlotinib as shown in (Table 2 of Supplementary Material), where the oxygen of the hydroxyl group showed 2 additional hydrogen bonds with LYS692 and VAL693, oxygen of hydroxamic carbonyl group and hydrogen of hydroxamic amine showed extra hydrogen bonds with LYS692 and LEU694, respectively. Moreover, the methyl group of (20d) showed additional hydrophobic interaction with VAL702, and chlorine showed additional hydrophobic interaction with ALA719 and LYS721. Compound (12c) bearing urea moiety exhibited significant EGFR inhibitory activity (66%, IC50 = 1.14 µM), but it showed lower inhibition than other compounds bearing aniline moiety. Compound (12c) revealed a missing hydrogen bond with the essential MET796 although having better binding affinity than erlotinib (-CDOCKER energy = 43.31 Kcal/mol). Compound (12c) showed five extra hydrogen bonds, where the oxygen of carbonyl urea showed hydrogen bond with THR766, oxygen of hydroxamic carbonyl group showed hydrogen bond with CYS773, oxygen of hydroxyl group showed hydrogen bond with PHE771, F of terminal phenyl ring showed hydrogen bond with PHE832 and hydrogen of ureaamine showed hydrogen bond with ASP831. Additionally, the three fluorine atoms of CF3 showed three halogen interactions with CYS751 and one halogen interaction with ARG752. Pyrimidine ring of (12c) forms hydrophobic interactions with LEU694 and LEU820, thiophene ring forms hydrophobic interaction with LEU694, phenyl ring forms hydrophobic interaction with VAL702 and LEU820 and terminal phenyl ring shows four hydrophobic interactions with extra hydrophobic pocket including LEU753, LEU764, LEU834, and MET742. Moreover, a methyl group and C of CF3 of (12c) showed additional hydrophobic interactions with LEU694 and CYS751, respectively.
Figure 6.
(a) 3D image of binding interaction of Erlotinib (-CDOCKER energy = 27.18) with EGFR (with PDB code: 4HJO) using Accelrys Discovery Studio 4.5. (b) 2D image of binding interaction of Erlotinib with EGFR. Hydrogen bond interaction is described as green dotted line, hydrophobic interactions are described as purple-dotted lines and electrostatic interactions are described as orange-dotted lines.
Figure 7.
(a) 3D images of binding interaction of compound (15c) (-CDOCKER energy = 28.34) with EGFR (with PDB code: 4HJO), (b) 3D images of binding interaction of compounds (20d) (-CDOCKER energy = 42.52) with EGFR (c) 3D images of binding interaction of compound (12c) (-CDOCKER energy = 43.31) with EGFR. Hydrogen bond interaction is described as a green-dotted line, hydrophobic interactions are described as purple dotted lines and halogen bonds are described as cyan-dotted lines.
Design strategy of hybrid multitarget (EGFR, VEGFR, and HDAC inhibitors). Colour codes show essential pharmacophore for HDAC inhibitors composed of surface recognition moiety (purple), linker (red), and zinc-binding group (blue). The cyan oval shape represents the area binding to the hinge region of the ATP binding site of VEGFR2 and EGFR and the green square represents the area binding to the hydrophobic pocket of VEGFR2 and EGFR, while the rest of the structure represents the area fitting in the solvent-accessible area of VEGFR2 and EGFR.(a) 3D image of binding interaction of Erlotinib (-CDOCKER energy = 27.18) with EGFR (with PDB code: 4HJO) using Accelrys Discovery Studio 4.5. (b) 2D image of binding interaction of Erlotinib with EGFR. Hydrogen bond interaction is described as green dotted line, hydrophobic interactions are described as purple-dotted lines and electrostatic interactions are described as orange-dotted lines.(a) 3D images of binding interaction of compound (15c) (-CDOCKER energy = 28.34) with EGFR (with PDB code: 4HJO), (b) 3D images of binding interaction of compounds (20d) (-CDOCKER energy = 42.52) with EGFR (c) 3D images of binding interaction of compound (12c) (-CDOCKER energy = 43.31) with EGFR. Hydrogen bond interaction is described as a green-dotted line, hydrophobic interactions are described as purple dotted lines and halogen bonds are described as cyan-dotted lines.
Docking study on VEGFR-2
The second docking study was done on VEGFR-2 (PDB code: 3VHE) co-crystalised with the reference compound; (pyrrolo[3,2-d]pyrimidine derivative;). RMSD between the re-docked and the co-crystallised conformer of the reference is 0.5 Å, which proves the validity of the C-DOCKER algorithm. Noticeably, redocked pyrrolo[3,2-d]pyrimidine reference compound (-C-DOCKER energy = 42.91 Kcal/mol) retrieved the reported binding mode into the active site of VEGFR-2 (pdb code: 3VHE) as depicted in (Figure 8). Where nitrogen of pyrimidine ring showed hydrogen bond with CYS919 in the ATP-binding site (hinge region), the oxygen of urea carbonyl group formed hydrogen bond with ASP1046, while the two NH of urea moiety showed bifurcate hydrogen bonds with GLU885. Additionally, the pyrrole ring showed hydrophobic interactions with CYS919, LEU840, ALA866, VAL848 and LEU1035, the pyrimidine ring formed hydrophobic interactions with CYS919, ALA868, VAL848 and LEU1035. Furthermore, the phenyl ring showed hydrophobic interactions with VAL848, VAL899, LYS868, CYS1045 and pi-pi interaction with PHE1047. Finally, the terminal phenyl ring formed hydrophobic interaction with LEU889. Docking of all synthesised derivatives was done on VEGFR-2 (PDB code: 3VHE), where the binding mode of representative thieno[2,3-d]pyrimidine derivatives was explored to justify the potent or significant VEGFR-2 inhibition exhibited by some of the synthesised compounds (Table 3 of Supplementary Material).
Figure 8.
(a) 2D image of binding interaction of reference with VEGFR-2 (PDB code: 3VHE), (-CDOCKER energy = 42.91). (b) 3D image of binding interaction of reference with VEGFR-2 (with PDB code: 3VHE). Hydrogen bond interaction is described as a green-dotted line, hydrophobic interactions are described as purple dotted lines, Pi-Pi interaction is described as a pink-dotted line.
(a) 2D image of binding interaction of reference with VEGFR-2 (PDB code: 3VHE), (-CDOCKER energy = 42.91). (b) 3D image of binding interaction of reference with VEGFR-2 (with PDB code: 3VHE). Hydrogen bond interaction is described as a green-dotted line, hydrophobic interactions are described as purple dotted lines, Pi-Pi interaction is described as a pink-dotted line.In general, docking results revealed that most of the thieno[2,3-d]pyrimidine derivatives with urea moiety (7c and 12c) showed higher binding affinities than the reference compound and compounds with aniline moiety (15c) which is consistent with VEGFR-2 inhibition assay. In which, compounds bearing urea moiety, such as (12c) with (IC50 = 185 nM), and (7c) with (IC50 = 191 nM) exhibited potent VEGFR-2 inhibition compared to compound (15c) bearing aniline moiety with (IC50 = 5.58 µM). Where (12c) and (7c) retrieved all the essential hydrogen bond and hydrophobic interactions of the reference compound, while compound (15c) missed the essential hydrogen bond interactions with GLU885 and ASP1046 although its pyrimidine ring retrieved the essential hydrogen bond with CYS919 (Table 3 of Supplementary Material). The binding interactions and energies of two representatives docked compounds are presented in (Figure 9). Additionally, we found that compounds (12c) and (7c); both bearing urea moiety with CF3 substituent and compound (15c) bearing 4-chloro aniline moiety showed higher binding affinities and made extra hydrophobic bonds than the reference compound. Where CF3 in compound (12c) made extra hydrophobic bonding with ILE892, LEU1019 and HIS1026, F of CF3 in compound (7c) made an extra halogen hydrogen bond with HIS1026 and 4-Cl of the terminal phenyl ring in compound (15c) forms additional hydrophobic interaction with LEU889, VAL899 and LYS868. According to previously reported SAR, this area is an allosteric hydrophobic binding pocket revealed for type II VEGFR inhibitors, in which the phenylalanine residue of the DFG loop flips out of its lipophilic pocket to increase the binding affinity of inhibitors to VEGFR receptors and increase residence time,.
Figure 9.
(a) 3D image of binding interactions of compound (7c) (-CDOCKER energy = 51.35) with VEGFR-2 (PDB code: 3VHE). (b) 3D image of binding interaction of compound (15c) (-CDOCKER energy =30.81) with VEGFR-2. Hydrogen bond interaction is described as green-dotted line.
(a) 3D image of binding interactions of compound (7c) (-CDOCKER energy = 51.35) with VEGFR-2 (PDB code: 3VHE). (b) 3D image of binding interaction of compound (15c) (-CDOCKER energy =30.81) with VEGFR-2. Hydrogen bond interaction is described as green-dotted line.
Docking study on HDAC6
The final docking study was done on HDAC6 (PDB code: 5G0H) co-crystallised with the reference compound (Trichostatin A; TSA). RMSD is 0.44 Å, which signifies the validity of CDOCKER algorithm. Moreover, the reported binding mode of TSA was retrieved into the X-ray crystal structure of HDAC6 (-CDOCKER energy = 25.49) as shown in (Figure 10), where the oxygen of hydroxamic hydroxyl group forms a hydrogen bond with HIS573 and oxygen of both hydroxamic carbonyl and hydroxyl groups forms coordinate bonds with zinc metal in a bidentate fashion with distances 2.24 Å and 2.20 Å, respectively. Additionally, the phenyl ring of TSA shows pi-pi stacking with PHE643 and the aliphatic chain of the linker showed hydrophobic interaction with PHE583 and PHE643 (Figure 10). Docking of all synthesised compounds was performed on HDAC6 (PDB code: 5G0H), where the binding mode of representative thieno[2,3-d]pyrimidine derivatives were explored to justify the significant inhibition of compound (20 b) (56% inhibition) against HDAC6 among other derivatives that showed remarkably lower % inhibition against HDAC6 (Table 4 of Supplementary Material).
Figure 10.
(a) The 2D image of binding interactions of Trichostatin A with HDAC6 (PDB code: 5G0H). (b) 3D image of binding interaction of redocked TSA with HDAC6 (-CDOCKER ENERGY= 25.54), where hydrogen bond interaction is described as green-dotted line, hydrophobic interactions are described as purple-dotted lines, pi-pi stacking is described as pink-dotted line and coordinate bond with zinc metal is described as grey-dotted line (2 D image) and red-dotted line (3 D image).
(a) The 2D image of binding interactions of Trichostatin A with HDAC6 (PDB code: 5G0H). (b) 3D image of binding interaction of redocked TSA with HDAC6 (-CDOCKER ENERGY= 25.54), where hydrogen bond interaction is described as green-dotted line, hydrophobic interactions are described as purple-dotted lines, pi-pi stacking is described as pink-dotted line and coordinate bond with zinc metal is described as grey-dotted line (2 D image) and red-dotted line (3 D image).The binding interactions and energies of three representatives docked compounds (12c, 19 b and 20 b) are presented in (Figure 11). Compound (20 b) with unsubstituted aniline moiety and three carbons in the linker region exhibited the most significant HDAC6 inhibition among all the synthesised compounds (56% inhibition) and it revealed higher binding affinity than TSA (-CDOCKER energy = 34.32). This significant inhibitory activity for (20 b) was justified by its consistent binding mode with the reported binding mode of TSA against HDAC6 (PDB: 5G0H), where the oxygen of hydroxamic carbonyl and hydroxyl groups forms coordinate bonds with zinc metal in a bidentate fashion with distances 2.24 Å and 2.36 Å, respectively. The oxygen of the hydroxamic hydroxyl group forms a hydrogen bond with HIS573 and oxygen of hydroxamic and amidic carbonyl group forms two additional hydrogen bonds with TYR745 and HIS614, respectively. Also, the thiophene ring of (20 b) shows pi-pi stacking with PHE643. While, compound (19 b) with unsubstituted aniline moiety and two carbons in the linker region exhibited insignificant inhibitory activity (10%), although it revealed higher binding affinity (-CDOCKER energy = 37.68) than (20 b) (Figure 11). By comparing the binding mode of compounds (20 b) and (19 b), it was found that the insignificant inhibitory activity of the compound (19 b) is certainly accounted for by having a shorter linker between the hydroxamate moiety and the connecting amide group (two carbons) than compound (20 b) (3 carbons) since they both bear the same hydrophobic cap group. The shorter two carbons linker of compound (19 b) contributes to its failure to form a coordinate bond with zinc metal in a bidentate fashion, where the oxygen of hydroxamic carbonyl group forms only one coordinate bond with zinc metal (Figure 11). Whereas compound (12c) with urea moiety exhibited lower HDAC6 inhibition (23%) than (20 b), although it revealed higher binding affinity (-CDOCKER energy = 37.58) than (20 b) and TSA. The decrease in HDAC6 inhibition by compound (12c) is justified by its binding mode as described in (Figure 11), in which the oxygen of hydroxamic carbonyl group forms only one coordinate bond with zinc metal in a monodentate fashion, despite fulfilling all other key interactions with reported essential amino acids (Table 4 of Supplementary Material). Unexpectedly, compound (15 b) with unsubstituted aniline moiety, exhibited considerate HDAC6 inhibition (30%), although (i) it has no linear hydrocarbon linker between the thieno[2,3-d]pyrimidine fragment and its hydroxamate moiety. (ii) it has a low binding affinity to HDAC6 of (-CDOCKER energy = 15.78). (iii) it could not form a coordinate bond between the hydroxamate moiety and the zinc metal. (iv) it also missed the essential hydrogen bond with HIS573 as shown in (Table 4 of Supplementary Material). Finally, future directions towards optimisation of the identified lead inhibitor (20 b) involve: (i) hydrocarbon linker from five to seven carbons to enhance the potency of novel thieno[2,3-d]pyrimidine derivatives against HDACs, to be able to capture the zinc metal in bidentate fashion with shorter distances as TSA. (ii) Hydrophobic cap group of thieno[2,3-d]pyrimidine derivatives bearing aniline moiety instead of the larger thieno[2,3-d]pyrimidine derivatives with biphenyl urea or thiourea moiety, in which it has a better binding affinity and showed higher inhibition to HDAC6.
Figure 11.
(a) 3D image of binding interactions of (20 b) (-CDOCKER ENERGY= 34.32) with HDAC6 (PDB code: 5G0H). (b) 3D image of binding interaction of (19 b) (-CDOCKER ENERGY= 37.68) with HDAC6. (c) 3D image of binding interaction of (12c) (-CDOCKER ENERGY= 37.58) with HDAC6. Hydrogen bond interaction is described as a green-dotted line, Pi-Pi interaction is described as a pink-dotted line, halogen interaction is described as cyan dotted lines and coordinate bond with zinc metal is described as a red dotted line.
(a) 3D image of binding interactions of (20 b) (-CDOCKER ENERGY= 34.32) with HDAC6 (PDB code: 5G0H). (b) 3D image of binding interaction of (19 b) (-CDOCKER ENERGY= 37.68) with HDAC6. (c) 3D image of binding interaction of (12c) (-CDOCKER ENERGY= 37.58) with HDAC6. Hydrogen bond interaction is described as a green-dotted line, Pi-Pi interaction is described as a pink-dotted line, halogen interaction is described as cyan dotted lines and coordinate bond with zinc metal is described as a red dotted line.
In silico ADMET predictive study
A computer-aided ADMET study was performed using Accelrys Discovery Studio 4.5 software protocols to investigate pharmacokinetics properties of newly synthesised thieno[2,3-d] pyrimidine hydroxamic acid derivatives and their ester derivatives, such as absorption, distribution, metabolism, excretion and toxicity. This study was done to investigate the unexpected moderate in vitro antiproliferative activity of most of the synthesised hydroxamic acid compounds despite their potent in vitro EGFR and VEGFR-2 enzyme inhibition. Aqueous solubility, blood brain barrier (BBB) penetration, CYP2D6 binding, human intestinal absorption (HIA), toxicity and plasma protein binding descriptors are predicted. 2 D ADMET plot for thieno[2,3-d] pyrimidine hydroxamic acid derivatives is drawn between ADMET_PSA_2D versus ADMET_AlogP98 (Figure 1 of Supplementary Material). Two analogous 95% and 99% confidence ellipses are shown in the biplot corresponding to HIA and BBB models. The results of the calculated representative parameters for thieno[2,3-d] pyrimidine hydroxamic acid derivatives are tabulated in (Table 5 of Supplementary Material). There is an inverse relationship between polar surface area (PSA) with percent HIA and thus cell membrane permeability. ALogP is used to indicate the lipophilicity of the investigated compounds, thus calculating ALogP with PSA is used to bring information about hydrogen bond characteristics. Optimum cell membrane permeability is predicted if the compound has (PSA < 140˚A2 and AlogP98 < 5). In this study, all hydroxamic acid compounds showed ALogP98 < 5, 14 compounds showed high accepted values of PSA: 140˚A2 < PSA > 100˚A2, indicating considerable cell membrane permeability but not highly permeable due to its high value of PSA (>100˚A2), which is inversely proportional with cell membrane permeability. Whereas, the rest of the compounds bearing urea moiety and with linker between the hydroxamic acid moiety and the amide connecting unit (11a-c, 12a-c) showed PSA > 140, which indicated poor cell membrane permeability and absorption.By comparing the AlogP98 and the PSA of the synthesised thieno[2,3-d] pyrimidine hydroxamic acids and their esters (Table 6 of Supplementary Material), it was found that almost all of the prepared thieno[2,3-d]pyrimidine esters showed lower PSA (PSA < 140) and higher AlogP than most of the synthesised target thieno[2,3-d]pyrimidine hydroxamic acids. These results coincide with the unsatisfactory poor to moderate in vitro antiproliferative activities of the hydroxamic acid derivatives which could be attributed to their poor cell membrane permeability. Also, it was found that the synthesised hydroxamic acid compounds bearing aniline moiety were predicted to show better intestinal absorption and cell membrane permeability than compounds bearing urea moieties. While compounds (7a-d) bearing urea moiety with no linear hydrocarbon linker between the hydroxamic acid moiety and the thieno[2,3-d]pyrimidine fragment were predicted to have acceptable cell membrane permeability as an exception.
Conclusion
In the current study, series of hybrids of thieno[2,3-d]pyrimidine fragment bearing arylamino and diaryl urea moieties at 4-position linked to hydroxamic acid moiety were designed and synthesised as inhibitors for EGFR/VEGFR-2/HDAC6. Compounds (12c, 15 b and 20 b) may be considered as potential multitarget inhibitors against EGFR/VEGFR-2/HDAC6. Additionally, other synthesised novel hydroxamic acid compounds may be considered as EGFR/VEGFR-2 dual inhibitors, such as (7a, 15c and 15d), in which the hydroxamic acid moiety is directly attached to the thieno[2,3-d]pyrimidine fragment, contributing to maximum inhibitory dual activity. Compound (15c) could be the most potent dual EGFR/VEGFR-2 inhibitor in this study. Moreover, promising selective VEGFR-2 (7c and 11c) and EGFR (20d) inhibitors were identified. The results of the enzyme inhibitory assays were found to coincide with the molecular docking studies performed on the three enzymes EGFR, VEGFR-2 and HDAC6.These results supported the rationale of the design strategy, where the thieno[2,3-d]pyrimidine fragment bearing urea moiety or aniline moiety may be considered as a promising capping group. Unexpectedly, it is observed that most of the compounds showed common inhibition (weak to moderate inhibition) towards Renal UO-31 and CNS cancer SNB-75 cell lines (20–29% inhibition) although they showed potent enzyme inhibitory profiles. Results of the ADMET study revealed poor cell membrane permeability, explaining their unsatisfactory in vitro antiproliferative activities. Future structural modifications will be studied by implementing specific structural modifications to optimise the hydrophobicity of the designed compounds for enhancing their cell membrane permeability.
Experimental protocols
Reagents and starting materials were obtained from Alfa-Aeser organics and Sigma Aldrich without further purification. Melting points were uncorrected and recorded using BUCHI B-540 apparatus. Thin layer chromatography (TLC) is performed on 0.255 mm silica gel plates and is purchased from Merck to monitor reactions under U.V. light (254 nm). Mass spectrum was recorded in Thermo Scientific GCMS model on Direct Inlet part to mass analyser ISQ at the Regional Centre for Mycology and Biotechnology, Al-Azhar University. 1H and 13 C NMR spectra were recorded on a Bruker Avance III HD FT-high resolution- NMR 400 MHz and scaled as δ in ppm at the Centre for Drug Discovery Research and Development, Faculty of Pharmacy, Ain Shams University. Description of coupling patterns is as follows: s, singlet; br s, singlet broad; d, doublet; t, triplet; m, multiplet; and 1H, 2H, 3H, etc. The coupling constants are described as J, which was rounded off to one decimal place. IR spectra were done at the Micro analytical centre, Cairo University, and recorded on a 4100 Jasco spectrophotometer. Elemental analysis was performed at the Regional Centre for Mycology and Biotechnology, AL Azhar University.
A solution of ethyl cyanoacetate (0.01 mol), ethyl acetoacetate (0.01 mol), sulphur (0.01 mol) and morpholine (0.01 mol) was stirred and heated (70 °C) in absolute ethanol (30 ml) for 4 h. The mixture was then left for 24 h at 0 °C. The formed solid was filtered, rinsed with ethanol (20 ml), dried and then crystallised from absolute ethanol to yield the titled compound as yellow crystals (1) (1.85 g, 72%), m.p. 107–109 °C (as reported),.
A solution of (1) (2 g, 7.8 mmol) in formamide (16 ml) and acetic acid (0.5 ml) was heated under reflux at 150 °C for 40 h, cooled to 70 °C and then 90 ml water was added. The suspension was then cooled to room temperature with stirring. The formed solid was collected by filtration, washed with water (2 × 10 ml) and dried to give the titled compound as pale yellow solid (2) (1.8 g, 90%), m.p. 246–248 °C (as reported),.
A mixture of compound (2) (3.5 g, 14.7 mmol; 1 equiv) and phosphorus oxychloride (29 ml, 278 mmol; 18.9 equiv.) was refluxed at 90 °C for 4 h. The mixture was then slowly poured on iced cold water with continuous stirring. The formed solid was filtered off immediately and dried thoroughly over anhydrous sodium sulphate to afford the titled compound (3) (3.46 g, 92%) as light brown crystals, m.p. 114–115 °C (as reported).
General procedure for preparation of compounds (4a-d)
The appropriate isocyanate (7 mmol;1 equiv.) was added to a solution of p-nitroaniline (7 mmol;1 equiv.) in dry methylene chloride (20 ml) and stirred at room temperature for 24 h. The formed solid was filtered, stirred again with dry methylene chloride then collected by filtration and dried. Crystallisation from ethanol yielded the titled compounds (4a-d).
1–(4-nitrophenyl)-3-phenylurea (4a)
Yield 55% as yellowish white crystals, m.p. 226–228 °C (as reported).
1–(3-methoxyphenyl)-3–(4-nitrophenyl)urea (4 b)
Yield 50% as yellow crystals, m.p. 210–212 °C (as reported).
Yield 69% as greenish yellow crystals, m.p. 258–260 °C (as reported).
1–(4-nitrophenyl)-3-phenylthiourea (4d)
Yield 59% as yellow crystals, m.p. 147–149 °C (as reported).
General procedure for preparation of compounds (5a-d)
A mixture of the appropriate nitro phenyl urea or thiourea derivatives (4a-d) (4 mmol) and Pd-C (0.1 g, 10%) in methanol (100 ml) was stirred under H2 at 60 bars, at room temperature for 6 h. Then Pd-C was completely removed through filtration over celite. The filtrate was evaporated under vacuum and dried to afford titled compounds (5a-d) which were then recrystallized from its suitable solvent.
1–(4-aminophenyl)-3-phenylurea (5a)
Yield 85%, crystallised from ethyl acetate/hexane as white crystals, m.p. 297–300 °C (as reported).
1–(3-methoxyphenyl)-3–(4-aminophenyl)urea (5 b)
Yield 70%, crystallised from ethanol as white crystals, m.p. 164–168 °C (as reported).
Yield 80%, crystalised from ethanol as greyish white crystals, m.p. 148–150 °C (as reported).
1–(4-aminophenyl)-3-phenylthiourea (5d)
Crystallised from ethanol as dark purple liquid, b.p. 410 °C (as reported).
General procedure for preparation of 5-aminobenzimidazole (13a)
A mixture of 5-nitro benzimidazole (1 g, 30.65 mmol) and 10% Pd/C (0.1 g) in methanol (40 ml) was stirred under H2 at 60 bars, at room temperature for 6 h. The Pd/C was removed completely through filtration over celite, then the filtrate was evaporated under vacuum and dried to afford brown crystals of (13a) (0.84 g, 85%), m.p. 165–166 °C (as reported).
General procedure for preparation of compounds (6a-d) and (14a-d)
To a solution of the ethyl (4-chloro-5-methylthieno[2,3-d]pyrimidine)-6-carboxylate (3) (0.256 g, 1 mmol; 1 equiv.) in ethanol (15 ml), the appropriate amino phenyl urea and thiourea derivatives (5a-d)/5 amino benzimidazole derivative (13a)/appropriate aniline derivatives (13 b-d) (1 mmol; 1 equiv.) and triethylamine TEA (0.3 ml, 2 mmol; 2equiv.) were added. The mixture was heated for 24–48 h under reflux. The formed solid was filtered, washed with hot ethanol, dried, and crystallised from appropriate solvent to afford the titled compounds (6a-d) and (14a-d).
Yield 97% as faint yellow crystals, crystallised from ethanol, m.p 154–156 °C (as reported).
General procedure for preparation of compounds (7a-d) and (15a-d)
A cooled solution of sodium (27 mmol; 10 equiv.) in absolute ethanol (20 ml) was added to a solution of hydroxylamine hydrochloride (54 mmol; 20 equiv.) in absolute ethanol cooled in an ice bath. The mixture was stirred for 15 min. The formed precipitate was filtered, and the free hydroxylamine solution was prepared. The freshly prepared hydroxylamine solution was placed in a rounded bottom flask and cooled in an ice bath. Compounds (6a-d) and (14a-d) (2.7 mmol; 1 equiv.) were added to the solution and the mixture was stirred overnight. The reaction is monitored by TLC till the complete disappearance of starting material using an appropriate eluting system [DCM: MeOH (9.5:0.5)]. The solvent was removed under vacuum to afford crude products, washed with water and crystalised from ethanol to afford target compounds (7a-d) and (15a-d),.
Yield 91% as buff solid, m.p 320 °C. 1.72 (s, 1H, OH, exchangeable by D2O), 3.07 (s, 3H, CH3), 7.39–7.44 (t, 1H, aromatic H5′, J = 8 Hz), 7.65–7.69 (m, 1H, aromatic H2′), 7.96–7.99 (dd, 1H, aromatic H6′, J = 4, 12 Hz) , 8.41 (s, NH, exchangeable by D2O), 8.42 (s, 1H, pyrimidine H), 8.50 (s, 1H, NH, exchangeable by D2O). MS: 354 (M++2, 20%), 353 (M++1, 10%), 352 (M+, 100%). Anal. Calcd. For CC 47.67, H 2.86, N 15.88, S 9.09; Found: C 47.80, H 3.12, N 16.15, S 9.23. FT-IR (ú max, cm: 3441 (OH/NH broad), 2985 (CH aromatic), 2916–2769 (CH aliphatic), 1612 (C = O), 1573 (NH bend).
General procedure for preparation of compounds (8a-c, 16b-d)
An aqueous solution of lithium hydroxide (30 mmol) in water (50 ml) was added to a solution of compounds (6a-c, 14b-d) (10 mmol) in THF (30 ml) and ethanol (9 ml). The mixture is stirred for 7 h at 50–60 °C. Complete hydrolysis is monitored by TLC, followed by the addition of ethyl acetate (50 ml), 1 N HCl (10 ml) and water (20 ml). The aqueous layer was separated from the organic layer and washed with ethyl acetate (2 × 50 ml). Then, the combined organic solutions were dried over anhydrous sodium sulphate (20 g), filtered and removed under vacuum to afford titled compounds (8a-c, 16 b-d).
Yield 84% as white crystals, m.p 305–307 °C (as reported).
General procedure for preparation of compounds (9a-c, 17b-d)
To a solution of beta-alanine ethyl ester hydrochloride (8.94 mmol), thieno[2,3-d]pyrimidine-6-carboxylic acid derivatives (8a-c, 16b-d) (8.94 mmol) and 1‐hydroxybenzotriazole (1.32 g, 9.18 mmol) in dry dimethylformamide (20 ml), N‐methyl morpholine (4.51 g, 44.65 mmol) was added at 0 °C and stirred for 30 min. Then EDCI.HCl (1.8 g, 13.12 mmol) was added and stirred for another 30 min maintaining the same temperature. The reaction mixture was stirred for 8 h at room temperature and then followed by the addition of water (20 ml) and ethyl acetate (2 × 20 ml) for extraction. Organic extracts were combined, washed with water (2 × 10 ml) and aqueous sodium bicarbonate solution (2 × 10 ml). Then organic solution was dried with anhydrous Na2SO4, and evaporated under vacuum to afford the titled compounds (9a-c, 17 b-d).
Yield 92% as buff solid, m.p 184 °C. 1.30–1.34 (t, 3H, OCH2CH3, J = 8.8 Hz), 2.67–2.71 (t, 2H, CH2CO, J = 8 Hz), 3.04 (s, 3H, CH3), 3.73–3.77 (q, 2H, CH2NH, J = 8 Hz), 4.19–4.25 (q, 2H, OCH2CH3, J = 8.8 Hz), 6.73 ((s, 1H, NH, exchangeable by D2O), 7.19–7.23 (t, 1H, aromatic H5′, J = 8 Hz), 7.51 (s, 1H, aromatic H2′), 7.87–7.89 (d, 1H, aromatic H6′, J = 8 Hz), 8.60 (s, 1H, pyrimidine H). MS: 437 (M++1, 18%), 436 (M+, 26%), 285 (100%). Anal. Calcd. For CC 52.23, H 4.15, N 12.82, S 7.34; Found: C 52.40, H 4.28, N 13.09, S 7.46. FT-IR (ú max, cm: 3444, 3367 (2 × NH stretch), 3120 (CH aromatic), 2981–2858 (CH aliphatic), 1728 (C = O ester), 1639 (C = O amide), 1612 (NH bend).
General procedure for preparation of compounds (10a-c, 18 b-d)
To a solution of ethyl 4-amino butyrate hydrochloride (8.94 mmol), prepared thieno[2,3-d]pyrimidine-6-carboxylic acid derivatives (8a-c, 16 b-d) (8.94 mmol) and 1‐hydroxybenzotriazole (1.32 g, 9.18 mmol) in dry dimethylformamide (20 ml), N‐methyl morpholine (4.51 g, 44.65 mmol) was added at 0 °C and stirred for 30 min. Then, EDCI.HCl (1.8 g, 13.12 mmol) was added and stirred for another 30 min maintaining the same temperature. The reaction mixture was stirred for 8 h at room temperature and then followed by the addition of water (20 ml) and ethyl acetate (2 × 20 ml) for extraction. Organic extracts were combined, washed with water (2 × 10 ml) and aqueous sodium bicarbonate solution (2 × 10 ml). Then the organic solution was dried with anhydrous Na2SO4, and evaporated under vacuum to afford the titled compounds (10a-c, 18 b-d).
Yield 90% as faint yellow solid, m.p 131–133 °C. 1.18–1.21 (t, 3H, OCH2CH3, J = 8 Hz), 1.77–1.84 (m, 2H, CH2CH2CH2), 2.36–2.40 (t, 2H, CH2CO, J = 8 Hz), 2.87 (s, 3H, CH3), 3.26–3.32 (q, 2H, CH2NH, J = 8 Hz), 4.05–4.10 (q, 2H, OCH2CH3, J = 8 Hz), 7.43–7.47 (t, 1H, aromatic H5′, J = 8 Hz), 7.63–7.67 (m, 1H, aromatic H2′), 7.92–7.94 (dd, 1H, aromatic H6′, J = 2.8, 6.8 Hz) , 8.48 (s, 1H, NH, exchangeable by D2O), 8.50 (s, 1H, pyrimidine H), 8.60 (s, 1H, NH, exchangeable by D2O). MS: 450 (M+, 30%), 434 (M++2, 35%), 110 (100%), 66 (85%). Anal. Calcd. For C: C 53.27, H 4.47, N 12.43, S 7.11; Found: C 53.49, H 4.61, N 12.69, S 7.17.
General procedure for preparation of compounds (11a-c, 19 b-d, 12a-c and 20b-d)
A cooled solution of sodium (27 mmol; 10 equiv.) in absolute ethanol (20 ml) was added to a solution of hydroxylamine hydrochloride (54 mmol; 20 equiv.) in absolute ethanol cooled in an ice bath. The mixture was stirred for 15 min. The formed precipitate was filtered, and the free hydroxylamine solution was prepared. The freshly prepared hydroxylamine solution was placed in a rounded bottom flask and cooled in an ice bath. Compounds (9a-c, 17 b-d, 10a-c, 18 b-d) (2.7 mmol; 1 equiv.) were added to the solution and stirred overnight. The reaction is monitored by TLC till the complete disappearance of starting compounds using an appropriate eluting system [DCM: MeOH (9.5:0.5)]. The solvent was removed under vacuum to afford crude products, washed with water, and crystallised from ethanol to yield target compounds (11a-c, 19 b-d, 12a-c, 20b-d),.
Yield 80% as brown solid, m.p 244–246 °C. 1.71–1.76 (m, 2H, CH2CH2CH2), 2.09–2.13 (t, 2H, CH2CO, J = 8 Hz), 2.87 (s, 3H, CH3), 3.22–3.28 (q, 2H, CH2NH, J = 8 Hz), 7.40–7.44 (t, 1H, aromatic H5′, J = 8 Hz), 7.63–7.65 (m, 1H, aromatic H2′), 7.91–7.93 (d, 1H, aromatic H6′, J = 8 Hz) , 8.47 (s, 1H, pyrimidine H), 9.43 (s, 1H, NH, exchangeable by D2O). MS: 439 (M++2, 26%), 437 (M+, 100%), 436 (M+-1, 84%). Anal. Calcd. For CC 49.37, H 3.91, N 15.99, S 7.32; Found: C 49.60, H 4.04, N 16.12, S 7.30.
In vitro EGFR and VEGFR-2 tyrosine kinase inhibitory activity
In vitro enzyme inhibition determination for the synthesised compounds was carried out in Thermo Fischer Scientific, Madison, WI, USA (SelectScreenServices@thermofisher.com). The EGFR and VEGFR-2 tyrosine kinase % inhibitory activity at a single-dose concentration of 10 µM was performed initially then representative compounds were evaluated at single-dose concentrations of (0.1 nM, 1 nM, 10 nM, 100 nM, 1 µM) using ZYLTE technology. A fluorescence-based, coupled-enzyme format is employed in the biochemical assay and depends on the sensitivity of proteolytic cleavage to non-phosphorylated and phosphorylated peptides. Then IC50 values for selected compounds were determined by GraphPad Prism at Misr international university.
In vitro HDAC inhibitory activity
Materials and methods
The HDAC inhibitory activity was carried out at BPS Bioscience (www.bpsbioscience.com). All of the synthesised compounds were tested at 10 μM for their inhibitory percent against HDAC6. HDAC6 (BPS#50006), HDAC substrate 3 (BPS#50037), HDAC Assay Buffer (BPS catalogue number 50031) and HDAC Assay Developer (BPS catalogue number 50030) was used in testing. TSN is also used as a reference compound purchased from Selleckcom.com (#S1045).
HDAC6 inhibitory activity assay protocols
All compounds are dissolved in DMSO. Serial dilutions of the compounds were carried out with 10% DMSO in HDAC assay buffer and 5 µl of the dilution was mixed to a 50 µl reaction so that the final concentration of DMSO is 1% in all of the reactions. Duplicate enzymatic reactions were conducted at 37 °C in a 50 µl mixture containing HDAC assay buffer, 5 µg BSA, HDAC substrate 3, HDAC6 (catalogue# 50006), and a test compound (10 µM) for 30 to 60 min. 50 μl of 2 × HDAC Developer was added to each well after enzymatic reaction, and the plate was incubated at room temperature for another 15 min. Fluorescence intensity was measured using a Tecan Infinite M1000 microplate reader at an excitation of 360 nm and an emission of 460 nm
Data analysis
Duplicates of HDAC activity assays were performed at each concentration. Analysis of the fluorescent intensity data was performed using the Graphpad Prism computer software., the % activity was calculated according to the equation: % activity = (F-Fb)/(Ft-Fb), where F = the fluorescent intensity in each data set in the presence of the compound, Ft = the fluorescent intensity in each data set, defined as 100% activity in the absence of the compound, Fb = the fluorescent intensity in each data set, defined as 0% activity in the absence of HDAC.
In vitro anticancer activity
All of the synthesised compounds (20) were submitted to NCI and selected for in vitro anticancer assay against a panel of 60 cancer cell lines. The methodology of assay and calculations is as reported [54,55].Click here for additional data file.
Authors: Mahender B Dewal; Amit S Wani; Celine Vidaillac; David Oupický; Michael J Rybak; Steven M Firestine Journal: Eur J Med Chem Date: 2012-02-25 Impact factor: 6.514
Authors: M C Alley; D A Scudiero; A Monks; M L Hursey; M J Czerwinski; D L Fine; B J Abbott; J G Mayo; R H Shoemaker; M R Boyd Journal: Cancer Res Date: 1988-02-01 Impact factor: 12.701
Authors: Stephen R Wedge; Donald J Ogilvie; Michael Dukes; Jane Kendrew; Rosemary Chester; Janet A Jackson; Sarah J Boffey; Paula J Valentine; Jon O Curwen; Helen L Musgrove; George A Graham; Gareth D Hughes; Andrew P Thomas; Elaine S E Stokes; Brenda Curry; Graham H P Richmond; Peter F Wadsworth; Alison L Bigley; Laurent F Hennequin Journal: Cancer Res Date: 2002-08-15 Impact factor: 12.701