Next Article in Journal
Pathogens: Journal Section Readjustment
Next Article in Special Issue
Neuropathological Changes in Nakalanga Syndrome—A Case Report
Previous Article in Journal
Enlisting the Ixodes scapularis Embryonic ISE6 Cell Line to Investigate the Neuronal Basis of Tick—Pathogen Interactions
Previous Article in Special Issue
Effect of Ivermectin Treatment on the Frequency of Seizures in Persons with Epilepsy Infected with Onchocerca volvulus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Integrated Approach to Identify New Anti-Filarial Leads to Treat River Blindness, a Neglected Tropical Disease

1
Division of Infectious Diseases, Department of Medicine, Washington University School of Medicine, 4523 Clayton Ave., St. Louis, MO 63110, USA
2
Department of Pharmaceutical Chemistry, University of California San Francisco, 1700 4th Street, San Francisco, CA 94158, USA
3
ANDI Centre of Excellence for Onchocerciasis Drug Research, Biotechnology Unit, Faculty of Science, University of Buea, Buea CM-00237, Cameroon
4
Division of Neurocritical Care and Hospital Neurology, Department of Neurology, Mayo Clinic, 200 First Street SW, Rochester, MN 55905, USA
5
Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, CA 92093, USA
6
Calibr, a Division of The Scripps Research Institute, 11119 Torrey Pines Road, La Jolla, CA 92037, USA
7
Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, 660 S. Euclid Ave., St. Louis, MO 63110, USA
8
Lindsley F. Kimball Research Institute, New York City, NY 10065, USA
9
McDonnell Genome Institute, Washington University School of Medicine, 4444 Forest Park Ave., St. Louis, MO 63108, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 17 December 2020 / Revised: 5 January 2021 / Accepted: 11 January 2021 / Published: 14 January 2021
(This article belongs to the Special Issue Onchocerciasis and River Epilepsy)

Abstract

:
Filarial worms cause multiple debilitating diseases in millions of people worldwide, including river blindness. Currently available drugs reduce transmission by killing larvae (microfilariae), but there are no effective cures targeting the adult parasites (macrofilaricides) which survive and reproduce in the host for very long periods. To identify effective macrofilaricides, we carried out phenotypic screening of a library of 2121 approved drugs for clinical use against adult Brugia pahangi and prioritized the hits for further studies by integrating those results with a computational prioritization of drugs and associated targets. This resulted in the identification of 18 hits with anti-macrofilaricidal activity, of which two classes, azoles and aspartic protease inhibitors, were further expanded upon. Follow up screening against Onchocerca spp. (adult Onchocerca ochengi and pre-adult O. volvulus) confirmed activity for 13 drugs (the majority having IC50 < 10 μM), and a counter screen of a subset against L. loa microfilariae showed the potential to identify selective drugs that prevent adverse events when co-infected individuals are treated. Stage specific activity was also observed. Many of these drugs are amenable to structural optimization, and also have known canonical targets, making them promising candidates for further optimization that can lead to identifying and characterizing novel anti-macrofilarial drugs.

1. Introduction

River blindness (onchocerciasis) and lymphatic filariasis (LF) are two major neglected tropical diseases (NTD) caused by parasitic nematodes that, together, affect millions of people worldwide in mostly poor, developing countries [1]. Approximately 20 million people are infected with Onchocerca volvulus; 14.6 million of the infected people have skin disease, 1.2 million people are visually impaired, and 270,000 are blind by river blindness, a chronic disease caused by the first larval stage, microfilariae (mf). Mf are released from female worms residing in subcutaneous tissues and migrate throughout the skin causing severe itchiness, as well as inflammatory responses in the skin or eyes. When the inflammatory responses occurs in the eye they may eventually lead to impaired vision and ultimately to blindness. LF (elephantiasis) is caused by damage to lymphatic tissues by adult Wuchereria bancrofti, Brugia malayi and B. timori worms, and is characterized by pain and severe lymphedema, often involving the extremities leading to great economic losses as well as social stigma [1].
To date, there are no vaccines to prevent these diseases, and no drugs that directly kill the adult stages (macrofilaricidal drugs) [2,3,4] and can be used in mass drug administration (MDA). A promising triple drug regimen for LF that has some macrofilaricidal effects is currently being evaluated against onchocerciasis [5,6,7]. International control programs attempt to interrupt transmission of infection with annual or biannual MDA using microfilaricidal drugs (ivermectin since 1989 and more recently also moxidectin [8] for onchocerciasis; albendazole and ivermectin or diethylcarbamazine for LF) that kill mf over the lifetime of the adult worms (10–14 years for O. volvulus, 6–8 years for Wuchereria and Brugia spp.) [2,3,4,9,10,11,12,13,14,15,16,17]. These drugs target different critical processes in the mf or make them more susceptible to immune system, e.g., targeting microtubule polymerization (albendazole), glutamate-gated chloride channels and other transporters (ivermectin and moxidectin), and sensitizing microfilariae to phagocytosis by host immune cells (diethylcarbamazine).
Given the longevity and high fecundity of these worms and the current lack of macrofilaricidal drugs, it is unlikely that the WHO goal of eliminating LF and onchocerciasis by 2030 [18,19] will be met using only microfilaricidal drugs [11,20,21]. According to the 2013 Global Burden of Disease Study, only a 31% reduction in onchocerciasis was achieved since MDA with ivermectin began in the 1990s [22]. Indeed, the African Programme for Onchocerciasis Control estimated that elimination of onchocerciasis would require some 1.15 billion treatments with ivermectin (IVM), with MDA efforts continuing until 2045 [11,23,24,25]. It is believed that moxidectin, a potent microfilaricide approved by the FDA to treat human onchocerciasis in 2018 [26], could substantially reduce the time it will take to eliminate onchocerciasis [8], as moxidectin has a superior clinical performance compared to ivermectin [27,28]. This would be especially helpful in places where resource limitation prevents a biannual ivermectin strategy, since biannual ivermectin and annual moxidectin treatment has been shown to achieve similar reductions in program duration [29].
Moreover, MDA with IVM is also confounded in Africa by the fact that it cannot be distributed in areas co-endemic for Loa loa (another filarial nematode), due to the risk of severe adverse events, especially toxic encephalopathy when infected individuals have high loads of L. loa mf [23,24,30,31]. Presently, treatment in these areas requires a test-and-treat approach, which is more resource-intensive and may result in incomplete treatment for O. volvulus [32]. The prohibitive use of IVM for 12 million people, in 11 Loa-affected central African countries, impedes elimination efforts and creates a reservoir of Onchocerca infections which can re-infect neighboring communities [30]. In addition, LF and onchocerciasis elimination programs in sub-Saharan Africa do not implement MDA in hypoendemic areas (low prevalence of infections), also leading to concerns of putative spreading of reinfections [11] in areas that might have controlled transmission. Furthermore, the potential emergence of IVM-resistant O. volvulus limits the long-term effectiveness of present MDA with IVM [4,33,34,35,36], and in time may undermine gains achieved by the MDA programs. The restrictions on MDA in children adds to the complexity of elimination efforts. IVM is not administered to children under 5 years of age, and the only proven indirect macrofilaricidal drug, doxycycline cannot be given to children under 9 years because of drug contraindications. Children thus remain vulnerable and serve as reservoirs of transmission [31]. Moreover, doxycycline requires long treatment periods of 4–6 weeks, which is not feasible for MDA programs [37].
Thus, it is still critical to identify and develop novel, effective and safe macrofilaricidal drugs for use in integrated anti-filarial MDA programs. A few approaches have been studied, including targeting respiratory enzymes and using antibiotics (e.g., doxycycline) effective against Wolbachia – an endosymbiont essential to survival of many filarial worms [38]. Ultimately, macrofilaricidal drugs will also have the potential to shorten the time to successfully eliminate onchocerciasis. Unfortunately, there is no animal model to facilitate harvesting of adult O. volvulus worms, and even the pre-adult developmental stage (OvL5) which can be used in motility and viability in vitro assays [39] is prohibitively expensive for large scale screenings. Therefore, in this study, we decided on a new stepwise screening approach, with initial large-scale screenings against adult B. pahangi worms, followed by secondary screenings of prioritized hits with adult O. ochengi worms (a closely related and readily accessible surrogate of O. volvulus and a clinical model of human onchocerciasis), and then with pre-adult L5 worms of O. volvulus.
In order to facilitate novel macrofilaricidal drug discovery, we also undertook an integrated multidisciplinary study that leveraged our recent progress in the field of nematode genomics [40] and performed a systematic search for genes essential for the survival of filarial nematodes at a level not previously possible. Computational prioritization interfaced with experimental identification of repurposed drugs that are active against filarial nematodes resulted in a short list of prioritized targets and drug pairs. This was accomplished by first undertaking a target class repurposing approach and testing a library of compounds approved for clinical use (n = 2121, purchased by the Small Molecule Discovery Center, University of California San Francisco, CA) for their activity in vitro on adult B. pahangi female worms. Because one of the major challenges faced by investigators working in the field of discovery of drugs with activity against adult worms is the inability to screen large compound libraries and to follow up on the many new leads, we next intersected the actives with the omics-driven, computationally identified and prioritized target:drug pairs (n = 4442), which resulted in a short list of 25 prioritized drug candidates. We subsequently screened a subset of these prioritized drugs with adult worms the cattle filarial nematode O. ochengi. Finally, to demonstrate proof-of-concept, i.e., our ability to expand the chemical space and obtain supportive evidence for the identified potential target(s), we followed up on actives associated with two very different parasite protein target classes, thus establishing a solid platform for informed rational medicinal chemistry approaches in future studies. Those studies will lead to target validation and drug optimization, followed by in vivo activity confirmation using human equivalent dosages predicted based on jird pharmacokinetics data and human safety and clinical information. This successful approach can now also be used for identifying novel drugs and corresponding targets essential for the survival of other parasites, leveraging the extensive omics datasets for the human host.

2. Results and Discussion

2.1. Experimental Screening to Identify Active Macrofilaricidal Compounds Against Adult B. pahangi

A library of 2121 drugs approved for clinical use (obtained from the UCSF Small Molecule Discovery Center) was screened in a phenotypic assay against adult female B. pahangi at 10 µM (Table S1) using the WormAssay software (University of California San Francisco, San Francisco, CA, USA) and dark-field plate imaging system [41]. A total of 124 drugs showed ≥50% inhibition of motility on day 3 (Figure 1A; Table S2) and a subset of 62 among these active drugs (based on motility inhibition and commercial availability) was screened in a dose response assay to determine their IC50 values (Figure 1B). The drugs showing higher levels of inhibition of motility tended to also have lower IC50 values (Figure 1C, correlation coefficient = −0.82 for the top hits).

2.2. Integrating Active Drugs with Computational Prioritization

