New polycyclic heterocycles were synthesised and evaluated as potential inhibitors of thymidine phosphorylase (TP). Inspired by the pharmacophoric pyrimidinedione core of the natural substrate, four series have been designed in order to interact with large empty pockets of the active site: pyrimidoquinoline-2,4-diones (series A), pyrimidinedione linked to a pyrroloquinoline-1,3-diones (series B and C), the polycyclic heterocycle has been replaced by a pyrimidopyridopyrrolidinetetraone (series D). In each series, the tricyclic nitrogen heterocyclic moiety has been synthesised by a one-pot multicomponent reaction. Compared to 7-DX used as control, 2d, 2l, 2p (series A), 28a (series D), and the open intermediate 30 showed modest to good activities. A kinetic study confirmed that the most active compounds 2d, 2p are competitive inhibitors. Molecular docking analysis confirmed the interaction of these new compounds at the active binding site of TP and highlighted a plausible specific interaction in a pocket that had not yet been explored.
New polycyclic heterocycles were synthesised and evaluated as potential inhibitors of thymidine phosphorylase (TP). Inspired by the pharmacophoric pyrimidinedione core of the natural substrate, four series have been designed in order to interact with large empty pockets of the active site: pyrimidoquinoline-2,4-diones (series A), pyrimidinedione linked to a pyrroloquinoline-1,3-diones (series B and C), the polycyclic heterocycle has been replaced by a pyrimidopyridopyrrolidinetetraone (series D). In each series, the tricyclic nitrogen heterocyclic moiety has been synthesised by a one-pot multicomponent reaction. Compared to 7-DX used as control, 2d, 2l, 2p (series A), 28a (series D), and the open intermediate 30 showed modest to good activities. A kinetic study confirmed that the most active compounds 2d, 2p are competitive inhibitors. Molecular docking analysis confirmed the interaction of these new compounds at the active binding site of TP and highlighted a plausible specific interaction in a pocket that had not yet been explored.
Angiogenesis – the formation of new blood vessels from pre-existing vasculature – has been validated as a target for several tumours and has been shown to promote tumour growth and metastasis. Given the complexity of this process, there is a need for new compounds targeting multiple pro-angiogenic factors. In this context, thymidine phosphorylase (TP) seems to be an interesting and understudied therapeutic target.The physiological role of TP is to catalyse the reversible phosphorolysis of thymidine into thymine and 2-deoxyribose-1-phosphate (2dR1P) that is metabolised into 2-deoxyribose (2dR) (Figure 1). TP has been shown to be up-regulated in the hypoxic regions of many solid tumours (stomach, pancreas, kidney, oesophagus, breast, ovary, lung, colon, bladder, uterus, kidney …),,. Usually, this over-expression is highly associated with tumour micro vessel level, infiltration, metastasis and also correlates with the aggressiveness and invasiveness of the cancer. Via 2dR, TP stimulates the secretion of several pro-angiogenic factors such as VEGF, MMP-1, IL-8 by both malignant cells and stromal cells located in the tumour microenvironment,. Consequently, TP inhibitors seem to be promising agents.
Figure 1.
TP and 2dR1P roles in tumour development.
TP and 2dR1P roles in tumour development.Various crystal structures of Escherichia coli TP and human TP (hTP) complexes with thymine or thymidine analogues have been published. Briefly, the binding of thymidine is ensured by three hydrogen bonds between the two carbonyls and the NH function with an arginine, a serine and a lysine. These studies revealed a large empty space (pocket 1) facing the C-5 and C-6 of the pyrimidinedione and another (pocket 2) close to the 1-position (Figure 2).
Figure 2.
Schematic representation of the pyridimidinedione binding site of hTP and chemical structures of 7-DX and TPI.
Schematic representation of the pyridimidinedione binding site of hTP and chemical structures of 7-DX and TPI.Various TP inhibitors have been reported in the literature,,. For a recent review see Sajid et al.. In order to interact with the thymidine binding site, most of them are pyrimidine-2,4-dione derivatives. 7-Deazaxanthine (7-DX) and 5-chloro-6-[1-(2-iminopyrrolidinyl)methyl]uracil hydrochloride (TPI) have emerged (Figure 2),. Up to now, the most potent inhibitor is TPI (IC50 hTP = 0.035 µM), whose anti-angiogenic activity was attested through several in vitro and in vivo models. TPI has been co-crystallized with hTP. In combination with trifluridine, TPI has been recently approved (TAS-102, Lonsurf®) for the treatment of metastatic colorectal cancer. The role of TPI is dual, it prevents the degradation of trifluridine by TP and exerts anti-angiogenic effect.Only a few polycyclic TP inhibitors have been reported (for a review see Bera et al.). Crystallographic studies data allowed us to design inhibitors in which the natural ligand feature is linked to a polycyclic aromatic nitrogen heterocycle either by ring annelation or via various linkers.In the present study, we describe the synthesis and the biological in vitro evaluation as TP inhibitors of four series of new original aromatic derivatives. To elucidate the mechanism of enzyme inhibition of the hits, a brief kinetic study was attempted. In silico molecular docking studies have also been performed to explore the binding site and possible interactions mode of these new derivatives with TP.
Results and discussion
Chemistry
In a first time, we have designed three series of new heterocycles as potential TP inhibitors. The general structures are presented in Figure 3. In series A, the pyrimidinedione was annelated to a quinoline giving a rigid tricyclic heterocycle. The substituents on the 5 to 9-positions could modulate the interaction with pocket 1. Compounds of series B and C were more flexible. Thus, a pyrimidinedione substituted by a methylenic chain either on the 1-position (series B) or on the 6-position (series C) was linked to the 2-position of a quinolopyrrolidinedione. In a second time, the results of a first docking study allowed us to design compounds of series D in which the quinolopyrrolidinedione has been replaced by a pyrimidopyridopyrrolidinetetraone.
Figure 3.
General structures of the designed derivatives in series A–D.
General structures of the designed derivatives in series A–D.Our strategies for synthesising the tricyclic heterocycles are based on one-pot multicomponent reactions, developed in our group, that involved an aniline, an aldehyde, and a 1,3-dicarbonyl derivative.
Pyrimido[4,5-b]quinoline-2,4(1H, 3H)-diones (series A)
We have described a very simple and convenient one-pot reaction for the synthesis of pyrimido[4,5-b]quinoline-2,4(1H,3H)-diones involving an aldehyde, an aniline, and barbituric acid 1 as the 1,3-dicarbonyl reagent (Table 1). This straightforward method circumvents the preparation of unstable substituted 2-aminobenzaldehydes that limits the scope of previously described syntheses. The use of commercially available anilines allowed the facile syntheses of pyrimido[4,5-b]quinolinediones substituted on the 6 to 9-positions with electron donor or electron-withdrawing groups. Access to the 5-substituted derivatives is also possible starting from aliphatic or aromatic aldehydes. We have previously described the synthesis of compounds 2a–d, 2f–h, 2j–k, 2m, 2p–r which have been obtained in 25% to 88% yield. According to this method, six new derivatives (2e, 2i, 2l, 2n, 2o, 2s) have been synthesised in modest-to-good yields (Table 1).
Table 1.
Synthesis of the pyrimido[4,5-b]quinoline-2,4(1H,3H)-diones 2.
Compound
R5
R6
R7
R8
R9
Method
Reaction time (h)
Yield (%)
2aa
H
H
CH3
H
H
D
72
45b
2ba
H
H
H
CH3
H
C
2
56b
2ca
H
H
H
H
CH3
A
1e
25d
2da
H
H
H
OCH3
H
C
2
70b
2e
H
H
H
OPh
H
C
6
57d
2fa
H
H
H
OCH2Ph
H
C
8
66b
2ga
H
H
O–CH2–O
H
A
19
75b
2ha
H
OCH3
OCH3
OCH3
H
A
2
73b
2i
H
H
H
CH2Ph
H
C
3
41d
2ja
H
H
H
Cl
H
A
100
50b
2ka
H
H
H
CF3
H
A
1e
51b
2l
H
H
Cl
Cl
H
A
5
17b
2ma
H
H
H
A
48
65b
2n
H
H
H
A
20
32b
2o
H
H
H
B
3
67d
2pa
CH3
H
H
OCH3
H
A
24
77b
2qa
CH3
H
O-CH2-O
H
A
4.5
44b
2ra
Arc
H
O-CH2-O
H
A
67
88b
aPreviously described compound. bYield of pure crude product. cAr = 3,4,5-trimethoxyphenyl. dYield of purified product. eIncomplete reaction.
Synthesis of the pyrimido[4,5-b]quinoline-2,4(1H,3H)-diones 2.aPreviously described compound. bYield of pure crude product. cAr = 3,4,5-trimethoxyphenyl. dYield of purified product. eIncomplete reaction.It is worth noting that the presence of two electron-withdrawing groups on anilines is allowed but provided compound 2l in low yield (17%). Tetra and pentacyclic derivatives have been obtained respectively using 1-naphthylamine and 1-aminoanthracene (compounds 2n and 2o). In the last case, it was necessary to dilute the reaction mixture (330 mL of AcOH/mmol instead of 100 mL) to avoid the precipitation of the poorly soluble dihydro intermediate. The angular structures of compounds 2n, and 2o have been unambiguously determined by NMR.The preparation of compound 2s has required another reaction scheme. In the 1,3-dimethyl series, the literature described the synthesis of 5-dimethylaminopyrimidoquinolines using 6-arylamino-1,3-dimethyluracils as starting materials. When the synthesis was performed by nucleophilic substitution on an intermediate 5-thiomethyl heterocycle, the overall yields were very low due to the preparation of the thiomethylheterocycle. The use of phosgene iminium chloride (Viehe’s salt) furnished pyrimidoquinolines substituted either on the 5-position by a 5-dimethylamino group (in the presence of triethylamine) or by a chloride (without triethylamine). The dimethylpyrimidoquinoline 2s has been prepared using this second procedure, so 6-arylaminouracil 5 was obtained by reaction of 6-chlorouracil 3 with 3,4-methylenedioxyaniline 4, under microwave irradiation, according to the method described by Fang et al.. Then, the reaction of the Viehe’s salt on 5 gave an iminium intermediate that cyclized into 2s (Scheme 1). Surprisingly, the reaction leads to the dimethylamino compound even in the absence of triethylamine. The low solubility of 2s, which precipitates in the reaction mixture, probably prevents the nucleophilic attack by chloride ions observed in other series.
Scheme 1.
Synthesis of compound 2s.
Synthesis of compound 2s.
2H-pyrroloquinoline-1,3-diones/pyrimidinediones (series B and C)
Compounds of series B et C have been obtained by N2-alkylation of the 6,7-methylenedioxy-2H-pyrrolo[3,4-b]quinolines (10) using a pyrimidinedione substituted by a halogenoalkyl chain.Although many publications have already reported the synthesis of N2-substituted 2H-pyrrolo[3,4-b]quinolines, only three have been covered regarding the N2-unsubstituted heterocycle. The reported methods are multistep syntheses. Starting from the corresponding 2-aminoacetophenones, the 9-methyl derivatives were obtained as a guanidinium salt in two steps and 38% overall yield and the 9-phenyl analogue in 5 steps and 12% overall yield. The ultimate step of the third chemical scheme involved the condensation of urea on the quinoline-2,3-dicarboxylic diacid to give the unsubstituted heterocycle.In a first time, we intended to apply this last method to the synthesis of the 6,7-methylenedioxy derivative 10 using a 6,7-methylenedioxyquinoline-2,3-dicarboxylate. In the literature, an analogous dimethyl ester has been already described either by reaction of the 2-aminobenzaldehyde with dimethyl acetylenedicarboxylate and subsequent acidic treatment or by reaction of the methylenedioxyaniline with the dimethyl acetylenedicarboxylate and subsequent treatment with the Vilsmeier reagent. In both cases, the yield of the second step was not specified.We thought that it would be possible to prepare the diethyl ester 8 in a one-pot reaction. Indeed, 8 was obtained in 40% yield by refluxing in EtOH a mixture of aniline 4, diethyl acetylenedicarboxylate (7) and formaldehyde (6) (Scheme 2). Unfortunately, the condensation of 8 with urea afforded 10 in only 12% yield. Assays with the corresponding diacid were unsuccessful.
Scheme 2.
Synthesis of compound 10.
Synthesis of compound 10.Finally, we have synthesised 10 in 39% yield according to an adjustment of our three-component one-pot method replacing the cyclic 1,3-diketone by bromomaleimide (9) which is more easily accessible than the hydroxymaleimide (Scheme 2). With this example, we have shown that the scope of the one-pot reaction could be extended to 1,3-diketone analogues such as a cyclic β-halogenoketone.It is well known that pyrimidinediones could be alkylated on the N1 and/or on the N3. As an example, the reaction of thymine with the ω-dibromobutane and the ω-dibromohexane lead to a mixture of N1-substituted derivatives and N1,N1′-dimers. With the dibromopentyl chain in the same conditions, only the N1,N3-disubstituted pyrimidinedione has been isolated. In several papers, the N1-alkylation was achieved in two steps via the formation of the 2,4-O-disilyl pyrimidine in order to avoid the formation of side products. We have prepared compounds 12a,b and 15 according to the described procedure but without isolation of the unstable disilyl intermediates (Scheme 3).
Scheme 3.
Synthesis of compounds 13a, 13b, 16a, 16b, and 17.
Synthesis of compounds 13a, 13b, 16a, 16b, and 17.Final compounds 13a and 13b were obtained by alkylation of 10 by 12a and 12b in DMF in the presence of K2CO3. Unfortunately, in the bromo-uracil series, it was not possible to isolate compound 16b which reacted again with another molecule of 15 to give 17 (Scheme 3).Access to the 5′-bromo derivatives required the use of a benzoyl protecting group on the N3 position of 5-bromouracil (14). After alkylation of the N3-benzoylbromouracil (18), the reaction of the resulting compounds 19a,b with 10 and subsequent deprotection of the pyrimidinedione in acidic medium afforded compounds 16a and 16b (Scheme 3). A NOE NMR experiment realised on compound 19a confirmed the N1-alkylation of the pyrimidinedione from bromouracil (14).For the synthesis of series C, it was necessary to protect the 6-chloromethyluracil (21) by bis(trimethylsilyl)acetamide (BSA) before the condensation with 10. The reaction of the resulting 5-unsubstituted compound 22 with N-halogenosuccinimides lead to the 5′-halogeno compounds 23a–c (Scheme 4).
Scheme 4.
Synthesis of compounds 22 and 23a–c.
Synthesis of compounds 22 and 23a–c.Biological results and molecular modelling have prompted us to synthesise the N2-benzyl substituted compound 24 (Table 3) in which the pyrimidinedione moiety of derivative 22 is replaced by a phenyl nucleus. Thus, the reaction of benzyl bromide with 10 in the presence of K2CO3 in DMF furnished 24 in 65% yield.
Table 3.
Escherichia coli TP inhibition of series B and C derivatives.
Derivatives of series D have been synthesised from the unsubstituted tricyclic heterocycle 27 (Scheme 5). We have previously described easy access to this heterocyclic skeleton via a new synthon, the 4-formyl-3-hydroxy-2,5-dioxo-2,5-dihydro-1H-pyrrole (26). The functionalization of 27 was achieved by reaction with two equivalents of the appropriate amine in refluxed DMF, and cyclisation of the resulting diamide using two equivalents of p-toluenesulfonic acid (PTSA). Four new derivatives 28c, 28e, 28f, and 28h have been synthesised for the present study.
Scheme 5.
Synthesis of compounds 28a–h and 30.
Synthesis of compounds 28a–h and 30.In an attempt to prepare 28i, the reaction of the 6-aminomethyluracil hydrochloride (29), with 27 afforded only the diamide 30 which could not be cyclized into 28i whatever the acidic conditions used (PTSA, TFA, HCl, H3PO4, B(OH)3) (Scheme 5).
In vitro thymidine phosphorylase inhibiting activity and molecular docking study
Inhibitory activities of the new heterocyclic derivatives in series A-D have been tested in vitro against recombinant E. coli TP. The adopted protocol is a modification of the original method developed by Khan. TP activity and inhibition assays were based on the phosphorolysis of thymidine to thymine, this conversion results in a significant spectrophotometric decrease. At 290 nm and pH 7.4, thymidine has a higher extinction coefficient than thymine (Δε = −480 M−1 cm−1). The absorbance at 290 nm was measured every 2 min for 30 min in 96-well plates at 25 °C.A preliminary screening has been carried out at 50 µM or at the maximum concentration allowed by the solubility of the tested compound in the incubation medium (aqueous buffer containing 2.5% of DMSO). Results are given in percentage of inhibition (average of three measurements). For the most promising compounds (inhibition > 30% and solubility > 50 µM), IC50 values have been determined (Tables 2–4). 7-DX, which presents the activity of the same order of magnitude (IC50
TP = 40 µM) was used as a reference.Escherichia coli TP inhibition of series A derivatives.aMaximum solubility.Escherichia coli TP inhibition of series B and C derivatives.aMaximum solubility.Escherichia coli TP inhibition of series D derivatives.aMaximum solubility.To help understand structure–activity relationships, the geometry optimised structures of the synthesised compounds were docked into the active site of hTP using GOLD 5.1. In order to validate the docking protocol, an X-ray crystal of TP complexed with TPI was retrieved, and the ligand was redocked into the active site of the enzyme. The predicted binding mode and interaction pattern of TPI generated in our study were found to be in accordance with the reported crystallographic study with an RMS deviation of 0.654 Å.Finally, a brief kinetic study was performed to elucidate the mechanism of inhibition of the two most interesting derivatives.
Series A
Series A has been designed as a conformationally constrained series in which the pyrimidinedione moiety should interact with the thymine binding site and the other nuclei with pocket 1 facing the C-5 and C-6 of the thymine (Figure 2). In order to modulate the interactions and/or the solubility, various substituents have been introduced on the 5 to 9-positions.Seven derivatives exhibited inhibiting activity towards E. coli TP (compounds 2c, 2d, 2e, 2i, 2l, 2m, 2p), three of them (compounds 2d, 2l, 2p) are interesting in terms of activity and solubility as compared to 7-DX (Table 2). This preliminary enzyme inhibition study suggested the following structure-activity relationships for series A compounds:
Table 2.
Escherichia coli TP inhibition of series A derivatives.
Compound
Structure
TP inhibition (%) at 50 µM or at maximum solubilityaIC50 (µM)
2a
0% at 50 µM
2b
0% at 50 µM
2c
14% at 50 µM
2d
89% at 50 µMIC50 = 28 ± 1 µM
2e
30% at 50 µMa
2f
0% at 50 µMa
2g
0% at 15 µMa
2h
0% at 15 µMa
2i
30% at 50 µMa
2j
0% at 50 µMa
2k
0% at 50 µM
2l
71% at 50 µMIC50 = 42 ± 7 µM
2m
29% at 50 µM
2n
0% at 15 µMa
2o
0% at 25 µMa
2p
67% at 50 µMIC50 = 26 ± 3 µM
2q
0% at 50 µM
2r
0% at 25 µMa
2s
0% at 15 µMa
7-DX
60% at 50 µMIC50 = 28 ± 6 µM
aMaximum solubility.
A methyl on the benzenic ring was unfavourable except for the 9-methyl derivative 2c that was weakly active.Monosubstitution on the 8-position by an ether group seemed to be interesting in the case of a methoxy group (compound 2d) but the activity decreased with the substituent size (compounds 2e and 2f). The replacement of the oxygen by methylene seemed to be possible (compare compounds 2e–2i). Compounds 2g and 2h were sparingly soluble; consequently, it is not possible to conclude about the effect of the presence of two ether groups.The presence of an electron-withdrawing substituent on the 8-position such as a chlorine atom (compound 2j) or a trifluoromethyl group (compound 2k) was unfavourable. However, the 7,8-dichloro derivatives 2l exhibited an IC50 in the same order of magnitude as the 8-methoxy derivative 2d.The enzyme seemed to accept the presence of an additional angular benzene ring in compound 2m. It is not possible to conclude concerning the geometrical isomer 2n or the homologue 2o which were less soluble.The introduction of methyl on the 5-position of the 8-methoxy derivative 2d was possible (compound 2p). In the case of compound 2g, the 5-substituted derivatives remained inactive whatever the substituent sizes (compounds 2q, 2r, 2s).In this series, results were found to be optimum for the two 8-methoxypyrimido[4,5-b]quinoline-2,4(1H,3H)-diones 2d and 2p, and the 7,8-dichloro derivative 2l.As expected, preliminary docking study showed that the pyrimidinedione moiety of the series A derivatives interacted with the active site of hTP in a similar binding mode as the 5-chlorouracil fragment of TPI, that is, by establishing hydrogen bonds with arginine 202, serine 217 and lysine 221 of the α domain and histidine 116 of the α/β domain. As an example, docking of compound 2d superposed on TPI (in green) is presented in Figure 4. Results were the same for the other active derivatives 2p and 2l (Supporting information, Figure S1).
Figure 4.
Docking superposition in the hTP active site of compounds 2d (in blue), 2e (in pink), 2i (in orange), and TPI (in green).
Docking superposition in the hTP active site of compounds 2d (in blue), 2e (in pink), 2i (in orange), and TPI (in green).The less marked activity of compounds 2c and 2e was also supported by the modelling results that highlighted two major binding modes corresponding in both 180° vertical and 180° horizontal rotations of the tricyclic moieties (Supporting information, Figure S2 for compound 2c as an example). It is worth noting that several substituents on the 8-position seemed to be located in the same area of the active site as showed in Figure 4: methoxy group of compound 2d, benzyl of compound 2i, and phenoxy for one of the conformations of compound 2e.
Series B and C
Series B and C have been designed as more flexible series than series A with the aim of interactions of the pyrimidinedione in the thymidine binding site of TP and of the methylenedioxypyrroloquinoline moiety in pocket 1 (see Figure 2).In terms of enzymatic inhibition, series B and C did not present really interesting activities in comparison with 7-DX. The unsubstituted derivative 22 and the 5-bromopyrimidinediones 16a, 16b, or 23b seemed to the more interesting but, unfortunately, they are the less soluble. The spacer size has no influence. In the case of a methylene linker (series C), the presence of a halogen on the 5-position was not favourable. More interestingly, the pyrimidinedione did not seem to be essential, it can be replaced by a phenyl nucleus (compound 24) (Table 3).Surprisingly, docking study showed that in both series, the dioxolane of the pyrroloquinoline moiety is positioned in the same place as the pyrimidinedione of TPI suggesting interactions of the methylenedioxy group with the active site; the pyrimidinedione moiety interacting with a large pocket which has never been explored (superposition of compound 16b and TPI on Figure 5).
Figure 5.
Docking superposition in the hTP active site of TPI (in green), compounds 16b (dark green cloud), and 23b (white cloud).
Docking superposition in the hTP active site of TPI (in green), compounds 16b (dark green cloud), and 23b (white cloud).However, as presented in Figure 5 for compound 16b (dark green cloud) and compound 23b (white cloud), two differences between series B and C were noticed. On the one hand, the tricyclic moieties were not superposed (180° vertical rotation). On the other hand, the aliphatic chains and pyrimidinediones of the B-series molecules were positioned in the upper part of the pockets, whereas pyrimidinediones of series C derivatives interacted with the lower part of this pocket. Interestingly, the phenyl nucleus of compound 2i (series A) was also located in this white cloud (Supporting information, Figure S3).
Series D
The B and C series modelling results led us to design the D-series in which the methylenedioxyphenyl rings were replaced by a pyrimidinedione with the aim to optimise interactions with the thymidine binding site of the natural ligand. Furthermore, various substituents have been introduced on the imide nitrogen in an attempt to explore the lower part of the “new” pocket.In a physicochemical point of view, series D compounds are more soluble than the previous series. However, they were uninteresting in terms of enzymatic inhibition activity except for the N-benzyl derivative 28a (Table 4).
Table 4.
Escherichia coli TP inhibition of series D derivatives.
Compound
Structure
TP inhibition (%) at 50 µM or maximum solubilitya
28a
44% at 50 µMIC50 = 62 ± 7 µM
28b
0% at 50 µM
28c
0% at 50 µM
28d
0% at 30µMa
28e
0% at 50 µM
28f
0% at 50 µM
28g
0% at 30 µMa
28h
0% at 50 µM
27
0% at 50 µM
30
31% at 50 µMIC50 = 87 ± 4 µM
7-DX
60% at 50 µMIC50 = 28 ± 6 µM
aMaximum solubility.
As expected, the tricyclic moiety of compound 28a was well superposed with compound 2i, the benzyl analog in series A. The benzyl nuclei were positioned in the same area but not in the same position (Figure 6).
Figure 6.
Docking superposition in the hTP active site of compounds 2i (in orange), 28a (in purple), and 30 (in blue).
Docking superposition in the hTP active site of compounds 2i (in orange), 28a (in purple), and 30 (in blue).Surprisingly, the flexible open intermediate derivative 30 was as active as compound 28a. Interestingly, one 6-pyrimidinomethylamino arm was positioned in the same area as previously described compounds (Figure 6), the other one was located at the entry of the pocket.
Mechanism of enzyme inhibition
The synthesised compounds have been designed as TP inhibitors interacting with the thymidine fixation site. However, TP being a two substrates enzyme, it was necessary to verify that new inhibitors do not bind to the phosphate site. This study has been realised with one of the most active and/or the most soluble compounds in each series, that is, compounds 2p, 23b, 28a, and 30. E. coli TP activity was determined in the presence of a saturating concentration of thymidine (equivalent to three times the value of the Michaelis-Menten constant [km]) and variable amounts of phosphate (2–30 mM). Tested compounds have been used at a concentration corresponding to the IC50. For each tested derivative, phosphate concentration did not have any significant influence on the percentage of inhibition, indicating that the selected inhibitors did not bind to the phosphate site.To explore whether 2d and 2p, the most active compounds which showed the same activities as 7-DX (TP CI50: 26–28 µM), acted as a competitive, uncompetitive, or mixed TP inhibitor, a brief kinetic study was performed according to previously described method. The percentage of TP inhibition was determined for different concentrations of the inhibitors and thymidine. Lineweaver-Burk plots showed that all straight lines converged at the same point on the positive side of the y-axis (Figure 7). Consequently, 2d and 2p exhibited competitive inhibition kinetic on TP with thymidine as substrate. This was further supported by the fact that the Km values increased in the presence of 2d and 2p, while the Vmax values did not change significantly as the inhibitor concentration increased.
Figure 7.
Lineweaver–Burk plots of E. coli TP inhibition by 2d and 2p, in the presence of variable concentrations of thymidine demonstrating competitive-type enzyme inhibition. Results are presented as means ± SD; SD denoted by error bars (experiments carried out in triplicate).
Lineweaver–Burk plots of E. coli TP inhibition by 2d and 2p, in the presence of variable concentrations of thymidine demonstrating competitive-type enzyme inhibition. Results are presented as means ± SD; SD denoted by error bars (experiments carried out in triplicate).
Cell proliferation assay (MTT)
The two most active compounds against TP (compounds 2d and 2p) were evaluated for growth-inhibiting properties in two cells lines namely human umbilical vein endothelial cells (HUVEC) and epidermal carcinoma cells (A431) which are well known to overexpress epidermal growth factor receptor (EGFR), a pathway involved in the cell proliferation and angiogenesis.The inhibitory effect on cell proliferation was assessed using MTT assay after 72 h of treatment through dose-response assays performed in the 100 µM to 0.1 µM concentration range. The results were expressed as percentages of growth inhibition.Compound 2d exhibited better activity against A431 cell line (57% at 50 µM) than against HUVEC cell line (56% at 100 µM) but antiproliferative effect remains low. Compound 2p can be considerate as less cytotoxic against the two cell lines (50% and 31% at 100 µM against respectively A431 and HUVEC cells).
Conclusion
In conclusion, a small library of 38 derivatives has been synthesised and evaluated for its TP inhibiting activity. Around the pharmacophoric pyrimidinedione core of the natural substrate thymidine, four series have been designed in order to interact with wide empty pockets of the active site.The natural ligand has been annealed to a quinoline in series A (pyrimidoquinoline-2,4-diones), substituted via a methylenic chain by a quinolopyrrolidinedione in the more flexible series B and C, and the polycyclic heterocycle has been replaced by a pyrimidopyridopyrrolidinetetraone in series D. The tricyclic moieties of these new polycyclic nitrogen heterocycles have been synthesised by one-pot multicomponent reactions that involved an aniline, an aldehyde, and a 1,3-dicarbonyl derivative or an analogue.The biological evaluation identified several structurally distinct TP inhibitors in series A and D. Among them, compounds 2d, 2l, 2p, 28a, and surprisingly the open intermediate 30 showed a modest to good TP inhibition (IC50 values ranging from 26 to 87 µM) when compared to 7-DX used as a positive control. The two most active compounds 2d and 2p were shown to interact with the thymidine fixation site and to exhibit a competitive mode of inhibition towards TP.Molecular docking analysis confirmed the interaction of these newly synthesised compounds at the active binding site of TP. Moreover, docking studies highlighted a plausible specific interaction in a wide pocket that had not been yet explored.For the first time, our study showed that it is possible to inhibit TP with tricyclic heterocycles. It is worth noting that the active compounds have in common a pyrido[2,3-d]pyrimidinedione nucleus and possess a chain interacting with the same part of this pocket. Interestingly, the open intermediate 30 becomes a starting point for a novel series that will be further exploited in order to improve the activity against TP.
Experimental section
General procedures: Commercial reagents were used as received without further purification. Microwave irradiation reaction was performed with an InitiatorTM 2.0 device, Biotage. Reactions were followed with thin-layer chromatography (TLC) (using 0.20 mm silica or alumina gel 60 F254 aluminium plates, Merck) and visualisation was achieved with UV light (254 and 365 nm). Purifications were achieved through recrystallization or flash chromatography (using 40–63 µM silica, Merck). 1H and 13C NMR spectra were recorded on a Bruker AC 300 or 400 spectrometer by using DMSO-d6 as the solvent and internal standard. Chemical shifts are reported in ppm and coupling constants (J) are given in Hertz. Spin multiplicities are reported as follows: s = singlet, d = doublet, dd = doublet of doublet, m = multiplet, q = quadruplet, quint = quintuplet, t = triplet. Melting points were measured on a Stuart SMP3 melting point apparatus and are uncorrected. IR spectra were obtained on Perkin-Elmer 1600 spectrophotometer. Elemental analyses were performed at the CNRS Analysis Laboratory, Gif-sur-Yvette, France.
Pyrimido[4,5-b]quinoline-2,4 (1H, 3H) -diones (series A)
General procedure for the synthesis of compounds 2
A suspension of the requisite aniline (1.00 mmol, 1.0 eq.), barbituric acid 1 (128 mg, 1.00 mmol, 1.0 eq.), and paraformaldehyde (30.0 mg, 1.00 mmol, 1.0 eq.) was refluxed in AcOH (100 mL, method A; 300 mL, method B) or heated at 120 °C in a 1:1 mixture of AcOH/DMF (10 mL, method C). Reaction times and yields are given in Table 1.
Method A: To a solution of Na (23 mg, 1.00 mmol, 2.8 eq.) in EtOH (0.5 mL) were added 8 (115 mg, 0.36 mmol, 1.0 eq.) and urea (34 mg, 0.56 mmol, 1.6 eq.). The reaction mixture was heated under reflux for 1.5 h. The resulting brown solid was isolated by filtration, washed with a 1:1 AcOH/H2O mixture (5 mL) to afford 10 as a yellow solid (10 mg, 12% yield). Method B: 3,4-Methylenedioxyaniline (4) (411 mg, 3.00 mmol, 1.0 eq.), paraformaldehyde (270 mg, 9.00 mmol, 3.0 eq.) and 3-bromomaleimide 9 (528 mg, 3.00 mmol, 1.0 eq.) were stirred in EtOH (30 mL) at r.t. for 1 night. The reaction mixture was then heated at 60 °C for 6 h. The resulting suspension was filtered off and the solid was washed with EtOH and H2O to afford crude 10 (283 mg, 39%) which was used without further purification. A small amount was recrystallized from DMF to give an analytical sample (yellow solid). mp: 366–368 °C (decomposition). 1H NMR (400 MHz, DMSO-d6) δ: 11.50 (s, 1H), 9.01 (s, 1H), 7.88 (s, 1H), 7.57 (s, 1H), 6.35 (s, 2H) ppm. 13C NMR (100 MHz, DMSO-d6) δ: 170.2, 169.3, 153.6, 151.2, 151.0, 141.3, 134.6, 123.2, 119.0, 106.3, 103.8, 98.7 ppm. IR ν: 3449, 2956, 2723, 1771, 1720, 1623, 1573, 1471, 1411, 1330, 1265, 1100, 1030, 935, 872, 790, 748 cm−1. Anal. Calcd. for C12H6N2O4 · 0.25 H2O (246.69): C, 58.42; H, 2.66; N, 11.36. Found: C, 58.35; H, 2.49; N, 11.35%.
General method for preparation of compounds 12 and 15
To a suspension of thymine (11) or 5-bromouracile (14) (4.00 mmol, 1.0 eq.) in MeCN (6 mL) was added BSA (2.5 mL, 10.00 mmol, 2.5 eq.) at r.t. When the reaction mixture became clear (15 min), the ω-dibromoalkyl derivative (6.00 mmol, 1.5 eq.) and I2 (cat. amount) were added. The mixture was refluxed for 2 h and then either maintained at r.t. for 1 night (12a,b) or refluxed for 2 days (15). After elimination of the solvent under vacuum, the residue was stirred with water (10 mL), isolated by filtration, and then purified by flash column chromatography.
To a suspension 21 (482 mg, 3.00 mmol, 1.5 eq.) in anhydrous MeCN (10 mL) was added BSA (1.5 mL, 6.00 mmol, 3.0 eq.) at r.t. When the reaction mixture became clear (15 min), a suspension of imide 10 (484 mg, 2.00 mmol, 1.0 eq.) and K2CO3 (415 mg, 3.00 mmol, 1.5 eq.) in anhydrous MeCN (10 mL) was added and the mixture was refluxed for 5 days. Then, another amount of 21 (241 mg, 1.50 mmol, 0.75 eq.) and BSA (0.85 mL, 3.48 mmol, 1.74 eq.) in MeCN (5 mL) were added and the reflux was maintained for 2 more days. After evaporation of the solvent under vacuum, the residue was stirred with H2O (10 mL) and an aqueous 1 N HCl solution was added until a pH = 3 was reached. The resulting solid was isolated by filtration, washed with H2O and EtOH. Recrystallization from DMF afforded 22 as a yellow solid (579 mg, 79% yield). mp: 407–409 °C (decomposition). 1H NMR (400 MHz, DMSO-d6) δ: 11.07 (s, 1H), 11.03 (s, 1H), 9.10 (s, 1H), 7.89 (s, 1H), 7.60 (s, 1H), 6.37 (s, 2H), 5.55 (s, 1H), 4.51 (s, 2H) ppm. 13C NMR (100 MHz, DMSO-d6) δ: 167.8, 167.0, 163.8, 153.4, 151.3, 150.9, 150.8, 150.7, 140.8, 133.8, 122.1, 118.5, 105.9, 103.4, 98.2, 97.1, 37.4 ppm. IR ν: 3185, 3106, 3056, 3001, 2817, 1772, 1702, 1650, 1494, 1455, 1414, 1392, 1341, 1322, 1307, 1264, 1230, 1177, 1104, 1030, 1014, 930, 866, 847, 828, 794, 752 cm−1. Anal. Calcd. for C17H10N4O6 · 0.5 H2O (375.29): C, 54.41; H, 2.95; N, 14.93. Found: C, 54.10; H, 3.05; N, 14.94%.
To a suspension of 27 (70.0 mg, 0.30 mmol, 1.0 eq.) in DMF (444 µL), the requisite amine (0.60 mmol, 2.0 eq.) was added. The reaction mixture was refluxed for a time T1. To the resulting solution, PTSA (115 mg, 0.60 mmol, 2.0 eq.) was added and reflux was maintained during a time T2. The resulting precipitate was filtered off and washed to give a crude pure compound.
To a suspension of 27 (232 mg, 1.00 mmol) in DMF (10 mL), Et3N (0.697 mL, 5.00 mmol, 5 eq.) and 29, (533 mg, 3.00 mmol, 3 eq.) were added and the mixture was heated at 90 °C overnight. Once cooled at r.t., the reaction mixture was poured on Et2O (10 mL) and maintained at 4 °C for 2 h. To a solution of the resulting solid in H2O (12 mL), a 1 N HCl solution was added until a pH < 3 was reached. The solid was filtrated, washed successively with H2O and EtOH to afford pure 30 as a white solid (410 mg, 82% yield). mp: 259–261 °C (decomposition). 1H NMR (400 MHz, DMSO-d6) δ: 12.10 (s, 1H), 11.70 (s, 1H), 10.97 (s, 2H), 10.90 (s, 1H), 10.82 (s, 1H), 9.10 (t, J = 5.8 Hz, 1H), 8.96 (t, J = 6.1 Hz, 1H), 8.58 (s, 1H), 5.53 (s, 1H), 5.48 (s, 1H), 4.18–4.07 (m, 4H) ppm. 13C NMR (100 MHz, DMSO-d6) δ: 166.0, 165.3, 164.0, 161.6, 156.9, 153.4, 153.3, 152.4, 151.4, 151.3, 150.3, 136.8, 124.6, 109.8, 96.8, 96.6, 39.3, 38.7 ppm. IR ν: 3563, 3489, 3247, 3118, 2988, 2819, 1689, 1657, 1609, 1583, 1527, 1502, 1419, 1357, 1321, 1284, 1233, 1020, 841, 819, 766, 743, 703 cm−1. Anal. Calcd. for C19H15N9O8 · 3.5 H2O (560.43): C, 40.72; H, 3.96; N, 22.49. Found: C, 40.97; H, 3.73; N, 22.02%.
Biologicals
A spectrophotometric assay method was adopted to evaluate in vitro E. coli TP inhibiting activity of all the synthesised compounds. The conversion of thymidine to thymine was recorded at 290 nm. Thymidine was dissolved in buffer (3 mM). DMSO stock solutions of the tested compounds were diluted in buffer, the incubation medium contained 2.5% of DMSO. All experiments were conducted in triplicate.
In vitro thymidine phosphorylase essay
Initially, all synthesised compounds were tested against E. coli TP at 50 µM or at the maximum concentration allowed by their solubility.The method reported by Khan et al. has been adapted to a screening in 96-well plates (200 µL well, UV-Star® 96-well microplate, Greiner). Each well contained 160 µL of 50 mM KH2PO4 buffer (pH 7.4), 0.006 U of E. coli TP (T2807-1KU, Sigma Aldrich), and 20 µL of the tested compound solution. The reaction was initiated by the addition of 20 µL thymidine (i.e. a 300 µM final concentration corresponding to the Km). The decrease in absorbance due to conversion of thymidine to thymine was followed during 30 min (PowerWave HT Microplate Reader, Biotek) at 25 °C, with one measure every 2 min and a 30 s shaking before each measure. Data were processed using Biotek KC4 software and the results were expressed in percentage of inhibition at 50 µM or at the maximum concentration allowed by the solubility of the tested compound (Tables 2–4). For the most active and soluble compounds 2d, 2l, 2p, 28a, and 30 (inhibition > 30% and solubility > 50 µM), IC50 values have been determined.
Phosphate competition study
For this study, the reaction mixture was composed of 160 µL of 10 mM Tris buffer (pH 7.4), 150 mM NaCl, 2 mM EDTA, 0.006 U of E. coli TP, different concentrations of KH2PO4 (2, 5, 10 and 30 mM). Compounds were tested at a concentration close to their IC50. The reaction was initiated by the addition of 20 µL of thymidine (3 mM and 10 mM solutions).
Enzyme inhibition kinetic study
The TP inhibiting activity at differents concentrations of compound 2d (0, 15, 20, 30 µM) and 2p (0, 4, 6, 8 µM) was evaluated in the presence of different concentrations of thymidine (200, 250, 300, 500, 750, 1000 mM) in KH2PO4 buffer (pH 7.4).
In vitro anti-proliferative assay
Cell viability was determined on HUVEC and A431 cell lines using 3–(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT; Sigma-Aldrich) assay after treatment with compounds 2d and 2p. A total of 5 × 103 cells/well (A431 cell line) or 104 cells/well (HUVEC cell line) were plated in a 96-well plate and treated, in duplicates, with compounds 2d and 2p at five concentrations (0.1 µM, 1 μM, 10 μM, 50 μM, 100 μM) for 72 h. Absorbance was read at 570 nm using Labsystem Multiskan MS microplate reader. The results were expressed as percentages of growth inhibition compared to positive control cell growth (100%).
Molecular docking study
Molecular modelling studies were carried out with the GOLD software (version 5.1). The ligands were constructed using the standard fragments of the Sybyl library 6.9.1. Their geometry was optimised using the Tripos force field by assigning partial loads calculated by the Gasteiger–Hückel method to a gradient of 0.0001 Kcal/mol/Å.The conformations obtained for each compound were classified by a consensus scoring function and then the conformation groups were observed visually to evaluate their consistency and their complementarity with the active site. Only the most representative conformations of each group were selected.In order to check the validity of the results, the TPI reference inhibitor was relocated to the hTP binding site (PDB code: 1UOU) following the same protocol as for the inhibitors. The best solution has an RMS on heavy atoms of 0.654 Å with respect to the 1UOU crystallographic structure. Results were visualised by Chimaera and Maestro.Supplementary data includes docking superposition of compounds 2d, 2p, and 2l in the hTP active site; two docked conformations of compound 2c in the hTP active site; docking superposition of compounds 2i, and 23b in the hTP active site; NMR spectra of the final compounds.Click here for additional data file.
Authors: F K Brown; P J Brown; D M Bickett; C L Chambers; H G Davies; D N Deaton; D Drewry; M Foley; A B McElroy; M Gregson Journal: J Med Chem Date: 1994-03-04 Impact factor: 7.446
Authors: J Balzarini; A E Gamboa; R Esnouf; S Liekens; J Neyts; E De Clercq; M J Camarasa; M J Pérez-Pérez Journal: FEBS Lett Date: 1998-10-30 Impact factor: 4.124
Authors: I V Bijnsdorp; F Capriotti; F A E Kruyt; N Losekoot; M Fukushima; A W Griffioen; V L Thijssen; G J Peters Journal: Br J Cancer Date: 2011-03-08 Impact factor: 7.640
Authors: Eirini Mitsiki; Anastassios C Papageorgiou; Shalini Iyer; Nethaji Thiyagarajan; Steven H Prior; Darrell Sleep; Chris Finnis; K Ravi Acharya Journal: Biochem Biophys Res Commun Date: 2009-06-23 Impact factor: 3.575