ACS Medicinal Chemistry Letters
Letter
structure with those reported in previous literature gave further
support to the reliability of the computational docking
protocol. The latter was then applied also to thiophene and
pyrazolo-pyrimidine hit compounds in the attempt to predict
their binding mode and significant interactions with GLI1.
The thiophene derivative SST0682 was accommodated in
such a way that the phenyl ring at the amide portion showed a
cation−π interaction with the side chain of Arg223 (Figure 5),
thus accounting for the importance of a phenyl ring at this
position. Moreover, a similar interaction occurred between the
phenyl ring on the thiophene cycle (corresponding to the
aromatic ring feature of the pharmacophore) and the side
chain of Lys209. Further stabilization came by a hydrogen
bond between the ester side chain (that corresponds to an
HBA feature of the model) and the His220 ring. The
trichloromethyl group was embedded in a pocket comprised
within the side chains of Thr224 and Thr243, where also the
gem-dimethyl group of SST0704 was located (Figure 6).
The same binding pocket occupied by both GlaB and
SST0682 also accommodated the pyrazolo-pyrimidine deriv-
ative SST0704 (Figure 6). In particular, the pyrazole ring gave
a cation−π interaction with the terminal ammonium group of
Lys209, and the amide oxygen of the ligand (that corresponds
to one of the HBA features of the pharmacophore) interacted
by hydrogen bond with the His220 heterocycle. Hydrophobic
contacts involving both the chromene and the cyclopropyl
rings further stabilized the complex. In particular, the latter
pointed toward Cys207 that constituted the metal binding
pocket for zinc ion of the zinc finger 4.
To further confirm the effects of these compounds on the
GLI transcription factors, two pyrazolo-pyrimidines and two
thiophene derivatives were selected to perform a gene
expression analysis after treatment for 48 h at 5 μM.
Quantitative real-time PCR (qPCR) showed that all four
compounds are able to reduce the expression of GLI1, PTCH1,
HIP1, GLI2, and Cyclin D1 mRNA (Figure 7). Similar results
were obtained using GlaB. These results reinforce the effects of
our compounds as GLI inhibitors.
To confirm the specificity of our compounds for GLI
transcription factors, the effects of SST0683 and SST0790 on
melanoma cell viability in GLI1- or GLI2-silenced A375 cells
were also tested. This assay showed that depletion of GLI1 or
GLI2 greatly attenuated the effects of our putative GLI
inhibitors on cell viability, suggesting the specificity of our
compounds for GLI1 and GLI2 (Figure 8). Finally, to exclude
nonspecific cytotoxic effects, we demonstrated that the four
selected compounds had no effect on cell viability in a non-
neoplastic mammary epithelial cell line (MCF10A, Supporting
Figure 6. Best-ranked docking pose of SST0704 within GLI1 protein.
Important interaction sites are represented by the pyrazole ring that
gave a π−cation interaction with Lys209, the carbonyl group that is at
hydrogen bond distance from His220 of the zinc binding pocket, the
cyclopropyl ring that points toward the zinc binding pocket, and the
gem-dimethyl moiety that is allocated between the side chains of
Arg223 and Thr243. The structure of SST0704 is represented by
thick sticks, amino acid side chains in thin sticks, the zinc ion as a blue
sphere, hydrogen bonds as dotted yellow lines, and π−cation
interaction as dotted green lines.
associated to an almost complete inhibition of GLI1 protein
expression (Figure 3).
The hydrophobic character of pyrazolo-pyrimidines ap-
peared to be a crucial parameter that affected compound
activity. LogP values (calculated by MarvinSketch software,
inactive compounds, while activity occurred in compounds
with higher values, with the only exception SST0896 that was
inactive with a logP value of 4.84.
The ability of such compounds to affect GLI1 protein level
and inhibit luciferase assay suggested that they might interact
with GLI1 itself. To support this hypothesis, molecular
docking simulations (Glide software)26 were performed on
the X-ray structure of the five-finger GLI1/DNA complex (2gli
within the protein data bank).27 Best scored poses of GlaB
showed that it accommodates in the region between zinc finger
4 and 5. Moreover, Lys209, Lys219, and Arg223 (we adopt
amino acid numbering of the protein data bank even if a list to
convert it into the full length numbering is reported in
described.18 The agreement between the binding pose of
docked GlaB and its major interactions within the GLI1
Figure 7. Gene expression analysis of Hedgehog pathway components upon treatment with compounds and GlaB. Quantitative PCR (qPCR)
analysis of GLI1, PTCH1, HIP1, GLI2, and Cyclin D1 in A375 cells treated with DMSO, SST0682, SST0683, SST0704, and SST0790 or GlaB (5
μM) for 48 h. The y axis represents expression ratio of gene/(HPRT and EIF2α average), with DMSO control equated to 1. Data shown as mean
SD.
E
ACS Med. Chem. Lett. XXXX, XXX, XXX−XXX