Y.-P. Chen et al. / Bioorg. Med. Chem. Lett. 24 (2014) 349–352
351
Figure 2. The representative conformational snapshots of compounds 5, 20 and 21 in the catalytic domain of SIRT1 produced by 0.5
ls explicit solvent MD simulations.
activities. Next, the activity data for Group C including the weaker
activity for compound 12 with methyl ester than that of compound
5 and the lack of activity for compounds 13–16 indicated that the
carboxylic acid appeared to be required for inhibition activity. It is
interesting that sulfonic acid 14 was inactive despite the similarity
between the functional groups. In Group D, compound 17, in which
the dimethylamino group was replaced by a methyl group, showed
diminished activity whereas compound 18, in which the dimethyl-
amino group was replaced by a similarly sized and shaped isopro-
pyl group, showed similar inhibitory activity to that of compound
5. The results indicated that the size of the substituent is very
important for inhibitory activity. Therefore, we evaluated the two
derivatives in Group E, compounds 19 and 20 with trimethylam-
monium and tert-butyl groups as respective substituents.
depicted in Figure 2. The more favorable binding energy for com-
pound 21 can be explained by its tighter conformational ensemble
observed for this ligand in the binding domain than those of com-
pounds 5 and 20. For compound 5, MD simulations produced two
distinct conformational families, as clearly seen by the position of
the carbonyl group. One conformational family was similar to
those of compounds 20 and 21, in which the polar functional
groups positioned towards the pocket entrance. In the second con-
formational family, compound 5 rotated by approximately 90 de-
grees. This rotation placed the carbonyl group adjacent to the
nicotinamide group of the NAD+, resulting in slight displacement
of the NAD+ from the active site and thus creating a possible escape
route from the binding pocket. This behavior was observed after
first 100 ns of simulation, and could also support weaker binding
of 5 than 20. This conformational flexibility resulted in the less
favorable binding energy of this ligand as compared to 20 and
21. The smaller size of compound 20 (as compared to 21) allowed
it more freedom to move inside the binding pocket than compound
21. Based on this conformational analysis, it is possible to propose
that increase in conformational flexibility of the compound in the
SIRT1 binding pocket leads to decrease in binding affinity. The re-
sults obtained from the MD calculations suggest that bulkier aro-
matic acid derivatives may be better SIRT1 inhibitors by
complementing the size and hydrophobic environment of the en-
zyme binding pocket, and we plan to test such derivatives in the
near future.
Compound 20 showed the best activity (MIC 50 lM) among the
derivatives tested, which supported the importance of the size of
the substituent as suggested above. The lack of activity seen
in compound 19 suggested that bulky but not charged groups
are required for enhanced inhibitory activity for benzoic acid
derivatives.
In light of their Sir2p inhibitory activities, compounds 5 and 20
were further evaluated for SIRTs inhibitory activities. Compounds 5
inhibited SIRT1 and SIRT2 with 25.3% and 30.3% at 1.6 mM whereas
compound 20 inhibited SIRT1 and SIRT2 with 54.8% and 28.0% at
1.6 mM, respectively. These results suggested that compound 20
was a weak but selective SIRT1 inhibitor (IC50 1.0 mM). Further-
more, the SIRT1 inhibitory activity was enhanced twofold by
replacing the dimethylamino group with tert-butyl group, which
supported the SAR pattern described above.
The crystal structure of the SIRT1 catalytic domain with NAD+
and an EX-527 analogue [(S)-2-chloro-5,6,7,8,9,10-hexahydrocy-
clohepta[b]indole-6-carboxamide (21)] has recently been
reported.21 Based on the SIRT1/NAD+/21 coordinates, we constructed
SIRT1/NAD+/5 and SIRT1/NAD+/20 complexes to evaluate the
In this study, 4-tert-butylbenzoic acid (20) was discovered as a new
Sir2p inhibitor based on the SAR study on the naturally occurring
weak Sir2p inhibitor, 4-dimethylaminobenzoic acid (5), isolated
from the Streptomyces sp. CP27-53. Compound 20 also showed a
weak but selective inhibitory activity against SIRT1. It is quite inter-
esting that the structure of 5 was identical to the capping group of
the potent class I/II HDAC inhibitor trichostatin A.29 This study also
demonstrated a reasonable correlation between the calculated
binding energy and potency of SIRT1 inhibition activity, suggesting
that it would be possible to establish a SIRT1 virtual screening
method by collecting more data points. The SAR study and MD
binding free energies (
DGbind) of compounds 5 and 20 relative to
compound 21. The Gbind energies were calculated by using the
D
MM-GBSA approach for post processing of explicit solvent molec-
ular dynamics (MD) trajectories. This approach has been proven
to provide absolute free energies in good agreement with experi-
mental data.22–28 The potent and selective SIRT1 inhibitor 21
(IC50 60–100 nM17) showed a binding energy of À37.5 kcal/mol
whereas compounds 5 and 20 were found to have binding energies
of À26.8 and À31.6 kcal/mol, respectively. Conformational ensem-
bles for compounds 5, 20 and 21 complexes produced by MD are