Deborah A Smithen1,2, Leo M H Leung1,2, Mairi Challinor1, Rae Lawrence1, HaoRan Tang3, Dan Niculescu-Duvaz1,2, Simon P Pearce4, Robert Mcleary1,2, Filipa Lopes1,2, Mohammed Aljarah1,2, Michael Brown1,2, Louise Johnson2, Graeme Thomson1, Richard Marais3, Caroline Springer1,2. 1. Drug Discovery Unit, Cancer Research UK Manchester Institute, University of Manchester, Alderley Park, Macclesfield SK10 4TG, United Kingdom. 2. Cancer Research UK Centre for Cancer Therapeutics, The Institute of Cancer Research, 15 Cotswold Road, London SM2 5NG, United Kingdom. 3. Molecular Oncology Team, Cancer Research UK Manchester Institute, University of Manchester, Alderley Park, Macclesfield SK10 4TG, United Kingdom. 4. Clinical and Experimental Pharmacology, Cancer Research UK Manchester Institute, University of Manchester, Alderley Park, Macclesfield SK10 4TG, United Kingdom.
Abstract
The lysyl oxidase (LOX) family of extracellular proteins plays a vital role in catalyzing the formation of cross-links in fibrillar elastin and collagens leading to extracellular matrix (ECM) stabilization. These enzymes have also been implicated in tumor progression and metastatic disease and have thus become an attractive therapeutic target for many types of invasive cancers. Following our recently published work on the discovery of aminomethylenethiophenes (AMTs) as potent, orally bioavailable LOX/LOXL2 inhibitors, we report herein the discovery of a series of dual LOX/LOXL2 inhibitors, as well as a subseries of LOXL2-selective inhibitors, bearing an aminomethylenethiazole (AMTz) scaffold. Incorporation of a thiazole core leads to improved potency toward LOXL2 inhibition via an irreversible binding mode of inhibition. SAR studies have enabled the discovery of a predictive 3DQSAR model. Lead AMTz inhibitors exhibit improved pharmacokinetic properties and excellent antitumor efficacy, with significantly reduced tumor growth in a spontaneous breast cancer genetically engineered mouse model.
The lysyl oxidase (LOX) family of extracellular proteins plays a vital role in catalyzing the formation of cross-links in fibrillar elastin and collagens leading to extracellular matrix (ECM) stabilization. These enzymes have also been implicated in tumor progression and metastatic disease and have thus become an attractive therapeutic target for many types of invasive cancers. Following our recently published work on the discovery of aminomethylenethiophenes (AMTs) as potent, orally bioavailable LOX/LOXL2 inhibitors, we report herein the discovery of a series of dual LOX/LOXL2 inhibitors, as well as a subseries of LOXL2-selective inhibitors, bearing an aminomethylenethiazole (AMTz) scaffold. Incorporation of a thiazole core leads to improved potency toward LOXL2 inhibition via an irreversible binding mode of inhibition. SAR studies have enabled the discovery of a predictive 3DQSAR model. Lead AMTz inhibitors exhibit improved pharmacokinetic properties and excellent antitumor efficacy, with significantly reduced tumor growth in a spontaneous breast cancer genetically engineered mouse model.
The
lysyl oxidase (LOX) family of copper-dependent extracellular
proteins comprises the founder member enzyme, LOX, and four LOX-like
enzymes (LOXL1–4).[1−6] While there is greater than 50% sequence identity between the isoforms,
which includes a conserved C-terminal catalytic domain across the
family containing the copper binding site and the lysine tyrosylquinone
(LTQ) cofactor, the enzymes can be divided into two subgroups based
on differences to the N-terminal structure. Indeed, LOX and LOXL1
contain a variable N-terminal propeptide that undergoes proteolytic
cleavage to form the active enzyme extracellularly. LOXL2–4
differ in that they do not possess this propeptide region and instead
contain four scavenger receptor cysteine-rich (SRCR) domains at the
N-terminus, which are thought to mediate protein–protein interactions
in the ECM. In the case of LOXL2, the protein undergoes proteolytic
cleavage of the first two SRCR domains upon secretion; however, unlike
LOX, processing is not required for catalytic activation.[7,8]The most widely studied function of the LOX enzymes is their
ability
to form cross-links in fibrillar elastin and collagens through oxidative
deamination of specific lysyl residues, thus stabilizing the ECM.[3] However, recent reports suggest that these enzymes
have a multitude of biological functions, which include cell proliferation
and epithelial–mesenchymal transition (EMT).[3,5,9] Consequently, the more widely studied LOX
and LOXL2 isoforms have been implicated in tumor progression, where
they are highly expressed and actively involved in remodeling the
tumor microenvironment.[10−16]The LOX family has thus become an attractive therapeutic target
for the treatment of many types of invasive cancers, particularly
those with poor patient outcomes. Targeting LOXs with small molecule
inhibitors is very challenging owing to the lack of crystal structures
useful for drug design for any of the isoforms (the only reported
LOXL2 structure is a precursor state without cofactor formed)[17] and difficulties associated with isolating several
of the enzymes in an active form. Nevertheless, in recent years several
LOXL2-selective inhibitors have been reported, including haloallylamine-based
inhibitors PXS-S1A and the highly potent PXS-S2A (full structures
not disclosed),[18] as well as dual LOXL2/LOXL3
inhibitors PXS-5153A (1)[19] and aminomethylenepyridine 2 (Figure ).[20] Our dual LOX/LOXL2 inhibitor CCT365623 (3a)[21,22] is an orally efficacious aminomethylenethiophene (AMT) based
inhibitor, which has been used to help elucidate mechanisms by which
LOX drives tumor progression. These novel inhibitors offer significant
advantages over the prototypical pan-LOX inhibitor β-aminopropionitrile
(BAPN),[23,24] whose lack of sites amenable for chemical
modification precludes preclinical optimization.
Figure 1
Small molecule inhibitors
of LOX family enzymes.
Small molecule inhibitors
of LOX family enzymes.We recently reported
the discovery of AMT inhibitor CCT365623 (3a) following
a significant medicinal chemistry campaign to
elucidate the structure–activity relationship (SAR) of this
class of compound with respect to LOX inhibition.[22] Systematic modifications were made to a hit compound identified
following a high-throughput screen (HTS), leading to development of
submicromolar half maximal inhibitory concentration (IC50) inhibitors possessing desirable selectivity and pharmacokinetic
(PK) properties.During the course of these studies, the 2,5-substituted
thiophene
core was replaced with various other five-membered heterocyclic rings
to ascertain the importance of this moiety on activity. Of those assessed,
only a 2-aminomethylene-5-sulfonyl thiazole core retains activity,
with naphthalenesulfonyl-substituted thiazole 6 showing comparable levels of LOX inhibition to the analogous thiophene
compound 4 (Figure ).[22] By contrast, thiazole
regioisomer 5 is relatively inactive. Interestingly,
in the case of the active thiazole compound (6) we also
observe a modest increase in potency toward LOXL2 inhibition, with 6 proving equipotent against both isoforms. These observations
prompted parallel investigations into the development of 2-aminomethylene-5-sulfonylthiazoles
(AMTz) as dual LOX/LOXL2 inhibitors, and given the commercial availability
of purified LOXL2 enzyme, these studies were carried out using LOXL2
enzyme, both as a target and as a surrogate to assess LOX activity.
Figure 2
LOX and
LOXL2 activity of analogous AMT and AMTz inhibitors.
LOX and
LOXL2 activity of analogous AMT and AMTz inhibitors.
Results and Discussion
Time-Dependent Inhibition
Initial
studies concerned
the enzyme–inhibitor preincubation time in our biochemical
assay, whereby we assessed whether our compounds inhibit enzyme activity
in a time-dependent manner. Dual LOX/LOXL2 inhibitors 4 and 6, along with the less active 5-aminomethylene-2-sulfonylthiazole regioisomer 5, were assessed using longer preincubation
times of 1 and 3 h, and the activity was compared to previously obtained
20 min preincubation data (Table ). The results of this study demonstrate that longer
preincubation times result in increased levels of enzyme inhibition
with the greatest difference in effect observed when the time is increased
from 20 min to 1 h, upon which up to a 5-fold increase in potency
is observed. Further increasing the time to 3 h has a smaller positive
effect on activity while remaining within 2-fold of the 1 h data.
On the basis of these findings, we decided to employ a 1 h preincubation
time for the remainder of the studies described herein, which provides
an optimal determination of enzyme inhibition at a physiologically
relevant time point.
Table 1
Effect of Preincubation
Time on LOXL2
Potency
entry
compd
X
Y
LOXL2, 20 min, IC50 (μM)a,b
LOXL2, 1 h, IC50 (μM)a,c
LOXL2, 3 h, IC50 (μM)a,d
1
4
CH
CH
3.40
0.81
0.47
4.76
0.85
0.59
2
5
CH
N
26.8
6.41
3.54
27.0
6.02
3.19
3
6
N
CH
0.86
0.28
0.12
0.69
[0.22, 0.36]
0.12
Reported IC50 values
were determined in at least two separate experiments (n ≥ 2). When n = 2, individual IC50 values are shown. When n > 2, the values are
reported
as the geometric mean with the corresponding 95% confidence interval
in square brackets.
Compound
was preincubated with enzyme
for 20 min.
Compound was
preincubated with enzyme
for 1 h.
Compound was preincubated
with enzyme
for 3 h.
Reported IC50 values
were determined in at least two separate experiments (n ≥ 2). When n = 2, individual IC50 values are shown. When n > 2, the values are
reported
as the geometric mean with the corresponding 95% confidence interval
in square brackets.Compound
was preincubated with enzyme
for 20 min.Compound was
preincubated with enzyme
for 1 h.Compound was preincubated
with enzyme
for 3 h.
Thiophene vs Thiazole Core
Optimization studies began
with our recently published AMT inhibitor (3a),[22] whereby we modified the core to the analogous
1,3-thiazole compound (7a) and assessed LOXL2 activity
inhibition (Table , compare entries 1 and 4). Pleasingly, this gives a modest increase
in potency, with an IC50 of 0.086 μM. We then assessed
whether this increase in activity, attributed to the presence of a
nitrogen atom in the ring, could be mimicked through the incorporation
of an electronegative halogen atom on the thiophene ring, which could
also engender the formation of an intramolecular hydrogen bond to
improve ligand binding:[25] 3-chlorothiophene
(entry 2) was found to be less potent by >10-fold compared to the
parent compound (entry 1), while a fluorine substituent was 4-fold
less active (entry 3). A second matched pair was synthesized to confirm
this trend, and again the potency achieved with the thiazole analogue
was slightly greater than that of the thiophene compound (compare
entries 6 and 7). This suggests that the presence of a nitrogen atom
in the heterocyclic core is advantageous, providing additional stabilization
to the protein–inhibitor complex through either resonance mechanisms
or noncovalent interactions.
Table 2
Effects of Thiophene
vs Thiazole Cores
on LOXL2 Potency
entry
compd
X
R1
R2
SO2Me position
LOXL2, IC50 (μM)a
1
3a
CH
H
H
3
0.176
[0.105, 0.295]
2
3b
CCl
H
H
3
2.488
2.164
3
3c
CF
H
H
3
0.671
0.714
4
7a
N
H
H
3
0.086
[0.061, 0.119]
5
7b
N
CH3
H
3
68.26
[43.77, 106.5]
6
8
CH
H
CH3
4
0.686
[0.456, 1.031]
7
9
N
H
CH3
4
0.423
[0.302, 0.593]
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.In order to confirm that the aminomethylene group
was still essential
for inhibition, we assessed the effect of an N-methyl
substituent group (entry 5). As expected, this results in significant
loss of activity, which is consistent with previous conclusions that
the aminomethylene group is required to form a stable Schiff base
between the inhibitor and the LTQ cofactor that is a feature of both
LOX and LOXL2.
Modification of Thiazole C-5 Group
We next assessed
the C-5 substituent group of the thiazole-based inhibitors to determine
whether the same SAR trends were observed as in the thiophene series
and whether the lead scaffold remained optimal. An aryl sulfonyl group
is preferred to either an alkyl or cycloalkyl group (Table , compare entries 1–3).
Focusing on the aryl substituents, we find that monosubstitution with
either a methane sulfonyl (entry 4) or a phenyl group (entry 6) is
tolerated since these are equipotent with the unsubstituted example
(entry 3), but a disubstituted aryl group (entry 8) remains preferred
for anti-LOXL2 activity.
Table 3
Effect of Thiazole
C-5 Linker (L)
and Substituent (R) Groups on LOXL2 Potency
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.We assessed the effect of the oxidation level of the
sulfonyl linker
and find that use of a sulfoxide or sulfide linker results in less
active compounds (compare entries 8–10). This is expected based
on previous studies; however, it is noted that the impact resulting
from a decrease in oxidation state is less significant in the AMTz
series herein than it was in the AMT series.[22] Removal of the sulfonyl linker is tolerated but results in partial
loss of activity, with the most significant impact observed in the
biaryl systems (compare entries 4 and 5; 6 and 7; 8 and 11). As explored
later (Table ), this
increased tolerance to a range of linkers appears to be a feature
of the thiazole compared to thiophene core, whereby the presence of
the nitrogen atom in the ring increases the electron-withdrawing properties,[26,27] thus mitigating the need for an electron-withdrawing linker.
Table 6
Activity of AMT and AMTz Inhibitors
against LOX, LOXL2, and LOXL3
IC50 (μM)a
compd
LOXL2
LOX
LOXL3
BAPN
0.665
3.14
0.31
[0.57, 0.77]
[1.79, 5.49]
[0.22, 0.45]
3a
0.176
0.898
13.49
[0.105, 0.295]
0.728
18.54
7a
0.086
1.77
0.92
[0.061, 0.12]
1.49
0.86
7d
0.233
1.31
4.70
[0.19, 0.28]
1.13
3.85
7e
0.38
7.17
6.81
[0.27, 0.54]
6.91
6.78
17
0.652
51.96
9.31
0.716
46.49
9.52
21a
0.109
3.42
1.22
[0.066, 0.18]
2.65
1.51
21b
0.151
3.3
1.25
[0.12, 0.19]
1.86
1.43
22b
0.302
>100
3.86
[0.22, 0.42]
>100
4.26
22c
0.37
>100
7.56
[0.23, 0.596]
>100
6.69
22d
0.425
>100
24.24
[0.35, 0.51]
>100
36.86
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.
Variation
of the Phenylsulfonyl Ring Substituents
Incorporation
of a bis-sulfonyl biarylC-5 group is advantageous for activity and
was previously found to improve oral in vivo PK exposure in the AMT
series of LOX inhibitors.[22] We subsequently
looked to ascertain if there was scope to optimize the substituent
effects further (Table ). Varying the methanesulfonyl group (R2) to an ethyl
or isopropyl group is tolerated, albeit resulting in a drop in potency;
however, a larger tert-butyl group proves detrimental
to activity (compare entries 1–4). A methoxy substituent is
similarly well tolerated (entry 5), while a methylamino group appears
to be less favorable (entry 6).
Table 4
Effects of 3,5-Disubstituted
Phenylsulfonyl
Groups on LOXL2 Potency
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.Small modifications to the aryl substituent (R1) do
not result in a significant change in potency. Inhibitors bearing
either an electron-donating p-methyl group or an
electron-withdrawing p-fluoro or p-trifluoromethyl substituent demonstrate comparable levels of activity
to the parent compound (compare entries 7–9 to entry 1). Replacing
the phenyl group with an N-methyl pyrazolyl group
is found to be favorable in conjunction with an ethyl group (compare
entries 2 and 11), though it does not appear to confer additional
potency when R2 is a methanesulfonyl group (compare entries
1 and 10). Further modification of the ethyl group of 21b to an electron-withdrawing trifluoromethyl or chloro group has a
slightly negative effect on activity, though a fluoro substituent
is well tolerated (compare entries 11–14). Overall, from these
results we conclude that the parent inhibitor 7a is well
optimized while gaining further understanding of the SAR and identifying
additional scaffolds for further study.
Activity of AMTz Regioisomers
Final SAR studies concerned
the substitution pattern around the thiazole core, whereby we compared
the anti-LOXL2 activity of selected 2,5-substituted AMTz inhibitors
to that of their 2,4-substituted regioisomers (Table ). Although there is a trend favoring 5-substitution
vs 4-substitution, it is again less pronounced than that previously
observed in the thiophene series against LOX[22] (compare entries 1 and 2; 3 and 4; 5 and 6; 7 and 8).
Table 5
Effects of 2,4- vs 2,5-Substitution
of Thiazole Core on LOXL2 Potency
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.
Selectivity Studies
We assessed the selectivity profiles
of our AMTz inhibitors against LOX and LOXL3 isoforms (Table ) and against common amine oxidases and the potassium ion
channel hERG (Table ).
Table 7
Potency and Selectivity of AMTz Inhibitors
over Common Amine Oxidases and hERG
IC50 (μM)a
compd
LOXL2
MAO-A
MAO-B
DAO
SSAO
hERG
6
0.28
>100
>100
>100
47.2
51
[0.22, 0.36]
>100
>100
>100
47.1
7a
0.086
63.8
>100
>100
3.2
10
[0.061, 0.12]
59.3
>100
>100
1.2
7e
0.38
29.2
64.3
>100
5.2
36
[0.27, 0.54]
31.4
46.2
>100
3.6
18
0.103
100
16.4
>100
2.7
3.2
[0.083, 0.13]
88.3
20
>100
1.3
21a
0.109
>100
>100
>100
13.3
nd
[0.066, 0.18]
>100
>100
>100
17.2
21b
0.151
>100
>100
>100
28.6
47
[0.12, 0.19]
>100
>100
>100
19.6
21e
0.143
61.7
>100
21.3
36.3
66
[0.057, 0.36]
66.8
>100
22.5
40.3
22b
0.302
>100
>100
>100
>100
nd
[0.22, 0.42]
>100
>100
>100
>100
22c
0.37
>100
>100
>100
>100
nd
[0.23, 0.596]
>100
>100
>100
>100
22d
0.425
>100
>100
>100
102.2
29
[0.35, 0.51]
>100
>100
>100
74.3
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets. nd: not determined.
Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets.Reported IC50 values
were determined using 1 h preincubation in at least two separate experiments
(n ≥ 2). When n = 2, individual
IC50 values are shown. When n > 2,
the
values are reported as the geometric mean with the corresponding 95%
confidence interval in square brackets. nd: not determined.
Activity Profiles against LOX and LOXL3 Isoforms
While
use of LOXL2 as a surrogate for LOX was a practical and informative
approach, it was necessary to assess the activity of a range of AMTz
inhibitors against LOX to confirm our belief that these are in fact
dual LOX/LOXL2 inhibitors. A diverse set of compounds were selected,
encompassing 2-aminomethylene-5-sulfonyl- (7a, 7e, 21a, and 21b) and 2-aminomethylene-4-sulfonyl-
(22b, 22c, and 22d) AMTz regioisomers,
sulfide-linked 7d, and 17 with no linker,
as well as thiophene 3a and BAPN. From this study we
were able to confirm that 2-aminomethylene-5-sulfonyl thiazoles 7a, 21a, and 21b, thiophene 3a, and BAPN exhibit good anti-LOX activity that is in accordance
with previously obtained data.[21,22,28] A slight decrease in potency toward LOX inhibition is observed with
compound 7e, which is consistent with LOXL2 activity.
Removal of the sulfonyl linker (17) is detrimental to
activity; however, a sulfide linker (7d) is well tolerated
which, as discussed previously, can be attributed to the presence
of a nitrogen atom in the thiazole core which presumably increases
the ability of these compounds to form a stable covalent bond upon
binding. Interestingly, 2,4-AMTz regioisomers (22b, 22c, and 22d) are found to be inactive against
LOX, which follows a similar trend to the thiophene series in which
2,4-regioisomers were demonstrated to be 15-fold less potent.[22] As seen previously, all AMTz inhibitors demonstrate
increased potency toward LOXL2 inhibition versus LOX.As discussed
previously, literature compounds 1 and 2 are in fact dual LOXL2/LOXL3 inhibitors,[19,20] while selectivity data concerning LOXL isoforms have not been reported
for other literature inhibitors. As a member of the LOX family, LOXL3
is known to modulate the ECM,[29−31] though there is tissue expression
variance compared to other LOX proteins, and it has been shown to
play a significant role in muscular, skeletal, and lung development
in mice.[32−34] More recently, studies have demonstrated an involvement
of LOXL3 in cancer and metastasis, suggesting it to be a potential
therapeutic target for malignant disease.[35,36]With regards to LOXL3 inhibition, it is interesting to observe
that our AMT and AMTz compounds exhibit moderate to high selectivity
toward LOXL2 in all cases (Table ), unlike the nonselective LOX-family inhibitor BAPN
and reported literature compounds 1 and 2 (activity within 3-fold for LOXL2 and LOXL3 in all cases).[19,20] 2,5-AMTz inhibitors 7a, 7e, 21a, and 21b demonstrate ≥10-fold selectivity toward
LOXL2, with IC50 values of approximately 1 μM; compound 7e is consistently less potent against all isoforms. Modification
of the sulfonyl linker to a sulfide (7d) or direct aryl-linked
compound (17) results in a decrease in potency, comparable
to that observed against LOXL2. Increased selectivity is observed
with 2,4-AMTz regioisomers (22b–d), in particular compound 22d which does not possess
a bis-sulfonyl group. AMT inhibitor 3a is found to be
a weak inhibitor of LOXL3, providing selectivity in excess of 100-fold.
This study demonstrates that 2-aminomethylene-5-sulfonyl thiazoles
are potent inhibitors of three LOX isoforms, while selectivity toward
LOXL2, particularly 2,4-AMTz regioisomers, supports their use as valuable
tool compounds to study the biology and functions of LOXL2. Owing
to the lack of availability and difficulties involved in obtaining
other LOXL enzymes in an active form,[20] we have been unable to assess selectivity of our compounds over
other LOX-family members.
Selectivity over Common Amine Oxidases and
hERG
A selection
of compounds including our most potent AMTz inhibitors were assessed
for their selectivity over the flavin-containing monoamine oxidases
(MAO) A and B, copper-containing diamine oxidase (DAO) and semicarbazide-sensitive
amine oxidase (SSAO), and the hERG channel (Table ). In general, our inhibitors show excellent
selectivity over MAO-A and -B and DAO. Bis-sulfonyl compounds 7a, 18, and 21a, along with biphenyl
compound 7e, demonstrate moderate SSAO inhibition. N-Methyl pyrazolyl inhibitors 21b and 21e show improved selectivity, and further studies carried
out with 21b indicate that it is not a substrate of SSAO,
unlike 3a.[22] Bis-sulfonyl
compounds 7a and 18 are also found to be
moderate inhibitors of the hERG channel, while others compounds assessed
show good selectivity. On the basis of the in vitro activity and selectivity
profiles, 2-aminomethylene-5-sulfonylthiazole inhibitors 6, 7e, 21b, and 21e, along
with 7a for direct comparison with thiophene 3a and 22d as an exemplar of the 2,4-AMTz subseries, were
advanced to metabolic stability assessment and in vivo mouse PK studies.
Pharmacokinetic Evaluation
All compounds assessed in
vivo demonstrate good metabolic stability against mouse liver microsomes
(MLM) and inhibitor exposure following oral administration in mouse
PK studies (Table ). Naphthalenesulfonyl compound 6 shows low to moderate
plasma exposure (AUC = 1–10 μM), as does 3-ethyl-5-phenyl
inhibitor 7e and 21e bearing N-methyl pyrazolyl and fluoro substituents. Pleasingly, a number of
compounds assessed demonstrate desirable PK profiles, achieving greater
plasma exposure levels than that of our previously published inhibitor
(3a),[22] and excellent oral
bioavailability. Indeed, compounds 7a, 21b, and 22d have AUCs above 18 μM·h, achieving Cmax concentrations of up to 32 μM. These
compounds also exhibit good permeability in the Caco-2 assay, used
to model adsorption of orally administered drugs in the small intestine,
and present lower efflux levels than those seen previously for thiophene 3a. Given our observations that these compounds demonstrate
time-dependent inhibition (see Tables and 9), it is likely that efficacy
would be driven by Cmax.
Table 8
In Vitro Mouse Liver Microsome (MLM)
Stability, in Vivo PK Properties, and Caco-2 Permeability of AMTz
Inhibitors
Caco-2 Pappe,f (10–6 cm/s)
compd
MLM stability (%)a
Cmax(PO) (μM)c
AUC(PO) (μM·h)d
t1/2(PO) (h)
F (%)e
A → B
B → A
6
89
6.35
4.20
0.8
nd
nd
nd
21e
90
9.75
4.48
2.0
nd
nd
nd
7e
70
10.16
9.97
0.6
nd
nd
nd
3a
(i) mouse
(i) 64
(i) 16.83
(i) 14.89
(i) 0.45
(i) 45
8.5
35
(ii) rat
(ii) 35b
(ii) 0.34
(ii) 0.22
(ii) 0.5
(ii) 0.2
21b
(i) mouse
(i) 76
(i) 26.71
(i) 18.5
(i) 1.12
(i) 98
31
29
(ii) rat
(ii)
82b
(ii) 6.5
(ii) 18.5
(ii)
1.3
(ii) 68
22d
97
27.22
20.98
1.0
68
23
51
7a
100
32.03
21.61
0.89
nd
23
33
Mouse liver microsome
(MLM) stability
values represent the percentage of compound remaining after 30 min.
Rat liver microsome (RLM) stability
values represent the percentage of compound remaining after 30 min.
Mouse plasma PK parameters were determined following a single po dose
at 50 mg/kg or iv dose at 10 mg/kg. Rat plasma PK parameters were
determined following a single po dose at 20 mg/kg or iv dose at 4
mg/kg.
Cmax: maximum concentration.
AUC: area under curve.
nd: not determined.
Papp: permeability coefficient.
Table 9
Time-Dependent Activity
of Lead AMTz
Inhibitors
compd
LOXL2, 20 min, IC50 (μM)a,b
LOXL2, 1 h, IC50 (μM)a,c
LOXL2, 3 h, IC50 (μM)a,d
BAPN
4.26
0.665
0.372
[3.63, 5.00]
[0.57, 0.77]
[0.31, 0.45]
3a
0.439
0.176
0.157
0.46
[0.105, 0.295]
0.126
7a
0.587
0.086
0.069
0.648
[0.061, 0.12]
0.075
21a
0.544
0.109
0.072
0.331
[0.066, 0.18]
0.069
21b
0.71
0.151
0.078
0.687
[0.12, 0.19]
0.079
22b
1.252
0.302
0.09
0.53
[0.22, 0.42]
0.085
22c
0.795
0.37
0.146
0.787
[0.23, 0.596]
0.157
22d
1.238
0.425
0.199
1.373
[0.35, 0.51]
0.176
Reported IC50 values
were determined in at least two separate experiments (n ≥ 2). When n = 2, individual IC50 values are shown. When n > 2, the values are
reported
as the geometric mean with the corresponding 95% confidence interval
in square brackets.
Compound
was preincubated with enzyme
for 20 min.
Compound was
preincubated with enzyme
for 1 h.
Compound was preincubated
with enzyme
for 3 h.
Mouse liver microsome
(MLM) stability
values represent the percentage of compound remaining after 30 min.Rat liver microsome (RLM) stability
values represent the percentage of compound remaining after 30 min.
Mouse plasma PK parameters were determined following a single po dose
at 50 mg/kg or iv dose at 10 mg/kg. Rat plasma PK parameters were
determined following a single po dose at 20 mg/kg or iv dose at 4
mg/kg.Cmax: maximum concentration.AUC: area under curve.nd: not determined.Papp: permeability coefficient.Reported IC50 values
were determined in at least two separate experiments (n ≥ 2). When n = 2, individual IC50 values are shown. When n > 2, the values are
reported
as the geometric mean with the corresponding 95% confidence interval
in square brackets.Compound
was preincubated with enzyme
for 20 min.Compound was
preincubated with enzyme
for 1 h.Compound was preincubated
with enzyme
for 3 h.Further PK studies
were carried out using AMT compound 3a and AMTz 21b in a rat model. Metabolic stability against
rat liver microsomes (RLM) is moderate in 3a and good
in 21b; however, a significant difference in oral bioavailability
is observed between the compounds, with 3a demonstrating
very poor levels of plasma exposure. In contrast, 21b exhibits comparable AUCs of 18 μM·h between the species
and maintains a good Cmax concentration
(6.5 μM) and oral bioavailability (68%). On the basis of the
balance of potency, selectivity profile, PK, and permeability, compound 21b is determined to have the best profile overall and thus
was chosen for further in vivo antitumor efficacy evaluation.
Evaluation
of Anti-Tumor Efficacy
In vivo efficacy
studies involving compound 21b were carried out using
a genetically engineered mouse model (GEMM) that functions as a LOX-driven
spontaneous breast cancer model.[21] Mice
were dosed daily via oral gavage (70 mg/kg) from around 60 days after
birth, once primary tumors become palpable. Inhibitor 21b was very well tolerated, with no observed body weight loss.Pleasingly, we observe a delay in primary tumor development and a
significant reduction in tumor growth rate in the 21b-treated group compared to that of the controls (Figure a), with no deaths due to tumor
volume reaching ethical size limits necessitated in the inhibitor-treated
group during the course of this study (Figure b).
Figure 3
(a) Antitumor efficacy of compound 21b in MMTV-PyMT
GEMM model: control animals (n = 5; gray) or treated
with compound 21b at 70 mg/kg q.d. (n = 3; green); day 91 *p = 0.0367, Welch’s t test. (b) Kaplan–Meier survival analysis.
(a) Antitumor efficacy of compound 21b in MMTV-PyMT
GEMM model: control animals (n = 5; gray) or treated
with compound 21b at 70 mg/kg q.d. (n = 3; green); day 91 *p = 0.0367, Welch’s t test. (b) Kaplan–Meier survival analysis.
Potency and Mode of Inhibition (MOI) of AMT
and AMTz Inhibitors
We assessed the time-dependent inhibitory
activity of select AMT
and AMTz compounds (Table ) and investigated the mode of enzyme inhibition for key series
examples (Figure ).
Figure 4
LOXL2
jump dilution assay to assess MOI of key LOXL2 inhibitors.
LOXL2
jump dilution assay to assess MOI of key LOXL2 inhibitors.
Time-Dependent Inhibition of Lead Compounds
We previously
demonstrated that series exemplars of AMT and AMTz compounds inhibit
LOXL2 in a time-dependent manner (Table ); we now sought to ascertain if this trend
remained true for the series in general. As such, a range of AMTz
inhibitors including potent 2-aminomethylene-5-sulfonyl compounds 7a, 21a, and 21b, 2,4-regioisomers 22b, 22c, and 22d, along with AMT
compound 3a and BAPN, were assessed for anti-LOXL2 activity
following different enzyme preincubation times (Table ). All compounds were found to exhibit time-dependent
inhibition, with the most significant effect again observed on increasing
the time from 20 min to 1 h, whereupon up to a 7-fold increase in
activity is observed; further increasing the preincubation time to
3 h provides a small additional increase in activity, with IC50 < 0.1 μM attained for all 2-aminomethylene-5-sulfonylthiazoles
(7a, 21a, and 21b) assessed.
MOI Studies
Given the time-dependent inhibition observed
with these compounds, we wanted to clarify the mechanism of LOXL2
inhibition; as such, we set up a jump dilution assay, which readily
distinguishes between reversible and irreversible modes of enzyme
inhibition (Figure ).[37] Leading AMT and AMTz inhibitors 3a and 21b, respectively, along with BAPN were
assessed, whereby the enzyme was preincubated for 1 h with 10 ×
IC50 of these compounds, at which concentration we see
complete inhibition of the enzyme in all cases. The enzyme/inhibitor
mixture was then diluted 100-fold into a solution containing all enzyme
reaction components, and the activity of the enzyme was assessed.The resulting curve obtained for 3a displays about 80%
recovery of activity following dilution, as compared to the DMSO control,
which suggests that this compound behaves as a reversible inhibitor
under these assay conditions, unlike BAPN which can be characterized
as an irreversible inhibitor. In the case of 21b a more
modest regain in activity of around 30% is observed following the
jump dilution. This is suggestive of either a slowly reversible compound
or an irreversible inhibitor whereby a residual amount of enzyme has
not committed to forming a stable covalent bond during the two-step
inhibition process. The results obtained from this study suggest that
the electronic nature of the heterocyclic core affects the binding
mechanism of these inhibitors, with the AMTz inhibitors better able
to form stable, irreversible Schiff base intermediates following initial
reversible binding to the LTQ cofactor.[20,38] One possible
explanation for this is that the initial Schiff base formed is more
susceptible to hydrolysis in the case of 3a compared
with 21b, which is better able to rearrange to form a
more stable intermediate as a result of either increased resonance
stabilization or noncovalent electrostatic interactions. Alternatively,
this could also be the result of differences in enzyme–inhibitor
kinetics, with initial rates of binding affecting the ability of these
inhibitors to form covalent bonds.
Pharmacophore and QSAR
Modeling
Binding Mode/Pharmacophore Hypothesis Generation
With
the dearth of protein–ligand crystal structures for LOXL2 in
its active form, we decided to embark upon a ligand-based approach
to propose a legitimate binding conformation for SAR analysis and
modeling. FieldTemplater[39,40] was used to derive
a pharmacophore model by comparing conformational ensembles of molecules
using their electrostatic and hydrophobic characters to identify common
motifs. The field point pattern for a conformer of 7a is shown in Figure .
Figure 5
Sample conformation of 7a showing electrostatic fields
and 3D field point pattern, which provides a concise 3D pharmacophore.
Sample conformation of 7a showing electrostatic fields
and 3D field point pattern, which provides a concise 3D pharmacophore.The field-based alignment can be independent of
chemical structure,
allowing alignment of molecules from different series. The compounds
reported in this article are thus augmented with a set of LOXL2 inhibitors
obtained from the literature,[20] along with
in-house inhibitors not explicitly described in this article, with
only those compounds with pIC50 > 5 included (see Supporting Information). The collection of 54
molecules (with IC50 data) were visually inspected, with
the four most active, most structurally diverse compounds selected
for pharmacophore modeling (7a, pIC50 = 7.07; 21e, pIC50 = 6.84; JMC2017-31,[20] pIC50 = 6.55; and JMC2017-33,[20] pIC50 = 6.51).The FieldTemplater experiment
was run using Normal (large molecules)
conformation hunt settings, with a group constraint placed on the
cationic NH3+ groups to force their alignment
in the templated result, since these compounds are all believed to
follow the same mechanism of inhibition via covalent binding of the
aminomethylene group to the LTQ cofactor. The best scoring template
consisting of all four compounds was taken forward to Forge,[40,41] and the 7a and 21e conformations were
used as references for the field-based alignment of the other molecules
in the data set.
Structure–Activity Relationship Modeling
We
attempted to calculate statistically relevant mathematical models
to predict activities of new compounds. Pleasingly, we were rewarded
with two predictive, complementary models as described below.
Field-Based
3DQSAR[40,41]
Cresset’s
approach to 3DQSAR is similar to traditional CoMFA;[42,43] however, there are some striking differences around how sampling
points are selected and the use of irregular grids, such that calculation
speed is greatly improved. The 54 compound data set was randomly partitioned
to put 15% in the test set (8 molecules), leaving 46 in the training
set, upon which the model was built; tested with 50 y-scrambles, followed
by leave-one-out (LOO) cross-validation. The resulting three-component
model features an r2 of 0.856, q2 of 0.665, RMSE of 0.147, and a Kendall’s
tau[44] value of 0.760, a model we believe
to be statistically relevant[45] and able
to predict rank order of activity.Forge can be used to visualize
the field/steric contributions to predicted activities, as demonstrated
in Figure c; these
plots were helpful in rationalizing the SAR observed for the reported
compounds. The details of this field-based 3DQSAR model are shown
in Figure a and Figure b, with the model
details included in the Supporting Information.
Figure 6
(a) Plot of predicted vs actual activity for training, training
cross-validation, and test sets. (b) 7a with the electrostatic
and steric coefficient positions/values displayed. (c) Example of
the field/steric contributions to predicted activity for 7a.
(a) Plot of predicted vs actual activity for training, training
cross-validation, and test sets. (b) 7a with the electrostatic
and steric coefficient positions/values displayed. (c) Example of
the field/steric contributions to predicted activity for 7a.As a follow-up and to build confidence
in the model and binding
mode hypothesis, an alternative set of machine-learning algorithms,
k-nearest neighbor (kNN), random forest (RF), support vector machine
(SVM), and relevance vector machine (RVM) regression models were used
to model the data.[41] The statistics for
RVM for the full model on the training set were r2 = 0.852, RMSE = 0.150, and Kendall’s tau = 0.757.
For the cross-validation, the model statistics were q2 = 0.649, RMSE = 0.230, Kendall’s tau = 0.638.
This model is similarly predictive to the field-based 3DQSAR model
described earlier, and the calculation of multiple predictive models
from a set of molecules aligned to a common binding mode hypothesis
lends support to our having derived a sensible binding mode in the
absence of protein–ligand crystallographic information. While
not reported here, it is noteworthy that the other machine learning
models also have statistics that suggest that they are predictive.
Activity Atlas Analysis: Activity Cliff Summary
In
tandem with the predictive 3DQSAR model, Activity Atlas/Activity Miner[40,41] was used to conduct a qualitative assessment of the SAR for the
data set. The summary plots of the calculated electrostatic and steric
activity cliffs are shown in Figure and compare favorably to the 3DQSAR model visualizations.
This information is built up from doing a pairwise analysis of all
molecules in their aligned conformation and automatically examining
activity cliffs, highly similar pairs where there is a large change
in activity. The most potent and least potent compounds are shown
in the context of the activity cliff summary.
Figure 7
Activity Atlas/activity
cliff summary plots for LOXL2 activity
shown with (a) 7a around thiophene/thiazole core, (b) 7a, 7e, and 7g to illustrate steric
and electrostatic contributions (using field points, as described
above) that rationalize the observed SAR. All isosurfaces are shown
at ≥2.0 confidence level.
Activity Atlas/activity
cliff summary plots for LOXL2 activity
shown with (a) 7a around thiophene/thiazole core, (b) 7a, 7e, and 7g to illustrate steric
and electrostatic contributions (using field points, as described
above) that rationalize the observed SAR. All isosurfaces are shown
at ≥2.0 confidence level.The predicted SAR around the thiazole moiety (Figure a) is consistent with the observed
decrease in potency when the thiazolenitrogen is replaced by the
larger CX groups (3b, 3c), suggesting that
this position is sterically restricted and requires negative electrostatics.
In addition, the requirement for negative electrostatics is consistent
with the moderate decrease in potency of the thiophene analogues (3a) and the 2,4 substituted thiazoles (22b).
Examination of the SAR around the thiazole C-5 substituent using Activity
Atlas suggests that the meta-position of the phenyl
group (R1) is predicted to favor large groups; in 21b, replacement of phenyl with N-methyl
pyrazolyl substituent meets both the desired steric and electrostatic
conditions near the 5-position. With regard to the 3-position, the
activity cliff summary (Figure b) can be used to rationalize the drop in potency on varying
the methanesulfonyl group of 7a to an ethyl (7e) or tert-butyl (7g) group, which increases
the number of cliff violations (i.e., a mismatching of field points
with activity cliff summary): In 7a, the 3-sulfonyl group
has favorable electrostatics and is able to thread the needle of a
sterically unfavorable region; a larger tert-butyl
group has unfavorable electrostatics and multiple steric clashes that
result in diminished activity. Further, in the context of the model,
the central phenyl group is predicted to have a less electron-rich
π-cloud, in line with the observed improvement in potency when
small electron-withdrawing groups are present in the 3-position.We have developed predictive field-based 3DQSAR and visual qualitative
models of activity based on the LOXL2 inhibition data acquired during
the development of the aminomethylenethiazole inhibitors. This
will be a useful tool to aid further development of the series as
well as to design new chemical inhibitors in the future.
Synthetic
Chemistry
Sulfone-linked AMTz analogues were synthesized
from N-Boc-protected (5-bromothiazol-2-yl)methanamine 23 and
the appropriate thiol, using palladium-catalyzed coupling methods,
as shown in Scheme . In the case of aryl thiols, tris(dibenzylideneacetone)dipalladium
catalyst and XantPhos ligand were used in conjunction with sodium tert-butoxide base in a 4:1 solvent mixture of toluene/tert-butanol. In the case of alkyl thiols, DIPEA was found
to be preferable as a base and toluene was used as the solvent. Oxidation
of the sulfide was achieved using m-CPBA to afford
the desired sulfone, and subsequent Boc-removal using 4 M HCl in dioxane
gave rise to the target AMTz inhibitors.
Scheme 1
General Route to
Sulfur-Linked 2,5-AMTz Analogues
Reagents and conditions:
(a)
Pd2dba3, XantPhos, NaOBu, PhMe, BuOH, 110 °C; (b)
Pd2dba3, XantPhos, DIPEA, PhMe, 110 °C;
(c) m-CPBA, DCM, rt; (d) 4 M HCl in dioxane, rt.
General Route to
Sulfur-Linked 2,5-AMTz Analogues
Reagents and conditions:
(a)
Pd2dba3, XantPhos, NaOBu, PhMe, BuOH, 110 °C; (b)
Pd2dba3, XantPhos, DIPEA, PhMe, 110 °C;
(c) m-CPBA, DCM, rt; (d) 4 M HCl in dioxane, rt.Substituted aryl thiols that were not commercially
available were
synthesized according to the method shown in Scheme .[46] Palladium-catalyzed
coupling of an aryl halide with 2-ethylhexyl 3-mercaptopropanoate
afforded a thiol surrogate that also functions as a thiol protecting
group for further chemical modification. The thiol surrogate could
then be used directly in a palladium-catalyzed coupling reaction,
whereby deprotection was achieved in situ, or the protecting group
can be removed using sodium ethoxide to afford the desired aryl thiol.
Scheme 2
General Synthesis of Noncommercially Available Thiols
Reagents
and conditions: (a)
2-ethylhexyl 3-mercaptopropanoate, Pd2dba3,
XantPhos, DIPEA, PhMe, 110 °C; (b) NaOEt, PhMe, EtOH, rt; then
acid; (c) KOBu, THF, rt; then acid.
General Synthesis of Noncommercially Available Thiols
Reagents
and conditions: (a)
2-ethylhexyl 3-mercaptopropanoate, Pd2dba3,
XantPhos, DIPEA, PhMe, 110 °C; (b) NaOEt, PhMe, EtOH, rt; then
acid; (c) KOBu, THF, rt; then acid.AMTz analogues bearing aryl groups directly attached
to the thiazole
ring were synthesized according to the method shown in Scheme . Suzuki reaction of N-Boc-protected (5-bromothiazol-2-yl)methanamine (23) with a boronic acid or ester, followed by acid-mediated
deprotection, afforded the desired phenyl-linked target compounds.
Scheme 3
General Route to Phenyl-Linked AMTz Analogues
Reagents
and conditions: (a)
Pd(PPh3)4, K2CO3, 1,2-DME,
H2O, 95 °C; (b) 4 M HCl in dioxane, rt.
General Route to Phenyl-Linked AMTz Analogues
Reagents
and conditions: (a)
Pd(PPh3)4, K2CO3, 1,2-DME,
H2O, 95 °C; (b) 4 M HCl in dioxane, rt.Synthesis of the thiol surrogate used in the synthesis
of 7a and close analogues was achieved starting from
1,3-dibromo-5-(methylsulfonyl)benzene, Scheme . Palladium cross-coupling
with 2-ethylhexyl 3-mercaptopropanoate installed the protected thiol
group and Suzuki coupling with phenyl boronic acid at the remaining
bromo-position completed the synthesis of the required thiol reactant.
In situ deprotection followed by coupling with a 5-bromo-2-AMTz intermediate
(23 or 24) resulted in the corresponding
sulfide-linked compound. Direct N-Boc-deprotection
provided the sulfide-linked final compound 7d, or oxidation
using either 1 or 2 equiv of m-CPBA prior to deprotection
resulted in the sulfoxide or sulfone-linked inhibitors, 7c and 7a or 7b, respectively.
Scheme 4
AMTz Sulfide,
Sulfoxide, and Aminomethylene Modifications
AMTz Sulfide,
Sulfoxide, and Aminomethylene Modifications
Reagents
and conditions: (a)
2-ethylhexyl 3-mercaptopropanoate, Pd2dba3,
XantPhos, DIPEA, PhMe, 110 °C; (b) PhB(OH)2, Pd(PPh3)4, K2CO3, 1,2-DME, H2O, 95 °C; (c) Pd2dba3, XantPhos,
NaOBu, PhMe, BuOH, 110 °C; (d) 4 M HCl in dioxane, rt; (e) m-CPBA (1 equiv), DCM, 0 °C; (f) m-CPBA (2.2
equiv), DCM, rt.Thiophene analogues 3b and 3c, bearing
a halogen atom in the 3-position, were synthesized from the corresponding
3-halo-5-bromo-AMT precursors (25 and 26), which were synthesized according to literature procedures,[47,48]Scheme . Palladium
coupling was then carried out with the thiol surrogate, and m-CPBA oxidation of the resulting bis-sulfide followed by
treatment with HCl afforded the desired halogenated thiophene analogues.
Scheme 5
General Route to Sulfur-Linked AMT Analogues
Reagents
and conditions: (a)
NaOEt, PhMe, EtOH, rt; then acid; (b) KOBu, THF, rt; then acid; (c) Pd2dba3, XantPhos,
NaOBu, PhMe, BuOH, 110 °C; (d) m-CPBA, DCM, rt; (e) 4 M
HCl in dioxane, rt.
General Route to Sulfur-Linked AMT Analogues
Reagents
and conditions: (a)
NaOEt, PhMe, EtOH, rt; then acid; (b) KOBu, THF, rt; then acid; (c) Pd2dba3, XantPhos,
NaOBu, PhMe, BuOH, 110 °C; (d) m-CPBA, DCM, rt; (e) 4 M
HCl in dioxane, rt.Synthesis of the 2,4-thiazole
regioisomers was achieved in the
same manner as described previously, starting from N-Boc-protected (4-bromothiazol-2-yl)methanamine (27),[49]Scheme . Palladium-catalyzed coupling with the respective thiol followed
by subsequent oxidation using m-CPBA and Boc-deprotection
using HCl afforded the desired regioisomers.
Scheme 6
General Route to
Sulfur-Linked 2,4-AMTz Regioisomer Analogues
Reagents
and conditions: (a)
Pd2dba3, XantPhos, NaOBu, PhMe, BuOH, 110 °C; (b) m-CPBA, DCM, rt; (c) 4 M HCl in dioxane, rt.
General Route to
Sulfur-Linked 2,4-AMTz Regioisomer Analogues
Reagents
and conditions: (a)
Pd2dba3, XantPhos, NaOBu, PhMe, BuOH, 110 °C; (b) m-CPBA, DCM, rt; (c) 4 M HCl in dioxane, rt.
Conclusion
Following our recent discovery of a potent,
selective, and orally
bioavailable LOX inhibitor we observed that replacement of the thiophene
core with a 2-aminomethylene-5-sulfonyl thiazole core leads to potent,
irreversible LOXL2 inhibitors. SAR investigations revealed similar
trends as seen in the analogous AMT series, resulting in potencies
of <0.1 μM achieved, and enabled development of a predictive
LOXL2 3DQSAR model. Selectivity studies concerning LOX and LOXL3 isoforms
revealed a modest selectivity toward LOXL2 in our main series inhibitors,
while 2-aminomethylene-4-sulfonyl thiazole regioisomers exhibit excellent
selectivity for LOXL2 and thus have the potential to be used as probe
compounds. Further selectivity and ADME assessment leads to the discovery
of 21b, which has an improved PK profile and demonstrates
excellent antitumor efficacy in a LOX-driven GEMM breast cancer model.
Experimental Section
Synthesis of Inhibitors
All chemicals, reagents, and
solvents were purchased from commercial sources and were used as received.
Flash chromatography was performed on a Biotage Isolera or Combiflash
Rf+ UV–vis flash purification system using prepacked silica
gel cartridges (Biotage) with HPLC grade solvents. Thin layer chromatography
(TLC) analysis was performed using silica gel 60 F-254 thin layer
plates and visualized using UV light (254 nm) and/or developed with
vanillin stain. LCMS and HRMS analyses of chemical compounds were
performed on an Agilent 1200 series HPLC and diode array detector
coupled to a 6210 time-of-flight mass spectrometer with a multimode
ESI source; or a Waters Acquity UPLC or I-class UPLC with a diode
array detector coupled to a Waters G2 QToF, SQD, or QDa mass spectrometer
fitted with a multimode ESI/APCI source. 1H, 19F, and 13C NMR spectra were recorded using a Bruker Avance
500, 400, or 300 MHz spectrometer using an internal deuterium lock.
Chemical shifts are expressed in parts per million (ppm), and splitting
patterns are indicated as follows: br, broad; s, singlet; d, doublet;
t, triplet; q, quartet; p, pentet; h, hextet; m, multiplet. All coupling
constants (J) are reported in hertz (Hz). All final
inhibitors submitted for biological evaluation were at least 95% pure
by HPLC–MS. Synthesis of inhibitors 3a, 4, 5, 6, and 8 has
previously been described in the literature.[22] Below are a representative synthesis of compound 7a and analytical data for all final inhibitors. All tested inhibitors
have purity of >95% (LCMS/UV).
A solution of 1,3-dibromo-5-(methylsulfonyl)benzene (500
mg, 1.59 mmol) and DIPEA (0.55 mL, 3.18 mmol) in toluene (10 mL) was
degassed with nitrogen for 5 min. Pd2dba3 (36
mg, 2.5 mol %), XantPhos (46 mg, 5 mol %), and 2-ethylhexyl-3-mercaptopropionate
(0.36 mL, 1.59 mmol) were then added with stirring, and the solution
was bubbled with nitrogen for a further 5 min before sealing the flask
and heating to 110 °C, with stirring for 18 h. After cooling
to room temperature, the reaction mixture was diluted with ethyl acetate
(60 mL) and washed with water (50 mL) and brine (50 mL), dried over
anhydrous magnesium sulfate, filtered, and concentrated to give the
crude product, which was purified using flash column chromatography
(5–100% EtOAc in PE) to give the title compound (664 mg, 72%
purity) as a clear pale yellow oil, which is used in the subsequent
transformation as an impure mixture. 1H NMR (CDCl3, 300 MHz) δ 7.84 (t, 1H, J = 1.6 Hz), 7.76
(t, 1H, J = 1.6 Hz), 7.66 (t, 1H, J = 1.7 Hz), 4.04–4.01 (m, 2H), 3.26 (t, 2H, J = 7.2 Hz), 3.06 (s, 3H), 2.67 (t, 2H, J = 7.2 Hz),
1.61–1.51 (m, 1H), 1.39–1.25 (m, 8H), 0.91–0.85
(m, 6H) ppm; LCMS, did not ionize.
2-Ethylhexyl 3-((3-bromo-5-(methylsulfonyl)phenyl)thio)propanoate
(1.80 g, 72% purity, 2.88 mmol) and phenylboronic acid (422 mg, 3.46
mmol) were dissolved in 1,2-DME/H2O (5:1, 24 mL), and the
solution was bubbled with nitrogen for 5 min. Pd(PPh3)4 (333 mg, 10 mol %) and K2CO3 (796 mg,
5.76 mmol) were then added with stirring, and the mixture was bubbled
with nitrogen for a further 5 min before sealing the flask and heating
at 100 °C with stirring for 18 h. After cooling to room temperature
the reaction mixture was diluted with ethyl acetate (60 mL) and washed
with brine (50 mL), dried over anhydrous magnesium sulfate, filtered,
and concentrated to give the crude product, which was purified using
flash column chromatography (5–80% EtOAc in PE) to give the
title compound (998 mg, 72% purity) as a clear pale yellow oil, which
is used in the subsequent transformation as an impure mixture. 1H NMR (CDCl3, 500 MHz) δ 7.94 (t, 1H, J = 1.6 Hz), 7.84 (t, 1H, J = 1.7 Hz),
7.77 (t, 1H, J = 1.7 Hz), 7.60–7.58 (m, 2H),
7.48 (t, 2H, J = 7.4 Hz), 7.43 (t, 1H, J = 7.3 Hz), 4.06–3.99 (m, 2H), 3.31 (t, 2H, J = 7.2 Hz), 3.10 (s, 3H), 2.70 (t, 2H, J = 7.2 Hz),
1.59–1.55 (m, 1H), 1.37–1.31 (m, 2H), 1.30–1.25
(m, 6H), 0.88 (t, 6H, J = 7.5 Hz) ppm. LCMS, did
not ionize.
A solution of 2-ethylhexyl 3-((5-(methylsulfonyl)-[1,1′-biphenyl]-3-yl)thio)propanoate
(830 mg, 72% purity, 1.34 mmol), tert-butyl ((5-bromothiazol-2-yl)methyl)carbamate
(431 mg, 1.47 mmol) in toluene/BuOH (4:1,
15 mL) was degassed with nitrogen for 5 min. Pd2dba3 (128 mg, 10 mol %), XantPhos (162 mg, 20 mol %), and NaOBu (283 mg, 2.94 mmol) were then added with
stirring, and the solution was bubbled with nitrogen for a further
5 min before sealing the flask and heating to 110 °C, with stirring
for 18 h. After cooling to room temperature, the reaction mixture
was diluted with ethyl acetate (60 mL) and washed with water (50 mL)
and brine (50 mL), dried over anhydrous magnesium sulfate, filtered,
and concentrated to give the crude product, which was purified using
flash column chromatography (5–100% EtOAc in PE) followed by
reversed-phase chromatography (C18 silica, 5–95% CH3CN in water) to give the title compound (108 mg, 17%) as a yellow
oil. 1H NMR (CDCl3, 500 MHz) δ 7.94 (t,
1H, J = 1.6 Hz), 7.88 (s, 1H), 7.73 (t, 1H, J = 1.7 Hz), 7.63 (t, 1H, J = 1.6 Hz),
7.53–7.51 (m, 2H), 7.48–7.45 (m, 2H), 7.44–7.41
(m, 1H), 5.29 (br s, 1H), 4.62 (d, 2H, J = 6.0 Hz),
3.07 (s, 3H), 1.45 (s, 9H) ppm; LCMS m/z 421.0358 found (M – Bu + H)+, 421.0345 calculated for C18H17N2O4S3.
m-CPBA (77%, 41 mg, 0.184 mmol) was added over
2 min to a solution of tert-butyl ((5-((5-(methylsulfonyl)-[1,1′-biphenyl]-3-yl)thio)thiazol-2-yl)methyl)carbamate
(35 mg, 0.073 mmol) in dichloromethane (3 mL), with stirring at room
temperature under air for 18 h. The reaction mixture was then quenched
with aq NaHCO3 solution and extracted with dichloromethane
(2 × 20 mL). The combined organic extracts were washed with brine
(30 mL), dried over anhydrous sodium sulfate, filtered, and concentrated
to give the crude product, which was purified using flash column chromatography
(25–75% EtOAc in cyclohexane) to give the title compound (27
mg, 72% yield) as a clear pale yellow oil. 1H NMR (CDCl3, 500 MHz) δ 8.48 (t, 1H, J = 1.7 Hz),
8.41 (t, 1H, J = 1.7 Hz), 8.36 (t, 1H, J = 1.6 Hz), 8.28 (s, 1H), 7.63–7.61 (m, 2H), 7.55–7.49
(m, 3H), 5.27 (br s, 1H), 4.60 (d, 2H, J = 6.0 Hz),
3.14 (s, 3H), 1.44 (s, 9H) ppm; LCMS m/z 453.0214 found (M – Bu + H)+, 453.0243 calculated for C18H17N2O6S3.
All compounds were imported to Forge
and subjected to a field-based alignment to the reference structure
(the FieldTemplater model described in the article), using maximum
common substructure to guide the alignment; a similarity score was
calculated as the average of the field and shape similarity scores.[40] The calculated alignments were visually inspected
to ensure best alignment and adjusted as required.
PAINS Assessment
To identify reactive compounds that
might exhibit interference in biochemical assays, the PAINS filters
as described by Baell and Holloway[50] were
curated as SMARTS and scripted as a flagging protocol deployed in
Vortex (version 2018.09.76561.53-s, 2018, https://www.dotmatics.com/products/vortex) and Pipeline Pilot (Dassault Systèmes BIOVIA, BIOVIA Pipeline Pilot, release 2018, San Diego,
Dassault Systèmes, 2018). The 480 patterns were used to recognize
structures that may result in nonspecific binding to multiple biological
targets by virtue of comprising one or more fragments established
to be of concern. No LOXL2 inhibitor in this study showed any potential
PAINS liability when screened against this PAINS filter.
LOX Protein
Preparation and Enzyme Assays
LOX enzyme
was extracted from pig skin by the method of Shackleton and Hulmes.[51] LOXL2 and LOXL3 were purchased from R&D
Systems. LOX, LOXL2, and LOXL3 catalytic activity were determined
using the Promega ROS-Glo assay kit with cadaverine dihydrochloride
as substrate, BAPN as the reference inhibitor control, a preincubation
time of 20 min, 1 h, or 3 h, with nine dilutions from a top concentration
of 10 μM or 100 μM.
LOXL2 Jump Dilution Assay
LOXL2 catalytic activity
was determined using the Amplex Red hydrogen peroxide assay kit with
cadaverine dihydrochloride as substrate. LOXL2 at 100-fold final assay
concentration and compound at 10 × IC50 were preincubated
for 1 h. The enzyme/inhibitor mixture was then diluted 100-fold into
a solution containing substrate and detection reagents and read kinetically
every 5 min.
Amine Oxidase Assays
All amine oxidase
assays were
performed with concentrations as above. MAO-A and MAO-B enzymes were
purchased from Sigma. The catalytic activity of MAO-A and MAO-B was
determined using the Promega MAO-Glo assay kit (substrate included),
with clorgyline and deprenyl as reference inhibitor controls, respectively.
DAO was purchased from Sigma, and the catalytic activity was determined
using the Promega ROS-Glo assay kit, with aminoguanidine as the reference
inhibitor control. SSAO was purchased from Sigma. SSAO catalytic activity
was determined using the Promega MAO-Glo assay kit, with mofegiline
as the reference inhibitor control.
Assessment of Compound 21b as a Substrate for Amine
Oxidases
The catalytic activities of MAO-A, MAO-B, and SSAO
with compound 21b as a substrate were determined using
the respective enzymes described above, and the hydrogen peroxide
produced was quantified using an Amplex red monoamine oxidase assay
kit. p-Tyramine was used as the positive substrate
control for MAO-A and MAO-B, and benzylamine was used for SSAO.
MLM Stability Assay
Mouse liver microsomes (CD1 female;
M1500) and rat liver microsomes (Sprague Dawley female; R1500) were
purchased from Tebu-bio, and the assay was performed by methods previously
described.[21] Inhibitors at 10 μM
concentration incubated with the microsomes were assessed at 0, 15,
and 30 min. Control samples containing no microsomes and no cofactors
were also assessed at 0 and 30 min. Samples were extracted by protein
precipitation, and centrifugation for 20 min in a refrigerated centrifuge
(4 °C) at 3700 rpm. The supernatant was analyzed by LC–MS/MS
for % metabolized over time.
Animal Procedures
All procedures
involving animals
were performed under licenses PPL-70/7635, 70/7701, and PE3DF1A5B
and National Home Office regulations under the Animals (Scientific
Procedures) Act 1986. Procedures were approved by the Animal Welfare
and Ethical Review Bodies (AWERB) of the CRUK Manchester Institute
and the Institute of Cancer Research and reported in accordance with
ARRIVE guidelines. All mice and rats were maintained in pathogen-free,
ventilated cages in the Biological Resources Unit at Cancer Research
UK Manchester Institute and the Biological Services Unit at The Institute
of Cancer Research. All mice and rats were allowed free access to
irradiated food and autoclaved water in a 12 h light/dark cycle with
room temperature at 21 ± 2 °C. All cages contained wood
shavings, bedding, and a cardboard tube for environmental enrichment.
PyMT-driven breast cancer model mice (FVB background) were bred in
a specific pathogen-free facility at The University of Manchester
(U.K.) under a Home Office approved license.
In Vivo PK
Female
Balb/C or CD1nude mice (Charles
River Laboratories) at 6 weeks of age were used for the mouse PK analyses.
The mice were dosed orally by gavage (50 mg kg–1 in DMSO/water 1:19 v:v; n = 21) or intravenously
in the tail vein (10 mg kg–1 in DMSO/Tween20/saline
10:1:89 v:v:v; n = 24). Blood samples were taken
at seven (po) or eight (iv) time-points between 5 min and 24 h, after
one single dose of the inhibitor. Three mice were used per time-point
per route; average values are reported. They were placed under halothane
or isoflurane anesthesia, and blood for plasma preparation was taken
by terminal cardiac puncture into heparinized syringes. Female Sprague
Dawley (Charles River Laboratories) weighing between 170 g and 200
g were used for the rat PK analyses. The rats were dosed orally by
gavage (20 mg kg–1 in DMSO/water 1:19 v:v; n = 21) or intravenously in the tail vein (4 mg kg–1 in DMSO/Tween20/saline 10:1:89 v:v:v; n = 24).
Blood samples were taken at five (po and (iv) time-points between
5 min and 8 h, after 1 single dose of the inhibitor. One rat was used
per route with serial bleeds taken through the time points. They were
placed in a heated box for 10 min prior to sampling to increase vasodilation,
and blood for plasma preparation was taken by tail vain bleed into
heparinized tubes. Plasma samples, obtained after blood spun at 1300
rpm for 3 min, were pipetted into cryovials and immediately snap frozen
in liquid nitrogen and then stored at −80 °C prior to
analysis.
In Vivo Antitumor Efficacy
LOX inhibitor
treatment
was carried out in a genetically engineered MMTV-PyMT driven mousebreast cancer model where female mice were randomized as described
previously.[21] Mice were treated daily by
oral gavage with 70 mg/kg compound 21b (n = 3) in vehicle (5% DMSO/2.5% Tween20 in water), and controls (n = 5) received vehicle alone or were untreated. Oral administration
of 21b was initiated at day 57 (n =
1) and day 61 (n = 2) with all treatments continuing
for 34 days. The spontaneous breast tumors arising in the model were
measured twice weekly. All of the controls bar one were culled due
to large tumor size by day 95. All of the treated were culled at the
end of treatment at day 91 (n = 1) and day 95 (n = 2); none were culled due to reaching license limit tumor
volumes. Statistical significance was calculated using Welch’s t test on day 91 (*p = 0.0367), utilizing
the final measured tumor volume of the culled control mice and linear
interpolation for the two remaining control mice (between measurements
made on days 89 and 92). All animals allocated to the study were used.
Commercial ADME-T Services
hERG inhibition was determined
using the “hERGhuman potassium ion channel cell based antagonist
Qpatch assay” by Eurofins Ltd. Cell permeability was determined
using the “Caco-2 permeability assay” by Cyprotex Ltd.
Authors: J Molnar; K S K Fong; Q P He; K Hayashi; Y Kim; S F T Fong; B Fogelgren; K Molnarne Szauter; M Mink; K Csiszar Journal: Biochim Biophys Acta Date: 2003-04-11
Authors: Vincent M Crowley; Anuj Khandelwal; Sanket Mishra; Andrew R Stothert; Dustin J E Huard; Jinbo Zhao; Aaron Muth; Adam S Duerfeldt; James L Kizziah; Raquel L Lieberman; Chad A Dickey; Brian S J Blagg Journal: J Med Chem Date: 2016-04-04 Impact factor: 7.446