Julia
Fernández de Luco‡
a,
Alejandro I.
Recio-Balsells‡
a,
Diego G.
Ghiano‡
a,
Ana
Bortolotti
b,
Juán Manuel
Belardinelli§
b,
Nina
Liu
c,
Pascal
Hoffmann
d,
Christian
Lherbet
d,
Peter J.
Tonge
c,
Babu
Tekwani¶
e,
Héctor R.
Morbidoni
*bf and
Guillermo R.
Labadie
*ag
aInstituto de Química Rosario, UNR, CONICET, Suipacha 531, S2002LRK, Rosario, Argentina. E-mail: labadie@iquir-conicet.gov.ar; Fax: +54 341 4370477; Tel: +54 341 4370477
bLaboratorio de Microbiología Molecular, Facultad de Ciencias Médicas, Universidad Nacional de Rosario, Santa Fe 3100, S2002KTR, Rosario, Argentina. E-mail: morbiatny@yahoo.com
cDepartment of Chemistry, Institute of Chemical Biology & Drug Discovery, Stony Brook University, Stony Brook, NY 11794, USA
dLSPCMIB, UMR-CNRS 5068, Université Paul Sabatier-Toulouse III, 118 Route de Narbonne, 31062 Toulouse Cedex 9, France
eNational Center for Natural Products Research & Department of Biomolecular Sciences, School of Pharmacy, University of Mississippi, MS 38677, USA
fConsejo de Investigaciones, Universidad Nacional de Rosario, Argentina
gDepartamento de Química Orgánica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, Suipacha 531, S2002LRK, Rosario, Argentina
First published on 5th November 2020
Triclosan and isoniazid are known antitubercular compounds that have proven to be also active against Leishmania parasites. On these grounds, a collection of 37 diverse 1,2,3-triazoles based on the antitubercular molecules triclosan and 5-octyl-2-phenoxyphenol (8PP) were designed in search of novel structures with leishmanicidal activity and prepared using different alkynes and azides. The 37 compounds were assayed against Leishmania donovani, the etiological agent of leishmaniasis, yielding some analogs with activity at micromolar concentrations and against M. tuberculosis H37Rv resulting in scarce active compounds with an MIC of 20 μM. To study the mechanism of action of these catechols, we analyzed the inhibition activity of the library on the M. tuberculosis enoyl-ACP reductase (ENR) InhA, obtaining poor inhibition of the enzyme. The cytotoxicity against Vero cells was also tested, resulting in none of the compounds being cytotoxic at concentrations of up to 20 μM. Derivative 5f could be considered a valuable starting point for future antileishmanial drug development. The validation of a putative leishmanial InhA orthologue as a therapeutic target needs to be further investigated.
Leishmaniasis is one of the most concerning vector-borne NTDs, an expanding endemic disease caused by the protozoan parasite Leishmania spp. and transmitted by over 90 phlebotomine sand fly species. Visceral leishmaniasis (VL), or kala-azar, is the most severe form of the disease which, if not treated, has a mortality rate close to 100%, surpassing 30000 deaths per year. Unfortunately, there is no vaccine available and the regular treatment includes drugs that are toxic to the host, so they are used at low doses ending in the generation of drug resistant forms of the parasite.
The World Health Organization (WHO) has reported a rate of 1.3 million cases of leishmaniasis per year and 350 million people exposed to infections.1 This fact has alerted the scientific community and driven it to the urgent search for new effective and less toxic drugs with novel structures and different targets.
Bacterial fatty acid biosynthesis is a validated target for drug discovery;2 however, it has not been fully exploited. An interesting target to explore for the development of new antileishmanial drugs is the ortholog of the enzyme 2-trans-enoyl fatty acid reductase3 from the Mycobacterium tuberculosis FASII complex. This enzyme, called InhA, was discovered as a target for the drug isoniazid (INH) (Fig. 1), which is one of the current drugs for the treatment of tuberculosis.4 InhA catalyzes the reduction step that ends the process of a two carbon extension in the pathway of saturated fatty acid synthesis. Interestingly, antituberculosis drugs, such as triclosan (TRC) – an uncompetitive inhibitor of InhA – (Fig. 1), have also shown in vitro activity.5 The kinetoplastid type II FAS is found in mitochondrial organelles.6 Thus, since parasites have FASII enzymes they may be targeted by drugs active against bacterial fatty acid synthesis.
The M. tuberculosis FAS II complex is analogous to other bacterial FAS II, but is not capable of de novo fatty acid synthesis. However, FAS II from Mycobacterium is capable of elongating acyl-CoA generated by FAS I producing very long chain α-hydroxy, β-alkyl branched fatty acids known as mycolic acids. Even though mutations in the inhA gene are known to confer resistance to INH, the enzyme remains a validated and reliable target for drug discovery.
Triclosan7 has been described as an InhA targeting antitubercular drug which does not require KatG activation; however, its poor bioavailability prevented its use.8 In spite of that, TRC was used to inspire the development of other scaffolds such as 5-alkyl-diphenyl-ethers, one of which, 5-octyl-2-phenoxyphenol (8PP, Fig. 1),9 has good anti-tubercular activity and higher solubility. Therefore, 8PP and similar molecules are good starting points for structure-based drug discovery.
One approach to enhance the activity of a drug candidate is to incorporate additional heterocycles as building blocks.10 A recently introduced and popular strategy is to use 1,2,3-triazoles to link structurally diverse fragments with different pharmacophores.11–13 A straightforward way to prepare 1,2,3-triazoles is the well-known click chemistry reaction which involves azides and alkynes as starting materials. This kind of approach is a very efficient way to build new compound libraries as drug candidates.14 Also, when used as linkers, 1,2,3- triazoles are inert to degradation, soluble and prone to bind to biomolecular targets.13
Stec et al. have already taken advantage of this approach by preparing a collection of derivatives exchanging one of the TRC phenyl rings by 1,2,3-triazoles, introducing ketones, amides and isoxazoles.15 One of those analogs with a 1,2,3-triazole on position 4′ and a butyl substituent in position 4 of the triazole showed micromolar activity against MTB and 98% inhibition of pure InhA (Fig. 1).
We have previously reported bioactive 1,2,3-triazolyl derivatives against Leishmania parasites (amino sterols) and M. tuberculosis (fatty acid derivatives).16,17 In this opportunity, our aim was to develop novel structures based on the structures of TRC and 8PP, to find new chemical entities as candidates for antileishmanial drug discovery. The presence of putative enoyl-ACP reductase encoding genes in Leishmania prompted us to prepare and assay a library of 1,4 disubstituted 1,2,3-triazolyl catechols. The collection was initially assayed against Leishmania donovani promastigotes and M. tuberculosis H37Rv to determine IC50 and MIC. Cytotoxicity against Vero cells was also evaluated and afterwards, to understand the mechanism of action of these catechols and to test our hypothesis of an InhA ortholog being the target of these molecules, the most promising compounds were tested against purified Mtb InhA.
The search in the TDR Targets database21 for enoyl-[acyl-carrier-protein] reductases (ENRs) provided an ORF for an InhA-like ortholog present in the Leishmania major genome (TDR Targets ID: 27659; GeneDB ID: LmjF.34.0610). Moreover, analysis of the L. donovani genome by Ravooru et al.22 revealed that the E9B7Z4 encoded ORF is closely associated with the enzyme of the class trans-2-enoyl-CoA reductase (1.3.1.38), which was also well conserved along Leishmania spp. The activity of TRC against the intracellular and free form of L. donovani has been recently linked to ENRs by an in silico modelling based on the GenBank CBZ37546.1 sequence.3 Based on these precedents we performed a multiple alignment of M. tuberculosis, E. coli ENRs and different Leishmania spp. putative ENR sequences. The result showed highly conserved residues on the substrate binding and active site loops. Also, as was shown before, putative ENR proteins of different Leishmania species were highly conserved (Fig. S1, ESI†). These precedents led us to design a collection of compounds based on TRC and 8PP that could potentially target leishmanial ENRs (Fig. 1).
As was previously mentioned, the 1,2,3-triazole motif is widely used in medicinal chemistry, being present in compounds with anticancer, anti-inflammatory, antituberculosis, antitrypanosomal, antibacterial and antiviral activities. Most of the reported bioactive 1,2,3-triazoles are exclusively 1,4-disubstituted, with little information on the results of expanding the chemical space by obtaining their 1,5-isomers.23 Particularly, there are only a couple of examples of libraries of 1,4,5-trisubstituted 1,2,3-triazoles which have been assayed for antileishmanial or antituberculosis activity.24 One of those libraries was composed of triazolopyridopyrimidines which have shown a moderate activity against different Leishmania species.24
Aiming to expand the SAR studies of the TRC scaffold, a new set of 1,2,3-triazolyl analogs was synthesized using thermal and Cu(I) catalyzed-cycloaddition (CuAAC) reactions. The CuAAC reaction is completely selective being extremely useful for the synthesis of 1,4-disubstituted-1,2,3-triazoles. The thermal cycloaddition provides a mixture of 1,5- and 1,4-regioisomers, being useful to expand our SAR studies (Scheme 1). Two series of 1,2,3-triazole analogs of TRC were proposed (Scheme 1, series 1). Due to the critical role of the TRC hydroxy group, the designed compounds contained this group in different scaffolds.25 Additionally, the chlorine was removed since it is known that it diminishes the inhibitory activity and interferes with the bioavailability.26 The series 1 scaffold has a phenol connected by a linker to a substituted 1,2,3-triazole that was introduced by a CuAAC reaction. In this series, the linker can be either a methyleneoxy group preserving the catechol moiety of TRC or a methylene group (Scheme 1, series 1). The methyl ether derivatives were also included for the sake of comparison with the free hydroxy derivatives.
Scheme 1 1,4-Disubstituted- and 1,4,5-trisubstituted-1,2,3-triazoles generated by CuAAC reaction and thermal cycloaddition. |
To increase the diversity of the collection expanding the explored chemical space, 1,4,5-trisubstituted-1,2,3-triazoles were also prepared (Scheme 1, series 2). These 1,4,5-trisubstituted-1,2,3-triazoles were prepared from different azides and internal alkynes by thermal cycloaddition, resulting in mixtures of isomers.
The synthesis started with the phenol etherification of catechol 1 and guaiacol 2 according to Williamson's reaction, producing the ether formation with propargyl bromide and potassium carbonate in methanol (Scheme 2). O-Propargyl-catechol 3 and O-propargyl-2-methoxyphenol 4 were obtained in 46% and 40% yield, respectively. The reaction with catechol 3 also provided the di-propargylated product in 13% yield as a by-product. In both reactions, some of the starting material was recovered and could be reused. Once the alkyne intermediates were obtained, the CuAAC reaction was performed with different azides containing aromatic and carbocyclic rings, and aliphatic, prenyl and alkyl ethyl ester chains (5a–i, 6a–e, Scheme 2). The average yield for scaffold 5 was 74% and 77% for 6 (Table 1).
Comp. | Yield (%) | R | Comp. | Yield (%) | R |
---|---|---|---|---|---|
5a | 76 | CH2Ph | 6a | 70 | CH2Ph |
5b | 81 | (CH2)3Ph | 6b | 82 | (CH2)3Ph |
5c | 71 | Cinnamyl | 6c | 73 | Cyclohexyl |
5d | 68 | Cyclohexyl | 6d | 79 | C8H17 |
5e | 65 | C8H17 | 6e | 80 | (CH2)COOEt |
5f | 82 | C13H27 | 9a | 76 | C3H7 |
5gZ | 70 | Neryl | 9b | 81 | C5H11 |
5gE | 70 | Geranyl | 9c | 71 | C8H17 |
5h | 79 | (CH2)COOEt | 9d | 68 | Ph |
5i | 78 | (CH2)4COOEt | 9e | 65 | (CH2)2Ph |
The other members of series 1 were prepared with salicyl azide 8. This key intermediate was prepared from commercially available salicyl alcohol 7 by reaction with sodium azide and triphenylphosphine and in DMF and carbon tetrachloride.27 The CuAAC reaction was performed with commercially available terminal alkynes. The alkynes used were 1-propyne, 1-pentyne, 1-octyne, ethynylbenzene and pent-4-yn-1-ylbenzene (Scheme 2). The average yield of the reaction was 65% (Table 1).
Finally, a thermal cycloaddition reaction in toluene was used to prepare series 2 with the aim of mimicking the phenol aromatic ring by a hydroxylated triazole motif. Three different azides and three internal alkynes were employed obtaining a mixture of 1,4,5-substituted-1,2,3-triazoles. Benzyl azide 10, 3-phenylpropyl azide 12 and 2-(azidomethyl)naphthalene 14 were combined with the internal alkynes oct-2-yn-1-ol, 3-phenylprop-2-yn-1-ol, and hex-3-yn-2-ol to afford the products 11, 13 and 15 (Scheme 3). Phenylpropyl azide 12, with a longer linker, is more flexible than benzyl azide 10, and the naphthyl derivative is bulkier and has an extended pi-stacking surface.
The yields were calculated including both isolated isomers, which were obtained as a 1:1 mixture of both regioisomers. The ratio was calculated by 1H NMR. The average yield for each scaffold was 61%, 51% and 77% for triazoles 11a–e, 13a–f and 15a–f, respectively (Table 2).
Comp | R1 | R2 | Comp | R1 | R2 | Yld (%) |
---|---|---|---|---|---|---|
11a | CH2OH | C5H11 | 11b | C5H11 | CH2OH | 52 |
11c | Ph | CH2OH | 76 | |||
11d | CH(OH)CH3 | C2H5 | 11e | C2H5 | CH(OH)CH3 | 56 |
13a | CH2OH | C5H11 | 13b | C5H11 | CH2OH | 46 |
13c | CH2OH | Ph | 13d | Ph | CH2OH | 66 |
13e | CH(OH)CH3 | C2H5 | 13f | C2H5 | CH(OH)CH3 | 42 |
15a | CH2OH | C5H11 | 15b | C5H11 | CH2OH | 69 |
15c | CH2OH | Ph | 15d | Ph | CH2OH | 94 |
15e | CH(OH)CH3 | C2H5 | 15f | C2H5 | CH(OH)CH3 | 75 |
The collection of 37 triazolyl TRC analogs were characterized using 1H and 13C NMR and mass spectrometry. The regiochemistry of the series 2 products was determined using nuclear Overhauser effect (NOE) experiments.
For example, for product 11a the signal of the methylene of the methylenehydroxy group at 4.53 ppm showed a NOE with the signal of the methylene of the benzyl group at 5.59 ppm and with the methylene from the pentyl group at 2.56 ppm (Fig. 2). For isomer 11b, irradiation on the same methylene only showed a NOE with the methylene of the pentyl chain (Fig. 2).
Comp | Ld IC50 (μM) | Mt H37Rvb MIC (μM) | Vero IC50 (μM) | Mt InhAc inh (%) |
---|---|---|---|---|
TRC = triclosan, PM = pentamidine, AMB = amphotericin B. a Ld: L. donovani promastigotes, evaluated by the Alamar blue assay. b Mt H37Rv: Mycobacterium tuberculosis H37Rv, evaluated by a serial microdilution colorimetric assay using MTT as a viability indicator. c Mt InhA: remaining M. tuberculosis InhA activity. d N.I.: no inhibition. | ||||
5f | 13.4 | >100 | >20 | 63 |
5gZ | 55.0 | >100 | >20 | 85 |
5gE | 51.9 | >100 | >20 | N.I. |
9a | >60 | 20.0 | >20 | 92 |
9b | >60 | 20.0 | >20 | N.I.d |
9d | >60 | 20.0 | >20 | N.I. |
15b | 31.4 | >100 | >20 | N.I. |
TRC | 43.0 | |||
PM | 4.40 | |||
AMB | 0.11 |
Only four compounds of the series showed IC50 below 60 μM against L. donovani, with 5f being the most active compound showing an IC50 of 13.4 μM (Table 3). Regarding the M. tuberculosis assays, the results showed that only analogs 9a, 9b and 9d have MIC below 100 μM. These 3 compounds belong to series 1 and have MIC of 20.0 μM (Table 3).
Having completed the screening stage, the next step was validating the target of the compounds. To do this, the most promising analogs were assayed on the purified mycobacterial InhA.19
The initial assay was performed at a fixed compound concentration, 40 μM for derivatives 5f, 5gZ, and 5gE and 50 μM for derivatives 9a, 9b, 9d and 15b. The reaction velocity in the absence (v)0 and presence of the inhibitor (v)I was determined. The results were disappointing since the tested compounds proved to be extremely poor in inhibitory efficacy with only marginal inhibition shown by 3 of the most promising compounds compared to TRC which showed an IC50 of 1.0 μM.9 The remaining enzyme activities were 63, 85 and 92% for compounds 5f, 5gZ and 9a, respectively (Table 3).
When discussing the most active compounds against L. donovani, the disappointing results could be explained with the fact that they were assayed on the mycobacterial enzyme and not on the actual L. donovani orthologue. There are substantial differences between Leishmania spp. and M. tuberculosis fatty acid synthesis end products. While the parasite fatty acids are composed of C14 to C22 fatty acids,29 the mycobacterial FASII could elongate fatty acids up to ∼40 carbons.30 These variations would probably justify the sequence differences between the M. tuberculosis and Leishmania spp. ENRs that will ultimately may result in the lack of inhibition of InhA. Consequently, we should explore this in detail to obtain a definitive conclusion of whether minor but structurally important differences between the different ENRs may condition the inhibitory activity.
Also, there is low or null inhibition of the most promising compounds against M. tuberculosis. The compounds that display some activity, analogs 5f, 5gZ and 9a, share the critical phenolic hydroxyl on their structure with TRC and 8PP (Table 3), nevertheless they did not show an inhibitory effect on InhA. Those results could be explained by the presence of different target/s, a hypothesis that could be further investigated.
As regards L. donovani promastigotes, two populations of active compounds were identified. The first group comprises analogs from series 1 (5f, 5gZ, 5gE; IC50 = 13.4, 55.0 and 51.9 μM respectively); while the second consists of only compound 15b (IC50 = 31.4 μM) from series 2 (Fig. 4). The active compounds contain long lipophilic substituents on the triazole (pentyl, octyl, geranyl, neryl and tridecyl) and their activity seems to be related to a triazole side chain longer than 5 carbons. The restricted conformation produced by the double-bond on the isoprenyl derivatives seems to be detrimental to the activity. Comparison of analogs 15b and 5f clearly showed that the longer aliphatic saturated chain improved the activity. This structural feature resembles the 5-alkyl-diphenyl ether moiety present in 8PP where the saturated lipophilic chain plays a critical role in the binding to InhA.9 Three out of the four active compounds have a phenol in their structure, and the fourth compound has a methylenehydroxy group attached to the triazole. None of the methylated phenol derivatives was active, demonstrating the importance of the hydroxy group in the activity. The fact that three active compounds are from series 1 confirms that the incorporation of the ether on the structure mimics TRC and 8PP.
On the other hand, the results on M. tuberculosis show that compounds 9a and 9b have a lipophilic saturated aliphatic chain on the 1,2,3-triazole (Fig. 5) being like 8PP, but notably less active (8PP MIC99 = 6.4 μM). The other active analog is 9d, which has a phenyl ring decorating the 1,2,3-triazole, a fact that shows that the flexible lipophilic chain can be replaced by that non-polar bulkier ring. Surprisingly, derivative 9c that holds an octyl chain, the derivative that is structurally the most closely related to 8PP, was inactive. The activity vs. aliphatic tail relationship for the C5 derivative is less significant than for the C8 derivative when comparing the diphenyl ethers.9 Ethers 5e and 6d, which also resemble 8PP, were also inactive. The rest of the library does not show activity against M. tuberculosis H37Rv with MICs higher than 32 μM. These compounds are considerably less active than the previous collection based on 8PP reported by our group, which reached low micromolar MIC on H37Rv.16,31 Nevertheless, the active members of the collection are twice more active than TRC that has a MIC of 43.0 μM.
A close comparison with the triazolyl-diphenyl ethers reported by Spagnuolo et al.32 showed that not only the substitution pattern on the phenol ring, but mostly the presence of the phenoxy moiety, plays a critical role in the antitubercular activity (Fig. 5). This is clearly seen for compounds 9a, 9b, and 9d that are 1.8, 5.1 and 9.9 times less active than the diphenyl ethers PT163, PT509 and PT510, respectively (Fig. 5).
Lipinski's rule of five is widely used for determination of good pharmacokinetic properties. The parameters considered were the molecular weight (MW < 500 Da), the partition coefficient between water and octanol (0 < logP < 5) and the number of hydrogen bond acceptors (nNHOH < 10) and donors (nON < 5). If there were two or more violations of these rules, the compounds would present poor oral administration, but there are exceptions to these rules. Another useful parameter is the topological polar surface area (TPSA) that considers the electron density distribution in the molecule and gives us an idea of the capacity to penetrate biological membranes. This parameter is commonly used to establish if the compounds could penetrate the blood brain barrier. All the active compounds comply with the rule of five and TPSA parameters. Only compound 5f violates one rule having a logP of 5.67. Interestingly, this compound is the most active in the collection against L. donovani promastigotes. This result is not discouraging since the physicochemical requirements could be taken in a lighter way in order to obtain antiparasitic lead compounds susceptible of further improvement.34
The assay was performed in 96-well tissue culture-treated plates as described earlier.38 The cells were seeded to the wells of the plate (25000 cells per well) and incubated for 24 h. Samples were added and the plates were again incubated for 48 h. The number of viable cells was determined by the neutral red assay. IC50 values were determined from logarithmic graphs of growth inhibition versus concentration. Doxorubicin was used as a positive control (IC50 = 14 mM, Vero cells), while DMSO was used as a vehicle control.
In summary, compound 5f could be considered a valuable starting point for future antileishmanial drug development, based on its simple preparation and promising activity. The possibility that the activity could be linked to a putative leishmanial InhA orthologue should be a matter of future investigation, together with the exploration of new structural modifications to enhance the activity. For this purpose, the purification of the putative Leishmania ORF and its assay for ENR activity would lead the priority list. This will be followed by testing the analogues against this potential InhA orthologue. In parallel, the most promising candidates will be assayed against other Leishmania species including the intracellular stage of the parasite.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0md00291g |
‡ These authors contributed equally to this work. |
§ Present address: Department of Microbiology, Immunology & Pathology, Colorado State University, 1682 Campus Delivery, Fort Collins, CO 80523, USA. |
¶ Present address: Department of Infectious Diseases, Division of Drug Discovery, Southern Research, Birmingham, AL 35205, USA. |
This journal is © The Royal Society of Chemistry 2021 |