1. Introduction
Parasitic diseases represent an overwhelming health problem for impoverished populations living in developing countries with poor sanitary conditions. Previous reports have shown that the three most important protozoan diseases have an estimated disability-adjusted life year (DALY) of approximately one million [
1]. These diseases include malaria, leishmaniasis and trypanosomiasis, which are caused by vector-borne protozoan parasites, including
Plasmodium,
Leishmania and
Trypanosoma species, respectively. In malaria, for instance, the convergence of the COVID-19 pandemic and escalating medication resistance has propelled this pathological condition into the realm of the world’s deadliest diseases [
2], with an estimated 249 million cases and 608 000 malaria deaths in 85 countries in 2022. In 2022, the high rates of malaria burden was endorsed by the World Health Organization (WHO)’s African region, with 94% of malaria cases (233 million) and 95% (580 000) of malaria deaths [
3]. Core activities for the management of malaria include vector control and treatment of patients with appropriate antimalarial drugs, all impaired by the endless development of mosquito and
Plasmodium spp. resistant strains and vaccine administration [
2]. However, the vaccine approach has shown little success as the world’s first RTS, and the S vaccine offers only 30% protection [
4]. Leishmaniasis and trypanosomiasis are neglected tropical diseases (NTDs) that are fatal if left untreated [
5]. Leishmaniasis is endemic in Africa, Asia, the Americas, and the Mediterranean region [
6], where it can develop into 3 main clinical forms according to the involved
Leishmania species. These include cutaneous and mucocutaneous (
Leishmania major,
L. tropica and
L. braziliensis, etc.) and the most severe visceral (
L. donovani, and
L. infantum, etc.) forms [
7]. Annually, 700 000 to 1 million people are newly infected with leishmaniasis, with 20 000 to 30 000 deaths [
7]. The treatment options for leishmaniasis include pentavalent antimonials, sodium stibogluconate (pentostam), and meglumine antimoniate (glucantime) as first-line drugs and pentamidine, amphotericin B, paromomycin and miltefosine as second-line treatments [
8]. African trypanosomiasis, which affects both humans (human African trypanosomiasis) and animals (animal African trypanosomiasis) is a health concern in endemic areas. Human African trypanosomiasis (HAT) threatens millions of people in 36 sub-Saharan African countries, with 55 million people at risk of being infected [
9]. With the multiple control strategies initiated, an important decline in the number of new cases from approximately 40 000 in 1998 to 663 in 2020 has been recorded. However, important efforts should be made to completely fulfil the WHO’s initiative toward the eradication of sleeping sickness by 2030 [
9].
Trypanosoma brucei gambiense and
T. brucei rhodesiense are the main pathogens responsible for HAT. Animal African trypanosomiasis (AAT) is among the most important diseases in cattle with severe economic consequences [
10]. Indeed, the disease was reported to cause 3 million deaths in cattle [
11]. AAT is mainly caused by three species of the genus
Trypanosoma, including
T. brucei gambiense,
T. vivax, and
T. congolense. Pentamidine, eflornithine, nifurtimox, fexinidazole, suramin, and melarsoprol are among the current treatments of HAT, whereas homidium bromide, diminazene aceturate, and suramin are the main drugs for AAT [
9]. However, these drugs reveal poor efficacy, unacceptable toxicity, and drug resistance as the main limitations that make them inappropriate treatments [
10,
12]. The currently developed acoziborole drug, which showed 95% success rate and effectiveness during phases 2 and 3 trials against HAT is noteworthy; however, recent development of
Trypanosoma resistance to this drug has been identified [
13]. Thus, there is a crucial need to search for alternative treatments against
Trypanosoma infections. In spite of the recent advances in research to control these infectious diseases, they remain prevalent, supporting the need to search for effective and safe treatments against malaria, leishmaniasis and trypanosomiasis. The natural origin of antimalarial drugs is undeniable as quinine and artemisinin are some of the foremost examples that were respectively identified from cinchona bark [
14] and
Artemisia annua [
15]. Moreover, synthesis and structural modifications of natural product scaffolds have provided a number of active principles for antileishmanial (miltefosine) [
16,
17] and antitrypanosomal (fexinidazole) [
9,
18] drug development. The mechanistic basis of antiprotozoal action of these active principles revealed inhibition of a number of enzymes, such as dihydrofolate reductase, which are crucial for the survival and virulence of
Plasmodium [
19],
Leishmania and
Trypanosoma [
20] species. According to reported studies, potassium channels are integral membrane proteins that are intricately involved in the maintenance of vital parameters, including cellular membrane potential and cell volume in malaria parasites [
21]. Noteworthy, the inhibitory effect of glibenclamide [
22] and a number of halogenated glucose analogs [
23] vis-à-vis K+ channels is reported in
Leishmania sp. and
Trypanosoma brucei, respectively. Thus, it is not unreasonable to speculate that compounds or drugs that inhibit dihydrofolate reductase or K+ channels can afford potentially active antiprotozoal leads against malaria, leishmaniasis, and African trypanosomiasis. Indeed, target-based drug discovery and drug repositioning can afford potentially active compounds that can serve as scaffolds for drug development, thereby reducing the overall cost and time frame used in the traditional drug discovery process [
24]. Based on the foregoing, this study aims to investigate the antiprotozoal activity of a series of 3 heterocyclic potential DHFR inhibitors (compound X, Y and Z) (
Figure 1) and 1 potassium channel blocker (E4031) (
Figure 1,
Table 1) against
Plasmodium,
Leishmania and
Trypanosoma parasites.
3. Results and discussion
Whole-cell phenotypic screening of compounds with known inhibitory effects against validated therapeutic targets is called target-based drug discovery. This approach includes target identification and validation using a number of tests, such as
in silico prediction, biochemical and genetic analyses to identify proteins or enzymes that are crucial for the survival and virulence of the parasites [
35]. In this regard, several targets, including dihydrofolate reductase have been validated as important components involved in the pathogenesis of
Plasmodium,
Leishmania and
Trypanosoma species [
36]. In fact, dihydrofolate reductase aids in the replication of several microorganisms by reducing dihydrofolate into tetrahydrofolate for DNA synthesis [
37,
38]. Moreover, the implication of K+ channels in the survival and virulence of these parasites is noteworthy. According to the literature, a number of anticancer drugs, such as methotrexate exert their activity through inhibition of DHFR [
39,
40]. Furthermore, the two antimalarial drugs cycloguanil and pyrimethamine, as well as a codified antimalarial compound (P218) were also found to inhibit DHFR enzyme [
41]. Unlike the other most popular antifolate agents, such as trimethoprim, cycloguanil, and pyrimethamine that exhibited weak inhibition of
Leishmania major DHFR, methotrexate inhibited this enzyme (
L. major DHFR) in a nanomolar range (IC
50: 5 nM). The whole cell phenotypic screening of methotrexate against
L. major revealed IC
50 value of 0.3 µM, confirming that this compound might have exerted antileishmanial activity through DHFR inhibition [
42,
43]. Recently, Dize et al. revealed the inhibitory potential of MMV675968 derivatives bearing the quinazoline scaffold on the
Trypanosoma brucei brucei DHFR together with their potent antitrypanosomal activity (IC
50 range: 45-60 nM) [
44]. On the other hand, potassium channels are among the proteins, which are found in the membrane of parasites, including
Plasmodium,
Leishmania and
Trypanosoma. A number of authors have described the importance of K+ channels in the survival of these parasites [
45,
46,
47,
48]. Several synthetic compounds, such as clofazimine derivatives that revealed antiprotozoal activity against
Leishmania and
Plasmodium species were found to inhibit K+ channels [
49]. Moreover, Waller et al. evaluated the antiplasmodial activity of several known K+ channel blockers, including quinidine, clotrimazole, haloperidol, charybdotoxin, bicuculline methiodide, tubocurarine chloride, trifluoperazine hydrochloride, and verruculogen, and the results showed inhibition of
Plasmodium falciparum 3D7 with IC
50 values stretching from 0.046 to 187.86 µM [
45].
In this study, the antiprotozoal activity of a series of DHFR inhibitors (compound X, Y, and Z) and the potassium channel blocker E4031 was evaluated against
Plasmodium,
Leishmania and
Trypanosoma species. As a result, compounds X, Y, and Z, and E4031 showed different degrees of antiprotozoal activity extending from poorly (IC
50S>10 µM) to highly active (IC
50S< 10 µM) [
50,
51]. Against
T. brucei, compounds X and Y exhibited IC
50 values of 6.49 and 0.81 µM, respectively, vs pentamidine (IC
50 value: 0.006 µM), whereas these compounds revealed IC
50 values of 0.0052 and 0.028 µM, when tested against
Plasmodium falciparum 3D7 with high selectivity in Raw cells (SI: 366.6 and 1179 for compounds X and Y, respectively), compared to the value obtained with artemisinin (IC
50: 0.03 µM) (
Table 2). In addition, compound Y showed IC
50 values of 12.47 and 4.28 µM (SI: >2.69; Vero and Raw cells), when tested against
L. donovani promastigotes and amastigotes, respectively, vs amphotericin B (IC
50 values: 20 and 247.81 µM, respectively). Furthermore, compound Z and E4031 showed IC
50 values > 10 µM when tested against
T. brucei, promastigotes of
L. donovani and
P. falciparum 3D7, and were predicted to be poorly active compounds. Although numerous authors have reported the anticancer and antifolate activities of compounds X, Y, and Z [
25,
27], the antiprotozoal activity of these compounds has not yet been unveiled. Compounds X, Y and Z are structurally similar compounds, in which one amino group of the pyrimidine moiety is replaced by an OH group in compound Z. This structural modification might have contributed for the decrease in the observed antiprotozoal activity of compound Z. Although there is reported evidence that the presence of OH groups tend to increase the antiprotozoal activity of bioactive compounds [
52], other reports highlight a different opinion about this matter [
53]. Compounds’ selection based on their positive drug-likeness scores increases the likelihood to identify potential drug candidates with favorable pharmacokinetics, efficacy, and safety. As realigning pharmacokinetic studies early in the discovery phase can assist in selecting an ideal drug candidate, active antiprotozoal compounds were subjected to
in silico screening using the pkCSM online tool. Regarding the
in silico examination of physicochemical properties, molecular weight (441.488, 442.472, 442.476 and 415.559) ClogP (lipophilicity; octanol-water partition coefficient) (1.741, 1.865, 1.221 and 3.289), rotatable bonds (10, 10, 10 and 8), hydrogen bond acceptors (HBA) (7, 7, 8 and 5), hydrogen bond donors (HBD) (5, 5, 6 and 1) and Topological Polar Surface Area (TPSA) (184.604, 184.058, 183.789, and 172.949) (
Table 3) of compounds X, Y, and Z and E4031 were respectively predicted. Notably, none of these physicochemical parameters were violated according to criteria by the Lipinski's rule of five, which states that a molecule has an increased chance of being directly bioavailable when it obeys to the conditions of having (i) no more than five hydrogen bond donors, (ii) ten hydrogen bond acceptors, (iii) a molecular weight under 500, and (iv) a LogP number under 5 [
54]. In fact, all tested compounds (X, Y, Z and E4031) qualify as orally active compounds as they abide by the "Lipinski's rule of five" criteria [lipophilicity (LogP) < 5, HBD ≤ 5, HBA ≤ 10, molecular weight (g/mol) ≤ 500 g/mol, and flexibility (rotatable bonds): < 10] [
55,
56]. Other rulesets for drug-likeness, including the Veber [flexibility (rotatable bonds): < 10; HBD ≤ 5; HBA ≤ 10; except for TPSA (140 Å2)], Ghose filter (-0.4<Log P<+ 5.6; HBD ≤ 5 and HBA ≤ 10; MW<500 g/mol), Muegge (HBD ≤ 5, except for E4031 and MW<500 g/mol) and Egan (-0.4<Log P<+5.6; MW<500 g/mol and HBA ≤ 10)] rules showed favorable drug-like characteristics [
55,
56]. Furthermore,
in silico tests on ADME (absorption, distribution, metabolism and excretion) revealed poor permeability of compounds X, Y and Z, and E4031 across the skin, blood brain barrier and central nervous system. Nevertheless,
in vitro and
in vivo pharmacokinetic studies need to be thoroughly investigated for the successful utilization of these compounds as scaffolds in antiprotozoal drug discovery.
In this study, the antiprotozoal activity of a series of DHFR inhibitors (compound X, Y, and Z) and the potassium channel blocker E4031 was evaluated against
Plasmodium,
Leishmania and
Trypanosoma species. As a result, compounds X, Y, and Z, and E4031 showed different degrees of antiprotozoal activity extending from poorly (IC
50S>10 µM) to highly active (IC
50S< 10 µM) [
50,
51]. Against
T. brucei, compounds X and Y exhibited IC
50 values of 6.49 and 0.81 µM, respectively, vs pentamidine (IC
50 value: 0.006 µM), whereas these compounds revealed IC
50 values of 0.0052 and 0.028 µM, when tested against
Plasmodium falciparum 3D7 with high selectivity in Raw cells (SI: 366.6 and 1179 for compounds X and Y, respectively), compared to the value obtained with artemisinin (IC
50: 0.03 µM) (
Table 2). In addition, compound Y showed IC
50 values of 12.47 and 4.28 µM (SI: >2.69; Vero and Raw cells), when tested against
L. donovani promastigotes and amastigotes, respectively, vs amphotericin B (IC
50 values: 20 and 247.81 µM, respectively). Furthermore, compound Z and E4031 showed IC
50 values > 10 µM when tested against
T. brucei, promastigotes of
L. donovani and
P. falciparum 3D7, and were predicted to be poorly active compounds. Although numerous authors have reported the anticancer and antifolate activities of compounds X, Y, and Z [
25,
27], the antiprotozoal activity of these compounds has not yet been unveiled. Compounds X, Y and Z are structurally similar compounds, in which one amino group of the pyrimidine moiety is replaced by an OH group in compound Z. This structural modification might have contributed for the decrease in the observed antiprotozoal activity of compound Z. Although there is reported evidence that the presence of OH groups tend to increase the antiprotozoal activity of bioactive compounds [
52], other reports highlight a different opinion about this matter [
53]. Compounds’ selection based on their positive drug-likeness scores increases the likelihood to identify potential drug candidates with favorable pharmacokinetics, efficacy, and safety. As realigning pharmacokinetic studies early in the discovery phase can assist in selecting an ideal drug candidate, active antiprotozoal compounds were subjected to
in silico screening using the pkCSM online tool. Regarding the
in silico examination of physicochemical properties, molecular weight (441.488, 442.472, 442.476 and 415.559) ClogP (lipophilicity; octanol-water partition coefficient) (1.741, 1.865, 1.221 and 3.289), rotatable bonds (10, 10, 10 and 8), hydrogen bond acceptors (HBA) (7, 7, 8 and 5), hydrogen bond donors (HBD) (5, 5, 6 and 1) and Topological Polar Surface Area (TPSA) (184.604, 184.058, 183.789, and 172.949) (
Table 3) of compounds X, Y, and Z and E4031 were respectively predicted. Notably, none of these physicochemical parameters were violated according to criteria by the Lipinski's rule of five, which states that a molecule has an increased chance of being directly bioavailable when it obeys to the conditions of having (i) no more than five hydrogen bond donors, (ii) ten hydrogen bond acceptors, (iii) a molecular weight under 500, and (iv) a LogP number under 5 [
54]. In fact, all tested compounds (X, Y, Z and E4031) qualify as orally active compounds as they abide by the "Lipinski's rule of five" criteria [lipophilicity (LogP) < 5, HBD ≤ 5, HBA ≤ 10, molecular weight (g/mol) ≤ 500 g/mol, and flexibility (rotatable bonds): < 10] [
55,
56]. Other rulesets for drug-likeness, including the Veber [flexibility (rotatable bonds): < 10; HBD ≤ 5; HBA ≤ 10; except for TPSA (140 Å2)], Ghose filter (-0.4<Log P<+ 5.6; HBD ≤ 5 and HBA ≤ 10; MW<500 g/mol), Muegge (HBD ≤ 5, except for E4031 and MW<500 g/mol) and Egan (-0.4<Log P<+5.6; MW<500 g/mol and HBA ≤ 10)] rules showed favorable drug-like characteristics [
55,
56]. Furthermore,
in silico tests on ADME (absorption, distribution, metabolism and excretion) revealed poor permeability of compounds X, Y and Z, and E4031 across the skin, blood brain barrier and central nervous system. Nevertheless,
in vitro and
in vivo pharmacokinetic studies need to be thoroughly investigated for the successful utilization of these compounds as scaffolds in antiprotozoal drug discovery.
Author Contributions
“Conceptualization, F.F.B.; methodology, D.D., M.B.T.T., C.A.N.N., K.H., and R.K.; software, D.D., M.B.T.T., C.A.N.N. and R.K.; validation, E.A.M.K., L.R.T.Y., B.P.K., P.V.T.F. and F.F.B; formal analysis, K.H., B.P.K., E.A.M.K., L.R.T.Y., and P.V.T.F.; investigation, K.H., D.D., M.B.T.T., C.A.N.N. and R.K.; resources, B.P.K., E.A.M.K., L.R.T.Y., P.V.T.F. and F.F.B; data curation, D.D., M.B.T.T., C.A.N.N. and R.K.; writing—original draft preparation, D.D., M.B.T.T., C.A.N.N., E.A.M.K., L.R.T.Y., K.H. and R.K.; writing—review and editing, E.A.M.K., L.R.T.Y., P.V.T.F., B.P.K. and F.F.B; visualization, E.A.M.K., L.R.T.Y., P.V.T.F. and B.P.K.; supervision, F.F.B.; project administration, F.F.B.; funding acquisition, F.F.B. All authors have read and agreed to the published version of the manuscript.”.