V. Liao et al. / Bioorg. Med. Chem. Lett. 22 (2012) 6200–6204
6201
acid, 2-phenylacetic acid or 2-phenoxyacetic acid; and b-alanine
ethyl ester hydrochloride. Following hydrolysis of the ethyl ester,
the hydroxamic acid group was installed from the reaction between
the mixed anhydride formed in situ from chloroethylformate and
N-methylmorpholine; and hydroxylamine (Scheme 1). After work-
up, when sprayed with an ethanol solution of Fe(III), the final spot
of the product on a TLC plate turned pale pink, as evidence of the
metal-coordinating ability of the hydroxamic acid group. The prod-
ucts (Table 1) were characterised by 1H and 13C NMR spectroscopy,
ESI-MS and microanalysis (Supplementary data). Relative to the free
carboxylic acid, the signals in the 1H NMR spectra assigned to the
coordinate bonds between the hydroxamic acid group and the
Zn(II) ion were re-installed. This method assumes that the biden-
tate hydroxamate coordination between 2 and 7 and Zn(II) will
not differ significantly from the 1-Zn(II) coordination observed in
the X-ray crystal structure. The method allows for the structure
of 2–7 distal to the active site Zn(II) to be optimized within the
constraints of the HDAC8 active site cavity. The veracity of the
method was supported by the close agreement between the con-
formation of native 1 and 1 which was minimized within the
HDAC8 active site as described for 2–7 (Fig. 2, far left). The overlaid
structures of minimized 1 and each of 2–7 are shown in Figure 2.
Despite the variation in the number of main-chain atoms in the
linker separating the hydroxamic acid head group from the
phenyl-substituted tail (six atoms: 2–4, 7; five atoms: 5, 6), the
phenyl group in most cases minimized to a position that coalesced
closely with the phenyl group of 1. The poorest overlay of phenyl
groups occurred between 1 and 4 (Fig. 2, middle). In this case,
the ethylene group proximal to the tolyl group was sufficiently
flexible to mine a local hydrophobic pocket in the protein, which
re-directed the orientation of the tolyl group. This local conforma-
tional flexibility was not available to other candidates, due to
unsaturation (2, 7), an insufficient number of methylene groups
(5, 6) or the presence of a methoxy group subject to electron
delocalization effects (3).
The linker region of 1 is more complex than 2–7, with two
methyl substituents, one chiral carbon (6R) and two unsaturated
carbon–carbon bonds (2E, 4E). In the X-ray crystal structure of 1
bound to HDAC8, the 40-methyl group is positioned in a sand-
wich-like fashion between the aromatic rings of F152 and F208
with the distance from each plane of the aromatic ring to the
methyl group carbon atom equal at 3.84 Å (Fig. 3). This hydropho-
bic interaction was not available to 2–7 (Fig. 3, shown for 2) and
may be a factor that contributed to their low potency as HDAC1
inhibitors. The structures of 2–7 as minimized in the HDAC8 ac-
tive site cavity revealed a close distance between the 50-carbonyl
oxygen atom and the sulfur atom of Met274. In 2 and 7, this dis-
tance was 3.5 Å, which is close to the sum of the sulfur and oxy-
gen van der Waals radii (3.3 Å)22 and raises the possibility of the
carbonyl oxygen atom-mediated oxidation of the Met274 sulfur
atom in HDAC8, as shown to occur in amyloid b-peptide.23,24
Inhibitor potency was determined against HDAC1, which instead
of the four amino acid sequence P273-M274-C275-S276 in
HDAC8, contains R273-L274-G275-C276.25 In silico mutation of
PMCS(HDAC8) ? RLGC(HDAC1) followed by minimization of 2–7
and RLGC only, with the remaining protein constrained, showed
that the 50-carboxyl oxygen atom of 2 was in close proximity to
the isobutyl side chain of Leu274, with one oxygen–hydrogen
(methyl) distance of 2.8 Å. It is proposed that the proximity of
the 50-carbonyl oxygen atom in 2-7 to Leu274 in HDAC1 is a con-
tributing factor to the low potency of this group of compounds.
Reduced binding to the active site would arise from steric repul-
sion and from repulsive dipole–dipole interactions, since the
a-methylene protons shifted upfield by about 0.2 ppm in the
corresponding hydroxamic acid.
The IC50 values of 2–7 against HDAC1 were determined using a
fluorometric assay and the anti-proliferative effects of the com-
pounds were determined against BE(2)-C neuroblastoma cells
using the Alamar blue assay (Fig. 1, Table 1).13–15 Compared to
the IC50 value of 1 against HDAC1 (IC50 12 nM), 2–7 had signifi-
cantly decreased potency (38–84
a maximum of 40% HDAC1 activity at 250
l
M). Compound 2 inhibited only
M, and was the least
l
potent compound. The IC50 value of SAHA against HDAC1 deter-
mined in the current work (110.8 1.1 nM) calibrated the results
to literature IC50 values for SAHA;16 and supported the veracity
of the data for 2–7. Given the similarity between 2–7 and 1, with
respect to the Zn(II) coordinating hydroxamic acid group; and be-
tween 2, 5 and 1, with respect to the 4-(dimethylamino)phenyl
group, the low potency against HDAC1 was unexpected. Most sur-
prising was the extraordinarily low potency of 2 against HDAC1.
The whole cell-based assay system gave an activity profile that dif-
fered from the profile established using isolated HDAC1. At 50 lM
the anti-proliferative activity of 2–7 against BE(2)-C neuroblas-
toma cells ranged between 57.0 and 88.6%, with 2 the most potent
compound. Other work has shown differences between the activi-
ties of compounds screened using isolated HDACs or whole-cell
assay systems, with belinostat and SAHA a notable example
(belinostat:SAHA; biochemical activity (ꢀ1:1), anti-proliferative
activity against a selection of cancer cell lines (ꢀ1.7–7:1)).17 This
is likely due to the higher complexity of the whole cell-based assay
system with regard to the involvement of HDACs in regulating the
acetylation of non-histone proteins (a-tubulin, tumor suppressor
p53) and the optimal activity of most HDACs requiring the forma-
tion of multiprotein complexes. The IC70 values of 1 or SAHA were
determined in BE(2)-C neuroblastoma cells as 0.05
lM or 2 lM,
respectively. Compound 7 has been reported to inhibit 19% of
HDAC activity in HeLa extract at 1
oped as a lead.18
lM, but was not further devel-
Molecular modelling (HyperChem 7.5, MM+ force field for en-
ergy minimization) was used to inform the IC50 results of 2–7
against HDAC1. Three X-ray crystal structures of human HDACs
with 119,20 and one structure of a prokaryotic HDAC homologue
with 1 have been solved.21 The X-ray crystal structure of 1 bound
to human HDAC819 was used as the template. The coordinate
bonds between 1 and Zn(II) were deleted and the isolated fragment
was used to build models of each of 2–7 with zero charge. The pro-
tein and the hydroxamic acid group (C(O)NH(OH)) of the model of
2–7 were frozen and the remaining structure of the free ligand was
minimized. At the completion of the minimization process, the
molecular dipole in leucine is orientated from the
a-carbon atom
(positive) along the isobutyl tail (negative). Due to the presence of
the electron donating dimethylamino substituent (rp –0.83) in 2
and the conjugated linker region, this compound could comprise
different contributing resonance structures (Fig. 4a). In the
Scheme 1. Synthesis of amide-linked derivatives of b-alanine hydroxamic acid: reagents and conditions (shown for 2). (i) HOBt, EDCꢁHCl, DIPEA, DMF, RT, 16 h; (ii) NaOH,
MeOH, RT, 30 min; (iii) EtOCOCl, NMM, NH2OH, MeOH, RT, 15 min.