3448
M. Nettekoven et al. / Bioorg. Med. Chem. Lett. 15 (2005) 3446–3449
Figure 1. Conformations corresponding to the minimum energies for the thiazole isomers 1 and 2. The arrows represent the torsional angles sampled
in both isomers.
compounds the general SAR trend is nicely
illustrated.
carbonyl oxygen. The potential energy difference be-
tween the bent minimum and the extended conforma-
tion of thiazole 2 is 10 kJ/mol with the CDIE model
and 7 kJ/mol with the GB/SA model, which is in agree-
ment with the 10-fold IC50 shift, corresponding to a free
energy difference of 5.7 kJ/mol. Of course, the change in
potential energies is only a fraction of the observed free
energy difference, which also includes the interaction en-
ergy and entropy. Nevertheless, this comparison clearly
shows that an extended conformation is preferred for
the binding to the NPY5 receptor. The results of this
analysis will be further pursued and integrated in the de-
sign of new NPY5 receptor ligands.
In general, thiazoles 1 show low nanomolar binding
affinities for the mouse receptor. The nature of the sub-
stituents influenced the inhibition potential. Electron
withdrawing and pushing (e.g., 1a–1c) substituents show
similar IC50 values and also the positioning of the sub-
stituents seemed to have a minor influence. The combi-
nation of two substituents led to thiazoles (e.g., 1d–1f)
with IC50 values in the picomolar range. In general, het-
eroaromatic substitution is allowed as well (e.g., 1g).
The isomeric counterparts 2 showed, in general, at least
a 10-fold lower binding affinity, however, still producing
compounds in the nanomolar range. The SAR trends
for thiazoles 1 and 2 are very similar (e.g., 2a–2g).
In conclusion, sets of isomeric thiazole derivatives 1 and
2 have been synthesised in a parallel iterative solution-
phase synthesis approach guided by the SAR analysis
derived from biological results and computer-aided de-
sign and analysis. This synergistic and streamlined
working procedure led to highly active isomeric ligands
for the NPY5 receptor. However, a 10-fold difference at
least in their respective binding affinities was consis-
tently found for all isomeric pairs 1 and 2. The analysis
of conformational differences due to heteroatom interac-
tions in 1 and 2 revealed a favourable C@OÁ Á ÁS interac-
tion in 1, whereas thiazoles 2 showed a repulsive
C@OÁ Á ÁN interaction. In the design of new NPY5 recep-
tor ligands, this result will be incorporated in the quest
towards new compounds with preferable biological
and pharmacokinetic profiles. The variation of the bis-
amino linker moiety needs to be further explored. This
is currently under investigation in our laboratories and
will be reported in due course.
To understand the origin of the 10-fold shift of the bind-
ing affinity between thiazoles 1 and 2, a conformational
analysis was performed using MacroModel 8.6.9 For
this analysis, thiazole isomers have been reduced to their
core scaffolds since the shift is consistent over the whole
sulfonamide SAR (see Fig. 1). Three torsional angles
have been sampled with 5000 steps of Monte Carlo
search (MCMM), followed by 1000 steps of energy
minimization (PRCG). The solvent effect has been mod-
elled either with a dielectric constant (CDIE) of 4 or
with a Generalized Born/Surface Area (GB/SA) water
model. The results are similar in both models.
Figure 1 displays the energy minima for the thiazole iso-
mers 1 and 2. They both have in common a rotation of
the phenyl ring with respect to the carbonyl function of
58° induced by the methyl substituent, but the two
structures differ in the torsional angle between the car-
bonyl function and the heterocycle. On the one hand,
the more potent thiazole 1 series displayed a favourable
interaction between the carbonyl oxygen and the sulfur
stabilising an extended conformation. This type of
C@OÁ Á ÁS interaction, with a distance significantly below
Acknowledgments
It is with real pleasure that we wish to thank all our col-
laborators whose contributions made the described
work possible and so enjoyable, especially Dr. F. Daut-
zenberg, Dr. A. Bourson-Sleight, Dr. M. Klug and Dr.
M. Stahl.
˚
the sum of van der Waals radii (<3.3 A), is highly
favourable as explored both statistically and theoreti-
cally by Iwaoka et al.10 This is further corroborated
by more than 100 molecules with a similar atom
arrangement, which can be found in the Cambridge
Structural Database.11 On the other hand, the thiazole
2 series favours a bent conformation due to the repulsive
interaction between the thiazole nitrogen and the
References and notes
1. Tatemoto, K.; Carlquist, M.; Mutt, V. Nature 1982, 296,
659.