The 2121 screened drugs were associated with 1961 potential targets in the ChEMBL database [42] (Figure 2A). Given the ‘many-to-many’ relationships between drugs and targets, we analyzed the 1961 potential ChEMBL targets to identify a subset with significant association with the 124 active drugs (as compared to the non-active screened drugs). This identified 31 targets showing significantly enriched association (Fisher’s exact test FDR-adjusted p ≤ 0.05) with the actives. These 31 targets showed functional association with multiple critical cellular processes (Figure 2B). These include receptor based signaling (e.g., adenosine, opioid, and chemokine receptors) [43,44,45], ubiquitin pathway (ubiquitin specific peptidase, ubiquitin conjugate enzyme) [46], transcription factors (e.g., p53, HIF1) [47,48] and oxidoreductase activity (e.g., ALOX12 and ALOX 15) [49]. These 31 targets were associated with 69 out of the 124 active drugs (Figure 2A).
Further prioritization among these 69 compounds was done by considering a set of potential anthelmintic targets previously reported by us, as a part of International Helminth Genomes Consortium [40]. Briefly, an omics-driven pipeline to identify the most promising targets from parasitic worms was implemented using 528,469 proteins from 33 parasitic worms (including all major filarial species of human importance), yielding identification of 3994 ChEMBL homologous targets (Figure 2A). Further reduction in redundancy, via recognition of members of the same orthologous protein families, reduced the 3994 targets to 1925 ChEMBL targets. Identifying the commercially available compounds associated with these top targets having desirable properties (e.g., Quantitative Estimate of Druglikeness (QED) score, FDA approved) and selecting a structurally diverse set among these resulted in 4442 drugs and drug-like compounds (Figure 2A).
Intersecting the two parallel analyses from both computational and experimental prioritizations, resulted in a list of 25 overlapping hit compounds (Figure 2A). This set of drugs covered wide structural space (median pairwise Tanimoto similarity = 0.096) and is comprised of different compound classes, including benzyl ether azoles (antifungals targeting membrane sterol synthesis [50]), dihydropyridines (calcium channel blockers, for hypertension [51]), quinoline derivatives (antimalarials and antimicrobials [52]), and phenylpropanes (vasoactive drugs [53]). This high priority set of hits was then manually curated and further modified (for details see Methods), resulting in 18 prioritized hits that were experimentally effective in inhibiting motility (>50% inhibition) of adult B. pahangi at 10 μM (Table 1), with most having known putative targets [54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79].
We noticed that many of the prioritized hits (6/18) were antifungal azoles. Henceforth we treated these azoles as belonging to a single group (Group A), separate from the remaining 12 compounds of diverse structures (Group B) (Figure 3). We further categorized the 12 Group B drugs based on the % motility inhibition into high (>80%), moderate (65–80%) and low motility inhibition (50–65%) group. Then, two drugs per each group were assayed to determine IC50. As expected, the lowest IC50 values were observed for the high % motility inhibition group, 2.5 and 3.7 μM (suloctidil and pimozide, respectively) and higher IC50 values for the other 2 groups (clemizole 6.5 μM, proroxan 8.5 μM, tafenoquine 14.1 μM and primaquine 5.1 μM). Compounds in Group B were classified into high priority because they are amenable to modifications, or low priority due to central nervous system (CNS) activity and due to unknown mode of action (Figure 3).

2.3. Screening of Prioritized Hits in Other Filarial Species

Our ultimate goal is to identify drugs with macrofilaricidal activity against the human parasite O. volvulus. One of the major challenges faced by investigators working in the field of drug discovery and neglected parasitic diseases is the ability to screen large compound libraries to identify new leads for the treatment of these highly prevalent, debilitating diseases. In this study we employed a multi-faceted screening funnel that utilizes several species of filarial nematodes and their relevant life history stages, in conjunction with a chemogenomic approach to identify new leads. Since the parasite that causes onchocerciasis only infects humans, several closely related surrogate species are used in the preclinical studies, including a primary screen with B. pahangi (laboratory species) and a secondary screen with the surrogate species O. ochengi. O. ochengi infects cow and has been used for various in vivo intervention studies [80]. As O. ochengi adult worms are readily accessible from West African infected cattle, we used this screen as a useful secondary screen that can establish proof-of-principle for potential therapeutics for human onchocerciasis. A tertiary assay with the target species O. volvulus (larval and pre-adult worms maintained in vitro) and a counter screen with Loa loa (from human volunteers) add to the robustness of the screening program.
Following our screening funnel, the Group A drugs and the top 3 Group B drugs (prioritized based on IC50 < 5.1 µM with adult female B. pahangi) were screened against adult female (viability) and male (motility) O. ochengi worms in vitro. None of the Group A compounds strongly affect (>70%) the viability of female O. ochengi (7-day assay), but 2 of the 3 Group B drugs (suloctidil and primaquine) were highly effective on female O. ochengi (>87% inhibition of viability); the IC50 values for these two Group B drugs were 4.1 and 1.3 µM, respectively (Table 2). In comparison, 3 out of the 6 Group A compounds (miconazole, sulconazole, and clotrimazole) were highly effective (> 99% motility inhibition) on male O. ochengi (5-day assay), and all 3 Group B compounds (suloctidil, pimozide and primaquine) also showed potent activity against male O. ochengi worms (>96% motility inhibition; Table 2), with IC50 values of 5.5 µM and 0.4 µM for suloctidil and primaquine, respectively. It is possible that this stage-specific difference is partly a result of differences in assay conditions between male and female worms, i.e., the treated female worms reside in the worm masses, while the males exit the worm masses and are exposed to the media, hence are possibly more directly exposed to the drugs.
The screening results showed that both suloctidil and primaquine were potent across the B. pahangi and O. ochengi species, including both female and male O. ochengi. Thus, the two compounds were also screened against an advanced developmental stage of pre-adult O. volvulus worms (L5 larval stage [39]), and counter-screened against L. loa microfilariae (Table 2). Both the compounds at 10 µM were effective on OvL5, showing 100% inhibition of motility already on day 21. The counter-screen with L. loa is crucial in order to avoid potential severe adverse effects that are known to result from drugs that target mf, such as ivermectin, when used by individuals co-infected with a high load of Loa loa mf [31,81]. It is therefore an important strategy to ensure that the potential macrofilaricidal drugs are inactive against L. Loa mf. While both the drugs retained high activity against O. volvulus L5, primaquine proved to be more inactive against L. loa mf (no detectable inhibition on day 1 and inhibition of motility slowly increased to 50% on day 5 at 10 µM, IC50 = 18.5 µM) as compared to suloctidil (100% inhibition already on day 1 at 10 µM, IC50 = 4.1 µM) (Table 2). A lack of commercial macrofilaricides makes it difficult to put these results in translational perspective for human use and further investigations to determine the pharmacokinetics and pharmacodynamics are also needed to assess the potency of these drugs in vivo. However, two drugs that are on the DNDi (The Drugs for Neglected Diseases initiative) portfolio as macrofilaricidal candidates, oxfendazole and emodepside, [82,83] have similar IC50 values against filarial species. The anthelmintic drug oxfendazole, which is currently used in veterinary medicine, was reported to have IC50 of 7.6 µM (day 14 in culture) and 28.6 µM (day 19 in culture) at inhibiting motility of O. volvulus L5 larvae [84] and emodepside, another veterinary anthelmintic, has been reported to have a submicromolar motility inhibition IC50 (0.8 µM) in adult female B. malayi [85].
In summary, our results provide a proof-of-concept that targeting filarial-conserved genes essential for adult worm survival is predictive of anti-filaricidal activity across species, and that differential selectivity against Brugia and Onchocerca compared to Loa can be detected, indicating that selective targeting can be achieved.

2.4. Expanding the List of Potential Anti-Macrofilarial Drugs by Focusing on Two Specific Classes

