We have explored the isoelectronic replacement of the C═C double bond found at the core of many nonsteroidal estrogen ligands with a simple Schiff base (C═N). Di- and triaryl-substituted imine derivatives were conveniently prepared by the condensation of benzophenones with various anilines without the need for phenolic hydroxy protection. Most of these imines demonstrated high affinity for the estrogen receptors, which, in some cases exceeded that of estradiol. In cell-based assays, these imines profiled as ERα agonists but as ERβ antagonists, showing preferential reliance on the N-terminal activation function (AF1), which is more active in ERα. X-ray analysis revealed that the triaryl-imines distort the ligand-binding pocket in a new way: by controlling the separation of helices 3 and 11, which appears to alter the C-terminal AF2 surface that binds transcriptional coactivators. This work suggests that C═N for C═C substitution might be more widely considered as a general strategy for preparing drug analogues.
We have explored the isoelectronic replacement of the C═C double bond found at the core of many nonsteroidal estrogen lin>an class="Chemical">gands with a simple Schiff base (C═N). Di- and triaryl-substituted imine derivatives were conveniently prepared by the condensation of benzophenones with various anilines without the need for phenolic hydroxy protection. Most of these imines demonstrated high affinity for the estrogen receptors, which, in some cases exceeded that of estradiol. In cell-based assays, these imines profiled as ERα agonists but as ERβ antagonists, showing preferential reliance on the N-terminal activation function (AF1), which is more active in ERα. X-ray analysis revealed that the triaryl-imines distort the ligand-binding pocket in a new way: by controlling the separation of helices 3 and 11, which appears to alter the C-terminal AF2 surface that binds transcriptional coactivators. This work suggests that C═N for C═C substitution might be more widely considered as a general strategy for preparing drug analogues.
Thne">ere are many motivations
for preparing analogues of pharmaceuticals
(e.g., impn>roving drug propn>an class="Gene">erties, reducing drug liabilities, seeking
unclaimed intellectual property space, and simplifying or improving
synthesis), and there are numerous approaches for the preparation
of such analogues (e.g., substituting heteroatoms and replacing peripheral
or structural core elements with sterically or electronically similar
entities). We have sought to expand the chemical diversity of ligands
for the estrogen receptor (ER) by replacing their internal scaffolding
with various heterocycles and other structurally related elements,
and in the process, we have obtained structurally novel compounds
that are generally easy to prepare.[1−7] Because estrogens, acting through their two receptors, ERα
and ERβ, regulate a wide range of physiological and pathological
processes and because various ER ligands can demonstrate marked tissue
selectivity (based on ER subtype selectivity[8,9] or
selective engagement of coregulator proteins, i.e., SERM selectivity[10,11]), it is not surprising that in some cases our structurally novel
estrogen analogues showed unusual patterns of estrogenic activity
and selectivity.[1,12]
Many nonsteroidal estrogens
are di- or triarylethylenes having
a C=C double bond core, and in a recent publication, we examined
the replacement of this C=C double bond with the isoelectronic
and isostructural B—N bond (Figure 1, right, middle).[13] To achieve hydrolytic
stability, the electrophilicity of the boron center had to be sterically
masked with a full array of flanking ortho methyl groups, and analogues
with p-OH groups on the B-phenyl groups were unstable.
Despite these restrictions, some of the anilino bis(2,6-dimethylphenyl)borane
derivatives that we prepared were very stable and demonstrated reasonably
high binding affinity and good cellular potency, being ERα agonists
and ERβ antagonists.[13] Those studies
defined the structural determinants of stability and cellular bioactivity
of a B—N for C=C substitution and provided a framework
for further exploration of “elemental isomerism” for
diversification of drug-like molecules.[13]
Figure 1
Triarylethylene
nonsteroidal estrogen as well as B—N and
C=N double bond for C=C bond isoelectronic replacements.
pan class="Chemical">Triarylethylenen>
nonpan class="Chemical">steroidal estrogen as well as B—N and
C=N double bond for C=C bond isoelectronic replacements.
Hne">ere, we have further explored
the chemical diversification of
ER ligands by another, simple isostructural and isoelectronic replacement,
substituting a core C=C double bond with a simple Schiff base
(C=N). As far as we are aware, there appears to be only limited
and rather distant precedent for use of a C=N bond as a core
element for estrogens.[14,15] Furthermore, on the basis of
a substructure search of the Merck Index, the C=N bond does
not appear to be well-recognized as a surrogate for a C=C bond
in the construction of bioactive molecules (except as a sub-element
in 2-aryl-benzimidazoles, benzodiazepines, and other related heterocycles).
As was the case earlier,[13] di- and triarylethylene
nonsteroidal estrogens offer a convenient C=C double bond core
(Figure 1, right, top) for replacement by the
isostructural C=N element (Figure 1,
bottom, right). These Schiff base or imine-core ligands are very easy
to synthesize, in most cases in one step (and much easier than their
C=C double bond analogues), and to this end, we have prepared
a series of triaryl-substituted (and some diaryl-substituted) Schiff
base derivatives by a simple condensation that proceeds without phenolic
hydroxy protection. Their binding affinities and cellular biological
activities showed distinctive structure–activity relationships
(SARs), and they profiled as potent ERα agonists and ERβ
antagonists.
Results
Synthesis
Representative
members having an imine-core structure were prepared in three classes:
4,4′-dihydroxybenzophenone and 2,4,4′-trihydroxybenzophenone
derivatives of various anilines and the corresponding diaryl imines.
On the basis of previous work,[16] we developed
a HCl(g)-catalyzed condensation reaction without phenolic hydroxy
protection for the synthesis of bisphenolic Schiff bases 2a–l in chlorobenzene with heating at 140–145
°C for 24 h (Scheme 1A). Although all
yields were moderate, this improved method was superior to the phenol-protection
strategy.[17] We were unable to prepare triphenolic
Schiff base 2n by the reaction of 4,4′-dihydroxybenzophenone 1c with 4-aminophenol or by demethylation of 2l. However, when 4,4′-dihydroxybenzophenone 1a was treated with 4-aminophenyl benzoate, it gave imine 2m in 74% yield, and deprotection with KOH proceeded at room temperature
to give desired triphenolic product 2n in 94% yield.
We also used this methodology for the synthesis of the corresponding
phenolic Schiff bases 4 (Scheme 1B). Benzophenone derivatives 1b and 1c were
prepared by a Friedel–Crafts acylation reaction and subsequent
condensation with various anilines then afforded imines 4a–m in good yield. Product 4n was
also prepared in 92% yield via a two-step procedure similar to that
used for 2n. In all cases (4a–l), the (E) form of the imines was the only diastereoisomer
formed, presumably because the intramolecular hydrogen bond, which
can only form in the (E) isomer, contributes to the imine stability.
These triaryl systems appeared to be hydrolytically stable indefinitely
under aqueous conditions. By contrast, diaryl imines prepared from
anilines and aldehydes or phenyl alkyl ketones displayed pronounced
hydrolytic lability, undergoing substantial hydrolysis in aqueous
MeOH at rt within 1–6 h. However, the diarylimines having an
internal hydrogen bond between the iminenitrogen and an ortho phenolic
OH group were considerably more stable.
Scheme 1
Synthesis of Triaryl-Substituted
Schiff Base Analogues 2 and 4 with the Improved
Method
Reagents and conditions: (a)
aniline derivatives (3 equiv), HCl(g), PhCl, 140–145 °C,
24 h; (b) KOH, MeOH, rt, 3 h; (c) resorcinol, AlCl3, sulfolane,
65–70 °C, 8 h.
Synthesis of Triaryl-Substituted
Schiff Base Analogues 2 and 4 with the Improved
Method
Reagents and conditions: (a)
aniline derivatives (3 equiv), HCl(g), PhCl, 140–145 °C,
24 h; (b) KOH, MeOH, rt, 3 h; (c) resorcinol, AlCl3, sulfolane,
65–70 °C, 8 h.
Estrogen Receptor Binding
Affinity Assays
The binding
affinity of Schiff base analogues 2 and 4–6 for both ERα and ERβ was determined
by a competitive radiometric assay, using methods that have been previously
described, and are reported in Table 1.[18] The affinities are represented here by relative
binding affinity (RBA) values, where estradiol has an affinity of
100 (absolute affinities for estradiol: Kd 0.2 nM on ERα and 0.5 nM on ERβ).
Table 1
Relative Binding Affinities (RBAs)
of Compounds 2 and 4–6 for Estrogen Receptor α and βa
Determined by a competitive radiometric
binding assay with [3H]estradiol; preparations of purified,
full-length human ERα and ERβ (Invitrogen) were used;
see the Experimental Section. Values are reported
as the mean ± the range or SD of two or more independent experiments;
the Kd for estradiol for ERα is
0.2 nM and for ERβ, 0.5 nM. Ki values
for the reported compounds can be readily calculated using the formula Ki = (Kd[estradiol]/RBA)100.
Detne">ermined by a competitive radiometric
binding assay with [n>an class="Chemical">3H]estradiol; preparations of purified,
full-length humanERα and ERβ (Invitrogen) were used;
see the Experimental Section. Values are reported
as the mean ± the range or SD of two or more independent experiments;
the Kd for estradiol for ERα is
0.2 nM and for ERβ, 0.5 nM. Ki values
for the reported compounds can be readily calculated using the formula Ki = (Kd[estradiol]/RBA)100.
As a global obsne">ervation, most
of the imines are gratifyingly high-affinity
ligands for both ERs, although the number and position of the hydroxyl
group in the phenyl ring of the benzophenone moiety as well as the
disposition and size of substituents on N-phenyl
group have a marked influence on their binding affinity and selectivity.
When assayed on the individual ER subtypes, ERα and ERβ,
most compounds show only a modest binding-affinity preference for
ERβ, at most 5-fold (compound 2d). Binding-affinity
comparisons suggest that one of the p-hydroxy groups
on the benzophenone moiety is playing the role of the phenolic hydroxyl
of estradiol, which is well-known to be the dominant functional group
ensuring high binding affinity.[19] (This
was confirmed by X-ray crystallography; see below.) Curiously, the
presence of a para hydroxyl group in the distal N-phenyl group, which is essential for the high-affinity binding of
a number of bisphenolic estrogens such as diethylstilbestrol and hexestrol
as well as certain phenylindenes,[20] actually
has a markedly detrimental effect on binding, reducing affinity by
about 7-fold (Table 1, compound 2a vs 2n and 4b vs 4n). Where
studied, the second para hydroxy group in the second ring of the benzophenone
unit has little effect on binding (Table 1,
compound 4a vs 4l), although the methyl
ether of this phenol is a good ERβ ligand (Table 1, 4o). Lastly, placement of a second hydroxyl
group at the ortho position of the phenyl ring of the benzophenone
moiety, which enables potential formation of an intramolecular hydrogen
bond as in the series 4 compounds, caused, in most cases,
a decrease in affinity for both ERα and ERβ; this was
somewhat unexpected but can be rationalized from conformational changes
noted in crystal structures (see Discussion and below). Parallel changes in the nature of the substituents on
the N-phenyl group were made within the two major
series, 2c–n and 4c–n, and in nearly every case, the binding affinities for both
ERα and ERβ changed in a coordinated fashion, although
with some bias for one or the otherER subtype. This is illustrated
graphically in Figure 2.
Figure 2
Graphical presentation
of RBA values for imines 2c–n (black,
H) and 4c–n (gray, OH)
Graphical n class="Chemical">presentation
of RBA values for pan class="Chemical">imines 2c–n (black,
H) and 4c–n (gray, OH)
In both sne">eries, the highest affinity analogues
were those having
somewhat bulky but nonpolar substituents, and where a substituent
(Me or Cl) was placed at the C-2, -3, and -4 positions, there was
a general decline in affinity with that progression of ring positions,
although this was more evident in the ERα series than in the
ERβ series. Among the various substituents, certain ones were
preferred at each of the three positions, but these differed between
the two series and with the two ER subtypes. Addition of a second
methyl group, as in 2,6-dimethyl compound 2b, resulted
in a dramatic, ca. 50-fold, reduction in affinity compared to 2-methyl
analogue 2c. Electron-withdrawing substituents (halogens
and CF3) were associated with higher affinity than alkyl
(methyl) or electron-donating substituents (OH or OMe). Overall, these
results illustrate that ER ligands having simple imine-core structures
can be readily prepared but that high-affinity binding, as one might
expect, requires an appropriate distribution of bulk, polarity, and
functionality.
Previously, in connection with studies of dimethylgallium
chelates
that mimic nonsteroidal estrogens, we prepared a number of smallerimine-core systems (Figure 3), and although
we published their synthesis and structures,[16] we had, until now, not determined their ER-binding affinities. More
recently, we also prepared a few additional diarylimines. In stark
contrast to the triaryl systems, all of these monoaryl and diarylimines have very low binding affinity (the data is given in Supporting InformationTable
S1). The highest affinity compounds are benzophenone imines
of alkyl or benzyl amines (9), which also have some ERβ-binding
preference. Even on ERβ, however, these compounds have RBA values
of less than 0.3%, and the affinities of the otherdiarylimines we
prepared are typically more than 10-fold lower (see Supporting InformationTable S1). As noted earlier, some of these diaryl systems were also rather
hydrolytically labile, which is not the case for all of the series 2 and 4 compounds.
Figure 3
Other mono and diaryl
imines (also see Supporting
InformationTable S1).
Othpan class="Gene">ern> mono and pan class="Chemical">diaryl
pan class="Chemical">imines (also see Supporting
InformationTable S1).
Transcriptional Activity
We detne">ermined
the effects
of these imines on ER transcriptional activity using an ER-responsive
luciferase reporter gene. Steroid-deprived HepG2 liver cells were
transfected with a widely used 3×ERE-luciferase reporter and
an ERα or ERβ expression plasmid for agonist activity
(% efficacy) and potency (EC50) determinations. These cells
were stimulated with increasing concentrations of 17β-estradiol
(E2) or compounds 2a, 2c–k, 2n, 4b–k, and 4n. For antagonist mode assays (% efficacy and IC50), cells were stimulated with a combination of estradiol (10 nM)
and an increasing concentration of the various compounds. The next
day, luciferase activity was measured. HepG2 cells are particularly
useful as a test system to distinguish the activity of estrogens through
the two activation functions of the ERs, the N-terminal AF1, in the
A/B domain, and AF2, which is in the ligand-binding domain.[21] SERMs, such as tamoxifen, show tissue-selective
agonist activity in some tissues, such as the uterus, and in the HepG2
cells via AF1 activity, but they work as antagonists in the breast
through structural mechanisms that are not understood.
In general,
these imines fully stimulated ER-mediated transcription in cells transfected
with wild-type ERα, indicating that they are potent and highly
efficacious ERα agonists (Table 2 and
Figure 4A). In most cases, the number and position
of the hydroxyl groups in the phenyl ring as well as the disposition
and size of substituents on the N-phenyl group had
no obvious effect on ERα-mediated transcription. However, compared
to E2, none of the imines fully stimulated ER-mediated transcription
in cells transfected with an ER construct that lacks AF1 because of
the deletion of the N-terminal AB domain (Figure 4B), indicating that these imines do not fully induce AF2-mediated
ER activity but also rely substantially on the AF1-mediated activity
of ER to drive transcription.
Table 2
Effects of Imines on the Transcriptional
Activities of Estrogen Receptor α and β
agonist
modea
antagonist modea
ERα
ERβ
ERα
ERβ
EC50 (nM)
eff (% E2)
EC50 (nM)
eff (% E2)
IC50 (nM)
eff (% E2)
IC50 (nM)
eff (% E2)
2a
3
110 ± 5
34 ± 12
126 ± 4
13
17 ± 3
2c
7
103 ± 6
8 ± 3
149 ± 6
6
27 ± 9
2d
5
109 ± 3
13 ± 9
93 ± 5
41
14 ± 0
2e
nd
nd
7 ± 9
nd
45
10 ± 2
2f
1
104 ± 3
37 ± 9
83 ± 9
19
30 ± 3
2g
10
102 ± 3
10 ± 3
109 ± 5
16
11 ± 2
2h
2
114 ± 4
9 ± 3
138 ± 8
9
6 ± 2
2i
23
89 ± 3
3 ± 2
104 ± 6
24
11 ± 5
2j
2
114 ± 4
10 ± 3
113 ± 4
7
8 ± 2
2k
1
105 ± 2
11 ± 1
128 ± 10
6
4 ± 3
2n
14
95 ± 2
332
16 ± 8
111 ± 4
300
33 ± 9
4b
7
102 ± 2
628
25 ± 8
134 ± 34
6
23 ± 12
4c
111 ± 4
2
31 ± 3
100 ± 5
16
29 ± 4
4d
14
95 ± 2
9 ± 2
101 ± 2
51
11 ± 3
4e
16
103 ± 3
10 860
20 ± 9
122 ± 34
241
21 ± 9
4f
0.1
113 ± 4
54 ± 1
136 ± 13
6
33 ± 6
4g
15
117 ± 6
15 ± 2
118 ± 23
67
15 ± 0
4h
4
114 ± 3
385
19 ± 5
123 ± 18
180
15 ± 1
4i
57
100 ± 4
11 ± 6
95 ± 9
28
16 ± 2
4j
10
108 ± 3
23 ± 8
124 ± 6
139
14 ± 2
4k
9
113 ± 4
12 460
24 ± 2
111 ± 5
100
13 ± 6
4l
17
102 ± 3
276
38 ± 13
102 ± 6
20 ± 24
4n
7
97 ± 3
51 ± 10
108 ± 4
417
25 ± 10
In the agonist
mode, ERE-luciferase
assays were performed with 12-point dose curves of the indicated ligands,
whereas in the antagonist mode, this was done in the presence of 10
nM E2.
Figure 4
Graphical comparison of the agonist-mode efficacies
of matched
compounds on full-length ERα (A), ERα with the N-terminal
AF1 deleted (B), and ERβ (C)
Graphical comparison of the agonist-mode efficacies
of matched
compounds on full-length ERα (A), n>an class="Gene">ERα with the N-terminal
AF1 deleted (B), and ERβ (C)
In the agonist
mode, pan class="Gene">ERE-lucifpan class="Gene">erase
assays were performed with 12-point dose curves of the indicated ligands,
whereas in the antagonist mode, this was done in the presence of 10
nM E2.
ne">ERβ exhibits
negligible AF1-mediated activity compared to
ERα, but the AF2-mediated activities of the ER subtypes are
similar in reporter assays,[22] suggesting
that despite their similar binding affinities for the two ER subtypes
(Figure 2) these imines would act as partial
ERβ agonists. Consistent with this supposition, all of the Schiff
bases tested failed to stimulate ER-mediated transcription fully in
HepG2 cells transfected with ERβ (Table 2 and Figure 4C). In fact, many of these compounds
acted as potent and nearly complete ERβ-selective antagonists,
suppressing the E2-induced activity of ERβ in the low nanomolar
range (Table 2) and thereby underscoring the
AF1-dependent ER subtype-selective properties of these compounds.
The N-phenyl substituents appear to affect the
potency of these compounds as ERβ antagonists; howevn>an class="Gene">er, whether
some of the obvious differences in IC50 reflect variation
in RBA toward ERβ is unclear. For example, an overview of the
compound 4 series, which in general had lower RBAs (Figure 2), suggests that they were also less potent than
the compound 2 series (Table 2); yet, in most cases, they were more efficacious on ERβ (Figure 4C). Overall, our results suggest that the imines
generally profile as potent, efficacious, subtype-selective ER ligands
that depend to a large extent on AF1 to induce ER-mediated transcription
fully, but they poorly stimulate the AF2-mediated activity of ERα
or ERβ. These differences underlie their ERα-selective
agonist and ERβ-selective antagonist properties.
Structural
Basis for the ER Subtype-Selective Profile of Triaryl-Substituted
Schiff Base Analogues
We obtained crystal structures of the
ERα ligand-binding domain complexed with the 2-Cl-substituted
analogues, compounds 2f and 4f, and compared
these new structures to previously reported full agonist- and antagonist-bound
ERα structures.[23−25] Full agonists, such as E2, fit into the ERα
ligand-binding pocket in an orientation that facilitates hydrogen
bonding of the phenolic OH to the side chains of helix 3 residue Glu353,
helix 6 residue Arg394, and, more variably, the D-ring 17β-OH
to helix 11 residue His524 (Figure 5A). This
binding orientation allows the switch helix, helix 12, to dock against
helices 3 and 11, where it forms one side of the coactivator binding
site on the surface that constitutes the functional core of AF2, that
is, it is a coregulator binding site.[23] In contrast, antagonists and SERMs, such as 4-hydroxytamoxifen,
contain an agonist-like core but have a bulky side group that protrudes
between helices 3 and 11 and directly displaces helix 12 from its
active conformation. As result, helix 12 docks along helix 3, thereby
occluding the AF2 surface.[24,25]
Figure 5
Imines induce a suboptimal
conformation of the ERα ligand-binding
domain. (A, B) Active and inactive ERα ligand-binding domain
conformations show a ∼1 Å difference in distance between
helices 3 and 11. The crystal structures of 17β-estradiol (E2;
PDB ID: 1GWR) and 4-hydroxytamoxifen (TAM; PDB ID: 3ERT) bound complexes are shown. (C, D) Crystal
structures of the ERα ligand-binding domain in complex with
compounds 2f and 4f show a TAM-like binding
orientation and increased h3–h11 distance compared to E2. (E,
F) Crystal structures of the imine-bound ERα complexes were
superposed. Compared to compound 2f (white), the additional
hydroxyl group of compound 4f (coral) leads to a subtle
distortion of the ligand-binding orientation.
Imines induce a suboptimal
conformation of the n>an class="Gene">ERα ligand-binding
domain. (A, B) Active and inactive ERα ligand-binding domain
conformations show a ∼1 Å difference in distance between
helices 3 and 11. The crystal structures of 17β-estradiol (E2;
PDB ID: 1GWR) and 4-hydroxytamoxifen (TAM; PDB ID: 3ERT) bound complexes are shown. (C, D) Crystal
structures of the ERα ligand-binding domain in complex with
compounds 2f and 4f show a TAM-like binding
orientation and increased h3–h11 distance compared to E2. (E,
F) Crystal structures of the imine-bound ERα complexes were
superposed. Compared to compound 2f (white), the additional
hydroxyl group of compound 4f (coral) leads to a subtle
distortion of the ligand-binding orientation.
Within the ligand-binding pocket, compn>ounds 2f and 4f mimicked the binding orientation of the n>an class="Chemical">4-hydroxytamoxifen
core without a protruding side chain (Figure 5C,D), consistent with their high binding affinities toward ER subtypes.
In addition, the N-phenyl groups were accommodated
between helices 8 and 11, contacting M421, H524, and L525. The additional
hydroxyl group of compound 4f was also accommodated easily,
as this ring was utilized as the A-ring mimetic that forms a hydrogen
bond with Glu353 (Figure 5E); however, this
hydroxyl substitution also led to a 0.8 Å rotation of the other
phenyl ring that forms a hydrogen bond with Thr347 (Figure 5F). Therefore, the different binding orientations
of these imines within the ligand-binding pocket account for the fact
that many of the matched 2 and 4 series
compounds have slightly different binding affinities (Figure 2) but similar ERα-mediated activity profiles
(Figure 4).
The structures of the imines
suggest that reduced AF2 activity
derives from a shift in the distance between helices 3 and 11. The
active helix 12 conformer docks across helices 3 and 11; therefore,
shifting helices 3 and 11 apart would undoubtedly destabilize helix
12 and lead to a suboptimal ligand-binding domain conformation that
would likely affect the binding of coregulators. In the full agonist-bound
conformation of ERα, the distance between helices 3 and 11,
measured from the α-carbon of Thr347 to the α-carbon of
Leu525, is a remarkably consistent 9 Å (Figure 5A).[23] ER ligands that are not full
agonists, including SERMs, increase this distance in a ligand-dependent
manner. When bound to 4-hydroxytamoxifen, this distance is increased
by about 1 Å (Figure 5B).[25] In contrast, compounds 2f and 4f increased this distance by 0.5 and 0.7 Å, respectively (Figure 5C,D), indicating that these compounds induce a suboptimal
conformation of the ERα ligand-binding domain, which explains
their inability to stimulate the AF2-mediated activity of ERα
fully (Figure 4B). Furthermore, the ER subtypes
have similar ligand-binding pockets; therefore, imines are likely
to also distort the active ERβ conformation through a similar
structural mechanism, which underlies their partial ERβ agonist/antagonist
phenotype (Table 2 and Figure 4C).
Discussion
The ne">estrogen receptors
are remarkable in binding and respn>onding
to lin>an class="Chemical">gands of great structural diversity,[8,26,27] and this eclectic acceptance of ligands
offers an opportunity to investigate chemically novel structures as
potential selective estrogen receptor modulators (SERMs)[10,11] or ER subtype-selective ligands.[8,9] For example,
recently, simple acyclic amide,[3] diphenylamine,[28] monoaryl- or diaryl-substituted salicylaldoxime,[29−33] and anthranylaldoxime[34,35] derivatives have been
reported as ER ligands. Thus, the development of new ER ligands remains
an important issue in medicinal chemistry because novel functions
of estrogen are still being found.[36,37]
In this
article, we have explored diversification of ligands for
the estrogen receptor by replacing the C=C double bond with
a simple Schiff Base or imine (C=N) core structure. A series
of triaryl-substituted Schiff base derivatives were conveniently prepared,
generally in one step, by the condensation of benzophenones with various
anilines without the need for phenolic hydroxy group protection. Many
of these compounds have very high binding affinities for ERs, rivaling
or exceeding that of estradiol (up to 97% RBA for ERα and 140%
for ERβ). Some of the compounds also show significant affinity
selectivity in favor of ERβ (4- to 5-fold), and in cell-based
assays for transcriptional activity, they profiled as ERα agonists
and ERβ antagonists, a form of ER-subtype selectivity that appears
to be based on their preferential activity through AF1, the N-terminal
activation function of the ERs that is more active in ERα than
ERβ. In addition, an unusual distortion of the ligand-binding
domain, revealed by X-ray analysis, suggests that the function of
AF2 is not fully engaged and highlights the diverse ways in which
ligands can regulate the conformation, and ultimately the activity,
of the ERs.Beyond our general interest in preparing ER ligands
of unusual
structure,[1] two other factors motivated
our investigation of imine-core systems for ER ligand design. Some
time ago, in an attempt to prepare gallium chelates that might have
good affinity for ER and could be used in Ga-67- or Ga-68-labeled
form for the imaging of ER in breast cancer using SPECT or PET, we
prepared some similar imines, notably, those having the 2-hydroxyl
group on the carbonyl component.[16] Unfortunately,
although we could prepare dimethylgallium complexes, engaging the
oxygen and nitrogen of the o-hydroxyphenyl imine
system as a bidentate chelate, and could even obtain their crystal
structures, the aqueous stability of these complexes was very low.[16] Nevertheless, at the time, we noted that a few
of the imines had some ER binding affinity. (The complete binding
data on these compounds is now given in Table
S1.) We also returned to the imines when we prepared the B—N
core compounds as part of our attempt to make the C=C and C=N
systems that were isoelectronic with the hydrolytically stable B—N
compounds; however, the one C=N analogue of these more sterically
encumbered systems that we were able to prepare lacked the critical
phenolic hydroxyl group and had other substituents that precluded
high-affinity binding.[13] By contrast, in
this investigation, we were gratified by how easy it was to prepare
a large series of triaryl-substituted imines and how many of them
had high binding affinity for both ERα and ERβ.
Tolerance of
Polar Groups in the Ligand Interior
Because
the intne">erior of the ER ligand binding pockets is lined strictly with
hydrophobic residues (except for the phenol–hydrogen bonding
glutamate and arginine residues and one threonine),[24] we expected that compounds in series 2 would
have lower affinity compared those in the series 4: the
isolated lone pair on the iminenitrogen in 2 was expected
to exact a large desolvation energy penalty, whereas in the series 4 imines, the intramolecular hydrogen bond would internally
solvate the iminenitrogen, thereby muting the impact of moving this
polar function from water into the binding pocket. Relevant to this
issue is the higher affinity of genistein than daidzein for both ERα
and ERβ: the ketone group in genistein is shielded by an intramolecular
hydrogen bond, whereas in daidzein, it is an isolated function (Figure 6).[38] A similar intramolecular
hydrogen-bonded system is present in the resorcylic acid lactone series
and macrolides with high affinity for the ERs exemplified by zeralanol[39] as well as in the salicylaldoximine and anthranylaldoxime
systems developed by Minutolo, where an intermolecular hydrogen bond
completes the formation of a pseudocycle that is thought to mimic
the phenolic ring of estradiol and is important for high affinity
binding (Figure 6).[34] (In addition, ER ligands with polar, hydrogen-bonding cores, such
as imidazoles and pyridazines, bind much less well than their less
polar analogues, pyrazoles and pyrazines, respectively.[1])
Figure 6
Structure and relative binding affinities (RBAs, estradiol
= 100)
of various nonsteroidal and seco-steroidal estrogens as well as torsional
angles of the triarylimine propellane conformation.
Structure and relative binding affinities (RBAs, n class="Chemical">estradioln>
= 100)
of various nonsteroidal and seco-steroidal estrogens as well as torsional
angles of the triarylimine propellane conformation.
Despite these precedents for the importance of
intramolecular hydrogen
bonding to shield a polar element in the intn>an class="Gene">erior of ER ligands, in
nearly every case, the imines of series 2, all of which
lack this hydrogen bond, have comparable or higher affinity than the
corresponding members in series 4 (Figure 2). This was particularly evident in their ERβ-binding
affinities, where the average ratio of RBA (series 2)/RBA
(series 4) is 1.88 ± 0.70, whereas it is only 1.16
± 0.78 for ERα.
A reexamination of the crystal structures
for compounds 2f and 4f shows that at the
imine side of the ligands
there is ample space between the ligands and the pocket residues,
with no evidence for interaction of the imine lone pair in compound 2f with any elements of the protein. Also, the additional
hydroxyl group in compound 4f does not engage in hydrogen
bonding with the protein and is nicely accommodated without any obvious
effects on the shape of the pocket at this side of the ligands. Perhaps
the only factor contributing to the differences between the two series
could be the increased twist angle of the N-phenyl
ring noted in the X-ray structure of the series 4 compound 4f versus the series 2 compound 2f (Figure 5C–F), which appears to interfere
with a productive interaction with the Thr347 residue. It is notable,
as well, that in a number of ligand series, ERβ appears more
tolerant to interior polar groups.[8]
Number
of Substituents on the C=N Core
Although
ligands having C=C or B—N core elements could be tetrasubstituted,
the imine-core ligands, for valency reasons, can, at most, be trisubstituted.
Nevertheless, triaryl ethylene ligands with a C=C core lacking
a fourth substituent often have very good ER-binding affinity (Figure 6),[41] as do our imines.
By contrast, the mono and diaryl imines (Table
S1) have uniformly low affinity, presumably because they are
too small. There are a number of ER ligands in which a C=C
core has been replaced by two nitrogens (i.e., an azo or N=N
group), and these, for valence reasons, can be only disubstituted.
These azophenol or azoresorcinol systems have low, but clearly measurable,
ER-binding affinity (Figure 6) and are known
to be estrogens, although of very low potency.[42,43] In these molecules, as with the imines, there are ortho-substituted
hydroxyl groups that can engage in intramolecular hydrogen bonds,
but again, these hydrogen bonds reduce binding (Figure 6).
Effect of Phenyl Ring Substituents
In our imines, addition
of a single ortho methyl group in the distal n>an class="Chemical">N-phenol
ring improved binding considerably (Table 1, compound 2c vs 2a), whereas addition
of a second ortho methyl group caused a precipitous drop in binding
affinity (Table 1, compound 2b vs 2c). We encountered the beneficial effect of single
ortho methyl substitution in A-CD estrogens (B-seco steroids lacking
an intact B-ring), where addition of a methyl group (or other small
substituent, Cl, CF3) ortho to the site of attachment of
the phenol to the C-ring increased affinity (Figure 6).[44] In these seco-estrogens, we
interpreted the increase in binding resulting from this single ortho
substitution as the dual effect of supplying bulk that was lost by
deletion of the B-ring as well as twisting the aryl ring relative
to the rest of the ligand, thereby increasing ligand volume, at least
on one side. Similarly, in otherER ligands that we have explored,
a single twist-inducing ortho substituent generally increased binding
affinity.[2,45,46]
Substituents
on the N-phenyl ring in our triarylimines are directed
to different regions of the ER ligand-binding pocket than those in
the B-secosteroids, and their enhancement of binding is most likely
just due to increased hydrophobic bulk. There is a general trend,
with binding affinity decreasing with the series C-2 > C-3 >
C-4,
which suggests that the ortho-disposed (C-2) groups might be causing
an increased twist of the N-phenyl group. Simple MM2 energy minimization,
however, shows that all of these triaryl imines adopt a propeller-like
conformation; even the unsubstituted systems, 2a and 4b, show a coordinated twist of all three rings, giving dihedral
angles of 160, 155, and 172° for τ1, τ2, and τ3, respectively (Figure 6). Notably, addition of the substituents shown in Table 1, even those at the C-2 position, changes these
torsions by only a few degrees. Even 2,6-dimethyl substitution in 2c causes only a 10° change in τ3, suggesting
that the marked drop in affinity results from a steric clash of the
second methyl group in the ligand-binding pocket.
Structural
Mechanisms for Ligand-Dependent Modulation of ER
Activity
ne">ER binds a diverse collection of ligands that modulate
its activity through distinct structural mechanisms. Full agonists
directly drive helix 12 to adopt an active conformation where it lies
across helices 3 and 11, positioned to form one side of the AF2 coactivator-binding
surface.[24] In contrast, the bulky side
chain of SERMs protrudes out of the pocket and directly displaces
helix 12 from this active conformation.[24,25] Yet, several
other ligands disrupt the active conformation of helix 12 indirectly
by distorting the C-terminal end of helix 11,[7,45] thereby
modulating ER activity. Here, we show that triaryl imine ligands modulate
ER activity through a new type of indirect mechanism that involves
a ligand-dependent increase in the distance separating helices 3 and
11 (Figure 5). Through this new mechanism,
these imine ligands are able to modulate the AF2-mediated but not
AF1-mediated activity of ER (Figure 4).
The ease with which these C=Nimine analogues can be prepared,
in some cases much more readily than their C=C double bond
countn>an class="Gene">erparts, and the high binding affinity that they have for the
ERs suggest that this form of analogue development might be worth
exploring more generally in drug discovery. Furthermore, the marked
pattern of ER-subtype differential cellular efficacy, based on differential
utilization of the two ER activation functions, displayed by the C=N
analogues appears to be based on a novel mode of conformational control
of the ER ligand-binding pocket that affects AF2 function, as revealed
by our X-ray structural analyses. This new ligand design paradigm
could also be more widely investigated for the development of estrogens
having a novel spectrum of activities.
Experimental
Section
Analytical Techniques and Instrumentation Used
Melting
points were determined on a X-4 Beijing Tech melting point apparatus. 1H and 13C NMR spectra were recorded on a Bruker
AV400 spectrometer (400 MHz, 1H NMR; 101 MHz, 13C NMR) at room temperature. NMR spectra were calibrated to the solvent
signals of CDCl3 (δ 7.26 and 77.00 ppm), acetone-d6 (δ 2.05 and 29.84 ppm, 206.26 ppm),
or DMSO-d6 (δ 2.50 and 39.43 ppm).
The chemical shifts are provided in ppm, and the coupling constants,
in Hz. The following abbreviations for multiplicities are used: s,
singlet; d, doublet; dd, double doublet; t, triplet; dt, double triplet;
q, quadruplet; m, multiplet; and br, broad. The purity of all compounds
for biological testing was determined by HPLC analysis in two different
solvent systems (normal and reversed phase), confirming >95% purity
(see the Supporting Information).
Chemical
Synthesis. General Procedure of the Improved Method
for the Synthesis of Compounds 2a–m and 4a–m
A mixture of 1 (1 mmol) and aniline derivatives (3 mmol, 3 equiv) was dissolved
in chlorobenzene (4 mL). Under a N2 atmosphere, the mixture
was briefly exposed to HCl(g) and then heated at 140–145 °C
for 24 h. The resulting solution was concentrated under reduced pressure,
and the crude product was purified by silica gel chromatography (petroleum
ether/ethyl acetate: Et3N = 4:1:0.5) to afford target molecules 2 and 4. Further purification was achieved by
recrystallization (ethyl acetate/petroleum ether).
4,4′-((Phenylimino)methylene)diphenol
(2a)
According to the general procedure of the
impn>roved method, 2a was obtained as a yellow solid (60%
yield) and was furthn>an class="Gene">er
purified by recrystallization from ethyl acetate/petroleum ether (mp
273–275 °C). 1H NMR (400 MHz, acetone-d6) δ 8.83 (br s, 2H), 7.66–7.60
(m, 2H), 7.15–7.08 (m, 2H), 7.00–6.94 (m, 2H), 6.92–6.83
(m, 3H), 6.77–6.71 (m, 2H), 6.66 (dt, J =
8.4 Hz, 1.6 Hz, 2H). 13C NMR (101 MHz, acetone-d6) δ 167.97, 160.74, 158.37, 153.24, 132.58,
131.96, 131.94, 129.18, 128.51, 123.11, 121.78, 115.68, 115.49. HRMS
(MALDI/DHB) calcd for C19H17NO2 (M
+ H)+m/z, 290.11811;
found, 290.11756.
According to the general procedure of the
impn>roved
method, 2m was obtained as a yellow solid (74% yield).
Then, 2m was dissolved in n>an class="Chemical">THF and MeOH, a solution of
KOH (1 equiv) in MeOH was added to the stirred mixture for 3h, the
mixture was concentrated under reduced pressure, and the residue was
isolated by silica gel chromatography (petroleum ether/ethyl acetate
= 2:1, including 0.5% Et3N) to afford the desired product 2n (94% yield), which was further purified by recrystallization
from ethyl acetate/petroleum ether (mp 297–301 °C). 1H NMR (400 MHz, acetone-d6) δ
7.46 (d, J = 8.7 Hz, 2H), 6.85–6.79 (m, 2H),
6.73 (d, J = 8.7 Hz, 2H), 6.64 (t, J = 8.1 Hz, 2H), 6.51–6.44 (m, 2H), 6.42–6.36 (m, 2H). 13C NMR (101 MHz, acetone-d6) δ167.29,
160.52, 158.30, 153.94, 145.10, 132.99, 131.95, 131.73, 128.93, 123.34,
116.04, 115.82, 115.59. HRMS (MALDI/DHB) calcd for C19H17NO3 (M + H)+m/z, 306.11302; found, 306.11247.
According to the general procedure of the
impn>roved method, 4m was obtained as a yellow solid (65%
yield). Then, 4m was dissolved in n>an class="Chemical">THF and MeOH, a solution
of KOH (1 equiv) in MeOH was added to the stirred mixture for 3 h,
and the mixture was concentrated under reduced pressure, and the residue
was isolated by silica gel chromatography (petroleum ether/ethyl acetate
= 2:1, including 0.5% Et3N) to afford desired product 4n (92% yield), which was further purified by recrystallization
from ethyl acetate/petroleum ether (mp 269–271 °C). 1H NMR (400 MHz, acetone-d6) δ
15.50 (s, 1H), 8.94 (s, 1H), 8.40 (s, 1H), 7.05 (d, J = 8.4 Hz, 2H), 6.96 (d, J = 8.8 Hz, 1H), 6.86 (d, J = 8.4 Hz, 2H), 6.65 (d, J = 8.7 Hz, 4H),
6.41 (d, J = 2.3 Hz, 1H), 6.27 (dd, J = 8.8 Hz, 2.3 Hz, 1H). 13C NMR (101 MHz, acetone-d6) δ 173.11, 162.84, 158.74, 144.22, 134.58,
132.02, 131.39, 129.65, 126.59, 124.95, 116.00, 115.93, 114.27, 107.23,
104.09. HRMS (MALDI/DHB) calcd for C19H17NO4 (M + H)+m/z, 322.10794; found, 322.10739.
Relative
binding affinities wne">ere determined by a competitive radiometric binding
assay, as previously described,[47,48] using 2 nM [3H]estradiol as tracer ([2,4,6,7-3H]-estra-1,3,5(10)-triene-3,17β-diol,
70–115 Ci/mmol, PerkinElmer, Waltham, MA) and purified full-length
humanERα and ERβ, which were purchased from PanVera/Invitrogen
(Carlsbad, CA). Incubations were for 18–24 h at 0 °C.
Hydroxyapatite (Bio-Rad, Hercules, CA) was used to absorb the receptor–ligand
complexes, and free ligand was washed away. The binding affinities
are expressed as relative binding affinity (RBA) values, with the
RBA of estradiol set to 100%. The values given are the average ±
range or SD of two to three independent determinations. Estradiol
binds to ERα with a Kd of 0.2 nM
and to ERβ with a Kd of 0.5 nM;
these values were determined by Scatchard analysis using the binding
assay protocol described previously.[47]
Gene Transcriptional Activity
Assays wne">ere pn>an class="Gene">erformed
as previously described.[45] HepG2 cells
cultured in Dulbecco’s minimum essential medium (DMEM) (Cellgro
by Mediatech, Inc. Manassas, VA) supplemented with 10% fetal bovine
serum (FBS) (Hyclone by Thermo Scientific, South Logan, UT), 1% nonessential
amino acids (Cellgro), penicillin–streptomycin–neomycin
antibiotic mixture, and Glutamax (Gibco by Invitrogen Corp. Carlsbad,
CA) were maintained at 37 °C and 5% CO2. HepG2 cells
were transfected with 10 μg of 3×ERE-luciferase reporter
plus 1.6 μg of ERα or ERβ expression vector per
10 cm dish using FugeneHD reagent (Roche Applied Sciences, Indianapolis,
IN). The next day, the cells were transferred to phenol red-free growth
media supplemented with 10% charcoal-dextran sulfate-stripped FBS
at a density of 20 000 cells/well, incubated in 384-well plates
overnight at 37 °C with 5% CO2, and assayed in dose
curves ranging from 10 μM to 100 pM for ERE luciferase assays
in HepG2 cells. Luciferase activity was measured after 24 h using
BriteLite reagent (PerkinElmer Inc., Shelton, CT) according to the
manufacturer’s protocol. Raw data measured as relative light
units were normalized for each plate using the average of DMSO-treated
samples as 0% and the average of the top of the E2 curve as 100%.
X-ray Crystallography
As previously described,[49] humanERα-Y537S ligand-binding domain
was expressed in BL21(DE3) Escherichia coli cells, purified, mixed with SRC2 peptide, and crystallized at room
temperature by the hanging-drop vapor-diffusion method. The ERα
crystals obtained were then soaked in the different imine ligands.
The X-ray diffraction data was scaled using HKL-2000 software.[50] The crystal structures were solved via molecular
replacement using the PHENIX software suite,[51,52] with the protein components of PDB 2B1V as a starting model.[53] The new structures were then completed upon ligand docking
and extensive combinatorial refinement.[54]
Authors: Robert W Hsieh; Shyamala S Rajan; Sanjay K Sharma; Yuee Guo; Eugene R DeSombre; Milan Mrksich; Geoffrey L Greene Journal: J Biol Chem Date: 2006-04-28 Impact factor: 5.157
Authors: Jian Min; Pengcheng Wang; Sathish Srinivasan; Jerome C Nwachukwu; Pu Guo; Minjian Huang; Kathryn E Carlson; John A Katzenellenbogen; Kendall W Nettles; Hai-Bing Zhou Journal: J Med Chem Date: 2013-04-15 Impact factor: 7.446
Authors: J C Gould; L S Leonard; S C Maness; B L Wagner; K Conner; T Zacharewski; S Safe; D P McDonnell; K W Gaido Journal: Mol Cell Endocrinol Date: 1998-07-25 Impact factor: 4.102
Authors: Hai-Bing Zhou; Kendall W Nettles; John B Bruning; Younchang Kim; Andrzej Joachimiak; Sanjay Sharma; Kathryn E Carlson; Fabio Stossi; Benita S Katzenellenbogen; Geoffrey L Greene; John A Katzenellenbogen Journal: Chem Biol Date: 2007-06
Authors: Jerome C Nwachukwu; Sathish Srinivasan; Yangfan Zheng; Song Wang; Jian Min; Chune Dong; Zongquan Liao; Jason Nowak; Nicholas J Wright; René Houtman; Kathryn E Carlson; Jatinder S Josan; Olivier Elemento; John A Katzenellenbogen; Hai-Bing Zhou; Kendall W Nettles Journal: Mol Syst Biol Date: 2016-04-22 Impact factor: 11.429
Authors: Sarah Preston; Junjie Luo; Yuezhou Zhang; Abdul Jabbar; Simon Crawford; Jonathan Baell; Andreas Hofmann; Min Hu; Hai-Bing Zhou; Robin B Gasser Journal: Parasit Vectors Date: 2016-06-16 Impact factor: 3.876