Christopher R M Asquith1,2, James M Bennett3, Lianyong Su4, Tuomo Laitinen5, Jonathan M Elkins3,6, Julie E Pickett1, Carrow I Wells1, Zengbiao Li4, Timothy M Willson1, William J Zuercher1,7. 1. Structural Genomics Consortium, UNC Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA. 2. Department of Pharmacology, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA. 3. Structural Genomics Consortium and Target Discovery Institute, Nuffield Department of Clinical Medicine, University of Oxford, Oxford OX3 7DQ, UK. 4. Drumetix Laboratories, Greensboro, NC 27409, USA. 5. School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, 70211 Kuopio, Finland. 6. Structural Genomics Consortium, Universidade Estadual de Campinas-UNICAMP, Campinas, São Paulo 13083-886, Brazil. 7. Lineberger Comprehensive Cancer Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA.
Abstract
SGC-GAK-1 (1) is a potent, selective, cell-active chemical probe for cyclin G-associated kinase (GAK). However, 1 was rapidly metabolized in mouse liver microsomes by cytochrome P450-mediated oxidation, displaying rapid clearance in liver microsomes and in mice, which limited its utility in in vivo studies. Chemical modifications of 1 that improved metabolic stability, generally resulted in decreased GAK potency. The best analog in terms of GAK activity in cells was 6-bromo-N-(1H-indazol-6-yl)quinolin-4-amine (35) (IC50 = 1.4 μM), showing improved stability in liver microsomes while still maintaining a narrow spectrum activity across the kinome. As an alternative to scaffold modifications we also explored the use of the broad-spectrum cytochrome P450 inhibitor 1-aminobenzotriazole (ABT) to decrease intrinsic clearance of aminoquinoline GAK inhibitors. Taken together, these approaches point towards the development of an in vivo chemical probe for the dark kinase GAK.
SGC-GAK-1 (1) is a potent, selective, cell-active chemical probe for cyclin G-associated kinase (GAK). However, 1 was rapidly metabolized in mouse liver microsomes by cytochrome P450-mediated oxidation, displaying rapid clearance in liver microsomes and in mice, which limited its utility in in vivo studies. Chemical modifications of 1 that improved metabolic stability, generally resulted in decreased GAK potency. The best analog in terms of GAK activity in cells was 6-bromo-N-(1H-indazol-6-yl)quinolin-4-amine (35) (IC50 = 1.4 μM), showing improved stability in liver microsomes while still maintaining a narrow spectrum activity across the kinome. As an alternative to scaffold modifications we also explored the use of the broad-spectrum cytochrome P450 inhibitor 1-aminobenzotriazole (ABT) to decrease intrinsic clearance of aminoquinolineGAK inhibitors. Taken together, these approaches point towards the development of an in vivo chemical probe for the dark kinase GAK.
Entities:
Keywords:
1-aminobenzotriazole (ABT); cyclin G associated kinase (GAK); in vivo kinase design; kinase inhibitor design
Cyclin G-associated kinase (GAK) is a 160 kDa member of the numb-associated kinase (NAK) family of serine/threonine kinases [1]. GAK was originally identified as a directly associated partner of cyclin G, and is ubiquitously expressed across tissues [2]. Within the cell, GAK localizes to the Golgi complex, cytoplasm, and nucleus [3]. GAK has been genetically associated with a diverse range of biological processes. Genome-wide association studies have identified single nucleotide polymorphisms in GAK that are associated with susceptibility to Parkinson’s disease [4]. Knock down of GAK has shown that it is required for the maintenance of centrosome maturation and progression through mitosis [5]. Conversely, GAK was found to be overexpressed in osteosarcoma cell lines and tissues, where it contributes to proliferation and survival [6].The biological implications of GAK depletion or over-expression prompt the question of whether there is therapeutic utility in targeting its kinase domain with small molecules. GAK is a collateral target of several clinical kinase inhibitors, including the epidermal growth factor receptor (EGFR) inhibitors erlotinib and gefitinib (Figure 1), but until recently potent selective cell-active inhibitors of GAK had not been identified [7] and thus the biological consequences of selective GAK inhibition remained untested. As part of an ongoing effort to generate chemical probes for dark kinases, we recently described SGC-GAK-1 (1) as a high-quality cell-active chemical probe for GAK (Figure 1) [8,9]. SGC-GAK-1 (1) showed IC50 = 48 nM in a GAK live cell target engagement assay and proved to be more than 50-fold selective for GAK across a panel of over 400 human kinases.
Figure 1
Examples of reported GAK inhibitors.
GAK is a potential therapeutic target for the treatment of androgen-independent prostate cancer. Expression of GAK is known to increase during prostate cancer progression to androgen independence and has been positively correlated with Gleason scores in resections from prostate cancerpatients [10]. We have demonstrated that 1 potently inhibits the viability of prostate cancer cells that contain constitutively-active splice variants of the androgen receptor [9]. These results prompted us to identify a chemical probe of GAK that is suitable for use in rodent models of prostate cancer. We now describe our efforts in the optimization of 4-aminoquinolines as potential in vivo chemical probes for GAK.
2. Results
2.1. Initial Investigation of P450 Metabolism of 1
The metabolic stability of 1 was evaluated in mouse liver microsomes (MLMs) to assess the potential for use as an in vivo chemical probe. SGC-GAK-1 (1) was incubated with MLMs in the presence or absence of added NADPH, and the compound was quantified at various time points over an hour. No degradation of 1 was observed in the absence of NADPH. However, 1 was rapidly consumed (T1/2 = 5.7 min; Clint = 990 mL/min/kg) in the presence of NADPH. This result suggests that 1 would be prone to rapid first pass metabolism when dosed in vivo.The necessity of NADPH for the metabolism of 1 is consistent with cytochrome P450-mediated oxidation. To further understand the metabolic fate of 1, we conducted a metabolite identification experiment. 10 µM of 1 was incubated for 5 min (ratio of 0.5 with respect to propranolol, n = 4) in MLMs at 37 °C in the presence of an excess of NADPH before quenching with methanol. After centrifugation, the supernatant was subjected to analysis by LCMS, and structures of metabolites were assigned on the basis of their molecular weight and fragmentation pattern. The primary metabolite 2 was identified as a loss of CH2 from the 3,4,5-trimethoxyaniline (Figure 2 and Figure S1). Minor metabolites were identified as an alternate loss of CH2 from the 3,4,5-trimethoxyaniline (3), the loss of 2 × CH2 (4), and the formation of a cyclic acetal (5). Formation of 5 likely occurs via the product of a mono-demethylation intercepting an intermediate oxocarbenium ion formed during a second oxidative demethylation. Finally, two minor species with quinoline ring oxidation (6 and 7) were observed at low abundance. The liver microsomal stability and the metabolite identification experiments identified cytochrome P450-mediated oxidation as the primary route of metabolism of 1, with the trimethoxyaniline as the major metabolic liability.
Figure 2
Metabolites of 1 identified (2–7) along with percentage abundance in respect to the parent compound. Compounds are ordered by HPLC retention time (see Supplementary Materials).
2.2. Synthesis of Analogs of 1
In an attempt to address the poor metabolic stability of 1, we prepared a focused array of analogs in which the trimethoxyaniline was replaced by various bioisosteres (8–43). Compounds 8–11, 13–31, and 33–43 were synthesized through nucleophilic aromatic displacement of the corresponding 4-chloroquinolines (Scheme 1) to furnish the products with good to excellent yields (54–89%) [8,9,11,12,13,14]. Additional analogs (23, 24, and 27) were synthesized by the same route with modest yields (12–38%) due to the reduced nucleophilicity of their respective anilines. An alternative route employing a Buchwald–Hartwig cross-coupling with 4-chloro-6-trifluoromethylquinoline enabled access to the 6-trifluoromethyl analogs (12) and (32) in 21% and 17% yields, respectively (Scheme 2) [8,9,12,13,14].
Scheme 1
General displacement synthetic procedure.
Scheme 2
General Buchwald-Hartwig cross-coupling procedure.
The metabolic stability of the 4-anilinoquinolines (8–13) were evaluated in MLMs (Table 1). Half-lives for metabolic clearance were normalized to propranolol in each experimental run to control for variations in the MLM preparations. Changing the substitution of the core quinoline heterocycle to 6,7-dimethoxy, as in 8, marginally increased stability relative to the 6-bromo substitution and maintained an aniline O-demethylation species as the primary metabolite. Modification of the aniline portion had a wide range effects on metabolic stability. The 2,3-methylenedioxy analog 9 showed similar stability to 8, with an aniline ring oxidation as the primary metabolite. Surprisingly, modification of the trimethoxy to the difluoro (10 and 11) yielded compounds with equal or less metabolic stability than 1. Metabolite identification experiments with the fluorinated analogs continued to show aniline oxidation as the primary site of metabolism: direct ring oxidation was observed with 10 and hydroxyl replacement of a fluorine with 11. However, replacement of the trimethoxyaniline with a perfluorinatedaniline resulted in a significantly more stable analog 12. Finally, the similar stability of 6-Br quinoline 1 and 6-CF3 quinoline 13 showed that switching of the 6-substituents had little effect on metabolism.
Table 1
Initial microsomal stability assessment (see Supplementary Materials for details).
Cmpd
R
X
T1/2 (min) a
T1/2 Ratio b
Primary Metabolism Site
Primary Metabolite Mass
8
3,4,5-(OMe)3
6,7-(OMe)2
12.2
1.6
Aniline
Parent − CH2
9
2,3-O-CH2-O-
6-Br
8.9
1.9
Aniline
Parent + O
10
2,5-F
6-Br
1.7
0.17
Aniline
Parent + O
11
2,4-F
6-Br
3.2
0.47
Aniline
Parent − F + OH
12
2,3,4,5,6-F5
6-CF3
45.6
7.0
nt
nt
13
3,4,5-(OMe)3
6-CF3
6.4
0.78
Aniline
Parent + O
nt: not tested; a mean average, n = 2; b T1/2 ratio is the normalization with respect to internal control (propranolol).
The analogs (8–13) were also assessed for GAK activity in a Time-resolved fluorescence energy transfer (TR-FRET) ligand-binding displacement assay and a live cell NanoBRET target engagement assay (Table 2) [8,9,12,13,14,15]. The 6,7-dimethoxyquinoline 8 displayed a high GAK potency biochemically and in cells (IC50 = 25 nM), as did the methylenedioxy compound 9 (IC50 = 22 nM). The two difluoroaniline isomers 10 and 11 had nanomolar Ki values in the biochemical assay but showed a larger drop-off in potency in the live cell NanoBRET assay, suggesting that they had poorer cell penetration. The metabolically stable pentafluoro aniline analog 12 had only weak GAK affinity in the biochemical assay (Ki = 880 nM), and did not progress to cellular evaluation. Finally, 13 showed a 4-fold loss in cellular GAK potency relative to 1. Considering both the MLM stability and cellular GAK activity of these initial analogs, we selected the methylenedioxy compound 9 for further rounds of optimization.
Table 2
GAK affinity measurements and microsomal clearance of initial quinoline series.
Cmpd
R
X
Biochemical aGAK Ki (nM)
Live Cell bGAK IC50 (nM) c
8
3,4,5-(OMe)3
6,7-(OMe)2
0.54 d
25 d
9
2,3-O-CH2-O-
6-Br
1.9
22
10
2,5-F
6-Br
15
870
11
2,4-F
6-Br
35
2700
12
2,3,4,5,6-F5
6-CF3
880
nt
13
3,4,5-(OMe)3
6-CF3
3.9 d
180 d
nt: not tested; a TR-FRET assay; b nanoBRET; c mean average, n = 2; d reference [12].
2.4. Design and Tesing of 4-Anilinoquinolines Analogs 14–43
Docking of 1 and 9 into a model of GAK derived from the co-crystal structure with gefitinib (Figure 3) was used to define key structural features for favorable binding using Schrödinger Maestro [16]. The compounds were minimized using LigPrep before docking into the ATP-binding site prepared from PDB structure 5Y7Z [17]. As expected, the most favorable docking pose placed the quinolinenitrogen in a position to accept a hydrogen bond from the backbone NH of the hinge region residue C126. The 5-membered ring of the methylenedioxy in 9 adopted a conformation where it could form a weak hydrophobic interaction with T123 [18]. Based on these results, we designed and prepared a series of isosteric replacements of the substituted aniline (14–43) that would maintain a favorable interaction with GAK and be less prone to oxidative metabolism. The analogs were assessed for GAK activity in the biochemical and live cell assays and for metabolic stability in MLM (Figure 4 and Table 3).
Figure 3
Docking of 1 (A) and 9 (B) into the ATP binding site of GAK (PDB: 5Y7Z). The protein is shown as a ribbon view with key amino acids as sticks with carbon atoms in grey, nitrogen in blue, oxygen in red, and sulphur in yellow. The GAK inhibitor is shown as sticks with carbon atoms in green, nitrogen in blue, and oxygen in red. The ATP-binding pocket is shown in light grey. Hydrogen bonds from the hinge residues to the inhibitor are highlighted as yellow dashed lines.
Figure 4
Aniline fragments incorporated into compounds 9 and 14–43 in Table 3.
Table 3
GAK affinity measurements and microsomal clearance of isosteric replacements of 9.
Cmpd a
X
Biochemical bGAK Ki (nM)
Live Cell cGAK IC50 (nM) d
T1/2 (Min) d
T1/2 Ratio e
9
Br
1.9
22
8.9
1.9
14
Br
1.4
120
1.8
0.27
15
Br
740
>5000
1.0
0.12
16
Br
460
>5000
2.0
0.24
17
Br
5.2
>5000
1.7
0.16
18
Br
29
2500
0.82
0.10
19
Br
1.9
180
1.7
0.16
20
Br
4.3
570
0.9
0.08
21
Br
37
1400
5.2
0.64
22
Br
19
>5000
1.4
0.17
23
Br
19
340
1.7
0.21
24
Br
11
220
5.3
0.68
25
Br
26
>5000
2.6
0.32
26
Br
23
3300
4.1
0.50
27
Br
4.1
>5000
8.1
1.0
28
Br
41
2500
1.2
0.15
29
Br
5.3
400
1.6
0.20
30
Br
2.8
>5000
0.73
0.09
31
Br
3.6
>5000
0.80
0.10
32
Br
14
760
2.5
0.32
33
CF3
720
>5000
5.9
0.90
34
CF3
34
2600
46.4
9.8
35
Br
20
1400
45.6
9.8
36
Br
nt
>5000
14.8
2.2
37
Br
47
>5000
7.5
0.92
38
Br
62
>5000
8.1
1.0
39
Br
360
>5000
>200
>25
40
Br
680
>5000
55.9
6.9
41
Br
110
>5000
27.8
3.4
42
Br
1600
>5000
33.6
4.2
43
Br
1700
>5000
15.6
1.9
nt: not tested; a R groups are shown in Figure 4; b TR-FRET assay; c nanoBRET; d mean average, n = 2; e T1/2 ratio is the normalization with respect to internal control (propranolol).
Results of the initial compound evaluation underscored the difficulty in the concurrent optimization of GAK potency and metabolic stability (Table 4). Consistent with the docking model supporting a favorable fit for bicyclic aniline fragments, a majority of the analogs (14–43) demonstrated potent GAK activity in a biochemical assay (Ki < 100 nM). Importantly, all GAK active analogs (14–43) maintained good NAK family selectivity in biochemical assays (Table S11). Unfortunately, potent biochemical GAK activity did not necessarily translate to potent GAK activity in cells. Most of the heterocycles showed lower GAK activity in the live cell NanoBRET assay compared to the methylenedioxy analog 9, although six compounds demonstrated submicromolar IC50 values in cells: benzofurans 19 and 20, benzoxazole 23, benzothiazole 24, benzothiadiazole 29, and indazole 32. Interestingly, we observed that a switch of the 2,3-substitution in the aniline heterocycle to the corresponding 3,4-substitution often led to greatly increased metabolic stability, with six compounds demonstrating T1/2 ratios ≥ 2-fold higher than 9: indazoles 34 and 35, benzimidazolone 39 and 42, and benzoxazolidinones 40 and 41. The benzimidazolone 39 was the most stable heterocycle in MLM with a T1/2 of >200 min and T1/2 ratio > 25. Unfortunately, the matched pair of indazoles 34 and 35 were the only compounds to maintain any measurable activity in the live cell GAK NanoBRET assay, with IC50 = 1.4 and 2.6 μM, respectively. However, this level of potency was 30- to 60-fold lower than SGC-GAK-1 (1) and did not meet our minimum criteria for use as an in vivo chemical probe for GAK.
Table 4
Pharmacokinetic profile of 1 and 8 with and without pre-administration of ABT.
Cmpd
T1/2 (min) a
T1/2 Ratio
Clint (mL/min/kg)
T1/2 (min) a
T1/2 Ratio b
Clint (mL/min/kg)
(−) ABT
(+) ABT
1
5.7
0.58
990
130
13.2
40
8
18
1.3
320
>200
14.7
<28
9
0.74
0.05
7610
1.8
0.13
3140
14
0.74
0.05
7550
2.5
0.18
2290
35
46
9.8
120
170
17.5
35
a mean average, n = 2; b T1/2 ratio following normalization with respect to an internal control (propranolol).
2.5. Metabolite ID Profiling and Kinome Scan of 35
Despite the limited utility of indazole 35 as an in vivo chemical probe for GAK we decided to perform a metabolite identification experiment to elaborate the origin of its improved stability. Following incubation in MLM, metabolites were identified that were nucleated around the aniline ring. Oxidation of the phenyl ring of the indazole yielded the major metabolites (44–45). Minor metabolites were identified as hydroxylation on an indazolenitrogen (46) (Figure 5). Notably, the quinoline core remained less metabolically labile than the indazole heterocycle.
Figure 5
Compound 35 and corresponding metabolites (44–46) identified along with percentage abundance ordered by retention time (see Supplementary Materials).
Compound 35 was also submitted to a KinomeScan assay to determine whether the indazole moiety had changed the selectivity profile across > 400 human kinases. The results (Figure 6) showed activity on only four kinases (GAK, MEK5, TGFBR2 & RIPK2) with a 1 µM concentration of 35. So, although activity on GAK was reduced compared to 1, the indazole 35 still maintained a narrow spectrum kinase inhibition profile in vitro.
Figure 6
Representation of the human kinome based on sequence identity-displaying KinomeScan results for 35 screened at 1 μM. Red dots represent hits differentiated by size; kinases inhibited at 1 μM GAK, 90%; MEK5, 89%; TGFBR2, 82%; RIPK2, 71% (see Supplementary Materials).
2.6. Electronic Fuikui Modelling Aniline of 4-Anilinoquinolines
In order to rationalize the observed structure activity relationships within the MLM results and gage propensity for reactive oxidation events to occur, we attempted to model the compound electronics using frontier molecular orbital theory and Fukui indices [19] (Figure 7). The negative values for f− correspond to regions that lose electron density when the molecule undergoes an electrophilic attack or when the molecule itself acts as a nucleophile [19]. Unfortunately, the results did not align with the structure activity relationships observed in the MLM.
Figure 7
Electronic mapping of the scaffolds showing Fuikui minus visualized with two isovalues, 0.002 (white) and 0.004 (green): (A) 9, (B) 17, (C) 19, (D) 20, (E) 14, (F) 32, (G) 35, (H) 37, (I) 39, (J) 40.
2.7. Modelling of Lead Compounds in XenoSite
We also analyzed the five lead compounds that retained GAK activity (1, 8, 9, 14, and 35) using the P450 modelling software XenoSite (https://swami.wustl.edu/xenosite/) [20]. The XenoSite models were trained on a large set of publicly available P450 metabolism data [20,21,22] to predict the sites of metabolism by CYP1A2, -2A6, -2B6, -2C19, -2C8, -2C9, -2D6, -2E1 and -3A4 (Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13) [20].
Figure 8
Compound 1 in XenoSite’s P450 Metabolism P450 1.0 algorithm.
Figure 9
Compound 8 in XenoSite’s P450 Metabolism P450 1.0 algorithm.
Figure 10
Compound 9 in XenoSite’s P450 Metabolism P450 1.0 algorithm.
Figure 11
Compound 14 in XenoSite’s P450 Metabolism P450 1.0 algorithm.
Figure 12
Compound 35 in XenoSite’s P450 Metabolism P450 1.0 algorithm.
Figure 13
Summary image highlighting the predicted metabolic profiles of 1, 8, 9, 14, and 35 in XenoSite’s P450 Metabolism P450 1.0 algorithm.
Several functional groups were predicted by XenoSite to be metabolically labile, including the trimethoxy functionality of 1 and 8 and the benzo[d][1,3]dioxol-4-yl of 9 and 14. The bromine substitution of 1 appeared to deactivate the core quinoline structure compared with 8, but this did not translate to a large increase in stability (Figure 8; Figure 9). This analysis highlights that the main metabolism is likely to occur on the aniline ring system. The methylenes of the benzo[d][1,3]dioxol-4-yl of 9 and 14 were also predicted to be a hot spot for metabolism (Figure 10; Figure 11). In contrast, the indazole of 35 showed less propensity for metabolism, a prediction that was consistent with the observed MLM experimental data (Figure 12).
2.8. 4-Anilinoquinolines Stablisation with ABT
Having modelled the P450 metabolism, we measured the stability of analogs 8, 9, 14, and 35 in MLM in comparison to 1. The experiments were run with and without the broad spectrum P450 inhibitor, 1-aminobenzotriazole (ABT) [23] (Table 4). In the absence of ABT, we observed increased stability of 8 and 35, with a 5- to 20-fold increase in the half-life in MLM compared to 1 (Table 5). However, the benzo[d][1,3]dioxol-4-yl analogs 9 and 14 showed higher intrinsic clearance and shorter half-lives. In the presence of ABT, 1, 8, and 35 showed greatly increased stability. Surprisingly, ABT did not, however, increase the stability of analogs 9 and 14. ABT is reported to inhibit 100% of CYP2A6 and -3A4 activity but residual activity (up to 20%) can be measured with CYP1A2, -2B6, -2C8, -2C9, -2C19, and -2D6 [24]. CYP2C9 was particularly resistant to inhibition by ABT, with 60% of the activity remaining. A separate study indicated that ABT caused a concentration-dependent inhibition of CYP1A2, -2B6, -2C9, -2C19, -2D6, and -3A4 [25]. Taken together, our results suggest that CYP2C9 or -2C8 may be responsible for the metabolism of the benzo[d][1,3]dioxol-4-yl group of 9 and 14 in MLM. XenoSite predicts a metabolism hotspot for these analogs with CYP2C9 or -2C8 that ABT may be unable to block in MLM.
Table 5
Pharmacokinetic profile of 1 and 8 with and without pre-administration of ABT.
Experiment
Tmax (h)
Cmax (ng/mL)
AUCinf (h·ng/mL)
10 mpk 1
0.5
97
98
2 h pretreat 50 mpk ABT 10 mpk 1
0.5
650
1900
10 mpk 8
1.0
19
71
2 h pretreat 50 mpk ABT 10 mpk 8
2.0
107
569
2.9. In Vivo Stabilisation of 1 and 8 with ABT
Having been unable to identify a potent (sub-µM) cell-active GAK inhibitor with improved metabolic stability in MLM, we decided to test an alternative strategy to enable dosing of a GAK chemical probe in vivo. Since in vitro inhibition of P450 enzymes by cotreatment with ABT in MLM had greatly increased the metabolic stability of 1 and 8, we explored the effect of co-dosing in vivo. A pharmacokinetic experiment was performed in which C57BL/6 mice were dosed orally with 10 mg/kg of 1 or 8 with or without 2 h pre-treatment with 50 mg/kg ABT (Table 5). As predicted by the microsomal stability experiment, orally dosed 1 or 8 were rapidly cleared in mice with a T1/2 < 1 h and a low Cmax (<100 ng/mL). In contrast, pretreatment with ABT resulted in a dramatically extended T1/2 with a >20-fold higher exposure for 1 and 8-fold higher exposure for 8. In the absence of ABT, the plasma concentrations fell below the cellular IC50 after approximately 1 h for 1 and after 2 h for 8. In contrast, pretreatment with ABT resulted in a sustained plasma concentration of 1 at 2-fold above the cellular IC50 for 6 h with a 10 mg/kg dose (Figure 14A,B). This shift was even more pronounced for 8 with 2-fold above the cellular IC50 maintained for >8 h with a 10 mg/kg dose (Figure 14C,D). These results demonstrate that the decreased metabolism observed upon co-dosing of ABT in MLM can translate to increased exposure of 4-anilinoquinolineGAK inhibitors in vivo.
Figure 14
Plasma concentration profiles of 1 (A,B) and 8 (C,D) in mice (n = 18) following a 10 mg/kg po dose with (green) or without (red) 2 h pretreatment with 50 mg/kg ABT.
3. Discussion
Chemical probes are potent and selective modulators of protein function that are useful as reagents to elucidate the biology of their targets [26,27]. The use of chemical probes in phenotypic models of disease can also provide powerful evidence to support target validation for drug discovery. Positive results in cellular experiments naturally lead to interest in using the chemical probe for target validation in animal models. However, many chemical probes are not amenable to use in vivo because they are unable to attain sufficient exposure to their target protein to support valid experiments. The pharmacokinetic profiles of probes can be addressed through careful chemical optimization of the series but at the cost of significant investments of time and resources.A common pharmacokinetic liability of compounds is that they are rapidly metabolized by cytochrome P450 enzymes in the liver and cleared from the animal, a process known as first pass metabolism. As a result, compounds are routinely evaluated for stability in MLM prior to progression to in vivo pharmacokinetic studies [28]. It has been well established that molecules with low MLM stability are prone to rapid first pass metabolism in vivo. The challenge of optimizing compounds to improve in vivo pharmacokinetic properties while maintaining their target potency and selectivity has vexed many medicinal chemistry programs, including our own exploration of the dark kinase GAK.The application of the machine learning algorithms such as XenoSite to develop predictive models of P450 metabolism has greatly improved the design of molecules with fewer metabolic liabilities. The XenoSite model accurately predicted the reduced metabolic liability of indazole 35, which showed improved stability in liver microsomes compared to 1. However, the modest cellular GAK activity of 35 (IC50 = 1.4 μM) meant that it was unsuitable for use as a chemical probe. Unfortunately, all of the other 4-aminoquinolines analogs with improved metabolic stability had even poorer GAK activity. We found that replacement of the 3,4,5-trimethoxy substitution with 5-membered heterocycles annulated at the aniline 2,3-positions yielded several compounds that maintained high GAK activity in biochemical assays, but without improvement in microsomal stability. On the other hand, compounds with 3,4-annulation of heterocycles in general proved to be significantly more metabolically stable, but at the cost of weakened GAK activity. The identification of 3,4,5-trimethoxyphenyl heterocyclic replacement groups with increased microsomal stability is a valuable learning tool. The 3,4,5-trimethoxyphenyl group is a commonly encountered substructure in ATP-competitive inhibitors of protein kinases and in a variety of other molecules with biological activity, such as the antibiotic trimethoprim, the antihypertensive reserpine, the anti-inflammatory colchicine, and the DHFR-inhibitor trimetrexate for the treatment of leiomyosarcoma [29]. Thus, other medicinal chemistry campaigns requiring isosteres of the 3,4,5-trimethoxyphenyl group may consider incorporation of select heterocyclic replacements identified here.Although our original goal of the chemical optimization of 1 to decrease intrinsic metabolism while retaining potent GAK activity remained unrealized, we were able to use co-dosing with ABT to reduce P450 metabolism and increase in vivo exposure. ABT is a potent inhibitor of CYP3A4, CYP2D6, and several other P450 enzymes, and has been successfully employed to enhance the pharmacokinetics of metabolically labile compounds [24,30]. Although not an approved drug for administration to patients, co-dosing with ABT to block oxidative metabolism in animal models is significantly less demanding of time and resources than chemical optimization, and may be broadly applicable for the development of in vivo chemical probes [31]. For example, co-dosing of ABT with the antifungal APX001 increased exposure of the compound and enhanced its efficacy in a murine model of invasive candidiasis [32]. Long-term dosing of ABT is also possible. ABT has displayed a well-behaved pharmacokinetic profile over five days when dosing with an osmotic pump [25]. Minimal side effects were observed for chronic ABT dosing, the most prominent of which was delayed gastric emptying [33]. ABT can also be utilized in species other than mice for in vivo studies. ABT is a competent cytochrome P450 inhibitor in rats, guinea pigs, and other species, and has been used to enhance pharmacokinetic exposure in these animals [30,34]. Our results further extend the utility of ABT by demonstrating an increase of in vivo exposure for certain 4-anilinoquinolineGAK inhibitors in mice.Not all 4-anilinoquinolines showed reduced MLM metabolism in the presence of ABT. The 4-anilinoquinolines 9 and 14 containing the benzo[d][1,3]dioxol-4-yl group had shorter half-lives than 1 in MLM and were still rapidly metabolized in the presence of ABT. Notably, the XenoSite model predicted that the benzo[d][1,3]dioxol-4-yl group was a hotspot for metabolism by CYP2C9 and CYP2C8, two P450s that are not strongly inhibited by ABT. For compounds such as these, it may be possible to use the co-dosing of other P450 inhibitors such as atipamezole to block CYP2C8/9 activity [35]. The development of a range of broad spectrum P450 inhibitors with overlapping profiles that can be used safely in animal models to block P450 metabolism would be a powerful addition to the medicinal chemist’s toolbox.
4. Materials and Methods
4.1. Computational Methods
Small molecules were prepared using the ligprep module of Schrödinger (Schrödinger Release 2018-4: LigPrep, Schrödinger, LLC, New York, NY, USA, 2018.) Structures were superimposed using flexible ligand alignment functionality (Schrödinger Release 2018-4: Maestro, Schrödinger, LLC, New York, NY, USA, 2018) using compound 8 as a template and employing the Bemis–Murcko method for common substructure recognition [14].
Compounds were incubated at a concentration of 10 µM with mouse liver microsomes in the presence of an excess of NADPH. The duplicate incubations were performed in 0.5 mL 96-well plates in a shaking water bath maintained at 37 °C for 30 min and quenched by the addition of an equal volume of methanol (see Supplementary Materials for full details and reports).
4.2.4. Pharmacokinetic Profiling
Pharmacokinetic profiling was performed by SAI Life Sciences (Hyderabad, India). All procedures of the present study were in accordance with the guidelines provided by the Committee for the Purpose of Control and Supervision of Experiments on Animals (CPCSEA) as published in the Gazette of India, 15 December 1998. Prior approval of the Institutional Animal Ethics Committee (IAEC) was obtained before initiation of the study. Healthy male C57BL/6 mice (8–12 weeks old) weighing between 20 and 35 g were procured from (Envigo Research Private Ltd., City, Hyderabad, India). Three mice were housed in each cage. Temperature and humidity were maintained at 22 ± 3 °C and 30–70%, respectively, and illumination was controlled to give a sequence of a 12 h light and 12 h dark cycle. Temperature and humidity were recorded by an automatically controlled data logger system. All animals were provided a laboratory rodent diet (Envigo Research Private Ltd., Hyderabad, India). Reverse osmosis water treated with ultraviolet light was provided ad libitum. Eighteen male mice were weighed and divided in to two groups (Group 1 and Group 2), with nine mice in each group. Animals in Group 1 were administered orally with a solution formulation of 1 or 8 at 10 mg/kg dose. Animals in Group 2 were administered orally with a solution formulation of ABT (2 h before dosing of 1 or 8) at 50 mg/kg dose. After 2 h of ABT administration, animals were administered orally with solution formulation of 1 or 8 at 10 mg/kg dose. Solution formulation of 1 or 8 was prepared in 5% NMP, 5% solutol HS-15 in normal saline. The dosing volume administered was 10 mL/kg for the oral route. Blood (~60 µL) samples were collected from a set of three mice at each time point in a labeled micro centrifuge tube containing K2EDTA solution as anticoagulant at pre-dose, 0.5, 1, 2, 4, 6, 8, 12, and 24 h. Plasma samples were separated by centrifugation at 4000 rpm for 10 min and stored below −70 °C until bioanalysis. Concentrations of 1 or 8 in mice plasma samples were determined by a LC-MS/MS method. The plasma concentration-time data for 1 or 8 was provided for data analysis by the bioanalytical group (Sai Life Sciences Ltd., Pune, India). The plasma concentration-time data were then used for the pharmacokinetic analysis. A non-compartmental analysis module in Phoenix WinNonlin (Version 6.3, Certara USA Inc., Raleigh, NC, USA) was used to assess the pharmacokinetic parameters. Maximum concentration (Cmax) and time to reach the maximum concentration (Tmax) were the observed values. The areas under the concentration time curve (AUClast and AUCinf) were calculated by linear trapezoidal rule.
4.2.5. KinomeScan Assay
KinomeScan Assay on 35 was performed as previously described at EuroFins DiscoverX (Fremont, CA, USA) [7].
4.3. Chemistry Methods and Compound Characterization
4.3.1. General Information
All reactions were performed using flame-dried round-bottomed flasks or reaction vessels. Where appropriate, reactions were carried out under an inert atmosphere of nitrogen with dry solvents, unless otherwise stated. Yields refer to chromatographically and spectroscopically pure isolated yields. Reagents were purchased at the highest commercial quality and used without further purification. Reactions were monitored by thin-layer chromatography carried out on 0.25-mm E. Merck silica gel plates (60F-254) using ultraviolet light as a visualizing agent. NMR spectra were recorded on a Varian Inova 400 (Varian, Palo Alto, CA, USA) or Inova 500 spectrometer (Varian, Palo Alto, CA, USA) and calibrated using residual undeuterated solvent as an internal reference (CDCl3: 1H NMR = 7.26, 13C NMR = 77.16). The following abbreviations or combinations thereof were used to explain the multiplicities observed: s = singlet, d = doublet, t = triplet, q = quartet, m = multiplet, br = broad. Liquid chromatography (LC) and high resolution mass spectra (HRMS) were recorded on a ThermoFisher hybrid LTQ FT (ICR 7T) (ThermoFisher, Waltham, MA USA). The University of Southampton (Southampton, UK) small molecule x-ray facility collected and analyzed all X-ray diffraction data.Compounds 1, 8, and 13 were synthesized as previously reported [8,9,12]. Compound spectra of 9–12, 14–43 are included in the Supplementary Materials. All compounds were >98% pure by 1H/13C NMR and LCMS.
4.3.2. General Procedure for the Synthesis of 4-anilinoquinolines (Scheme 1)
6-Bromo-4-chloroquinoline (1.0 equiv.) and 3,4,5-trimethoxyaniline (1.1 equiv.) were suspended in ethanol (10 mL) and refluxed for 18 h. The crude mixture was purified by flash chromatography using EtOAc:hexane followed by 1–5% methanol in EtOAc, solvent was removed under reduced pressure to afford the desired product (1, 8–11, 13–31, and 33–43).
4.3.3. General Procedure for the Synthesis of 4-anilinoquinolines (Scheme 2)
4-Chloro-6-(trifluoromethyl)quinoline (1.0 equiv.), aniline (1.1 equiv.), bis(dibenzylideneacetone)palladium(0) (Pd(dba)2) (0.1 equiv.), 2-dicyclohexylphosphino-2′,4′,6′-triisopropylbiphenyl (XPhos) (0.15 equiv.), and cesium carbonate (2.0 equiv.) were all suspended in dimethylformamide 15 mL and degassed for 5 min. The mixture was heated to 140 °C for 18 h. The crude mixture was passed through a plug of Celite 545 before been purified by flash chromatography 20–100% ethyl acetate/hexane followed by 1–5% methanol/ethyl acetate and solvent removed under reduced pressure to afford the desired product (12 and 32).
Samples were analyzed with a ThermoFisher Q Exactive HF-X (ThermoFisher, Bremen, Germany) mass spectrometer coupled with a Waters Acquity H-class liquid chromatograph system (Waters Corporation, Milford, MA, USA). Samples were introduced via a heated electrospray source (HESI) at a flow rate of 0.6 mL/min. Electrospray source conditions were set as: spray voltage 3.0 kV, sheath gas (nitrogen) 60 arb, auxiliary gas (nitrogen) 20 arb, sweep gas (nitrogen) 0 arb, nebulizer temperature 375 °C, capillary temperature 380 °C, RF funnel 45 V. The mass range was set to 150–2000 m/z. All measurements were recorded at a resolution setting of 120,000.Separations were conducted on a Waters Acquity UPLC BEH C18 column (2.1 × 50 mM, 1.7 μM particle size). LC conditions were set at 100% water with 0.1% formic acid (A) ramped linearly over 9.8 min to 95% acetonitrile with 0.1% formic acid (B) and held until 10.2 min. At 10.21 min the gradient was switched back to 100% A and allowed to re-equilibrate until 11.25 min. Injection volume for all samples was 3 μL.Xcalibur (ThermoFisher) was used to analyze the data. Solutions were analyzed at 0.1 mg/mL or less based on responsiveness to the ESI mechanism. Molecular formula assignments were determined with Molecular Formula Calculator. All observed species were singly charged, as verified by unit m/z separation between mass spectral peaks corresponding to the 12C and 13C12Cc-1 isotope for each elemental composition.cyclin G-associated kinase (GAK)mouse liver microsomes (MLM)nicotinamide adenine dinucleotide phosphate (NADPH)liquid chromatography mass spectra (LCMS)mitogen-activated protein kinase kinase 5 (MEK5/MAP2K5)transforming growth factor beta receptor 2 (TGFBR2)receptor-interacting serine/threonine-protein kinase 2 (RIPK2)1-aminobenzotriazole(ABT)milligrams per kilogram (MPK)N-methyl-2-pyrrolidoneolidone (NMP)Bis(dibenzylideneacetone)palladium(0) (Pd(dba)2)2-Dicyclohexylphosphino-2′,4′,6′-triisopropylbiphenyl (XPhos)
5. Conclusions
We found that optimization of a series of 4-anilinoquinolinesGAK inhibitors to decrease liver metabolism was not compatible with maintaining potent enzyme activity. Nonetheless, several heterocyclic replacements of the trimethoxyaniline moiety were identified that had greatly improved metabolic stability and may be useful as bioisosteres in other lead optimization projects. Co-dosing with the P450 inhibitor ABT proved to be a more general approach to blocking the metabolism of several 4-anilinoquinolinesGAK inhibitors in MLM and in vivo. The availability of additional non-toxic broad spectrum P450 inhibitors could further broaden the implementation of this approach to other chemical probes.
Authors: Miles A Fabian; William H Biggs; Daniel K Treiber; Corey E Atteridge; Mihai D Azimioara; Michael G Benedetti; Todd A Carter; Pietro Ciceri; Philip T Edeen; Mark Floyd; Julia M Ford; Margaret Galvin; Jay L Gerlach; Robert M Grotzfeld; Sanna Herrgard; Darren E Insko; Michael A Insko; Andiliy G Lai; Jean-Michel Lélias; Shamal A Mehta; Zdravko V Milanov; Anne Marie Velasco; Lisa M Wodicka; Hitesh K Patel; Patrick P Zarrinkar; David J Lockhart Journal: Nat Biotechnol Date: 2005-02-13 Impact factor: 54.908
Authors: Christopher R M Asquith; Neil Fleck; Chad D Torrice; Daniel J Crona; Christoph Grundner; William J Zuercher Journal: Bioorg Med Chem Lett Date: 2019-07-10 Impact factor: 2.823
Authors: Christopher R M Asquith; Tuomo Laitinen; James M Bennett; Paulo H Godoi; Michael P East; Graham J Tizzard; Lee M Graves; Gary L Johnson; Ronna E Dornsife; Carrow I Wells; Jonathan M Elkins; Timothy M Willson; William J Zuercher Journal: ChemMedChem Date: 2017-11-27 Impact factor: 3.466
Authors: Suresh K Balani; Ping Li; Joanne Nguyen; Kym Cardoza; Hang Zeng; Dun-Xue Mu; Jing-Tao Wu; Liang-Shang Gan; Frank W Lee Journal: Drug Metab Dispos Date: 2004-06-24 Impact factor: 3.922
Authors: Christopher R M Asquith; Louisa Temme; Michael P East; Tuomo Laitinen; Julie Pickett; Frank E Kwarcinski; Parvathi Sinha; Carrow I Wells; Gary L Johnson; Reena Zutshi; David H Drewry Journal: ChemMedChem Date: 2022-05-09 Impact factor: 3.540
Authors: Christopher R M Asquith; Tuomo Laitinen; Carrow I Wells; Graham J Tizzard; William J Zuercher Journal: Molecules Date: 2020-04-07 Impact factor: 4.411
Authors: Andrew A Bieberich; Tuomo Laitinen; Kaitlyn Maffuid; Raymond O Fatig; Chad D Torrice; David C Morris; Daniel J Crona; Christopher R M Asquith Journal: Sci Rep Date: 2022-07-27 Impact factor: 4.996