To further expand the list of hit drugs that could have macrofilaricidal activity, two new alternative approaches were undertaken: expanding the most successful compound class (azoles) and exploring a promising target class (aspartic proteases).
Azoles (Group A hit drugs), the most potent compound class in the B. pahangi adult assays, was expanded by screening additional compounds from the same class. Azoles had potent macrofilaricidal effects and have shared activities, i.e., they inhibit sterol demethylases (specifically 14-alpha demethylase) and are known antifungal agents [86] that act by inhibiting ergosterol synthesis which is needed for fungal membrane permeability. However, it is likely that these compounds have a different mechanism of action in Nematoda since they appear to lack an ortholog of this 14-alpha demethylase. Azoles have also been reported to inhibit other crucial genes like P-gp, multiple CYP proteins, certain ion channels and receptors, thromboxane synthase (an inflammation mediator), FYN (a Src family tyrosine kinase oncogene), heme oxygenase (heme metabolism; maintaining homeostasis under oxidative and other stresses), Indoleamine 2,3-dioxygenase (heme containing immunomodulator), etc., [87,88]. They have also been identified in our previous work as potential broad-spectrum anthelmintics [89] based on their known activity against malate dehydrogenase (MDH) [90,91]. One concern is that some of these drugs may exhibit poor water solubility and oral bioavailability and may also have unacceptable side-effects on oral administration. However, this is no longer an issue with second and third generation azoles (itraconazole, posaconazole, etc.) [92]. We decided to focus on this group and expand this set of compounds by including all active azoles from our primary screening. A total of 9 azole hits were identified (out of 49 screened; Table S3), all with % motility inhibition of >80%: 6 from our integrated hits, 1 from the primary screen (tioconazole), and 3 from the ReFRAME library (a library of repurposed drugs) [93] (fenticonazole, isoconazole and posaconazole) (Table 2). For 7 of the 10 azoles, IC50 values for adult female B. pahangi ranged from 1.1 to 5.5 µM. When tested in adult female O. ochengi, 1 of these 10 azoles inhibited female worm viability by over 50% (isoconazole), compared to 7 azoles inhibiting adult male worm motility greater than 75%. For female B. pahangi, the IC50 values for fenticonazole, isoconazole, and posaconazole were 1.2, 3.8 and 0.1 µM, respectively; for female O. ochengi, these compounds lacked high potency, with IC50 values of 95, 31, and >100 µM, respectively, whereas the potency in male O. ochengi was better with IC50 values of 10.6, 0.6, and 6.1 µM respectively (Table 2). In O. volvulus L5, the IC50 for the fenticonazole and isoconazole were 4.7 and 2.4 µM, respectively.
In addition to this class of drugs, we also focused on aspartic proteases (APs) as potential targets, using known aspartic protease inhibitors (APIs) [94,95]. This was based on our observation that one of our hit drugs, suloctidil, that showed consistently high potency across species and worm gender contains a hydroxyethylamine (HEA) functional group, and it has an overall substructure similar to ritonavir, a known API that is also among our top 18 candidate drugs (Figure 3B). Ritonavir is an inexpensive HEA group containing inhibitor of the aspartyl HIV-protease (HIV-PR) and is on the WHO list of essential medicines for its antiretroviral activity against HIV and AIDS [96]. Interestingly, levonordefrin, another hit drug (Table 1), also contains the HEA substructure. It is plausible that these three compounds might be targeting one or more of the Brugia aspartic proteases (Figure 4). APs are attractive targets due to their crucial functions in helminths [97], and they have been well studied in various parasitic nematode species [98,99,100], including filarial worms [101].
Protease inhibitors in general, and APIs in particular, have previously been developed as therapies for various diseases, e.g., hypertension, cancer, malaria, Alzheimer’s, and AIDS [102], with 5–10% of all drugs under development targeting proteases [103]. For APIs, this includes drugs designed using computational approaches, to specifically mimic the transition state, often via a HEA moiety or ‘warhead’ forming key H-bonding electrostatic contacts at the active site in the ‘oxyanion hole’, analogous to the substrate tetrahedral intermediate upon addition of water to the amide bond. APIs are especially amenable for in silico de novo design because APs form a tetrahedral substrate transition state with no covalent bonds (in contrast to serine or cysteine proteases) between the enzyme and the substrate and that can be mimicked by rationally designed molecules. Examples of such drugs include saquinavir and ritonavir, which were identified as potent inhibitors of HIV-PR. Interestingly, they have also been shown to be effective against other parasitic species, presumably targeting their homologous APs. For example, ritonavir is known to inhibit pathogenic Candida species [104], and was shown to dock in the Trypanosoma cruzi AP active site [105]. In some species, including HIV, the active site is reconstituted by homodimerization of the protein and requires two critical DTG or DSG active motifs, one from each monomer. We identified 11 APs in B. malayi, 8 in B. pahangi and 17 in O. volvulus using sequence-similarity against the MEROPS peptidase database [106] (BLASTP E-value <10−4). Among these, 4 in B. malayi and 3 in O. volvulus have conserved DTG motifs, suggesting that these are potential targets of ritonavir-like APIs (Figure 4). Unlike HIV, but similar to many other organisms (e.g., humans), the nematode aspartic proteases have an active site that appears to be part of the same protein and does not need dimerization. Multiple AP gene copies might result in functional redundancy (notably B. pahangi, the species used for primary screening, also has 3 orthologs—BPAG_0001244901, BPAG_0001278201 and BPAG_0000957001—for these 4 APs in B. malayi). However, the different stage-specific expression profiles (Figure 4) suggest that these may be crucial genes in different life cycle stages, including the developmental stages inhabiting the human host, e.g., mf, L3, L4, and adults. Specifically, Bm8660 and OVOC11635, which are the proteins most closely related to human Cathepsin D (CD). OVOC11635 shares 94% identity over the full length with previously characterized O. volvulus AP (cathepsin D [101]), and shows significant overexpression in L3 and adult male, respectively. CD is a more broadly active AP than HIV PR, and has multiple crucial physiological functions in vertebrates, e.g., metabolic degradation of intracellular proteins, activation of enzyme precursors, and regulation of apoptosis [107]. In Nematoda, CDs can function in invasion of host tissue, modulating host immune response and digestion of host proteins. In C. elegans, CD is the most significant enzyme for substrate macromolecular digestion in vitro [108].
In addition to the APIs included in our library of approved drugs for clinical use, 9 other drugs that inhibit HIV-PR are FDA-approved [109]. We were able to purchase 5 of these 9 APIs, and also purchased aliskiren (an FDA approved renin inhibitor) and pepstatin A (a highly potent general AP inhibitor not approved for clinical usage). The 7 APIs were screened against B. pahangi adult females, adult O. ochengi male and female worms, and O. volvulus L4 larvae using motility assays, to test their possible stage-specific efficacy (since some of the filarial aspartic proteases have stage-specific expression [110] (Figure 4B). Two of the HIV-PR inhibitors, nelfinavir and lopinavir (30 µM), inhibited the motility of female B. pahangi by >98% (Table 3), male O. ochengi by >88%, and O. volvulus L4 by >56%. In comparison, nelfinavir was highly effective against female O. ochengi viability (100%), while lopinavir was not effective at all against O. ochengi female worms, which might indicate a differential stage-specific activity. Pepstatin A showed species-specific activity, i.e., was only active against O. ochengi, with preferential activity against male (100% motility inhibition) compared to female (50% viability inhibition) worms. Higher transcriptional expression of Cathepsin D in Onchocerca male compared to female worms (Figure 4B) may partially explain this observation. In addition, cuticle permeability may be also determining the differential activity since suloctidil, which is one of our hits and was one of the factors in our decision to explore APs as targets of interest, along with 2 of our API hits (nelfinavir and lopinavir) have significantly higher cLogP values compared to the 6 inactive (6.0 vs. 3.8, P = 8.9 × 10−4, two-tailed T-test).
Ritonavir and the 7 APIs mentioned above were also assayed using a molting assay, observing the extent to which the of O. volvulus L3 to L4 larval molting, when incubated with the compounds is affected. All the APIs, except for ritonavir and atazanavir, showed moderate percentage molting inhibition on day 6 (≥40%) with 10 µM dosage. This means that some of the APIs, including darunavir, aliskiren and amprenavir, demonstrate a stage specific activity. This could potentially be related to developmental stage-based expression variation of some APs, the putative targets of these compounds (Figure 4B), with mid-to-high level expression (i.e., positive stagewise Z-scores) in the L3 larval stage with especially high expression seen in B. malayi L3. However, nelfinavir and lopinavir show activity in almost all the stages and species assayed (lopinavir was not effective at 30 µM on O. ochengi female worms) (Table 3). Interestingly, when we compared the IC50 values, the vasodilator suloctidil was more potent against adult O. ochengi (IC50 4.1 µM for female and 5.5 µM for male) than ritonavir (IC50 27.7 µM and 15.9 µM, in male and female worms, respectively), nelfinavir (IC50 19 µM in both male and female worms) and lopinavir (IC50 27 µM in male; not determined in female worms, >100 µM). Nonetheless, these IC50 values for these APIs are within the range of oxfendazole (veterinary anthelmintic in clinical trial for human onchocerciasis [82,83,84]) making them a great starting target class for optimization. It is notable, that even though these APIs consistently show only moderate activity in female O. ochengi worms, these worms are quite susceptible to other drugs, including our positive control for these experiments, Auranofin, a potent thioredoxin reductase inhibitor [111], with 100% loss of viability (10 µM) in all such assays. Nevertheless, these preliminary data spanning multiple filarial species and developmental stages suggest that APIs are effective in inhibiting the natural worm enzyme(s) and are excellent leads for expanding the chemical space as chemoprophylactic drugs as well as for macrofilaricidal drugs.
In summary, our primary and secondary phenotypic in vitro screens have identified 24 drugs with potential macrofilaricidal activity, of which 13 were hits in at least 2 filarial species, and 7 showed gender-specific activity. Significantly, one of the tested hit drugs was also active against the Brugia and Onchocerca adult stages but more inactive against the L. loa mf, which is required for drug use in L. loa co-endemic regions. While we followed up on two compound classes, the azoles and the APIs, other hits are equally intriguing. For example, pimozide, proroxan, and clemizole, which have high anthelminthic potency against B. pahangi (IC50 of 2.5, 8.5, and 6.5 µM, respectively) are structurally similar CNS-active drugs that target different GPCRs. Pimozide and clemizole are benzimidazoles known to be dopamine [56,57] and histamine antagonists [60], respectively, but also have other related targets [58,61], and proroxan is a known adrenergic blocker. Overall, the hit drugs identified in this study present a broad range of structural space that can be explored further (Figure 5).

3. Conclusions

Overall, our integrated computational and experimental in vitro screening approach has proven to be highly successful for identifying new early leads and their putative parasite targets, while also confirming their essentiality for filarial adult worm fitness in phenotypic in vitro assays. Phenotypic in vitro drug screening based on motility of adult filarial worms has been extensively used to assess the potency of various chemical classes on filarial worm viability [84,85,112,113,114,115]. This study, however, not only undertook multiple screening assays in different species, sex and developmental stages of filarial worms adding rigor to the identification of potential drugs, but also using the computational analyses has subsequently allowed us to prioritize and expand the set of promising target and drug paired compounds (azoles and APIs) to show broad pan-filarial anti-macrofilarial activity. With good anti-macrofilarial potency, high potential for structural optimization, and known canonical therapeutic targets these hits offer a promising starting point for identifying lead drug scaffolds and targets for optimizing and characterizing novel anti-filarial drugs. Once future studies validate these filarial parasite targets, these can be further used for selective drug optimization and followed by in vivo activity confirmation using human equivalent dosages predicted based on animal PK data and human safety and clinical information. We posit that this successful approach can now be also used for identifying novel drugs and corresponding targets essential for the survival of other parasites, leveraging the extensive omics datasets for the human host.

4. Materials and Methods

4.1. Ethics Statement

L. loa mf were collected from adult patient donors, aged 21 or older, with 2000+ mf per mL of blood living in the Edea Health District of the Littoral Region of Cameroon. Each patient provided written and signed informed consent, and ethical and administrative approval was obtained from the Cameroon National Ethics Committee (N°2013/11/371/L/CNERSH/SP) and the Cameroon Ministry of Health.
Adult female B. pahangi were collected from jirds (Meriones unguiculatus) at TRS labs Inc., Athens, Georgia and from the Univ. of Missouri–Columbia Institutional Animal Care Facility (IACUC approval #8623 and #9537).

4.2. Experimental Screening of a Library of Drugs Approved for Clinical Use Against B. pahangi

Individual adult female B. pahangi were placed in media (RPMI-1640 with 25 mM HEPES, 2.0 g/L NaHCO3, 5% heat inactivated fetal bovine serum (FBS), and 1X Antibiotic and Antimycotic solution) in 24-well plates. A Biomek FxP liquid handler was used to remove excess media so that each well contained one female worm in 500 µL media. Compounds dissolved in DMSO were added to each well at a concentration of 10 µM, with 4 replicates per compound and 1% DMSO was used as a vehicle control and 10 µM Auranofin as a positive control. Cultures were maintained in a 37 °C incubator over the course of the assay, and motility measurements were recorded using the WormAssay software [116], dark-field plate recording apparatus and 1080p digital camcorder on days 0–3. Compounds that caused 75% or greater inhibition of motility by day 3 were considered hits and were then tested in IC50 assays, using a 6-point serial dilution ranging from 30 µM to 0.1 µM. IC50 values were calculated using GraphPad Prism (GraphPad Software, San Diego, CA), and this data was used to prioritize compounds for further in vitro screening with Onchocerca spp. and L. loa microfilaria. Selected compounds (i.e., the aspartyl protease inhibitors) were screened at 30 µM and an additional motility reading was recorded on day 6.

4.3. Experimental Screening of Prioritized Drugs in O. ochengi Adults

Nodules containing adult O. ochengi were collected from the umbilical skin of infected cows from abattoirs in Douala, Cameroon. Worm masses containing one adult female and 0–7 adult males were excised from the nodules and incubated in 4 mL of complete Onchocerca culture medium (RPMI-1640, 5% newborn calf serum, 200 units/mL penicillin, 200 μg/mL streptomycin and 2.5 μg/mL amphotericin B (Sigma–Aldrich)) in 12-well culture plates. Worm masses were incubated at 37 °C, 5% CO2 overnight to allow adult males to egress from the mass into the surrounding media. Compounds were then added to wells at a final concentration of 10 µM for initial screening; 1% DMSO was used as a negative control. Testing was conducted as previously described [117] in quadruplicate in each assay, and assays were conducted twice on separate days to ensure reproducibility. On day 5, after compound addition, inhibition of male motility was scored according to the following criteria: 100% (complete inhibition of motility), 90% (only head or tail of worm moving or vibrating), 75% (worm very sluggish), 50% (worm sluggish), 25% (little change in motility), to 0% (no observable reduction in motility). On day 7, female viability was assessed using an MTT/formazan assay in which each worm mass was washed in PBS and then transferred to a 48-well plate with 500 μL of 0.5 mg/mL MTT (Sigma–Aldrich) in incomplete culture medium, and then incubated in the dark at 37o C protected from light for 30 minutes. After incubation, female worm masses were visually examined and scored to determine inhibition of formazan production, with higher percent inhibition of formazan indicating reduced worm viability as previously described [117]. To calculate the IC50, quadruplicate worm masses were incubated using a 7-point serial dilution ranging from 30 µM to 0.03 µM and assays were conducted as described above. GraphPad Prism was used to calculate IC50 values.

4.4. Counter Screening Compounds with L. loa Microfilariae in Vitro

Ten mL of venous blood in an EDTA tube provided by each consented donor was layered on a stepwise 46 and 43% Percoll (GE Healthcare, Piscataway, NJ, USA) gradient, and then centrifuged at 400 rcf for 20 minutes. L. loa mf were collected from the 43% layer and washed 3× in complete culture media. Loa mf were then cultured in 96-well plates (10–15 per well) containing a confluent layer of monkey kidney epithelial cells as previously described by Cho-Ngwa et al. (2010) [117]. Compounds were tested at 10 µM in duplicate wells and 10 μg/mL ivermectin was used as a positive control. Motility was scored daily for 5 days after compound addition using the same scale as the motility of O. ochengi adult males.

4.5. O. volvulus L5 Motility and Viability Assay

High priority compounds, i.e., those that had previously shown activity against Onchocerca spp. and Brugia spp., were tested with O. volvulus L5 in vitro. L3 larvae were cultured to L5 larvae as described by Voronin et al., 2019 [39]. Once reaching the pre-adult and L5 stage after at least 70 days in culture, L5 were placed in transwell inserts in a 24-well plate and treated with drugs, as previously described [39]. Briefly, each transwell insert was added to a well containing 500 µL OvL4-CMS media with a human umbilical vein endothelial cell (HUVEC) monolayer, and 8–10 L5 were cultured in each transwell insert. Compounds were tested in triplicate, with the concentration selected based on the potency found in B. pahangi in vitro assays. OvL5 were treated for 14 days and then in normal media as described [39]. Inhibition of motility was scored on a 0–100% scale continuously (every 2–3 days) throughout the 28-day assay. Beyond 28 days the control worms also start showing loss of motility, thus the 28 day was selected as an end point. For some (suloctidil and primaquine) drugs, the effectiveness was already highly evident by day 21, obviating the need to continue the assay beyond that point. Viability was measured at the end of the motility assay (day 29) by incubating the L5s with 0.1% MTT in PBS at 37 °C for 1 h, then visually examining the larvae for formazan production as indicated by blue coloration of the worm. If less than 50% of the worm was stained blue the worm was considered dead, and if greater than 50% of the worm was stained blue the worm was considered alive.

4.6. O. volvulus L4 Motility Assay

L4 (day 9 after molting) were cultured with 30 µM of the 7 APIs, and the motility of the L4 was monitored over 6 days. The anthelmintic flubendazole at 10 µM was used as a positive control.

4.7. Computational Identification and Prioritization of Targets and Inhibitors

Our starting step was to identify genes essential for filarial nematode survival by bioinformatically parsing the recently published database of genes essential for helminth survival [40]. That database of 81 helminth genomes was published by the International Helminth Genomes Consortium, and included a systematic search for genes essential for helminth survival [40]. In short, a target score was assigned to each of these worm proteins to quantify their likely quality as targets based on prior knowledge about their advantageous or disadvantageous properties: (i) the quality (E-value and target coverage) of the BLASTP match between the helminth and ChEMBL target proteins; (ii) whether the ChEMBL target had a close human BLAST match (since targeting a protein that lacks a human homolog is less likely to cause undesirable side-effects); (iii) whether the helminth gene had a C. elegans or D. melanogaster homolog with a severe phenotype (for example, a lethal or sterile phenotype; (iv) whether it was expressed in key life cycle stages (for example, adult); (v) whether it had homologs in most members of a major helminth clade (for example, in most nematodes); whether it lacked within-species paralogs; whether it belonged to a Compara family with a highly conserved alignment; (vi) whether the matching ChEMBL protein had a structure in the PDBe; and (vii) whether it was from a non-chordate animal.
The compounds linked to the helminth targets were ranked based on: (i) QED score [118] (a score that quantifies a compound’s drug-likeness by integrating some relevant molecular properties together and correlates with its likelihood to become a hit); (ii) at least phase III approval of the compound as a drug (since it is expected to be quicker and more efficient and economical to develop a drug which has already passed early phase clinical trials); (iii) whether the compound (if an approved drug) could be administered orally or topically; and (iv) lack of serious side effects or predicted toxicology targets. A total of 4442 compounds were associated with the top 15% (289) of the potential targets.
To identify potential targets of the 124 active compounds their association with putative targets in the ChEMBL database was analyzed as follows: The 2121 screened compounds included 1705 compounds that could be identified in the ChEMBL database using various forms or synonyms of the common names for these compounds. These 1705 ChEMBL compounds were associated with 1961 potential ChEMBL targets (i.e., had a ChEMBL assay recorded with pchembl value >5, corresponding to 10 µM target-compound half-maximal response metric—e.g., IC50, EC50, and Kd). Given the many-to-many relationships between compounds and targets, we analyzed the 1961 potential ChEMBL targets to identify a subset with significant association with the 124 active compounds. Out of the 1961 targets, 31 showed significantly enriched association (Fisher’s exact test FDR-adjusted p ≤ 0.05) with the active ChEMBL compounds. These 31 targets were associated (i.e., had ChEMBL records with pchembl >5) with 69 out of the 124 active compounds.
An intersection of these 69 drugs with anti-macrofilaricidal activity and the 4442 computational prioritized compounds resulted in 25 prioritized candidates for further evaluation. A subset of these compounds (13/25) was removed due to factors such as known issues with cytotoxicity and not being amenable to chemical modifications. Additionally, 6 other active drugs (from the 124 from the primary screen) were included in the final prioritized hits based on prior knowledge and literature support yielding the final list of 18 active drugs.

4.8. Expanding the List of Azoles as Specific Drug Class with Anti-Macrofilarial Potency

The azole class of compounds were initially identified from the library of drugs identified for clinical use obtained by the Small Molecule Discovery Center, University of California San Francisco, CA, as hits in the adult Brugia worm assay. Additional azoles were obtained from the Repurposing, Focused Rescue, and Accelerated Medchem (ReFRAME) library, which was generated by Calibr at Scripps Research, La Jolla, CA with support from the Bill and Melinda Gates Foundation. The ReFRAME library is a set of approximately 13,000 compounds that include FDA-approved or registered drugs, investigational drugs that are currently or have been tested previously in clinical phase trials [93].

4.9. Expanding the List of Aspartic Protease Inhibitors

The five FDA approved drugs that inhibit HIV-PR were identified from literature [109], aliskiren was included as an FDA approved renin inhibitor and pepstatin A was included due to being a highly potent general API not approved for clinical usage. APs in filarial worms were identified using orthology with known C. elegans APs, as annotated on Wormbase [119]. The orthologs were identified using Orthofinder [120]. Multiple sequence alignment was done using MAFFT [121], and visualized using Jalview [122]. The phylogenetic tree was estimated using PhyML 3.0 [123] and the node support values are calculated using an “aLRT SH-like” option. The developmental stage expression level values are obtained from [110].

4.10. Clustering of Hits Based on Structural Similarity

The clustering was based on (1—Tanimoto similarity measure) as distance metric, calculated using ChemmineR [124] package, and agglomerated using “complete linkage” method.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2076-0817/10/1/71/s1, Table S1: Motility inhibition percentage values and IC50 values for the primary screening against Brugia pahangi. Table S2: IC50 and motility inhibition percentage values for all species assayed., Table S3: Motility, viability and molting assays across multiple nematode parasite species for the expanded set of Azole compounds.

Author Contributions

Conceptualization, M.M. (Makedonka Mitreva), S.L., J.S.; phenotypic screening, C.A.B., C.F., F.C.-N., N.T., S.J., M.R.A.; computational omics analysis, R.T., M.M. (Makedonka Mitreva); chemogenomics, M.M. (Matthew Mahoney), J.W.J.; resources, J.S., J.H.M., M.R.A., C.W.M., C.M., M.M. (Makedonka Mitreva); writing—original draft preparation, R.T., C.A.B., S.L., J.S., M.M. (Makedonka Mitreva); writing—review and editing, S.L., J.W.J., J.S., M.R.A., C.W.M., C.M., and M.M. (Makedonka Mitreva); supervision, S.L., J.S. and M.M. (Makedonka Mitreva); funding acquisition, S.L., J.S., M.M. (Makedonka Mitreva); All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Bill and Melinda Gates Foundation (OPP1017584) and the National Institute of Health (AI081803). We thank TRS Labs, Athens, GA for supplying Brugia worms, the Small Molecule Discovery Center, University of California San Francisco, CA for providing the library of approved drugs and Calibr, a division of The Scripps Research Institute, La Jolla, CA for providing the Repurposing, Focused Rescue, and Accelerated Medchem (ReFRAME) library.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available in the manuscript main tables and Supplementary Tables.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the phenotypic screens, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Taylor, M.J.; Hoerauf, A.; Bockarie, M. Lymphatic filariasis and onchocerciasis. Lancet 2010, 376, 1175–1185. [Google Scholar] [CrossRef]
  2. Hoerauf, A.; Pfarr, K.; Mand, S.; Debrah, A.Y.; Specht, S. Filariasis in Africa-treatment challenges and prospects. Clin. Microbiol. Infect. 2011, 17, 977–985. [Google Scholar] [CrossRef] [Green Version]
  3. Lustigman, S.; Prichard, R.K.; Gazzinelli, A.; Grant, W.N.; Boatin, B.A.; McCarthy, J.S.; Basanez, M.G. A research agenda for helminth diseases of humans: The problem of helminthiases. PLoS Negl. Trop. Dis. 2012, 6, e1582. [Google Scholar] [CrossRef] [Green Version]
  4. Prichard, R.K.; Basanez, M.G.; Boatin, B.A.; McCarthy, J.S.; Garcia, H.H.; Yang, G.J.; Sripa, B.; Lustigman, S. A research agenda for helminth diseases of humans: Intervention for control and elimination. PLoS Negl. Trop. Dis. 2012, 6, e1549. [Google Scholar] [CrossRef] [Green Version]
  5. King, C.L.; Suamani, J.; Sanuku, N.; Cheng, Y.C.; Satofan, S.; Mancuso, B.; Goss, C.W.; Robinson, L.J.; Siba, P.M.; Weil, G.J.; et al. A Trial of a Triple-Drug Treatment for Lymphatic Filariasis. N. Engl. J. Med. 2018, 379, 1801–1810. [Google Scholar] [CrossRef]
  6. Weil, G.J.; Bogus, J.; Christian, M.; Dubray, C.; Djuardi, Y.; Fischer, P.U.; Goss, C.W.; Hardy, M.; Jambulingam, P.; King, C.L.; et al. The safety of double- and triple-drug community mass drug administration for lymphatic filariasis: A multicenter, open-label, cluster-randomized study. PLoS Med. 2019, 16, e1002839. [Google Scholar] [CrossRef] [Green Version]
  7. Irvine, M.A.; Stolk, W.A.; Smith, M.E.; Subramanian, S.; Singh, B.K.; Weil, G.J.; Michael, E.; Hollingsworth, T.D. Effectiveness of a triple-drug regimen for global elimination of lymphatic filariasis: A modelling study. Lancet Infect. Dis. 2017, 17, 451–458. [Google Scholar] [CrossRef] [Green Version]
  8. Milton, P.; Hamley, J.I.D.; Walker, M.; Basanez, M.G. Moxidectin: An oral treatment for human onchocerciasis. Expert Rev. Anti Infect. Ther. 2020, 18, 1067–1081. [Google Scholar] [CrossRef] [PubMed]
  9. Molyneux, D.H.; Bradley, M.; Hoerauf, A.; Kyelem, D.; Taylor, M.J. Mass drug treatment for lymphatic filariasis and onchocerciasis. Trends Parasitol. 2003, 19, 516–522. [Google Scholar] [CrossRef] [PubMed]
  10. Molyneux, D.H.; Taylor, M.J. Current status and future prospects of the Global Lymphatic Filariasis Programme. Curr. Opin. Infect. Dis. 2001, 14, 155–159. [Google Scholar] [CrossRef] [PubMed]
  11. Molyneux, D.H.; Hopkins, A.; Bradley, M.H.; Kelly-Hope, L.A. Multidimensional complexities of filariasis control in an era of large-scale mass drug administration programmes: A can of worms. Parasites Vectors 2014, 7, 363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Plaisier, A.P.; van Oortmarssen, G.J.; Remme, J.; Habbema, J.D. The reproductive lifespan of Onchocerca volvulus in West African savanna. Acta Trop. 1991, 48, 271–284. [Google Scholar] [CrossRef] [Green Version]
  13. WHO. Global programme to eliminate lymphatic filariasis. Wkly Epidemiol. Rec. 2010, 85, 365–372. [Google Scholar]
  14. Ottesen, E.A.; Hooper, P.J.; Bradley, M.; Biswas, G. The global programme to eliminate lymphatic filariasis: Health impact after 8 years. PLoS Negl. Trop. Dis. 2008, 2, e317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Chu, B.K.; Hooper, P.J.; Bradley, M.H.; McFarland, D.A.; Ottesen, E.A. The economic benefits resulting from the first 8 years of the Global Programme to Eliminate Lymphatic Filariasis (2000–2007). PLoS Negl. Trop. Dis. 2010, 4, e708. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Diawara, L.; Traore, M.O.; Badji, A.; Bissan, Y.; Doumbia, K.; Goita, S.F.; Konate, L.; Mounkoro, K.; Sarr, M.D.; Seck, A.F.; et al. Feasibility of onchocerciasis elimination with ivermectin treatment in endemic foci in Africa: First evidence from studies in Mali and Senegal. PLoS Negl. Trop. Dis. 2009, 3, e497. [Google Scholar] [CrossRef] [Green Version]
  17. Traore, M.O.; Sarr, M.D.; Badji, A.; Bissan, Y.; Diawara, L.; Doumbia, K.; Goita, S.F.; Konate, L.; Mounkoro, K.; Seck, A.F.; et al. Proof-of-principle of onchocerciasis elimination with ivermectin treatment in endemic foci in Africa: Final results of a study in Mali and Senegal. PLoS Negl. Trop. Dis. 2012, 6, e1825. [Google Scholar] [CrossRef] [Green Version]
  18. NTD Modelling Consortium Onchocerciasis Group. The World Health Organization 2030 goals for onchocerciasis: Insights and perspectives from mathematical modelling: NTD Modelling Consortium Onchocerciasis Group. Gates Open Res. 2019, 3, 1545. [Google Scholar] [CrossRef]
  19. NTD Modelling Consortium Lymphatic Filariasis Group. The roadmap towards elimination of lymphatic filariasis by 2030: Insights from quantitative and mathematical modelling. Gates Open Res. 2019, 3, 1538. [Google Scholar] [CrossRef] [Green Version]
  20. Turner, H.C.; Churcher, T.S.; Walker, M.; Osei-Atweneboana, M.Y.; Prichard, R.K.; Basanez, M.G. Uncertainty surrounding projections of the long-term impact of ivermectin treatment on human onchocerciasis. PLoS Negl. Trop. Dis. 2013, 7, e2169. [Google Scholar] [CrossRef] [Green Version]
  21. Turner, H.C.; Walker, M.; Churcher, T.S.; Osei-Atweneboana, M.Y.; Biritwum, N.K.; Hopkins, A.; Prichard, R.K.; Basanez, M.G. Reaching the London Declaration on Neglected Tropical Diseases Goals for Onchocerciasis: An Economic Evaluation of Increasing the Frequency of Ivermectin Treatment in Africa. Clin. Infect. Dis. Off. Publ. Infect. Dis. Soc. Am. 2014, 59, 923–932. [Google Scholar] [CrossRef] [PubMed]
  22. GBD 2015 Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990-2015: A systematic analysis for the Global Burden of Disease Study 2015. Lancet 2016, 388, 1545–1602. [Google Scholar] [CrossRef] [Green Version]
  23. Hopkins, A.D. Neglected tropical diseases in Africa: A new paradigm. Int. Health 2016, 8 (Suppl. 1), i28–i33. [Google Scholar] [CrossRef] [Green Version]
  24. Kim, Y.E.; Remme, J.H.; Steinmann, P.; Stolk, W.A.; Roungou, J.B.; Tediosi, F. Control, elimination, and eradication of river blindness: Scenarios, timelines, and ivermectin treatment needs in Africa. PLoS Negl. Trop. Dis. 2015, 9, e0003664. [Google Scholar] [CrossRef] [Green Version]
  25. African Programme for Onchocerciasis Control (APOC). Report of the Thirty-Eight Session of the Technical Consultative Committee: Ouagadougou; Document Number: DIR/COORD/APOC/REP/TCC38; WHO: Geneva, Switzerland, 2014. [Google Scholar]
  26. The U.S. Food and Drug Administration (FDA); The United States Department of Agriculture (USDA). Drug Approval Package: Moxidectin; The U.S. Food and Drug Administration (FDA): Silver Springs, MD, USA, 2018.
  27. Awadzi, K.; Opoku, N.O.; Attah, S.K.; Lazdins-Helds, J.; Kuesel, A.C. A randomized, single-ascending-dose, ivermectin-controlled, double-blind study of moxidectin in Onchocerca volvulus infection. PLoS Negl. Trop. Dis. 2014, 8, e2953. [Google Scholar] [CrossRef] [PubMed]
  28. Opoku, N.O.; Bakajika, D.K.; Kanza, E.M.; Howard, H.; Mambandu, G.L.; Nyathirombo, A.; Nigo, M.M.; Kasonia, K.; Masembe, S.L.; Mumbere, M.; et al. Single dose moxidectin versus ivermectin for Onchocerca volvulus infection in Ghana, Liberia, and the Democratic Republic of the Congo: A randomised, controlled, double-blind phase 3 trial. Lancet 2018, 392, 1207–1216. [Google Scholar] [CrossRef] [Green Version]
  29. Turner, H.C.; Walker, M.; Attah, S.K.; Opoku, N.O.; Awadzi, K.; Kuesel, A.C.; Basanez, M.G. The potential impact of moxidectin on onchocerciasis elimination in Africa: An economic evaluation based on the Phase II clinical trial data. Parasites Vectors 2015, 8, 167. [Google Scholar] [CrossRef] [Green Version]
  30. Kelly-Hope, L.A.; Cano, J.; Stanton, M.C.; Bockarie, M.J.; Molyneux, D.H. Innovative tools for assessing risks for severe adverse events in areas of overlapping Loa loa and other filarial distributions: The application of micro-stratification mapping. Parasites Vectors 2014, 7, 307. [Google Scholar] [CrossRef] [Green Version]
  31. Boussinesq, M.; Gardon, J.; Gardon-Wendel, N.; Chippaux, J.P. Clinical picture, epidemiology and outcome of Loa-associated serious adverse events related to mass ivermectin treatment of onchocerciasis in Cameroon. Filaria J. 2003, 2, S4. [Google Scholar] [CrossRef] [Green Version]
  32. D’Ambrosio, M.V.; Bakalar, M.; Bennuru, S.; Reber, C.; Skandarajah, A.; Nilsson, L.; Switz, N.; Kamgno, J.; Pion, S.; Boussinesq, M.; et al. Point-of-care quantification of blood-borne filarial parasites with a mobile phone microscope. Sci. Transl. Med. 2015, 7, 286re284. [Google Scholar] [CrossRef] [Green Version]
  33. Boatin, B.A.; Basanez, M.G.; Prichard, R.K.; Awadzi, K.; Barakat, R.M.; Garcia, H.H.; Gazzinelli, A.; Grant, W.N.; McCarthy, J.S.; N’Goran, E.K.; et al. A research agenda for helminth diseases of humans: Towards control and elimination. PLoS Negl. Trop. Dis. 2012, 6, e1547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Osei-Atweneboana, M.Y.; Awadzi, K.; Attah, S.K.; Boakye, D.A.; Gyapong, J.O.; Prichard, R.K. Phenotypic evidence of emerging ivermectin resistance in Onchocerca volvulus. PLoS Negl. Trop. Dis. 2011, 5, e998. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Doyle, S.R.; Bourguinat, C.; Nana-Djeunga, H.C.; Kengne-Ouafo, J.A.; Pion, S.D.S.; Bopda, J.; Kamgno, J.; Wanji, S.; Che, H.; Kuesel, A.C.; et al. Genome-wide analysis of ivermectin response by Onchocerca volvulus reveals that genetic drift and soft selective sweeps contribute to loss of drug sensitivity. PLoS Negl. Trop. Dis. 2017, 11, e0005816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lustigman, S.; McCarter, J.P. Ivermectin Resistance in Onchocerca volvulus: Toward a Genetic Basis. PLoS Negl. Trop. Dis. 2007, 1, e76. [Google Scholar] [CrossRef] [Green Version]
  37. Debrah, A.Y.; Specht, S.; Klarmann-Schulz, U.; Batsa, L.; Mand, S.; Marfo-Debrekyei, Y.; Fimmers, R.; Dubben, B.; Kwarteng, A.; Osei-Atweneboana, M.; et al. Doxycycline Leads to Sterility and Enhanced Killing of Female Onchocerca volvulus Worms in an Area with Persistent Microfilaridermia After Repeated Ivermectin Treatment: A Randomized, Placebo-Controlled, Double-Blind Trial. Clin. Infect. Dis. Off. Publ. Infect. Dis. Soc. Am. 2015, 61, 517–526. [Google Scholar] [CrossRef] [Green Version]
  38. Veale, C.G.L. Unpacking the Pathogen Box-An Open Source Tool for Fighting Neglected Tropical Disease. ChemMedChem 2019, 14, 386–453. [Google Scholar] [CrossRef]
  39. Voronin, D.; Tricoche, N.; Jawahar, S.; Shlossman, M.; Bulman, C.A.; Fischer, C.; Suderman, M.T.; Sakanari, J.A.; Lustigman, S. Development of a preliminary in vitro drug screening assay based on a newly established culturing system for pre-adult fifth-stage Onchocerca volvulus worms. PLoS Negl. Trop. Dis. 2019, 13, e0007108. [Google Scholar] [CrossRef] [Green Version]
  40. International Helminth Genomes Consortium. Comparative genomics of the major parasitic worms. Nat. Genet. 2019, 51, 163–174. [Google Scholar] [CrossRef] [Green Version]
  41. Marcellino, C.; Gut, J.; Lim, K.C.; Singh, R.; McKerrow, J.; Sakanari, J. WormAssay: A novel computer application for whole-plate motion-based screening of macroscopic parasites. PLoS Negl. Trop. Dis. 2012, 6, e1494. [Google Scholar] [CrossRef]
  42. Mendez, D.; Gaulton, A.; Bento, A.P.; Chambers, J.; De Veij, M.; Felix, E.; Magarinos, M.P.; Mosquera, J.F.; Mutowo, P.; Nowotka, M.; et al. ChEMBL: Towards direct deposition of bioassay data. Nucleic Acids Res. 2019, 47, D930–D940. [Google Scholar] [CrossRef]
  43. Sheth, S.; Brito, R.; Mukherjea, D.; Rybak, L.P.; Ramkumar, V. Adenosine receptors: Expression, function and regulation. Int. J. Mol. Sci. 2014, 15, 2024–2052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Waldhoer, M.; Bartlett, S.E.; Whistler, J.L. Opioid receptors. Annu. Rev. Biochem. 2004, 73, 953–990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Hughes, C.E.; Nibbs, R.J.B. A guide to chemokines and their receptors. FEBS J. 2018, 285, 2944–2971. [Google Scholar] [CrossRef] [PubMed]
  46. Kerscher, O.; Felberbaum, R.; Hochstrasser, M. Modification of proteins by ubiquitin and ubiquitin-like proteins. Annu. Rev. Cell Dev. Biol. 2006, 22, 159–180. [Google Scholar] [CrossRef] [Green Version]
  47. Kastenhuber, E.R.; Lowe, S.W. Putting p53 in Context. Cell 2017, 170, 1062–1078. [Google Scholar] [CrossRef] [Green Version]
  48. Balamurugan, K. HIF-1 at the crossroads of hypoxia, inflammation, and cancer. Int. J. Cancer 2016, 138, 1058–1066. [Google Scholar] [CrossRef]
  49. Mashima, R.; Okuyama, T. The role of lipoxygenases in pathophysiology; new insights and future perspectives. Redox Biol. 2015, 6, 297–310. [Google Scholar] [CrossRef] [Green Version]
  50. Sheehan, D.J.; Hitchcock, C.A.; Sibley, C.M. Current and emerging azole antifungal agents. Clin. Microbiol. Rev. 1999, 12, 40–79. [Google Scholar] [CrossRef] [Green Version]
  51. Epstein, B.J.; Vogel, K.; Palmer, B.F. Dihydropyridine calcium channel antagonists in the management of hypertension. Drugs 2007, 67, 1309–1327. [Google Scholar] [CrossRef]
  52. Marella, A.; Tanwar, O.P.; Saha, R.; Ali, M.R.; Srivastava, S.; Akhter, M.; Shaquiquzzaman, M.; Alam, M.M. Quinoline: A versatile heterocyclic. Saudi Pharm. J. 2013, 21, 1–12. [Google Scholar] [CrossRef] [Green Version]
  53. Deliorman, D.; Calis, I.; Ergun, F.; Dogan, B.S.; Buharalioglu, C.K.; Kanzik, I. Studies on the vascular effects of the fractions and phenolic compounds isolated from Viscum album ssp. album. J. Ethnopharmacol. 2000, 72, 323–329. [Google Scholar] [CrossRef]
  54. Malaisse, W.J. Calcium-antagonists and islet function X. Effect of suloctidie. Arch. Int. Pharmacodyn. Ther. 1977, 228, 339–344. [Google Scholar] [PubMed]
  55. Chatelain, P.; Reckinger, N.; Roncucci, R. Effect of suloctidil on Na+/K+ ATPase activity and on membrane fluidity in rat brain synaptosomes. Biochem. Pharmacol. 1979, 28, 3677–3680. [Google Scholar] [CrossRef]
  56. Chen, X.; Ji, Z.L.; Chen, Y.Z. TTD: Therapeutic Target Database. Nucleic Acids Res. 2002, 30, 412–415. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Wise, L.D.; Pattison, I.C.; Butler, D.E.; DeWald, H.A.; Lewis, E.P.; Lobbestael, S.J.; Nordin, I.C.; Poschel, B.P.; Coughenour, L.L. 1-[3-(Diarylamino)propyl]piperidines and related compounds, potential antipsychotic agents with low cataleptogenic profiles. J. Med. Chem. 1985, 28, 606–612. [Google Scholar] [CrossRef]
  58. Katagi, J.; Nakamura, Y.; Cao, X.; Ohara, H.; Honda, A.; Izumi-Nakaseko, H.; Ando, K.; Sugiyama, A. Why Can dl-Sotalol Prolong the QT Interval In Vivo Despite Its Weak Inhibitory Effect on hERG K(+) Channels In Vitro? Electrophysiological and Pharmacokinetic Analysis with the Halothane-Anesthetized Guinea Pig Model. Cardiovasc. Toxicol. 2016, 16, 138–146. [Google Scholar] [CrossRef]
  59. Bang, S.; Yang, T.J.; Yoo, S.; Heo, T.H.; Hwang, S.W. Inhibition of sensory neuronal TRPs contributes to anti-nociception by butamben. Neurosci. Lett. 2012, 506, 297–302. [Google Scholar] [CrossRef]
  60. Finkelstein, M.; Kromer, C.M.; Sweeney, S.A.; Delahunt, C.S. Some aspects of the pharmacology of clemizole hydrochloride. J. Am. Pharm. Assoc. 1960, 49, 18–22. [Google Scholar] [CrossRef]
  61. Richter, J.M.; Schaefer, M.; Hill, K. Clemizole hydrochloride is a novel and potent inhibitor of transient receptor potential channel TRPC5. Mol. Pharmacol. 2014, 86, 514–521. [Google Scholar] [CrossRef] [Green Version]
  62. Ivanov, I.V.; Bondarenko, R.A. Individual Approach to the Use of Medications to Normalize Situational Anxiety under the Psycho-Emotional Stress. Aviakosm. Ekolog. Med. 2015, 49, 32–35. [Google Scholar]
  63. Mahon, B.P.; Hendon, A.M.; Driscoll, J.M.; Rankin, G.M.; Poulsen, S.A.; Supuran, C.T.; McKenna, R. Saccharin: A lead compound for structure-based drug design of carbonic anhydrase IX inhibitors. Bioorg. Med. Chem. 2015, 23, 849–854. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Kempf, D.J.; Sham, H.L.; Marsh, K.C.; Flentge, C.A.; Betebenner, D.; Green, B.E.; McDonald, E.; Vasavanonda, S.; Saldivar, A.; Wideburg, N.E.; et al. Discovery of ritonavir, a potent inhibitor of HIV protease with high oral bioavailability and clinical efficacy. J. Med. Chem. 1998, 41, 602–617. [Google Scholar] [CrossRef] [PubMed]
  65. Ikezoe, T.; Hisatake, Y.; Takeuchi, T.; Ohtsuki, Y.; Yang, Y.; Said, J.W.; Taguchi, H.; Koeffler, H.P. HIV-1 protease inhibitor, ritonavir: A potent inhibitor of CYP3A4, enhanced the anticancer effects of docetaxel in androgen-independent prostate cancer cells in vitro and in vivo. Cancer Res. 2004, 64, 7426–7431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Gunthard, H.F.; Aberg, J.A.; Eron, J.J.; Hoy, J.F.; Telenti, A.; Benson, C.A.; Burger, D.M.; Cahn, P.; Gallant, J.E.; Glesby, M.J.; et al. Antiretroviral treatment of adult HIV infection: 2014 recommendations of the International Antiviral Society-USA Panel. JAMA 2014, 312, 410–425. [Google Scholar] [CrossRef]
  67. Schlessinger, A.; Geier, E.; Fan, H.; Irwin, J.J.; Shoichet, B.K.; Giacomini, K.M.; Sali, A. Structure-based discovery of prescription drugs that interact with the norepinephrine transporter, NET. Proc. Natl Acad. Sci. USA 2011, 108, 15810–15815. [Google Scholar] [CrossRef] [Green Version]
  68. Kim, S.; Thiessen, P.A.; Cheng, T.; Yu, B.; Shoemaker, B.A.; Wang, J.; Bolton, E.E.; Wang, Y.; Bryant, S.H. Literature information in PubChem: Associations between PubChem records and scientific articles. J. Cheminform. 2016, 8, 32. [Google Scholar] [CrossRef] [Green Version]
  69. Carvalho, L.; Luque-Ortega, J.R.; Manzano, J.I.; Castanys, S.; Rivas, L.; Gamarro, F. Tafenoquine, an antiplasmodial 8-aminoquinoline, targets leishmania respiratory complex III and induces apoptosis. Antimicrob. Agents Chemother. 2010, 54, 5344–5351. [Google Scholar] [CrossRef] [Green Version]
  70. Graves, P.R.; Kwiek, J.J.; Fadden, P.; Ray, R.; Hardeman, K.; Coley, A.M.; Foley, M.; Haystead, T.A. Discovery of novel targets of quinoline drugs in the human purine binding proteome. Mol. Pharmacol. 2002, 62, 1364–1372. [Google Scholar] [CrossRef] [Green Version]
  71. Chu, C.S.; Bancone, G.; Nosten, F.; White, N.J.; Luzzatto, L. Primaquine-induced haemolysis in females heterozygous for G6PD deficiency. Malar. J. 2018, 17, 101. [Google Scholar] [CrossRef]
  72. Anden, N.E.; Rubenson, A.; Fuxe, K.; Hokfelt, T. Evidence for dopamine receptor stimulation by apomorphine. J. Pharm. Pharmacol. 1967, 19, 627–629. [Google Scholar] [CrossRef]
  73. Anden, N.E.; Strombom, U. Adrenergic receptor blocking agents: Effects on central noradrenaline and dopamine receptors and on motor activity. Psychopharmacologia 1974, 38, 91–103. [Google Scholar] [CrossRef] [PubMed]
  74. Warrilow, A.G.; Parker, J.E.; Kelly, D.E.; Kelly, S.L. Azole affinity of sterol 14alpha-demethylase (CYP51) enzymes from Candida albicans and Homo sapiens. Antimicrob. Agents Chemother. 2013, 57, 1352–1360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Thomson, S.; Rice, C.A.; Zhang, T.; Edrada-Ebel, R.; Henriquez, F.L.; Roberts, C.W. Characterisation of sterol biosynthesis and validation of 14alpha-demethylase as a drug target in Acanthamoeba. Sci. Rep. 2017, 7, 8247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Carrillo-Munoz, A.J.; Tur-Tur, C.; Giusiano, G.; Marcos-Arias, C.; Eraso, E.; Jauregizar, N.; Quindos, G. Sertaconazole: An antifungal agent for the topical treatment of superficial candidiasis. Expert Rev. Anti Infect. Ther. 2013, 11, 347–358. [Google Scholar] [CrossRef] [PubMed]
  77. Warrilow, A.G.; Hull, C.M.; Rolley, N.J.; Parker, J.E.; Nes, W.D.; Smith, S.N.; Kelly, D.E.; Kelly, S.L. Clotrimazole as a potent agent for treating the oomycete fish pathogen Saprolegnia parasitica through inhibition of sterol 14alpha-demethylase (CYP51). Appl. Environ. Microbiol. 2014, 80, 6154–6166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Pace, J.R.; DeBerardinis, A.M.; Sail, V.; Tacheva-Grigorova, S.K.; Chan, K.A.; Tran, R.; Raccuia, D.S.; Wechsler-Reya, R.J.; Hadden, M.K. Repurposing the Clinically Efficacious Antifungal Agent Itraconazole as an Anticancer Chemotherapeutic. J. Med. Chem. 2016, 59, 3635–3649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Odds, F.C.; Brown, A.J.; Gow, N.A. Antifungal agents: Mechanisms of action. Trends Microbiol. 2003, 11, 272–279. [Google Scholar] [CrossRef]
  80. Makepeace, B.L.; Tanya, V.N. 25 Years of the Onchocerca ochengi Model. Trends Parasitol. 2016, 32, 966–978. [Google Scholar] [CrossRef] [Green Version]
  81. Gardon, J.; Gardon-Wendel, N.; Demanga, N.; Kamgno, J.; Chippaux, J.P.; Boussinesq, M. Serious reactions after mass treatment of onchocerciasis with ivermectin in an area endemic for Loa loa infection. Lancet 1997, 350, 18–22. [Google Scholar] [CrossRef]
  82. Drugs for Neglected Diseases Initiative. Filaria: River blindness, Oxfendazole. Available online: https://dndi.org/research-development/portfolio/oxfendazole/ (accessed on 11 November 2020).
  83. Drugs for Neglected Diseases Initiative. Filaria: River Blindness, Emodepside. Available online: https://dndi.org/research-development/portfolio/emodepside/ (accessed on 11 November 2020).
  84. Hubner, M.P.; Martin, C.; Specht, S.; Koschel, M.; Dubben, B.; Frohberger, S.J.; Ehrens, A.; Fendler, M.; Struever, D.; Mitre, E.; et al. Oxfendazole mediates macrofilaricidal efficacy against the filarial nematode Litomosoides sigmodontis in vivo and inhibits Onchocerca spec. motility in vitro. PLoS Negl. Trop. Dis. 2020, 14, e0008427. [Google Scholar] [CrossRef]
  85. Kashyap, S.S.; Verma, S.; Voronin, D.; Lustigman, S.; Kulke, D.; Robertson, A.P.; Martin, R.J. Emodepside has sex-dependent immobilizing effects on adult Brugia malayi due to a differentially spliced binding pocket in the RCK1 region of the SLO-1 K channel. PLoS Pathog. 2019, 15, e1008041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Maertens, J.A. History of the development of azole derivatives. Clin. Microbiol. Infect. 2004, 10 (Suppl. 1), 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Gaulton, A.; Hersey, A.; Nowotka, M.; Bento, A.P.; Chambers, J.; Mendez, D.; Mutowo, P.; Atkinson, F.; Bellis, L.J.; Cibrian-Uhalte, E.; et al. The ChEMBL database in 2017. Nucleic Acids Res. 2017, 45, D945–D954. [Google Scholar] [CrossRef]
  88. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef] [PubMed]
  89. Tyagi, R.; Elfawal, M.A.; Wildman, S.A.; Helander, J.; Bulman, C.A.; Sakanari, J.; Rosa, B.A.; Brindley, P.J.; Janetka, J.W.; Aroian, R.V.; et al. Identification of small molecule enzyme inhibitors as broad-spectrum anthelmintics. Sci. Rep. 2019, 9, 9085. [Google Scholar] [CrossRef] [Green Version]
  90. Seidler, J.; McGovern, S.L.; Doman, T.N.; Shoichet, B.K. Identification and prediction of promiscuous aggregating inhibitors among known drugs. J. Med. Chem. 2003, 46, 4477–4486. [Google Scholar] [CrossRef]
  91. Choi, S.R.; Pradhan, A.; Hammond, N.L.; Chittiboyina, A.G.; Tekwani, B.L.; Avery, M.A. Design, synthesis, and biological evaluation of Plasmodium falciparum lactate dehydrogenase inhibitors. J. Med. Chem. 2007, 50, 3841–3850. [Google Scholar] [CrossRef]
  92. Mast, N.; Zheng, W.; Stout, C.D.; Pikuleva, I.A. Antifungal Azoles: Structural Insights into Undesired Tight Binding to Cholesterol-Metabolizing CYP46A1. Mol. Pharmacol. 2013, 84, 86–94. [Google Scholar] [CrossRef] [Green Version]
  93. Janes, J.; Young, M.E.; Chen, E.; Rogers, N.H.; Burgstaller-Muehlbacher, S.; Hughes, L.D.; Love, M.S.; Hull, M.V.; Kuhen, K.L.; Woods, A.K.; et al. The ReFRAME library as a comprehensive drug repurposing library and its application to the treatment of cryptosporidiosis. Proc. Natl. Acad. Sci. USA 2018, 115, 10750–10755. [Google Scholar] [CrossRef] [Green Version]
  94. Eder, J.; Hommel, U.; Cumin, F.; Martoglio, B.; Gerhartz, B. Aspartic proteases in drug discovery. Curr. Pharm. Des. 2007, 13, 271–285. [Google Scholar] [CrossRef]
  95. Nguyen, J.T.; Hamada, Y.; Kimura, T.; Kiso, Y. Design of potent aspartic protease inhibitors to treat various diseases. Arch. Pharm. (Weinheim) 2008, 341, 523–535. [Google Scholar] [CrossRef] [PubMed]
  96. WHO. The Selection and Use of Essential Medicines: Report of the WHO Expert Committee on Selection and Use of Essential Medicines, 2019 (Including the 21st WHO Model List of Essential Medicines and the 7th WHO Model List of Essential Medicines for Children); WHO: Geneva, Switzerland, 2019. [Google Scholar]
  97. Tort, J.; Brindley, P.J.; Knox, D.; Wolfe, K.H.; Dalton, J.P. Proteinases and associated genes of parasitic helminths. Adv. Parasitol. 1999, 43, 161–266. [Google Scholar] [CrossRef] [PubMed]
  98. Xu, J.; Liu, R.D.; Bai, S.J.; Hao, H.N.; Yue, W.W.; Xu, Y.X.Y.; Long, S.R.; Cui, J.; Wang, Z.Q. Molecular characterization of a Trichinella spiralis aspartic protease and its facilitation role in larval invasion of host intestinal epithelial cells. PLoS Negl. Trop. Dis. 2020, 14, e0008269. [Google Scholar] [CrossRef] [PubMed]
  99. Park, J.N.; Park, S.K.; Cho, M.K.; Park, M.K.; Kang, S.A.; Kim, D.H.; Yu, H.S. Molecular characterization of 45 kDa aspartic protease of Trichinella spiralis. Vet. Parasitol. 2012, 190, 510–518. [Google Scholar] [CrossRef] [PubMed]
  100. Williamson, A.L.; Brindley, P.J.; Abbenante, G.; Datu, B.J.; Prociv, P.; Berry, C.; Girdwood, K.; Pritchard, D.I.; Fairlie, D.P.; Hotez, P.J.; et al. Hookworm aspartic protease, Na-APR-2, cleaves human hemoglobin and serum proteins in a host-specific fashion. J. Infect. Dis. 2003, 187, 484–494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Jolodar, A.; Fischer, P.; Buttner, D.W.; Miller, D.J.; Schmetz, C.; Brattig, N.W. Onchocerca volvulus: Expression and immunolocalization of a nematode cathepsin D-like lysosomal aspartic protease. Exp. Parasitol. 2004, 107, 145–156. [Google Scholar] [CrossRef] [PubMed]
  102. Hamada, Y.; Kiso, Y. New directions for protease inhibitors directed drug discovery. Biopolymers 2016, 106, 563–579. [Google Scholar] [CrossRef]
  103. Drag, M.; Salvesen, G.S. Emerging principles in protease-based drug discovery. Nat. Rev. Drug Discov. 2010, 9, 690–701. [Google Scholar] [CrossRef] [Green Version]
  104. Monika, S.; Malgorzata, B.; Zbigniew, O. Contribution of Aspartic Proteases in Candida Virulence. Protease Inhibitors against Candida Infections. Curr. Protein Pept. Sci. 2017, 18, 1050–1062. [Google Scholar] [CrossRef]
  105. Castilho, V.V.S.; Goncalves, K.C.S.; Rebello, K.M.; Baptista, L.P.R.; Sangenito, L.S.; Santos, H.L.C.; Branquinha, M.H.; Santos, A.L.S.; Menna-Barreto, R.F.S.; Guimaraes, A.C.; et al. Docking simulation between HIV peptidase inhibitors and Trypanosoma cruzi aspartyl peptidase. BMC Res. Notes 2018, 11, 825. [Google Scholar] [CrossRef] [Green Version]
  106. Rawlings, N.D.; Barrett, A.J.; Finn, R. Twenty years of the MEROPS database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2016, 44, D343–D350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Benes, P.; Vetvicka, V.; Fusek, M. Cathepsin D-many functions of one aspartic protease. Crit. Rev. Oncol. Hematol. 2008, 68, 12–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Sarkis, G.J.; Kurpiewski, M.R.; Ashcom, J.D.; Jen-Jacobson, L.; Jacobson, L.A. Proteases of the nematode Caenorhabditis elegans. Arch. Biochem. Biophys. 1988, 261, 80–90. [Google Scholar] [CrossRef]
  109. Lv, Z.; Chu, Y.; Wang, Y. HIV protease inhibitors: A review of molecular selectivity and toxicity. HIV AIDS (Auckl.) 2015, 7, 95–104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  110. Cotton, J.A.; Bennuru, S.; Grote, A.; Harsha, B.; Tracey, A.; Beech, R.; Doyle, S.R.; Dunn, M.; Hotopp, J.C.; Holroyd, N.; et al. The genome of Onchocerca volvulus, agent of river blindness. Nat. Microbiol. 2016, 2, 16216. [Google Scholar] [CrossRef] [PubMed]
  111. Rigobello, M.P.; Scutari, G.; Boscolo, R.; Bindoli, A. Induction of mitochondrial permeability transition by auranofin, a gold(I)-phosphine derivative. Br. J. Pharmacol. 2002, 136, 1162–1168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Njouendou, A.J.; Ritter, M.; Ndongmo, W.P.C.; Kien, C.A.; Narcisse, G.T.V.; Fombad, F.F.; Tayong, D.B.; Pfarr, K.; Layland, L.E.; Hoerauf, A.; et al. Successful long-term maintenance of Mansonella perstans in an in vitro culture system. Parasites Vectors 2017, 10, 563. [Google Scholar] [CrossRef] [Green Version]
  113. Verma, M.; Pathak, M.; Shahab, M.; Singh, K.; Mitra, K.; Misra-Bhattacharya, S. Moxidectin causes adult worm mortality of human lymphatic filarial parasite Brugia malayi in rodent models. Folia Parasitol. (Praha) 2014, 61, 561–570. [Google Scholar] [CrossRef]
  114. Verma, S.; Kashyap, S.S.; Robertson, A.P.; Martin, R.J. Functional genomics in Brugia malayi reveal diverse muscle nAChRs and differences between cholinergic anthelmintics. Proc. Natl. Acad. Sci. USA 2017, 114, 5539–5544. [Google Scholar] [CrossRef] [Green Version]
  115. Partridge, F.A.; Forman, R.; Bataille, C.J.R.; Wynne, G.M.; Nick, M.; Russell, A.J.; Else, K.J.; Sattelle, D.B. Anthelmintic drug discovery: Target identification, screening methods and the role of open science. Beilstein J. Org. Chem. 2020, 16, 1203–1224. [Google Scholar] [CrossRef]
  116. Marcellino, C. GitHub Repository. WormAssay. Available online: https://github.com/chrismarcellino/wormassay. (accessed on 23 September 2020).
  117. Bulman, C.A.; Bidlow, C.M.; Lustigman, S.; Cho-Ngwa, F.; Williams, D.; Rascon, A.A., Jr.; Tricoche, N.; Samje, M.; Bell, A.; Suzuki, B.; et al. Repurposing auranofin as a lead candidate for treatment of lymphatic filariasis and onchocerciasis. PLoS Negl. Trop. Dis. 2015, 9, e0003534. [Google Scholar] [CrossRef] [PubMed]
  118. Bickerton, G.R.; Paolini, G.V.; Besnard, J.; Muresan, S.; Hopkins, A.L. Quantifying the chemical beauty of drugs. Nat. Chem. 2012, 4, 90–98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  119. Harris, T.W.; Arnaboldi, V.; Cain, S.; Chan, J.; Chen, W.J.; Cho, J.; Davis, P.; Gao, S.; Grove, C.A.; Kishore, R.; et al. WormBase: A modern Model Organism Information Resource. Nucleic Acids Res. 2020, 48, D762–D767. [Google Scholar] [CrossRef] [PubMed]
  120. Emms, D.M.; Kelly, S. OrthoFinder: Phylogenetic orthology inference for comparative genomics. Genome Biol. 2019, 20, 238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  121. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  122. Waterhouse, A.M.; Procter, J.B.; Martin, D.M.; Clamp, M.; Barton, G.J. Jalview Version 2—A multiple sequence alignment editor and analysis workbench. Bioinformatics 2009, 25, 1189–1191. [Google Scholar] [CrossRef] [Green Version]
  123. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [Green Version]
  124. Cao, Y.; Charisi, A.; Cheng, L.C.; Jiang, T.; Girke, T. ChemmineR: A compound mining framework for R. Bioinformatics 2008, 24, 1733–1734. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Effect of 124 repurposed drugs on B. pahangi motility (3-day assay). (a) Distribution of % motility inhibition for 124 active drugs. (b) IC50 values for 62 of the active drugs. (c) % motility inhibition vs IC50.
Figure 1. Effect of 124 repurposed drugs on B. pahangi motility (3-day assay). (a) Distribution of % motility inhibition for 124 active drugs. (b) IC50 values for 62 of the active drugs. (c) % motility inhibition vs IC50.
Pathogens 10 00071 g001
Figure 2. Identifying top hit compounds and functions associated with their targets. (A) Integrating experimental and computational prioritization schemes to obtain the top 25 candidate compounds. (B) Molecular functions represented by the 31 ChEMBL target proteins associated with the 124 active drugs.
Figure 2. Identifying top hit compounds and functions associated with their targets. (A) Integrating experimental and computational prioritization schemes to obtain the top 25 candidate compounds. (B) Molecular functions represented by the 31 ChEMBL target proteins associated with the 124 active drugs.
Pathogens 10 00071 g002
Figure 3. Chemical structures of the 18 active drugs on adult female B. pahangi divided into (a) Group A, and (b) Group B compounds. Compounds in Group B are classified into high priority because they are amenable to modifications (Black), or low priority due to central nervous system activity and due to unknown mode of action (Grey). The HEA moieties in two Group B compounds are highlighted with green boxes.
Figure 3. Chemical structures of the 18 active drugs on adult female B. pahangi divided into (a) Group A, and (b) Group B compounds. Compounds in Group B are classified into high priority because they are amenable to modifications (Black), or low priority due to central nervous system activity and due to unknown mode of action (Grey). The HEA moieties in two Group B compounds are highlighted with green boxes.
Pathogens 10 00071 g003
Figure 4. (a) A maximum-likelihood phylogenetic tree based on sequence similarity of aspartic proteases in B. malayi, O. volvulus, HIV and two Candida species. (b) Developmental stage-specific (mf, L3, and male or female adult worms) gene expression levels of the filarial APs in B. malayi and O. volvulus, as indicated by their protein IDs starting with Bm and OVOC, respectively. (c) Amino acid sequence alignments of aspartic proteases (shown only around the DTG/DTG and DTG/DSG motifs), with the two aspartic acid residues in the active site indicated by red arrows. HIV-PR only uses a single DTG active site as it forms a homodimer to constitute its protease activity. The alignment is color coded using the Zappo color scheme of Jalview, colors the residues according to their physicochemical properties—aliphatic or hydrophobic (pink), orange (aromatic), red (negative), green (hydrophilic), purple (conformationally special). Only the residues with moderate or high conservation are colored (i.e., conserved in at least 5 out of 11 aligned residues).
Figure 4. (a) A maximum-likelihood phylogenetic tree based on sequence similarity of aspartic proteases in B. malayi, O. volvulus, HIV and two Candida species. (b) Developmental stage-specific (mf, L3, and male or female adult worms) gene expression levels of the filarial APs in B. malayi and O. volvulus, as indicated by their protein IDs starting with Bm and OVOC, respectively. (c) Amino acid sequence alignments of aspartic proteases (shown only around the DTG/DTG and DTG/DSG motifs), with the two aspartic acid residues in the active site indicated by red arrows. HIV-PR only uses a single DTG active site as it forms a homodimer to constitute its protease activity. The alignment is color coded using the Zappo color scheme of Jalview, colors the residues according to their physicochemical properties—aliphatic or hydrophobic (pink), orange (aromatic), red (negative), green (hydrophilic), purple (conformationally special). Only the residues with moderate or high conservation are colored (i.e., conserved in at least 5 out of 11 aligned residues).
Pathogens 10 00071 g004
Figure 5. Structural similarity of the hit drugs identified in the primary and secondary phenotypic in vitro screen. Hits: at least 50% motility inhibition at 10 μM or 30 μM drug treatment of either female B. pahangi, or viability (female) and motility (male) inhibition of O. ochengi, or motility inhibition of O. volvulus L4. Red font: Azoles (Group A from primary screening hits and active azoles from expanded set). Blue font: active aspartic protease inhibitors and HEA containing compounds from primary screening hits.
Figure 5. Structural similarity of the hit drugs identified in the primary and secondary phenotypic in vitro screen. Hits: at least 50% motility inhibition at 10 μM or 30 μM drug treatment of either female B. pahangi, or viability (female) and motility (male) inhibition of O. ochengi, or motility inhibition of O. volvulus L4. Red font: Azoles (Group A from primary screening hits and active azoles from expanded set). Blue font: active aspartic protease inhibitors and HEA containing compounds from primary screening hits.
Pathogens 10 00071 g005
Table 1. The 18 top hits (female Brugia pahangi, 3-day assay, >50% inhibition of motility and IC50 values from a subset of the hits).
Table 1. The 18 top hits (female Brugia pahangi, 3-day assay, >50% inhibition of motility and IC50 values from a subset of the hits).
Drug Common NameCompound ClassCompound SubclassLiterature Targets
(Not All are Listed)
Brugia pahangi (Adult Female)
% Motility Inhibition (10 μM)IC50 (μM)
Group AEconazole ImidazoleBenzyl etherSterol 14-demethylase/K+ VGIC [72]1004.0
Miconazole ImidazoleBenzyl etherSterol 14-demethylase/K+ VGIC [72]963.3
SulconazoleImidazoleThiobenzyl etherSterol 14-demethylase [73]984.6
SertaconazoleImidazoleBenzothiopheneSterol 14-demethylase [74]100-
ClotrimazoleImidazoleTriphenylmethaneSterol 14-demethylase/K+ VGIC [75]855.5
ItraconazoleTriazolePiperazinylXIAP/K+ VGIC [76,77]80-
Group BSuloctidilAmino alcoholPhenyl thioetherK+ VGIC/thyroid receptor [52,53]993.7
PimozideBenzimidazolonePiperidineDRD2/DRD3 (GPCRs)/KCNH2/Calmodulin [54,55,56]962.5
ButambenBenzoic esterAnilineVGICs [57]66-
Clemizole BenzimidazolePyrrolidineHRH1/K+ VGIC/NR ROR-gamma [58,59]786.5
ProroxanPhenyl ketonePhenyl pyrrolidineAdrenergic receptors [60]658.5
SaccharinBenzoisothiazoloneSulfonyl amidecarbonic anhydrase [61]57-
RitonavirAmino alcoholPeptidomimeticHIV protease/CYP3A/SLC47 [62,63,64]52-
NifenazoneDihydropyrazoleNicotinamide-50-
LevonordefrinAmino alcoholPhenolADRA2/HADH2/TF HIF1A/APE1 [65,66]54-
Tafenoquine QuinolineAlkyl aminecytochrome c reductase [67]6314.1
Primaquine QuinolineAlkyl aminequinone reductase [68,69]515.1
ApomorphineAporphineTetrahydroisoquinolineDopamine/serotonin/adrenergic receptor agonist [70,71]53-
Table 2. Expanded list of Azoles and 3 best Group B hit drugs against multiple filarial species.
Table 2. Expanded list of Azoles and 3 best Group B hit drugs against multiple filarial species.
Drug Common NameDay 3 Adult Female
B. pahangi
Day 7 Adult Female
O. ochengi
Day 5 Adult Male
O. ochengi
Day 28 O. volvulus L5Day 5 Loa loa mf
% Inhibition of Motility (10 μM)IC50 (μM)% Inhibition of Viability (10 μM)IC50 (μM)% Inhibition of Motility (10 μM)IC50 (μM)% Inhibition of Motility (10 μM)IC50 (μM)% Inhibition of Motility (10 μM)IC50 (μM)
Group AIsoconazole1003.85730.81000.6962.433.310
Fenticonazole1001.23294.77510.6834.75010
Sertaconazole 100-21-31-----
Econazole 1004.022-26-----
Tioconazole993.214-100-----
Sulconazole984.626-99-----
Miconazole 963.30-100-----
Clotrimazole855.58-100-----
Posaconazole820.138>100796.157 *---
Itraconazole80-13-21-----
Group BSuloctidil993.7874.11005.5100 **-1004.1
Pimozide962.513-96-----
Primaquine 515.11001.31000.4100 **-5018.5
* Screen was done at 1 µM; ** 100% effective on day 21 and 100% inhibition of viability as measured by MTT assay; -, not tested.
Table 3. Known aspartic protease inhibitors with anti-filarial activity.
Table 3. Known aspartic protease inhibitors with anti-filarial activity.
Aspartyl Protease Inhibitors (APIs)% Inhibition of Motility (30 µM)% Inhibition of Molting (L3 to L4)
B. pahangiO. ochengiO. volvulusO. volvulus (day 6)
Female (day 6)Female (day 7) **Male (day 5)L4 (day 7)3 μM10 μM
Nelfinavir991001007350.245.5
Lopinavir9808856074.1
Ritonavir52 *63100000
Pepstatin A0501000040
Darunavir030ND4.443.946.7
Aliskiren06ND030.451.7
Amprenavir240178043.5
Atazanavir1142ND3000
* Screen was done at 10 µM; ** % viability inhibition as measured by MTT assay; ND, could not be determined.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tyagi, R.; Bulman, C.A.; Cho-Ngwa, F.; Fischer, C.; Marcellino, C.; Arkin, M.R.; McKerrow, J.H.; McNamara, C.W.; Mahoney, M.; Tricoche, N.; et al. An Integrated Approach to Identify New Anti-Filarial Leads to Treat River Blindness, a Neglected Tropical Disease. Pathogens 2021, 10, 71. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010071

AMA Style

Tyagi R, Bulman CA, Cho-Ngwa F, Fischer C, Marcellino C, Arkin MR, McKerrow JH, McNamara CW, Mahoney M, Tricoche N, et al. An Integrated Approach to Identify New Anti-Filarial Leads to Treat River Blindness, a Neglected Tropical Disease. Pathogens. 2021; 10(1):71. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010071

Chicago/Turabian Style

Tyagi, Rahul, Christina A. Bulman, Fidelis Cho-Ngwa, Chelsea Fischer, Chris Marcellino, Michelle R. Arkin, James H. McKerrow, Case W. McNamara, Matthew Mahoney, Nancy Tricoche, and et al. 2021. "An Integrated Approach to Identify New Anti-Filarial Leads to Treat River Blindness, a Neglected Tropical Disease" Pathogens 10, no. 1: 71. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10010071

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop