Metabolic network analysis predicts efficacy of FDA-approved drugs targeting the causative agent of a neglected tropical disease
© Chavali et al; licensee BioMed Central Ltd. 2012
Received: 16 January 2012
Accepted: 27 April 2012
Published: 27 April 2012
Systems biology holds promise as a new approach to drug target identification and drug discovery against neglected tropical diseases. Genome-scale metabolic reconstructions, assembled from annotated genomes and a vast array of bioinformatics/biochemical resources, provide a framework for the interrogation of human pathogens and serve as a platform for generation of future experimental hypotheses. In this article, with the application of selection criteria for both Leishmania major targets (e.g. in silico gene lethality) and drugs (e.g. toxicity), a method (MetDP) to rationally focus on a subset of low-toxic Food and Drug Administration (FDA)-approved drugs is introduced.
This metabolic network-driven approach identified 15 L. major genes as high-priority targets, 8 high-priority synthetic lethal targets, and 254 FDA-approved drugs. Results were compared to previous literature findings and existing high-throughput screens. Halofantrine, an antimalarial agent that was prioritized using MetDP, showed noticeable antileishmanial activity when experimentally evaluated in vitro against L. major promastigotes. Furthermore, synthetic lethality predictions also aided in the prediction of superadditive drug combinations. For proof-of-concept, double-drug combinations were evaluated in vitro against L. major and four combinations involving the drug disulfiram that showed superadditivity are presented.
A direct metabolic network-driven method that incorporates single gene essentiality and synthetic lethality predictions is proposed that generates a set of high-priority L. major targets, which are in turn associated with a select number of FDA-approved drugs that are candidate antileishmanials. Additionally, selection of high-priority double-drug combinations might provide for an attractive and alternative avenue for drug discovery against leishmaniasis.
Over one billion people are infected by one or more neglected tropical diseases (NTDs) [1, 2]. These diseases comprise a group of parasitic and bacterial infections that affect some of the poorest and most marginalized populations around the world, many who live on less than 1.25 USD a day . Leishmaniasis is one such NTD that is endemic in 88 countries with a total of 350 million people at risk. The disease is associated with a global prevalence of 12 million cases, a yearly incidence of 1.5 to 2 million cases, and an annual mortality rate of over 59,000 deaths . Caused by Leishmania species, the disease can manifest itself in varying clinical pathologies including visceral, cutaneous or mucocutaneous forms [3, 4]. These parasites transition through different morphological stages - from flagellated promastigotes within a female phlebotomine sandfly (vector) gut to non-flagellated amastigotes in host macrophages . Leishmaniasis and other NTDs tend to be overshadowed by research that is focused on malaria, tuberculosis and HIV/AIDS . There may also be a lack of interest by pharmaceutical companies in pursuing drug development against NTDs due to a complicated profit motive . Issues with high toxicity, life-threatening side-effects, cost, parenteral administration and emergence of resistance are drawbacks common to existing antileishmanial drugs such as sodium stibogluconate, meglumine antimoniate, pentamidine and amphotericin B [4, 8]. Therefore, with limited viable treatment options and few alternatives on the pipeline, there exists an obvious need to develop radically new strategies and conceptual frameworks to identify compounds with novel antileishmanial activity.
Computational metabolic networks have been reconstructed for several pathogenic organisms [9, 10], including Leishmania major, an agent of cutaneous leishmaniasis . Implementing constraint-based modeling techniques, such as flux balance analysis (FBA), on network reconstructions has yielded valuable insight into gene essentiality and enzyme robustness under varying environmental conditions [11, 12]. Accordingly, analysis of pathogen intracellular processes can be used to predict critical protein targets that when perturbed singly or in combination have adverse effects on virulence and/or growth. Here, a direct metabolic network-driven method that incorporates single gene essentiality and synthetic lethality predictions is proposed that generates a set of high-priority L. major targets, which are in turn associated with a select number of Food and Drug Administration (FDA)-approved drugs that are candidate antileishmanials. Subsequently, candidate drugs, which are already in clinical use for other indications, can begin to be investigated for clinical use against leishmaniasis. In addition, the selection of high-priority double-drug combinations that demonstrate superadditivity might provide for an attractive and alternative avenue for drug discovery against leishmaniasis. By integrating publicly available resources using a direct and efficient method, the approach presented here offers significant implications to future drug discovery and drug repurposing strategies against leishmaniasis and other NTDs.
Conceptual platform of MetDP: A pipeline for prioritization of targets and FDA-approved drugs
The list of prioritized L. major targets
No. of drugs
urea cycle and metabolism of amino groups
> 30% (GR)
hydroxymethylglutaryl CoA reductase
fructose and mannose metabolism
beta-ketoacyl-acyl-carrier-protein synthase I
fatty acid biosynthesis
pentose phosphate pathway
The list of prioritized synthetic lethal targets
No. of drugs
No. of drugs
ribonucleoside-diphosphate reductase small chain
ribonucleoside-diphosphate reductase small chain
dihydrofolate reductase-thymidylate synthase
F-type H+-transporting ATPase alpha chain
F-type H+-transporting ATPase alpha chain
F-type H+-transporting ATPase beta chain
F-type H+-transporting ATPase beta chain
The final step in MetDP involved considerations of tolerance and toxicity in order to place emphasis on drugs that have potential clinical relevance against leishmaniasis. Drugs that are known common metabolites or chemical elements (e.g. NADH, adenine, ATP, amino acids, ethanol, iron and zinc) and present in the metabolic reconstruction were removed from further analysis. Also, any illicit drugs and controlled substances (Schedule I through V) were also removed. Second, lethal dose 50 (LD50) data were culled from DrugBank database and material safety data sheets available online, and a toxicity rating based on the Hodge and Sterner scale (see Additional file 1: Table S1) was applied. Following tolerance and toxicity analysis, the 15 single drug targets were associated with 240 drugs with a toxicity rating greater than or equal to 3 ('moderately toxic' on the Hodge and Sterner scale). Meanwhile, a total of 107 drugs were associated with the 8 non-trivial lethal gene combinations. Of these, 37 drugs were mapped to both genes involved in any one combination and 14 of them were unique when compared to the list of 240 drugs (see Figure 2B).
Consequently, upon implementation of MetDP by means of iterating through the pipeline (presented in Figures 2A and 2B), the prioritized list of single L. major targets included 15 metabolic genes (out of 560 genes in the metabolic reconstruction; 2.7%) while the prioritized list of synthetic lethal targets included 8 double-gene combinations (out of 156,520; 0.005%). Collectively, these lists of high-priority targets were associated with 254 FDA-approved drugs (out of 4329 drugs in DrugBank; 5.9%). The list of prioritized L. major targets is presented in Tables 1 and 2, the list of 254 FDA-approved drugs (henceforth referred to as 'Lm254') is presented in Additional file 1: Table S2, and the list of all drugs associated with the 8 synthetic lethal targets is presented in Additional file 1: Table S3.
Target validation: Comparison of target predictions with observations from previous literature
As an initial validation, several of the genes in Table 1 were identified as proven drug targets in the literature and some have previously identified inhibitors in Leishmania or other related species. A particularly interesting target is LmjF05.0350 encoding for trypanothione reductase. In trypanosomatids, trypanothione reductase takes the place of glutathione reductase for the maintenance of intracellular redox balance . Trypanothione reductase is required for the survival of Leishmania sp. and one important physiological role of this enzyme is defense against oxidative damage when parasites are resident inside macrophages . Host glutathione reductase and parasite trypanothione reductase have mutually exclusive substrate specificities conferring an advantage for the purposes of drug targeting against the parasite enzyme . Many inhibitors have been identified for trypanothione reductase, including trivalent antimony (Sb (III)) ions as recently demonstrated . Also, two nitrofuran derivative compounds were experimentally shown to be effective non-competitive inhibitors of trypanothione reductase in Trypanosoma cruzi, a trypanosomatid organism and causative agent for Chagas disease . MetDP predicted the efficacy of two nitrofurans associated with LmjF05.0350 namely, nitrofurazone and nitrofurantoin (see Figure 3). Other interesting targets in the prioritized list include LmjF12.0280 and LmjF04.0580 encoding for ornithine decarboxylase and spermidine synthase, respectively. Similarly to trypanothione reductase, both genes are essential to Leishmania sp. and model predictions are consistent with experimental observations . Ornithine decarboxylase is also an important target in Trypanosoma brucei, another trypanosomatid organism and causative agent for human African trypanosomiasis, with a well-established inhibitor (eflornithine) . See also a discussion on false negative target predictions in Additional file 1.
Drug validation: Comparison of drug predictions to literature findings and existing high-throughput screens
As additional preliminary validation, the literature was mined to consider drugs that have been evaluated clinically against leishmaniasis (see Additional file 1: Table S4). There were nine positives in Lm254: amphotericin B [17–20], ketoconazole [17, 21], fluconazole [17, 22], clotrimazole , itraconazole [19, 24], miconazole , terbinafine , metronidazole [19, 26] and allopurinol [17–19]. Note that for some of these drugs, the literature is controversial with regards to their efficacy in treating leishmaniasis. From the computational network analysis, seven of the positive candidates are predicted to exclusively target L. major enzymes involved in steroid biosynthesis, with exceptions being ketoconazole and allopurinol. The anti-fungal ketoconazole is predicted to be associated with four L. major enzymes: three involved with steroid biosynthesis and one participating in fructose and mannose metabolism. And, allopurinol is associated with an enzyme involved in methionine metabolism (see Additional file 1: Table S4; see also a note on 'Dependency on confidence of interactions in DrugBank and STITCH' in Additional file 1). Previous research has identified the azole drugs ketoconazole, itraconazole, miconazole and fluconazole to be sterol 14-demethylase inhibitors that have been evaluated in vitro against T. cruzi parasites or in mouse models of T. cruzi infection . Additionally, in the fungal species Candida albicans, clotrimazole and terbinafine have been shown to be inhibitors of sterol 14-demethylase and squalene monooxygenase, respectively [28, 29]. These references from other organisms serve as a check for some of the computationally predicted links in MetDP between L. major genes and drugs.
There were also four false negative results: imiquimod , paromomycin [17, 19, 20], pentamidine [17–20] and sodium stibogluconate [17–21]. Imiquimod is a modulator of the innate immune response , and paromomycin is an aminoglycoside antibiotic that binds to the aminoacyl decoding site of the ribosomal 16S RNA [31, 32]. Since both imiquimod and paromomycin do not primarily target L. major proteins, their exclusion from Lm254 was obvious. However, the reasons for exclusion of all four drugs from Lm254 were examined in further detail (see Additional file 1). Hence, nine out of thirteen (69.2%) clinically relevant drugs for leishmaniasis were present in Lm254. Three other drugs used clinically against leishmaniasis (meglumine antimoniate [17–20], miltefosine [20, 33] and sitamaquine [19, 20, 34]) were not considered in the analysis as they were not present in DrugBank.
Additionally, several high-throughput drugs screens against in vitro growth of Leishmania and Trypanosoma species have been recently published [35–38]. In particular, one study by Sharlow et al. screened 196,146 compounds at 10 μM against L. major promastigotes . A total of 187 FDA-approved drugs that were also in DrugBank database overlapped with the set of 196,146 compounds; 68 out of 187 were present in Lm254. Of these 68 drugs, seven (pimozide, furazolidone, perphenazine, bifonazole, disulfiram, clotrimazole and floxuridine) were active as primary hits in the high-throughput screen, and the remaining 61 were classified as inactive as they did not meet the 50% inhibition threshold when evaluated at 10 μM. Furthermore, during assay optimization and validation screening in the same study, two additional drugs not included in the list of 68, tamoxifen and mycophenolic acid, were also deemed to be primary hits . These two drugs were also included in Lm254. A second high-throughput screening (HTS) study for inhibitors cytotoxic to bloodstream form T. brucei tested 2,160 FDA-approved drugs, bioactive compounds and natural products and produced 35 hits when assayed at a concentration of 1 μM . One drug, paclitaxel, which was included in the list of 35 hits was also present in Lm254. Paclitaxel was not included in the set of 68 drugs evaluated in the Sharlow et al. study. In previous literature, paclitaxel along with IFNγ was shown to induce killing of L. major infected murine macrophages . In summary, from comparison to HTS data, 10 out of 71 (14.1%) drugs from Lm254 that were experimentally evaluated via HTS were classified to be potential antileishmanial hits.
Taken together, the inclusion of all clinically relevant candidates and existing HTS data yields a total of 18 drugs (7.1% of Lm254) as having potential for antileishmanial activity (see Additional file 1: Table S2 for the list of 18 drugs). The HTS data and data on drugs that have been used clinically against leishmaniasis account for only 83 out of 254 drugs. Any drugs in Lm254 not previously evaluated against L. major parasites (in vitro and/or in vivo) serve as prime candidates for future investigation into the discovery of antileishmanial compounds.
In vitroexperimental evaluation of halofantrine
Originally, from the DrugBank database, halofantrine was associated with vacuolar ATP synthase catalytic subunit A from Plasmodium falciparum. Upon implementation of the MetDP pipeline, halofantrine was associated with the following synthetic lethal pairs: LmjF05.0500 & LmjF05.0510, genes encoding for F-type H+-transporting ATPase alpha chain, and LmjF25.1170 & LmjF25.1180, genes encoding F-type H+-transporting ATPase beta chain. While the exact mechanism of action for halofantrine in L. major requires further investigation, the model predictions for halofantrine represent testable hypotheses. Towards that end, an ATP bioluminescence assay was performed to further investigate the effect of halofantrine on ATP levels in L. major promastigotes.
As Figure 4D demonstrates, parasites were incubated with or without 10 μM halofantrine in various conditions as follows: (a) media, (b) HEPES-buffered saline (HBS) for 'No inhibition', (c) HBS with sodium azide for 'Inhibition of F1-ATPase and cytochrome c oxidase', (d) HBS without glucose, but supplemented with 2-deoxy-D-glucose and sodium pyruvate for 'Inhibition of glycolysis', and (e) HBS without glucose, but supplemented with sodium azide, 2-deoxy-D-glucose and sodium pyruvate for 'Inhibition of both'. In the assay, sodium azide was used to inhibit mitochondrial oxidative ATP generation . Further, 2-deoxy-D-glucose (a glucose analog and competing substrate for hexokinase) was used along with sodium pyruvate in glucose-free buffer to inhibit glycolytic ATP synthesis . In Figure 4E, when the parasites were subject to 'No inhibition,' halofantrine caused a significant reduction in ATP levels. Figure 4E also shows that ATP levels in parasites under 'Inhibition of F1-ATPase and cytochrome c oxidase' dropped significantly when halofantrine was added. Likewise, ATP levels in parasites under 'Inhibition of glycolysis' dropped significantly when halofantrine was added. In both these instances, ATP levels dropped to a level comparable to parasites subject to the simultaneous inhibition of glycolysis and F1-ATPase/cytochrome c oxidase ('Inhibition of both'). Hence, the ATP bioluminescence assay suggests that halofantrine has an effect on either one or both of glycolytic and mitochondrial oxidative ATP generation mechanisms. Future experimental efforts can be directed at more precise characterization of the underlying mechanism of action of the drug and elucidating other potential metabolic targets of halofantrine in L. major. Correspondingly, absorbance was monitored to make certain that the effects that were seen with regard to ATP levels were not a function of variations in cell count across the various conditions (Figure 4F). Finally, although halofantrine is known to be associated with cardiotoxicity , its potential as an antileishmanial agent should be investigated in more detail.
Model-guided drug combinations as an alternative strategy against leishmaniasis
With the ability to computationally simulate synthetically lethal gene deletions, there exists an opportunity to predict and prioritize multiple combination drug therapies that may be superadditive/synergistic. A potential advantage of superadditive drug combinations is overcoming toxicity or side-effects linked to high doses of individual drugs needed to establish the same inhibitory effect as the combination . One strategy of prioritizing for clinically-relevant drug combinations was to selectively focus on one drug that demonstrated excellent antileishmanial activity in vitro and was associated with low-toxicity. All combinations involving that particular drug could be prioritized for experimental analysis. Selection of such high-priority drug combinations through the power of network analysis can augment single compound discovery strategies.
Since disulfiram was previously shown to have a sub-micromolar IC50  (see also dose response data for disulfiram in Additional file 1: Figure S7), the goal was to identify superadditive drug combinations that involve disulfiram. Marketed as Antabuse, disulfiram has an LD50 of 8600 mg/kg (rat; oral) and is used in the treatment of alcoholism. At a moderate STITCH confidence of greater than 0.4, none of the 8 high-priority synthetic lethal gene predictions were associated with disulfiram. Therefore, the search criterion was expanded to include gene-drug interactions at a lower STITCH confidence (greater than 0.15; see Methods for mapping between L. major genes and drugs; see also Additional file 2 for a list of initial gene-drug associations).
Two of the combinations that demonstrated superadditivity involved antibiotics - kanamycin is an aminoglycoside antibiotic and amoxicillin is a common β-lactam antibiotic. The other two combinations involved antipsychotics - clozapine and chlorpromazine are used in the treatment of symptoms associated with schizophrenia. Of the four combinations, disulfiram + amoxicillin had the best overall toxicity rating. Moreover, only clozapine and chlorpromazine demonstrated inhibitory effects on L. major promastigotes when tested individually (see Additional file 1: Figures S9 and S11). Kanamycin and amoxicillin did not demonstrate any inhibitory effect at the concentrations tested (see Additional file 1: Figures S8 and S10).
Interestingly, along with disulfiram, kanamycin and chlorpromazine are also present in Lm254. In the future, evaluating two-drug combinations where both drugs are predicted to be effective individually by acting on single lethal or growth-reducing targets may have certain advantages. Based solely on synthetic lethality analysis, in the event that the parasite develops resistance to one of the drugs (in a combination) due to mutations in an enzyme/gene product encoded by a single gene in a synthetic lethal deletion, the overall drug combination would be predicted to be ineffective. However, by considering only those drug combinations where the individual drugs involved are also independently effective antileishmanial agents, this limitation can be potentially overcome. If the organism becomes resistant to one of the drugs during infection, there is a potential compensatory mechanism in place for the other drug to inhibit growth. This redundancy may serve to prevent or delay the onset of resistance .
Recent literature has described the applicability of metabolic network analysis towards target identification and drug discovery in general. One study (targetTB) used a variety of network analyses and bioinformatics-based sequence/structural assessments to predict a list of targets against Mycobacterium tuberculosis . By prioritizing targets using a layered approach, proteins that did not pass sequential cut-offs were filtered out . Another study implemented network analysis and molecular docking simulations to identify small molecule inhibitors of type II fatty-acid biosynthesis enzymes in Escherichia coli and Staphylococcus aureus and experimentally evaluated the computational predictions . FBA was also used to simulate drug synergy effects by assigning synergy scores to combinations of enzymes that would inhibit growth of E. coli in various environmental conditions . Network reconstructions of P. falciparum and Vibrio vulnificus have also guided the prediction of drug targets, and investigational compounds were experimentally evaluated against these pathogens [9, 46].
Here, a novel pipeline for the prioritization of drugs and drug targets is presented that is in many ways distinct from approaches previously developed. Instead of restricting the search to compounds that may only target non-human proteins or have computed Tanimoto similarity to important metabolites, MetDP undertook a more expansive strategy in making use of protein-compound interactions present in both DrugBank and STITCH resources. Importantly, a clear majority of the approved drugs target human proteins as opposed to proteins in any other single organism (e.g. in DrugBank). Therefore, rather than avoiding the selection of pathogen targets that are similar to human proteins, MetDP explicitly sought out such targets in order to create a link between L. major genes and drugs that are already approved. The side-effects caused by any of the drugs in the prioritized list that act on related human targets are relatively well known.
In addition, by restricting the search to only FDA-approved drugs and eliminating investigational drugs from the selection process, major regulatory hurdles that novel compounds are subject to before approval for widespread clinical use are bypassed. For instance, on average, it can take more than 15 years and over 800 million USD to bring a single novel drug to market . Moreover, only 20 to 30 new drugs are approved by the FDA annually . Instead, drugs already approved for clinical use have known pharmacokinetics and toxicity/safety profiles (many approved drugs have met superior safety standards with phase IV post-market safety surveillance) . Therefore, by demonstrating novel antileishmanial activity for a drug already in clinical use for another indication, a candidate drug or perhaps combinations of drugs can begin to be effectively repurposed by being evaluated in phase II clinical trials and bypassing approximately 40% of the overall cost that it takes to bring a novel drug to market . Besides, an FDA-approved drug with demonstrated antileishmanial activity can also be prescribed off-label, thus making it accessible to patients urgently in need of an alternative treatment option.
Ultimately, metabolic reconstructions provide a framework for the interrogation of human pathogens and serve as a platform for generation of future experimental hypotheses. This metabolic network-driven approach identified 15 L. major genes as high-priority targets, 8 high-priority synthetic lethal targets, and 254 FDA-approved drugs as potential antileishmanial agents. As experimental validation, the antimalarial drug halofantrine was shown to have noticeable antileishmanial activity. Moreover, through the ATP bioluminescence assay, halofantrine affected either one or both of glycolytic and mitochondrial ATP generation mechanisms. Additionally, synthetic lethal predictions from the metabolic network aided in the selection of drug combinations with potential for superadditivity. For proof-of-concept, double-drug combinations were evaluated in vitro against L. major and four novel superadditive combinations involving the drug disulfiram were discovered. Selection of such high-priority double-drug combinations guided by metabolic network analysis might provide for an attractive and alternative avenue for drug discovery against infectious diseases like leishmaniasis.
A rational method of prioritizing drugs and drug targets with a metabolic network as an underlying framework also provides for a possible mechanism of action for many of the drugs selected (e.g., the network model enables the analysis of what biosynthetic pathways are inhibited when a particular enzyme's function is pharmacologically inhibited). Surely the mechanism of action is only predicted and needs to be validated with follow-on experimentation. However, given a drug and potential target, a starting experimental hypothesis is provided. A network-driven approach guiding the selection of prioritized drugs can be immensely advantageous in terms of cost and efficiency in the beginning phases of drug discovery and can offer significant implications to future drug repurposing strategies against a variety of NTDs.
The previously published metabolic reconstruction of L. major was the first of its kind for a protozoan. It accounted for 560 genes, 1112 reactions, and 1101 metabolites spanning eight unique sub-cellular localizations .
The DrugBank knowledgebase (version 2.5; http://www.drugbank.ca/) included 4774 drugs that were classified into several categories namely, (a) small molecule, (b) biotech, (c) approved, (d) experimental, (e) nutraceutical, (f) illicit, and (g) withdrawn. These drugs were linked to 4554 peptide/protein and non-peptide/non-protein drug targets.
The STITCH database (version 2.0; http://stitch.embl.de/) comprised over 74000 small molecules and 2.5 million+ proteins in 630 organisms. The STITCH database includes many potentially important off-targets to particular compounds (derived from experimental data, other curated databases or text mining) in addition to the primary known targets, and all protein-compound interactions in STITCH have an associated confidence metric. It should be noted that L. major was not present in the list of organisms whose targets are explicitly included in STITCH.
Druggability indices for L. major genes were downloaded from the TDR Targets database (version 3; http://tdrtargets.org/). The druggability index was calculated using a combination of approaches including sequence similarity to known biological targets associated with FDA-approved drugs from curated databases (e.g. DrugStore and StARLITe), and a sequence-based Bayesian learning algorithm, among other methods . The index ranges between 0 (not druggable) and 1 (highly druggable). In all, 1598 L. major genes have a corresponding druggability index in the TDR Targets database. Of these, 261 are accounted for in the metabolic reconstruction.
DEA controlled substances
The list of controlled substances by Schedule was downloaded from the Drug Enforcement Administration (DEA) Office of Diversion Control database (http://www.deadiversion.usdoj.gov/schedules).
Tools and programs used
A BLAST+ v.2.2.24 executable was downloaded from the NCBI FTP site (ftp://ftp.ncbi.nih.gov/blast/executables/blast+/). The program makeblastdb was implemented on both DrugBank and STITCH target protein sequences. The amino acid sequences for all 560 L. major genes from the metabolic reconstruction were obtained from GeneDB (version 2.1; http://www.genedb.org/genedb/leish). Subsequently, the program blastp was used to query 560 L. major sequences against both DrugBank and STITCH target databases.
Flux Balance Analysis (FBA) and Flux Variability Analysis (FVA) were implemented using the COBRA toolbox (version 1.3.3; http://opencobra.sourceforge.net/openCOBRA/Welcome.html) in MATLAB. FBA and FVA have been extensively implemented on metabolic networks of various organisms ranging from prokaryotes to eukaryotes .
Mapping between L. majorgenes and drugs in DrugBank and STITCH
First, a BLASTP search was performed with 560 L. major genes from the metabolic reconstruction (query sequences) against 4538 target proteins in DrugBank (target sequences). An E-value cutoff was set at 0.001, and repetitious links between L. major genes and DrugBank target proteins were removed. At this first cutoff, 440 L. major genes were linked to 1116 DrugBank target proteins and transitively to 1129 drugs. The list of drugs was refined to focus on only FDA-approved drugs by removing experimental and withdrawn drugs, which resulted in 313 L. major genes being linked to 257 FDA-approved drugs.
In parallel, a similar approach with some variations was undertaken with the STITCH database. In order to be as expansive as possible in making initial associations between L. major genes and drugs, both DrugBank and STITCH resources were used. A BLASTP search was performed between the 560 L. major genes and 2.5 million+ proteins (2,590,259) to associate L. major genes with the proteins in STITCH. Protein sequences were derived from a related STRING database (http://string-db.org/). Again, an E-value cutoff of 0.001 was chosen and repetitious associations between L. major genes and STITCH proteins were removed. At this cutoff, all 560 genes were associated with 131,250 protein targets from STITCH.
Subsequently, the STITCH compound-protein interaction dataset was refined to only consider compounds that were FDA-approved drugs (as classified in the DrugBank database). PubChem compound identifiers were used to map drugs in DrugBank with compounds in STITCH. Only 1337 of the 1466 approved drugs in DrugBank had readily available PubChem compound identifiers (and a few had similar PubChem IDs). In all, 1330 drugs with unique DrugBank IDs mapped to 1305 unique compounds in the STITCH dataset. Next, the STITCH database was further reduced to consider protein targets strictly present in the list of 131,250 unique protein targets that made the BLAST cut. Finally, the STITCH dataset was also constrained to only include medium confidence and higher interactions (STITCH scores greater than 0.4). The newly constrained STITCH dataset along with the BLAST results were used to link 538 L. major genes to 905 FDA-approved drugs. At this stage, results from DrugBank and STITCH were merged for further refinement and analysis. The merged list linked 538 L. major genes to 926 FDA-approved drugs (see Additional file 2).
Generating an FVA score for genes
Flux variability analysis (FVA) allows for ascertaining the complete range of numerical values for every flux in the biochemical network given an optimal value of the objective . FBA computes only one possible optimal reaction network state. Therefore, the flux distribution resulting from FBA (i.e. the enumeration of all possible fluxes in the network) is only one of many feasible solutions that can exist for the same optimal value of the cellular objective. Enzymatic reactions associated with small ranges of flux values may not necessarily be robust to various kinds of network perturbations. Therefore, an FVA score for genes was developed based on an assumption that genes linked to reactions with little flux variability might be good targets to disrupt in the metabolic network.
Step 1: Conferring weights on reactions based on FVA
This scoring method gives a high score to reactions with little variability in flux. For example, reactions that spanned the entire flux range [-999999, 999999] would receive a score of 0, while reactions with no change in flux would receive a score of 1.
Step 2: Evaluating the contributions of genes in gene-protein-reaction (GPR) relationships
N states refers to the number of possible states for a given GPR, which is calculated as 2n with n being the number of genes in a GPR. The GPR was evaluated to be TRUE or FALSE for every state. For a given reaction k with a GPR, the presence (P) or absence (A) of a particular gene i and the TRUE (T) or FALSE (F) outcome of the GPR were tabulated as matches (PT i and AF i ) and mis-matches (PF i and AT i ). Explicitly, for a given reaction k, 'PTi' refers to the number of gene states (out of 2n) where a particular gene i is present, and the GPR is evaluated to be TRUE. 'ATi' refers to the number of gene states where gene i is absent, and the GPR is evaluated to be TRUE (e.g. due to the presence of isozymes). 'PFi' refers to the number of gene states where a particular gene i is present, and the GPR is evaluated to be FALSE (e.g. subunits of a larger protein complex that need to be present together). Finally, 'AFi' refers to the number of gene states where a particular gene i is absent, and the GPR is evaluated to be FALSE. The weights (cg,r) were calculated by totaling up matches (counted positively) and mis-matches (counted negatively) and dividing by the number of possible states for a given GPR. Hence, a gene-reaction matrix (C) was constructed with elements of the matrix representing weights for genes linked to corresponding reactions.
Step 3: Conferring weights on genes based on FVA
Next, the FVA reaction scores (ri) were multiplied in with the gene weights from the GPR (cg,r) to yield an FVA-gene matrix with the same dimensions as C. The maximum value for genes across all reactions was considered as the FVA score for a gene (value resided between 0 and 1).
Calculations of additive results
Here, E1±δE1 refers to the effect of drug 1 and E2±δE2 refers to the effect of drug 2.
Here, and refer to the means of the samples, s1 and s2 refer to the standard deviations of the samples, and n1 and n2 refer to the number of samples. One-tail p-value was subsequently determined by referring to a t-table (degrees of freedom computed as follows: n1 + n2 - 2).
Black flat-bottom 96-well microtiter plates were purchased from Fisher Scientific (http://www.fishersci.com) and used in all alamarBlue experiments. White flat-bottom 96-well microtiter plates also purchased from Fisher Scientific were used in all bioluminescence experiments. alamarBlue was purchased from Invitrogen (http://www.invitrogen.com). CellTiter-Glo was purchased from Promega (http://www.promega.com). All compounds used in this study were purchased from Sigma-Aldrich (http://www.sigmaaldrich.com/). Compounds were solubilized in dimethyl sulfoxide (DMSO) or water.
Previously published protocols on culturing L. major  were adhered to in this study. L. major promastigotes and protocol for preparing media were kindly provided by Mary E. Wilson and Melissa A. Miller, University of Iowa. Parasites in complete HOMEM (see Additional file 1) were cultured in 25 cm2 plastic tissue culture flasks with sealed or vented caps and maintained at 26°C.
The assay was conducted in accordance with previously established protocols [51–53]. Briefly, promastigotes were diluted to 1 × 106 cells/mL, and in a black flat-bottom 96-well microtiter plate, 180 μL of suspension was incubated with varying concentrations of drugs (singly or in combination) in triplicate. Specifically, 160 μL of parasite samples were first seeded in triplicate. Next, sample wells were topped off with 20 μL of media + drug(s) (ratio altered to achieve specific concentrations of drug(s)) such that the total volume equaled 180 μL. Heat-killed parasite samples (incubated at 60°C for 20 minutes) prepared at 1 × 106 cells/mL were also seeded in triplicate (160 μL of sample + 20 μL of media) to serve as a positive control. Amphotericin B at 1 μM also served as another positive control. If DMSO was used to solubilize the drug(s), three wells with the highest relevant concentration of DMSO were included in the plate as a negative control. Additionally, three wells were seeded with 180 μL of media alone. The plate was incubated at 26°C for 24 hours at which time point 20 μL of alamarBlue dye was added to all control and experimental wells. Using a Gemini EM Microplate Spectrofluorometer, fluorescence was monitored at excitation/emission wavelengths of 544 nm/590 nm at 24 and 48 hours post addition of dye to wells. Calibration data for alamarBlue assay is provided in Additional file 1: Figures S4, S5 and S6.
The protocol for the bioluminescence assay was modified from . Parasites at 8 × 106 cells/mL were incubated in culture medium or various buffers for 2 hours at 26°C either alone or in the presence of 10 μM halofantrine. Mitochondrial oxidative ATP generation was inhibited by incubating the parasites in HBS buffer with glucose plus 20 mM sodium azide, an inhibitor of F1-ATPase and cytochrome c oxidase from complex IV . Glycolytic ATP generation was inhibited by incubating the parasites in glucose-free HBS buffer plus 5 mM 2-deoxy-D-glucose, a competitor with glucose for hexokinase binding, and 5 mM sodium pyruvate . In a white opaque flat-bottom 96-well microtiter plate, 25 μL of parasite samples from each condition were seeded in triplicate. Heat-killed parasite samples (incubated at 60°C for at least 20 minutes) prepared at 8 × 106 cells/mL were also seeded in triplicate. Additionally, three wells were seeded with 25 μL of media alone. Subsequently, 25 μL of CellTiter-Glo was added to all control and experimental wells. The plate was incubated in the dark at 26°C for 10 minutes. Luminescence was monitored using a FLUOstar Optima plate reader (BMG Labtech). For absorbance measurements, 100 μL of control and experimental samples were seeded in triplicate at the 18 hour time point. The plate was immediately transferred to a Tecan infinite200 Pro microplate reader, and absorbance was monitored at 600 nm. Calibration data for the bioluminescence assay is provided in Additional file 1: Figures S12 and S13.
AKC gratefully acknowledges and expresses appreciation to Kevin M. D'Auria for the many insightful and invaluable discussions during the course of the project. In addition, all authors wish to thank Erwin P. Gianchandani, Joseph Klembczyk and Matthew Oberhardt for their respective contributions. The authors acknowledge financial support from the National Science Foundation (CAREER program grant no. 0643548 to JP and GRFP to PAJ), the National Institutes of Health (grant no. GM088244 to JP), and Double Hoo research grant from the University of Virginia (awarded to AKC and Joseph Klembczyk).
- Hotez PJ, Fenwick A, Savioli L, Molyneux DH: Rescuing the bottom billion through control of neglected tropical diseases. Lancet. 2009, 373 (9674): 1570-1575. 10.1016/S0140-6736(09)60233-6.View ArticleGoogle Scholar
- Hotez PJ: Mass drug administration and integrated control for the world's high-prevalence neglected tropical diseases. Clin Pharmacol Ther. 2009, 85 (6): 659-664. 10.1038/clpt.2009.16.View ArticleGoogle Scholar
- Desjeux P: Leishmaniasis: current situation and new perspectives. Comp Immunol Microbiol Infect Dis. 2004, 27 (5): 305-318. 10.1016/j.cimid.2004.03.004.View ArticleGoogle Scholar
- Croft SL, Coombs GH: Leishmaniasis-current chemotherapy and recent advances in the search for novel drugs. Trends Parasitol. 2003, 19 (11): 502-508. 10.1016/j.pt.2003.09.008.View ArticleGoogle Scholar
- Molyneux D, Killick-Kendrick R: Morphology, ultrastructure and life cycles. The Leishmaniases in Biology and Medicine. Edited by: Peters W, Killick-Kendrick R. 1987, London: Academic Press Inc, 1: 121-176.Google Scholar
- Hotez PJ, Molyneux DH, Fenwick A, Savioli L, Takeuchi T: A Global Fund to Fight Neglected Tropical Diseases: is the G8 Hokkaido Toyako 2008 Summit ready?. PLoS Negl Trop Dis. 2008, 2 (3): e220-10.1371/journal.pntd.0000220.View ArticleGoogle Scholar
- Caffrey CR, Steverding D: Recent initiatives and strategies to developing new drugs for tropical parasitic diseases. Expert Opin on Drug Discov. 2008, 3 (2): 173-186. 10.1517/174604188.8.131.52.View ArticleGoogle Scholar
- Guerin PJ, Olliaro P, Sundar S, Boelaert M, Croft SL, Desjeux P, Wasunna MK, Bryceson AD: Visceral leishmaniasis: current status of control, diagnosis, and treatment, and a proposed research and development agenda. Lancet Infect Dis. 2002, 2 (8): 494-501. 10.1016/S1473-3099(02)00347-X.View ArticleGoogle Scholar
- Plata G, Hsiao TL, Olszewski KL, Llinas M, Vitkup D: Reconstruction and flux-balance analysis of the Plasmodium falciparum metabolic network. Mol Syst Biol. 2010, 6: 408-View ArticleGoogle Scholar
- Jamshidi N, Palsson BO: Investigating the metabolic capabilities of Mycobacterium tuberculosis H37Rv using the in silico strain iNJ661 and proposing alternative drug targets. BMC Syst Biol. 2007, 1: 26-10.1186/1752-0509-1-26.View ArticleGoogle Scholar
- Chavali AK, Whittemore JD, Eddy JA, Williams KT, Papin JA: Systems analysis of metabolism in the pathogenic trypanosomatid Leishmania major. Mol Syst Biol. 2008, 4: 177-View ArticleGoogle Scholar
- Oberhardt MA, Palsson BO, Papin JA: Applications of genome-scale metabolic reconstructions. Mol Syst Biol. 2009, 5: 320-View ArticleGoogle Scholar
- Dumas C, Ouellette M, Tovar J, Cunningham ML, Fairlamb AH, Tamar S, Olivier M, Papadopoulou B: Disruption of the trypanothione reductase gene of Leishmania decreases its ability to survive oxidative stress in macrophages. EMBO J. 1997, 16 (10): 2590-2598. 10.1093/emboj/16.10.2590.View ArticleGoogle Scholar
- Baiocco P, Colotti G, Franceschini S, Ilari A: Molecular basis of antimony treatment in leishmaniasis. J Med Chem. 2009, 52 (8): 2603-2612. 10.1021/jm900185q.View ArticleGoogle Scholar
- Jockers-Scherubl MC, Schirmer RH, Krauth-Siegel RL: Trypanothione reductase from Trypanosoma cruzi. Catalytic properties of the enzyme and inhibition studies with trypanocidal compounds. Eur J Biochem. 1989, 180 (2): 267-272. 10.1111/j.1432-1033.1989.tb14643.x.View ArticleGoogle Scholar
- Heby O, Persson L, Rentala M: Targeting the polyamine biosynthetic enzymes: a promising approach to therapy of African sleeping sickness, Chagas' disease, and leishmaniasis. Amino Acids. 2007, 33 (2): 359-366. 10.1007/s00726-007-0537-9.View ArticleGoogle Scholar
- Blum J, Desjeux P, Schwartz E, Beck B, Hatz C: Treatment of cutaneous leishmaniasis among travellers. J Antimicrob Chemother. 2004, 53 (2): 158-166. 10.1093/jac/dkh058.View ArticleGoogle Scholar
- Olliaro PL, Bryceson AD: Practical progress and new drugs for changing patterns of leishmaniasis. Parasitol Today. 1993, 9 (9): 323-328. 10.1016/0169-4758(93)90231-4.View ArticleGoogle Scholar
- Singh S, Sivakumar R: Challenges and new discoveries in the treatment of leishmaniasis. J Infect Chemother. 2004, 10 (6): 307-315. 10.1007/s10156-004-0348-9.View ArticleGoogle Scholar
- Mishra J, Saxena A, Singh S: Chemotherapy of leishmaniasis: past, present and future. Curr Med Chem. 2007, 14 (10): 1153-1169. 10.2174/092986707780362862.View ArticleGoogle Scholar
- Navin TR, Arana BA, Arana FE, Berman JD, Chajon JF: Placebo-controlled clinical trial of sodium stibogluconate (Pentostam) versus ketoconazole for treating cutaneous leishmaniasis in Guatemala. J Infect Dis. 1992, 165 (3): 528-534. 10.1093/infdis/165.3.528.View ArticleGoogle Scholar
- Alrajhi AA, Ibrahim EA, De Vol EB, Khairat M, Faris RM, Maguire JH: Fluconazole for the treatment of cutaneous leishmaniasis caused by Leishmania major. N Engl J Med. 2002, 346 (12): 891-895. 10.1056/NEJMoa011882.View ArticleGoogle Scholar
- Larbi EB, al-Khawajah A, al-Gindan Y, Jain S, Abahusain A, al-Zayer A: A randomized, double-blind, clinical trial of topical clotrimazole versus miconazole for treatment of cutaneous leishmaniasis in the eastern province of Saudi Arabia. AmJTrop Med Hyg. 1995, 52 (2): 166-168.Google Scholar
- Dogra J, Saxena VN: Itraconazole and leishmaniasis: a randomised double-blind trial in cutaneous disease. Int J Parasitol. 1996, 26 (12): 1413-1415. 10.1016/S0020-7519(96)00128-2.View ArticleGoogle Scholar
- Bahamdan KA, Tallab TM, Johargi H, Nourad MM, Ibrahim K, el Sherbini AH, Karkashan E, Khare AK, Nauri MM: Terbinafine in the treatment of cutaneous leishmaniasis: a pilot study. Int J Dermatol. 1997, 36 (1): 59-60. 10.1046/j.1365-4362.1997.00021.x.View ArticleGoogle Scholar
- Mishra M, Thakur BD, Choudhary M: Metronidazole and Indian kala-azar: results of a clinical trial. Br Med J (Clin Res Ed). 1985, 291 (6509): 1611-10.1136/bmj.291.6509.1611.View ArticleGoogle Scholar
- Buckner FS: Sterol 14-demethylase inhibitors for Trypanosoma cruzi infections. Adv Exp Med Biol. 2008, 625: 61-80. 10.1007/978-0-387-77570-8_6.View ArticleGoogle Scholar
- Ryder NS: Terbinafine: mode of action and properties of the squalene epoxidase inhibition. Br J Dermatol. 1992, 126 (Suppl 39): 2-7.View ArticleGoogle Scholar
- Marriott MS: Inhibition of sterol biosynthesis in Candida albicans by imidazole-containing antifungals. J Gen Microbiol. 1980, 117 (1): 253-255.Google Scholar
- Stanley MA: Imiquimod and the imidazoquinolones: mechanism of action and therapeutic potential. Clin Exp Dermatol. 2002, 27 (7): 571-577. 10.1046/j.1365-2230.2002.01151.x.View ArticleGoogle Scholar
- Vicens Q, Westhof E: Crystal structure of paromomycin docked into the eubacterial ribosomal decoding A site. Structure. 2001, 9 (8): 647-658. 10.1016/S0969-2126(01)00629-3.View ArticleGoogle Scholar
- Fernandez MM, Malchiodi EL, Algranati ID: Differential effects of paromomycin on ribosomes of Leishmania mexicana and mammalian cells. Antimicrob Agents Chemother. 2011, 55 (1): 86-93. 10.1128/AAC.00506-10.View ArticleGoogle Scholar
- Garnier T, Croft SL: Topical treatment for cutaneous leishmaniasis. Curr Opin Investig Drugs. 2002, 3 (4): 538-544.Google Scholar
- Wasunna MK, Rashid JR, Mbui J, Kirigi G, Kinoti D, Lodenyo H, Felton JM, Sabin AJ, Albert MJ, Horton J: A phase II dose-increasing study of sitamaquine for the treatment of visceral leishmaniasis in Kenya. AmJTrop Med Hyg. 2005, 73 (5): 871-876.Google Scholar
- Siqueira-Neto JL, Song OR, Oh H, Sohn JH, Yang G, Nam J, Jang J, Cechetto J, Lee CB, Moon S, et al: Antileishmanial high-throughput drug screening reveals drug candidates with new scaffolds. PLoS Negl Trop Dis. 2010, 4 (5): e675-10.1371/journal.pntd.0000675.View ArticleGoogle Scholar
- St George S, Bishop JV, Titus RG, Selitrennikoff CP: Novel compounds active against Leishmania major. Antimicrob Agents Chemother. 2006, 50 (2): 474-479. 10.1128/AAC.50.2.474-479.2006.View ArticleGoogle Scholar
- Sharlow ER, Close D, Shun T, Leimgruber S, Reed R, Mustata G, Wipf P, Johnson J, O'Neil M, Grogl M, et al: Identification of potent chemotypes targeting Leishmania major using a high-throughput, low-stringency, computationally enhanced, small molecule screen. PLoS Negl Trop Dis. 2009, 3 (11): e540-10.1371/journal.pntd.0000540.View ArticleGoogle Scholar
- Mackey ZB, Baca AM, Mallari JP, Apsel B, Shelat A, Hansell EJ, Chiang PK, Wolff B, Guy KR, Williams J, et al: Discovery of trypanocidal compounds by whole cell HTS of Trypanosoma brucei. Chem Biol Drug Des. 2006, 67 (5): 355-363. 10.1111/j.1747-0285.2006.00389.x.View ArticleGoogle Scholar
- Doherty TM, Sher A, Vogel SN: Paclitaxel (Taxol)-induced killing of Leishmania major in murine macrophages. Infect Immun. 1998, 66 (9): 4553-4556.Google Scholar
- Manzano JI, Carvalho L, Perez-Victoria JM, Castanys S, Gamarro F: Increased glycolytic ATP synthesis is associated with tafenoquine resistance in Leishmania major. Antimicrob Agents Chemother. 2011, 55 (3): 1045-1052. 10.1128/AAC.01545-10.View ArticleGoogle Scholar
- Bouchaud O, Imbert P, Touze JE, Dodoo AN, Danis M, Legros F: Fatal cardiotoxicity related to halofantrine: a review based on a worldwide safety data base. Malar J. 2009, 8: 289-10.1186/1475-2875-8-289.View ArticleGoogle Scholar
- Lehar J, Krueger AS, Avery W, Heilbut AM, Johansen LM, Price ER, Rickles RJ, Short GF, Staunton JE, Jin X, et al: Synergistic drug combinations tend to improve therapeutically relevant selectivity. Nat Biotechnol. 2009, 27 (7): 659-666. 10.1038/nbt.1549.View ArticleGoogle Scholar
- WHO: Guidelines for the treatment of malaria (2nd ed.). 2010, Geneva, Switzerland: WHO Press, World Health Organization, 2Google Scholar
- Raman K, Yeturu K, Chandra N: TargetTB: a target identification pipeline for Mycobacterium tuberculosis through an interactome, reactome and genome-scale structural analysis. BMC Syst Biol. 2008, 2: 109-10.1186/1752-0509-2-109.View ArticleGoogle Scholar
- Shen Y, Liu J, Estiu G, Isin B, Ahn YY, Lee DS, Barabasi AL, Kapatral V, Wiest O, Oltvai ZN: Blueprint for antimicrobial hit discovery targeting metabolic networks. Proc Natl Acad Sci USA. 2010, 107 (3): 1082-1087. 10.1073/pnas.0909181107.View ArticleGoogle Scholar
- Kim HU, Kim SY, Jeong H, Kim TY, Kim JJ, Choy HE, Yi KY, Rhee JH, Lee SY: Integrative genome-scale metabolic analysis of Vibrio vulnificus for drug targeting and discovery. Mol Syst Biol. 2011, 7: 460-View ArticleGoogle Scholar
- Chong CR, Sullivan DJ: New uses for old drugs. Nature. 2007, 448 (7154): 645-646. 10.1038/448645a.View ArticleGoogle Scholar
- Aguero F, Al-Lazikani B, Aslett M, Berriman M, Buckner FS, Campbell RK, Carmona S, Carruthers IM, Chan AW, Chen F, et al: Genomic-scale prioritization of drug targets: the TDR Targets database. Nat Rev Drug Discov. 2008, 7 (11): 900-907. 10.1038/nrd2684.View ArticleGoogle Scholar
- Mahadevan R, Schilling CH: The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab Eng. 2003, 5 (4): 264-276. 10.1016/j.ymben.2003.09.002.View ArticleGoogle Scholar
- Goldoni M, Johansson C: A mathematical approach to study combined effects of toxicants in vitro: evaluation of the Bliss independence criterion and the Loewe additivity model. Toxicol In Vitro. 2007, 21 (5): 759-769. 10.1016/j.tiv.2007.03.003.View ArticleGoogle Scholar
- de Oliveira-Silva F, de Morais-Teixeira E, Rabello A: Antileishmanial activity of azithromycin against Leishmania (Leishmania) amazonensis, Leishmania (Viannia) braziliensis, and Leishmania (Leishmania) chagasi. AmJTrop Med Hyg. 2008, 78 (5): 745-749.Google Scholar
- Mikus J, Steverding D: A simple colorimetric method to screen drug cytotoxicity against Leishmania using the dye Alamar Blue. Parasitol Int. 2000, 48 (3): 265-269. 10.1016/S1383-5769(99)00020-3.View ArticleGoogle Scholar
- Shimony O, Jaffe CL: Rapid fluorescent assay for screening drugs on Leishmania amastigotes. J Microbiol Methods. 2008, 75 (2): 196-200. 10.1016/j.mimet.2008.05.026.View ArticleGoogle Scholar
- Urbina JA, Concepcion JL, Rangel S, Visbal G, Lira R: Squalene synthase as a chemotherapeutic target in Trypanosoma cruzi and Leishmania mexicana. Mol Biochem Parasitol. 2002, 125 (1-2): 35-45. 10.1016/S0166-6851(02)00206-2.View ArticleGoogle Scholar