Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Highly Accessed Methodology article

Network-based assessment of the selectivity of metabolic drug targets in Plasmodium falciparum with respect to human liver metabolism

Susanna Bazzani*, Andreas Hoppe* and Hermann-Georg Holzhütter

Author affiliations

Institut für Biochemie, Charite-Universitätsmedizin, Reinickendorfer Str. 61, Haus 10, 13347 Berlin

For all author emails, please log on.

Citation and License

BMC Systems Biology 2012, 6:118  doi:10.1186/1752-0509-6-118


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/6/118


Received:28 March 2012
Accepted:23 July 2012
Published:31 August 2012

© 2012 Bazzani et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The search for new drug targets for antibiotics against Plasmodium falciparum, a major cause of human deaths, is a pressing scientific issue, as multiple resistance strains spread rapidly. Metabolic network-based analyses may help to identify those parasite’s essential enzymes whose homologous counterparts in the human host cells are either absent, non-essential or relatively less essential.

Results

Using the well-curated metabolic networks PlasmoNet of the parasite Plasmodium falciparum and HepatoNet1 of the human hepatocyte, the selectivity of 48 experimental antimalarial drug targets was analyzed. Applying in silico gene deletions, 24 of these drug targets were found to be perfectly selective, in that they were essential for the parasite but non-essential for the human cell. The selectivity of a subset of enzymes, that were essential in both models, was evaluated with the reduced fitness concept. It was, then, possible to quantify the reduction in functional fitness of the two networks under the progressive inhibition of the same enzymatic activity. Overall, this in silico analysis provided a selectivity ranking that was in line with numerous in vivo and in vitro observations.

Conclusions

Genome-scale models can be useful to depict and quantify the effects of enzymatic inhibitions on the impaired production of biomass components. From the perspective of a host-pathogen metabolic interaction, an estimation of the drug targets-induced consequences can be beneficial for the development of a selective anti-parasitic drug.

Keywords:
Plasmodium falciparum; Human hepatocyte; Drug targets; Drug selectivity; Genome-scale networks; Reduced fitness; Flux balance analysis

Background

An ideal drug should tackle the disease-causing processes in the most selective way, i.e. with no harm for the healthy cells. To our knowledge, such a perfectly selective drug does not exist. In fact, administrated chemicals can trigger at least secondary consequences (off-target effects) in the organism or in host cells [1]. In case of anti-parasitic drugs, side-effects on the metabolism of host cells may arise from the binding to homologous proteins that share a long evolutionary history with the parasite [2]. To minimize such side-effects, analyses on the performance of multiple networks and the consequences of enzymatic homologous inhibitions may be useful. Ideally, this analysis should be based on detailed mechanistic models of the metabolic networks of the parasite and the host cell pathways [3,4]. However, all the kinetic information needed for the creation of such models is rarely available. Therefore, alternative modeling frameworks, such as Monte Carlo simulations [5] and flux balance analysis (FBA) [6], have been developed for the simulation of the cellular metabolism. These methods allow the simulation of genome-scale networks, do not require kinetic enzymatic knowledge and are suited to test the outcome of gene deletions [6-8].

Genome-scale metabolic models of pathogens are extensively exploited to predict putative drug targets with FBA frameworks [9-11]. Furthermore, recent network-based analyses integrate host and pathogen models to describe more accurately the metabolic interactions and to improve the search of putative drug targets. For example, the topology of automatically inferred networks of Plasmodium falciparum and its human host are studied to identify essential enzymes [12]. The metabolic network of Mycobacterium tubercolosis is integrated with a human alveolar macrophage to describe three degrees of infection (latent, pulmonary, and meningeal tuberculosis) [13]. Similarly, a metabolic model of Plasmodium falciparum is embedded in a red blood cell model to simulate its intra-erythrocytic developmental stage [14]. Furthermore, the selectivity of enzymatic drug targets is already extensively predicted with large-scale metabolic networks of human cancer cells [15]. Thus, it is reasonable to assess the selectivity of enzymatic drug targets in host-pathogen metabolic interactions with genome-scale networks. The aim of this study was the prediction of selective enzymatic inhibitions with genome-scale networks of Plasmodium falciparum (PlasmoNet) and the human hepatocyte (HepatoNet1) [14,16]. Although the intra-erythrocytic plasmodial phase is still extensively investigated, the intra-hepatic phase of development is the first infection site and, thus, a promising stage of treatment [17]. To predict feasible metabolic phenotypes, for each model a set of metabolic objectives was assembled with extensive literature search. This set described the cellular composition that is specific to the modeled cell type. Then, computational methods were applied to predict the selectivity of antimalarials in absolute and in relative sense. To predict the enzymatic inhibitions in absolute terms, in silico gene deletions of homologous enzymes were performed. Alternatively, the concept of reduced fitness was applied to homologous essential enzymes to assess inhibitions in relative terms. While in silico gene deletions predicted scenarios where the enzyme is fully disrupted (e.g. with gene-excision experiments), the second approach better described the possibility of residual enzymatic activity. The latter is more appropriate for the administration of enzymatic inhibitors that only gradually reduce the enzyme activity.

The set of experimental antimalarial targets, that were tested in this work, were a merged and pruned list of published “gold standards”. The merged list contained 96 enzymes, out of which 48 were selected. These enzymes are known to be essential for Plasmodium falciparum metabolic homeostasis and growth [14,18,19]. To test the selectivity of these enzymes, in silico gene deletions of the enzymatic set were simulated. Twenty-four enzymes were found to be essential only for the parasite and did not compromise the host metabolic performance. Twelve enzymes were essential for both models and were investigated with the concept of reduced fitness [20], to identify which network could be more sensitive to an enzymatic impairment. Then, the fitness-based selectivity score predicted that 10 enzymes out of this last set were more sensitive for the host, as their in silico perturbations had a larger impact on the achievement of host metabolic objectives (Table 1).

Table 1. Outcome of the simulated drug target inhibitions

Results and discussion

The gene-deletion approach divided the set of drug targets in three groups (Table 1): (i) enzymes that are essential in both networks, (ii) enzymes that are essential only in PlasmoNet and (iii) enzymes that are dispensable in both networks. The largest set of experimental drug targets (24 out of 48) was not lethal for the host (due to alternative paths that bypassed the deleted reactions) and essential for the parasite network (at least one biomass flux was not achieved). These enzymes are, thus, predicted to be perfectly selective with respect to the Plasmodium falciparum network. Intriguingly, 12 drug targets were identified to be dispensable in both models. When the simulations included the experimental in vitro conditions that validated the drug targets, 9/12 enzymes turned out to be essential. This situation was defined as conditional essentiality, as the enzyme becomes essential when some substrates are not available in the medium. The remaining 3 glycolytic enzymes were, instead, predicted to be non-essential. This could be due to off-target effects, to extra-metabolic functions or to possible assembly of multi-enzymatic complexes.

Homologous drug targets essential in both networks

The gene-deletion approach predicted that 12 enzymes are essential in both networks (Table 1), thus indicating a possible metabolic hepatic impairment caused by the inhibition of homologous enzymes.

While this enzymatic set is already experimentally tested, a literature search was done to find any evidence of human cellular toxicity.

It was predicted that only acyl-CoA syntethase may be a selective target for antimalarial treatment (selectivity score=1.97). The essentiality of this enzyme is experimentally demonstrated in human lymphoblast-like cells [21].

The fitness approach predicted that thymidylate synthase was equally detrimental for both models (selectivity score=1). This enzyme is a possible therapeutic drug target that is efficient against a wide spectrum of human tumors [22,23] but one main drawback is the possibility of its enzymatic over-expression that may confer drug resistance.

Orotidine-5-phosphate decarboxylase and orotate phosphoribosyltransferase form a multienzymatic complex, but for this analysis the two enzymatic activities were assessed separately. Nevertheless, the results obtained with the fitness approach indicated the same degrees of biomass impairment. The cellular toxicity caused by the inhibition of orotidine-5-phosphate decarboxylase was demonstrated in human leukemia cell lines [24], while orotate phosphoribosyltransferase is an antitumoral target for human gastric carcinoma cells [25]. Equally to the effects caused by these two enzymes, also aspartate carbamoyltransferase was predicted to inhibit the same biomass components (Table 1). The essentiality of last enzyme is demonstrated in human hepatocytes [26].

Similarly to the cited enzymes, dihydroorotase can also impair the synthesis of the same biomass components. Its essentiality is demonstrated in CCRF-CEM leukemic cells [27].

The reduced fitness approach correctly predicted that the enzymatic inhibition of acetyl-CoA carboxylase causes an importing of plasma triglycerides and this phenomenon has been experimentally demonstrated in vitro[28]. Acetyl-CoA carboxylase is the rate-limiting step of the fatty acid anabolism and it is essential for human breast cancer cells [29].

Carbamoyl-phosphate synthase was predicted to be more detrimental for the hepatocyte than acetyl-CoA carboxylase. This enzyme catalyzes the formation of carbamoyl phosphate from carbon dioxide and ammonia (EC 6.3.4.1.6) or glutamine (EC 6.3.5.5). Here, it was not discriminated between the ammonia-donors (NH4+ or glutamine) because the inhibition of one reaction was able to be fully compensated by the other. The essentiality of carbamoyl phosphate synthetase (EC 6.3.5.5) was demonstrated in mouse liver [30].

The last enzyme of this set was glycerol-3-phosphate acyltransferase, that is essential for human hepatocyte cultures [31].

The gene deletions of this antimalarial drug target set were confirmed by the aforementioned experimental data. Moreover, the reduced fitness approach correctly predicted the specific case that was triggered by the inhibition of acetyl-CoA carboxylase: while in the reference state the triglycerides are part of the maintenance function, thus they are exported into the blood, under the perturbation there was an inversion of the boundary flux and they were instead imported.

The behavior of the predicted fitness functions, that were simulated for this enzymatic set, are depicted in Figure 1.

thumbnailFigure 1. HepatoNet1 and PlasmoNet fitness profiles.

Reduced fitness rank of essential drug targets

For the 12 essential enzymes, we applied the reduced fitness concept to assess in a more subtle way their impact on the two metabolic networks (see Table 1). To validate the outcome of the reduced fitness rank, three enzymes were shortlisted from Table 1: acyl-CoA synthetase (score: 1.97); aspartate carbamoyltransferase (score: 0.713); glycerol-3-phosphate acyltransferase (score: 0.37).

The impairment of acyl-CoA synthetase in HepatoNet1 caused a perturbation in the mitochondrial synthesis of cardiolipin (–50%), while in PlasmoNet the inhibition severely destroyed the production of sphingomyelin (–97.07%). HepatoNet1 mitochondrial cardiolipin was produced by the condensation of glycerol-3-phosphate and CDP-diacylglycerol, that was obtained by mitochondrial phosphatidate. This latter metabolite was synthesized from acyl-CoA mitochondrial pool (that collects many activated fatty acids as palmitoyl-CoA and oleoyl-CoA). Cytosolic acyl-CoA syntethase was responsible for the formation of these activated fatty acids, thus an enzymatic impairment can affect the cardiolipin synthesis. This is experimentally demonstrated in human tumoral cells [32]. Acyl-CoA synthetase inhibition impaired PlasmoNet sphingomyelin production: the main precursors of sphingomyelin are serine and palmitoyl-CoA, that was produced in the cytosol by the enzyme. The parasite shows a high activity of fatty acid anabolism, thus the enzyme is important for plasmodial growth [33]. During the intra-erythrocytic stage, the parasite synthesizes new sphingolipids which are necessary for the formation of the tubovesicular network [34]. The formation of this membrane structure connects the parasitic vacuole with the host membrane during invasion. A possible activation of sphingomyelinase, that degrades sphingomyelin pools, provokes the plasmodial death. This suggests that a certain amount of sphingolipids is essential for Plasmodium falciparum.

The second enzyme is aspartate carbamoyltransferase whose inhibition caused an impairment of UDP-glucose production in HepatoNet1 (–80.57%) and mRNA in PlasmoNet (–57.44%). It is reported that inhibitors of aspartate carbamoyltransferase cause a 10% reduction of the UTP intracellular pools in hepatoma cell culture [35].

The last enzyme is glycerol-3-phosphate acyltrans- ferase, that caused dramatic effects in the hepatocyte metabolism, hindering the production of phosphatidylethanolamine (–55.77%), phosphatidylinositol (–50%), phosphatidylcholine (–50%), cardiolipin (–50%) and triglycerides (–50%). Phosphatidate is a common precursor of these metabolites and is formed in the human in the linear chain from 1-acyl-glycerol-3-phosphate, a product of glycerol-3-phosphate acyltransferase. A specific inhibitor of this enzyme (FSG67) causes in obese rats the reduction of triglycerides and phosphatidylcholine [36] and this was also confirmed by our analysis. The in silico inhibition of the enzyme impaired the production of phosphatidylcholine (–94.55%) in PlasmoNet. The plasmodial gene sequence that encodes glycerol-3-phosphate acyltransferase is expressed in double yeast mutants and biochemically characterized [37]. The authors suggested that it is likely that glycerol-3-phosphate acyltransferase is essential for a growing Plasmodium falciparum, that requires high amount of phospholipids for membrane synthesis. This hypothesis is sustained by another work, whose biochemical analyses show that parasites in trophozoite and schizont stages have an high acyltransferase activity [38].

Although it was not possible to assess the quantitative aspect of the selectivity score with experimental evidence, there was a good agreement among the obtained results and the literature.

Homologous drug targets predicted non-essential in both networks

Twelve drug targets showed up to be non-essential in both networks (Table 1), although their essentiality is experimentally validated for the parasite. In 8 cases the conflicting outcomes could be sorted out with literature-based assessment of the in vitro medium compositions and restriction of indicated inbound fluxes (conditional essentiality).

The remaining discrepancies may be due to possible molecular interactions that were not considered in the network reconstruction, e.g. off-target effects. For example, Velanker and coworkers find that plasmodial glycolytic enzymes associate with membranes and cytoskeleton components and drain their substrates near to the invasion machineries, making contact with host microtubules [39]. The authors also suggest that the inhibition of the glycolysis is achieved with enzymatic inhibitors or, alternatively, with the disruption of the cytoskeleton assembly. It is then likely that plasmodial glycolytic enzymes form a multi-enzymatic complex that is associated to the cytoskeleton. This hypothesis is also supposed by an old theoretical work, that suggests that the glycolytic enzymes are not ”evenly distributed throughout the cytosol”, but are likely localized in restricted regions [40]. It is also likely that these cytoskeleton-associated enzymes have an important role during host invasion and, thus, their inhibitors may cause off-target effects. In this set, gene deletions predicted that 3 glycolytic enzymes were dispensable. These were aldolase, hexokinase and triose-phosphate isomerase. Parasitic aldolase, for example, binds AMA1 effector to initiate host invasion [41]. Furthermore, the parasite relies on the glycolysis for its own ATP production, that is invested for growth, replication, motility and invasion. Hexokinase is found to be essential for the intra-erythrocytic stage of the protozoan [42]. On the other hand, it was predicted that hexokinase in intra-hepatic stage was dispensable for the achievement of the biomass reactions. Then, to simulate the effects of the upper glycolytic branch disruption, a double knock-out of hexokinase and glucose-6-phosphate isomerase was simulated. In this case, the impaired biomass components were phosphatidylinositols. This impairment and its consequences were in agreement with experimental assays, where the inhibition of hexokinase provokes the total disruption of the synthesis of glycophosphatidylinositols [43]. It is then likely that hexokinase and glucose-6-phosphate isomerase belong to the same multienzymatic complex that was proposed by Huebscher [40].

The last enzyme of this little set is triose-phosphate isomerase, that is expressed on the membrane of infected erythrocytes where it triggers antibody selection and prolonged hemolytic anemia [39].

Because of its specific extra-metabolic functions (infection/immune system activation), it is very likely that this enzyme is connected to off-target effects.

Among the enzymes that are conditionally essential, ornithine decarboxylase was found to be essential under the restriction of polyamines (as putrescine and spermidine) and under the block of agmatinase. This enzyme synthesizes the first polyamine: putrescine. In PlasmoNet, a secondary bypass through agmatinase (EC 3.5.1.53) made this enzyme dispensable. Agmatinase was identified in Plasmodium berghei and hypothesized in Plasmodium falciparum but never characterized here, so this bypass may not be present. Here, a double enzymatic knock-down was simulated and both enzymes resulted to be essential for PlasmoNet. Recently, it has been reported that a full perturbation of ornithine decarboxylase triggers sophisticated compensatory mechanisms on the transcriptome, proteome and metabolome of the parasite [44]. It is then likely that the enzyme is per se essential for the parasite and that the “rescue” mechanism is a evolutionary survival strategy. The literature indicates that the same enzyme is dispensable for the human host, suggesting a similar underlying conserved mechanism of regulation [45]. In the set of conditional essential enzymes, phosphoethanolamine methyltransferase was found to become essential under depletion of external choline and phosphatidylcholine. This enzyme had a replenishing function in the phospholipid synthetic pathway, transferring 3 methyl groups on the ethanolamine to form the choline head (in case of choline depletion). In Plasmodium berghei it has two different substrates, ethanolamine phosphate and phosphatidylethanolamine [46]. In PlasmoNet only the first reaction was included (ethanolamine-phosphate + 3 methyl-donors → choline-phosphate) and it was not active since phosphatidylcholine and choline were available in the simulated medium. This bypass is important only if choline is missing in the external environment, so to assess its essentiality the choline/phosphatidylcholine transporters were blocked and the missing reaction was added in PlasmoNet. Under these conditions the enzyme was essential. The last enzyme that was detected as dispensable is NADH dehydrogenase. Its metabolic functions can be fulfilled alternatively by Complex bc1. Inhibitors of NADH dehydrogenase (as atovaquone) target the mitochondrial electron transport chain but also affect Complex bc1[47]. Thus, a double knockout of NADH dehydrogenase/Complex bc1was applied, whose outcome indicated that the synergistic activity of both enzymes was essential for the parasite.

Conclusions

To increase the selectivity of an anti-pathogenic drug that could target also host enzymes, a conventional way is to exploit structural differences between homologous proteins. On the other hand, network-based modeling methods are alternative and complementary strategies to assess enzymatic inhibitions. In fact, genome-scale networks can detect enzymes that are likely to be essential and selective. Ideally, these approaches identify potential drug targets that are specific to the parasite.

While this idea may be impressive in its simplicity, the amount of parasite-specific enzymes can be restricted. In fact, out of the merged list of 96 “gold standards” only 18 enzymes are specific to Plasmodium falciparum metabolism.

It should also be considered that Plasmodium falciparum, when exposed to drug-induced selective pressure, develops drug-resistance (e.g. pyrimethamine and chloroquine) [48]. In this context, it would be nice to have more predictive methods and this was the rationale behind this research. The question this work aimed to address was the assessment of the selectivity of antimalarial drug targets with FBA-based methods. A human hepatocyte metabolic network was chosen to represent the host metabolism since the liver is the first human infection site for malaria parasites [17]. Furthermore, the liver possesses the most versatile metabolism among human cells. In fact, it is likely that drug targets, that compromise the metabolism of any human cell, should also be identifiable in the hepatocyte. Furthermore, it is the preferred organ to investigate the drug-induced metabolic impairments, whose consequences may be not detectable in preclinical and clinical trials [49].

To apply FBA-based methods, the main requirement is the availability of curated genome-scale metabolic networks. Furthermore, to achieve a realistic flux distribution, each model should be simulated with a literature-based set of metabolic objectives. This set allows the simulation of an anabolic physiological state of the metabolic model, in similar way to the biomass objective function [50]. The applied metabolic objective sets were assembled, approximated and rescaled from a multitude of different sources. The relative concentration share of each biomass component was taken as an estimate of its production rate. A physiological assembly of the metabolic objectives was necessary to describe realistically the consequences of the inhibitions of chosen antimalarial drug targets.

Overall, this work predicted that a large drug target set was non-essential for the hepatocyte model (24/48) and that 12 enzymes were instead essential for both models. For the first set gene deletions were simulated, while for the second set the concept of reduced fitness was applied. This last method can analyze more in detail the drug-induced impairments and their consequences on the network performance.

Gene-deletion simulations are a mean to assess the importance of a given biochemical reaction for the functionality of the network. This type of simulations are comparable to experimental gene-excision methods: the enzyme sequence is a priori disrupted and the resulting metabolic consequences are then investigated. On the other hand, it is rare that an administrated drug can achieve 100% enzymatic inhibition. To simulate this last scenario, the concept of reduced fitness was applied. In this case, 12 enzymes, that gene deletions predicted to be essential in both models, were studied in the context of reduced fitness. This method allows to compare the relative enzymatic essentiality and to understand which model is more sensitive to a chosen enzymatic restriction. Unfortunately, the predicted selectivity score could not be validated with the available pharmacological data (e.g. drug binding constant for the target and cytotoxic IC50 index) for infected and non-infected hepatocyte cultures. Our work had clearly some limitations, that were due to lacking pharmacological data on antimalarials and approximations of the metabolic objectives. Then, genome-scale metabolic networks do not usually integrate regulatory feed-back loops and, thus, they may be not feasible to predict consequences of enzyme inhibition that are triggered by negative regulation. As mentioned above for the case of thymidylate synthase, enzymatic drug-induced inhibition may cause enzymatic over-expression and this aspect can not be exhaustively predicted with the methods that are here applied.

Despite the applied approximations and the missing data, the obtained results were in agreement with the available literature. Thus, this framework may be useful to detect putative selective drug targets, that gene-deletion simulations may discard. Further analyses on the selectivity of antimalarial targets (by means of RNA interference and covalent inhibitor assays) are therefore required to validate the predicted selectivity score. RNA interference, for example, allows to tune the degree of inhibition, avoiding off-target effects due to unspecific protein binding [51]. Alternatively, enzyme impairments by covalent inhibitor assay [52,53] could provide the experimental mirror of the computed fitness function. In this last case, drug binding constant and cytotoxic index for each inhibitor will be useful to assess the ”pure” network effects that were here predicted. In fact, cytotoxic index alone is not informative in this respect, as it can not discriminate among strong network effects under weak drug binding and weak network effects under strong drug binding. Finally, the last important aspect to consider is the choice of the strain of Plasmodium falciparum, as drug resistance and sensitivity may largely vary among different strains. This is the rationale why anti-pathogenic drug discovery is now focusing on drug combinations, that are effective against a wide spectrum of pathogenic strains in low doses [54].

Methods

Drug target selection

Three literature-based datasets of validated antimalarial targets [14,18,19] were collected, merged and pruned. The merged set contains 96 experimentally validated essential enzymes reported for the parasite Plasmodium falciparum. Three enzymes are targets of approved drugs and only two of them are common in all sets (dihydrofolate reductase, dihydropteroate synthase). For 18 targets, no homologous enzymes are present in HepatoNet1. These pathogen-specific enzymes are suitable as drug targets and require no further in silico flux-based investigation. Thirty other enzymes belong to genetic functions and are only remotely intertwined with metabolism (e.g histone deacetylase, telomerase). Thus, they were discarded, since they are not represented in the networks. The remaining 48 enzymatic drug targets were present in both models. In Figure 2 the overlapping of the three datasets is depicted; the full list is given in Additional file 1.

Additional file 1. Merged set of validated antimalarial targets.

Format: XLS Size: 19KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

thumbnailFigure 2. Antimalarial drug targets and their overlap. For each circle the name indicates the first author. Globally the scheme depicts the distribution of 96 antimalarial drug targets.

The metabolic models

HepatoNet1 is a primarily literature-based metabolic network of a human hepatocyte, whose reactions are individually curated and functionally tested [16]. It comprises 2539 reactions, 704 genes and 1149 metabolites. To simulate the pathogen, PlasmoNet model was chosen (Biomodels database ID: MODEL1111240000). PlasmoNet is a large metabolic network of Plasmodium falciparum, whose reconstruction is based on information from several databases and extensive literature search [14]. While the published version consists of 1622 metabolites, 1375 reactions and 579 genes, here the model was modified for the scope of the research and the current version (PlasmoNet v2.0, Biomodels database id: MODEL1206070000) contained 1394 reactions: 20 inbound reactions, measured with metabolomics assays [55], were included; one reaction of steroid hormone pathway was removed (KEGG id: R01836; EC 1.1.1.239), since it is likely to be present only in mammalian genitourinary system [56]. Previous versions of the KEGG database [57] indicated that this reaction belongs to the plasmodial metabolism. The correction was then introduced in KEGG version 57.0 and above. Furthermore, an irreversible directionality was added for the reactions catalyzed by S-adenosyl-L-methionine decarboxylase (EC 4.1.1.50) and phosphoenolpyruvate carboxylase (EC 4.1.1.31), as indicated in BRENDA database [58]. The PlasmoNet v2.0 is included as Additional file 2. Detailed information about the applied parameters (reaction directions, imported/exported metabolites) is available on request.

Additional file 2. Updated version of PlasmoNet. The network is available on Biomodels database under the ID MODEL1206070000.

Format: XML Size: 1.2MB Download fileOpen Data

Definition of metabolic objectives

The metabolic objectives of the two networks were formulated in terms of the so-called biomass reactions which: (i) yield building blocks for the cellular composition, (ii) remove potentially harmful metabolites (e.g. toxins) and (iii) are exported by the cell in the context of systemic physiological functions. For HepatoNet1, literature search collected concentration shares of 98 metabolites, that are cellular building blocks or that can be exported into the blood. The list of applied HepatoNet1 metabolic functions is included as Additional file 3. In some cases, human liver data were not available, thus murine information was applied and rescaled. Also several approximations and assumptions were applied (e.g. the amount of cardiolipin is measured in the periportal and perivenous liver areas [59] and in this case an arithmetic mean of the two concentrations was applied). The metabolites, that HepatoNet1 was allowed to import from the extra-cellular environment, were the human blood components and their inbound fluxes were left unconstrained [16]. The set of PlasmoNet biomass components consisted of 98 metabolites, whose concentration shares were retrieved from literature. An initial set of 57 values was already published along with PlasmoNet reconstruction [14], but for the scope of the research this set was enlarged to 98 metabolites. The environment in which PlasmoNet is embedded represented the host cytoplasm. To describe a feasible situation for a growing parasite, all host-pathogen exchange reactions were left unconstrained. The full list of PlasmoNet metabolic objectives is given as Additional file 4.

Additional file 3. HepatoNet1 applied metabolic objectives.

Format: TXT Size: 3KB Download fileOpen Data

Additional file 4. PlasmoNet applied metabolic objectives.

Format: TXT Size: 5KB Download fileOpen Data

Gene deletions

To test the essentiality of drug targets, gene deletions were performed under the flux minimization framework [60], in similar way to previous studies [6-8]. The optimization problem aimed (i) to minimize the sum of internal fluxes and (ii) the simultaneous achievements of all biomass reactions. To simulate gene-deletion of a drug target, the fluxes through the corresponding catalyzed reactions were constrained to zero. Under these conditions, a successful simulation would predict that the enzyme is non-essential and all biomass reactions can be fulfilled. If no feasible solutions can be reached, the enzyme is considered to be essential. A more detailed explanation of the flux minimization problem is given in Additional file 5: Appendix A.

Additional file 5. Appendix A. Flux minimization optimization framework.

Format: PDF Size: 94KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Reduced fitness approach

This approach has been applied to an erythrocyte model to simulate the impact of enzyme deficiencies on network performance [20]. The initial step is the calculation of a reference state without any impairment. Then, for each potential target enzyme, a metabolic scenario is simulated, where (i) the fluxes through the reactions catalyzed by this enzyme are progressively restricted and (ii) the deviation of the biomass fluxes from the reference state is minimized. The deviations express the model’s impairment under the enzymatic perturbation, the inverse is called fitness. A detailed explanation of this concept together with an exemplary application are given in Additional file 6: Appendix B.

Additional file 6. Appendix B. Reduced Fitness optimization framework and example.

Format: PDF Size: 410KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Enzymatic fitness estimation and selectivity score

To assess the network performance under an enzymatic impairment, the initial reference state is obtained with flux minimization method. This framework does not predict a unique solution but this aspect is common to many FBA methods that are based on optimization problems. To calculate the reduced fitness for an enzymatic knock-down, the enzyme-catalyzed null fluxes are blocked. The non-null fluxes, instead, are subjected to a progressive reduction of 1% of the reference flux values for each simulation run. For each decremental step of the non-null fluxes, a single value of reduced fitness is calculated. In this way, a single fitness curve is obtained for each enzymatic knock-down (Figure 1). To define selectivity in terms of the reduced fitness, the area under the fitness curves (AUC) of the networks is calculated and compared for each enzymatic inhibition. In case of a non-essential drug target, the obtained fitness curve is a straight horizontal line at fitness 1 and, applying the progressive decrement of 1% of the reference enzymatic activity, its AUC is 100 (as a per cent measure). The larger the area under the fitness curve (i.e. the smaller its deviation from 100), the less is the effect of the drug target on metabolic fitness. The AUC deviation from 100 defines the selectivity score (equation 1).

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/118/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/118/mathml/M1">View MathML</a>

(1)

Then, the selectivity score of a drug target is defined as

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/118/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/118/mathml/M2">View MathML</a>

(2)

Athe possible extra-metabolic role of the biomass components.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SB developed the original idea, carried out all computational analyses and drafted the manuscript. AH provided tools for the knock-out simulations and together with HGH and SB participated in the design and evaluation of the analyses. All authors contributed to and approved the final manuscript.

Acknowledgements

SB gratefully acknowledges the help provided by Dr. Bernd Binder for his help, the valuable suggestions and useful discussions. SB is also grateful to Dr. Wolfram Liebermeister for the useful discussions and to the anonymous referees for their constructive remarks and comments.

SB acknowledges the support by SFB618-Collaborative Research Center for Theoretical Biology. AH acknowledges the support by the BMBF Germany via the VirtualLiver Network funding initiative.

References

  1. Lin Z, Will Y: Evaluation of Drugs with Specific Organ Toxicities in Organ Specific Cell Lines.

    Toxicol Sci 2012, 126(1):114-27.

    [ http://view.ncbi.nlm.nih.gov/pubmed/22166485 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  2. Armstrong PB: Proteases and protease inhibitors: a balance of activities in host-pathogen interaction.

    Immunobiology 2006, 211(4):263-81.

    [ http://view.ncbi.nlm.nih.gov/pubmed/16697919 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  3. Singh VK, Ghosh I: Kinetic modeling of tricarboxylic acid cycle and glyoxylate bypass in Mycobacterium tuberculosis, and its application to assessment of drug targets.

    Theor Biol Med Model 2006, 3:27.

    [ http://view.ncbi.nlm.nih.gov/pubmed/16887020 webcite]

    PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  4. Hornberg JJ, Bruggeman FJ, Bakker BM, Westerhoff HV: Metabolic control analysis to identify optimal drug targets.

    Prog Drug Res 2007, 64(171):173-89.

    [ http://view.ncbi.nlm.nih.gov/pubmed/17195475 webcite]

    OpenURL

  5. Murabito E, Smallbone K, Swinton J, Westerhoff HV, Steuer R: A probabilistic approach to identify putative drug targets in biochemical networks.

    J R Soc Interface 2011, 8(59):880-95.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21123256 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Edwards JS, Palsson BO: Metabolic flux balance analysis and the in silico analysis of Escherichia coli K-12 gene deletions.

    BMC Bioinformatics 2000, 1:1.

    [ http://view.ncbi.nlm.nih.gov/pubmed/11001586 webcite]

    PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  7. Snitkin ES, Dudley AM, Janse DM, Wong K, Church GM, Segrè D: Model-driven analysis of experimentally determined growth phenotypes for 465 yeast gene deletion mutants under 16 different conditions.

    Genome Biol 2008, 9(9):R140.

    [ http://view.ncbi.nlm.nih.gov/pubmed/18808699 webcite]

    PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  8. Raman K, Rajagopalan P, Chandra N: Flux balance analysis of mycolic acid pathway: targets for anti-tubercular drugs.

    PLoS Comput Biol 2005, 1(5):e46.

    [ http://view.ncbi.nlm.nih.gov/pubmed/16261191 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Lee DS, Burd H, Liu J, Almaas E, Wiest O, Barabási AL, Oltvai ZN, Kapatral V: Comparative genome-scale metabolic reconstruction and flux balance analysis of multiple Staphylococcus aureus genomes identify novel antimicrobial drug targets.

    J Bacteriol 2009, 191(12):4015-24.

    [ http://view.ncbi.nlm.nih.gov/pubmed/19376871 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Mazumdar V, Snitkin ES, Amar S, Segrè D: Metabolic network model of a human oral pathogen.

    J Bacteriol 2009, 191:74-90.

    [ http://view.ncbi.nlm.nih.gov/pubmed/18931137 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Oberhardt MA, Puchałka J Fryer KE: Genome-scale metabolic network analysis of the opportunistic pathogen Pseudomonas aeruginosa PAO1.

    J Bacteriol 2008, 190(8):2790-803.

    [ http://view.ncbi.nlm.nih.gov/pubmed/18192387 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Fatumo S, Plaimas K, Adebiyi E, König R: Comparing metabolic network models based on genomic and automatically inferred enzyme information from Plasmodium and its human host to define drug targets in silico.

    Infect Genet Evol 2011, 11(4):708-15.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21515412 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  13. Bordbar A, Lewis NE, Schellenberger J, Palsson BO, Jamshidi N: Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions.

    Mol Syst Biol 2010, 6:422.

    [ http://view.ncbi.nlm.nih.gov/pubmed/20959820 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Huthmacher C, Hoppe A, Bulik S, Holzhütter HG: Antimalarial drug targets in Plasmodium falciparum predicted by stage-specific metabolic network analysis.

    BMC Syst Biol 2010, 4:120.

    [ http://view.ncbi.nlm.nih.gov/pubmed/20807400 webcite]

    PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  15. Folger O, Jerby L, Frezza C, Gottlieb E, Ruppin E, Shlomi T: Predicting selective drug targets in cancer through metabolic networks.

    Mol Syst Biol 2011, 7:501.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21694718 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Gille C, Bölling C, Hoppe A, Bulik S, Hoffmann S, Hübner K, Karlstädt A, Ganeshan R, König M, Rother K, Weidlich M, Behre J, Holzhütter HG: HepatoNet1: a comprehensive metabolic reconstruction of the human hepatocyte for the analysis of liver physiology.

    Mol Syst Biol 2010, 6:411.

    [ http://view.ncbi.nlm.nih.gov/pubmed/20823849 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Nussenzweig RS, Long CA: Malaria vaccines: multiple targets.

    Science 1994, 265(5177):1381-3.

    [ http://view.ncbi.nlm.nih.gov/pubmed/8073276 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  18. Fatumo S, Plaimas K, Mallm JP, Schramm G, Adebiyi E, Oswald M, Eils R, König R: Estimating novel potential drug targets of Plasmodium falciparum by analysing the metabolic network of knock-out strains in silico.

    Infect Genet Evol 2009, 9(3):351-8.

    [ http://view.ncbi.nlm.nih.gov/pubmed/18313365 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  19. Yeh I, Hanekamp T, Tsoka S, Karp PD, Altman RB: Computational analysis of Plasmodium falciparum metabolism: organizing genomic information to facilitate drug discovery.

    Genome Res 2004, 14(5):917-24.

    [ http://view.ncbi.nlm.nih.gov/pubmed/15078855 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Holzhütter HG: The generalized flux-minimization method and its application to metabolic networks affected by enzyme deficiencies.

    Biosystems 2006, 83(2-3):98-107.

    [ http://view.ncbi.nlm.nih.gov/pubmed/16229937 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  21. Tomoda H, Igarashi K, Cyong JC, Omura S: Evidence for an essential role of long chain acyl-CoA synthetase in animal cell proliferation. Inhibition of long chain acyl-CoA synthetase by triacsins caused inhibition of Raji cell proliferation.

    J Biol Chem 1991, 266(7):4214-9.

    [ http://view.ncbi.nlm.nih.gov/pubmed/1999415 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  22. Galvani E, Peters GJ, Giovannetti E: Thymidylate synthase inhibitors for non-small cell lung cancer.

    Expert Opin Investig Drugs 2011, 20(10):1343-56.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21905922 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  23. Wilson PM, Labonte MJ, Lenz HJD, Mack PC, Ladner RD: Inhibition of dUTPase induces synthetic lethality with thymidylate synthase-targeted therapies in non-small cell lung cancer.

    Mol Cancer Ther 2012, 11(3):616-28.

    [ http://view.ncbi.nlm.nih.gov/pubmed/22172489 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  24. Bello AM, Konforte D, Poduch E, Furlonger C, Wei L, Liu Y, Lewis M, Pai EF, Paige CJ, Kotra LP: Structure-activity relationships of orotidine-5’-monophosphate decarboxylase inhibitors as anticancer agents.

    J Med Chem 2009, 52(6):1648-58.

    [ http://view.ncbi.nlm.nih.gov/pubmed/19260677 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  25. Ochiai T, Nishimura K, Noguchi H, Kitajima M, Tsukada A, Watanabe E, Nagaoka I, Futagawa S: Prognostic impact of orotate phosphoribosyl transferase among 5-fluorouracil metabolic enzymes in resectable colorectal cancers treated by oral 5-fluorouracil-based adjuvant chemotherapy.

    Int J Cancer 2006, 118(12):3084-8.

    [ http://view.ncbi.nlm.nih.gov/pubmed/16425285 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  26. Pausch J, Rasenack J, Häussinger D, Gerok W: Hepatic carbamoyl phosphate metabolism. Role of cytosolic and mitochondrial carbamoyl phosphate in de novo pyrimidine synthesis.

    Eur J Biochem 1985, 150:189-94.

    [ http://view.ncbi.nlm.nih.gov/pubmed/4018077 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  27. Brooke J, Szabados E, Lyons SD, Goodridge RJ, Harsanyi MC, Poiner A, Christopherson RI: Cytotoxic effects of dihydroorotase inhibitors upon human CCRF-CEM leukemia.

    Cancer Res 1990, 50(24):7793-8.

    [ http://view.ncbi.nlm.nih.gov/pubmed/1979249 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  28. Keil S, Müller M, Zoller G, Haschke G, Schroeter K, Glien M, Ruf S, Focken I, Herling AW, Schmoll D: Identification and synthesis of novel inhibitors of acetyl-CoA carboxylase with in vitro and in vivo efficacy on fat oxidation.

    J Med Chem 2010, 53(24):8679-87.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21082864 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  29. Chajès V, Cambot M, Moreau K, Lenoir GM, Joulin V: Acetyl-CoA carboxylase alpha is essential to breast cancer cell survival.

    Cancer Res 2006, 66(10):5287-94.

    [ http://view.ncbi.nlm.nih.gov/pubmed/16707454 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  30. Marrero M, Prough RA, Putnam RS, Bennett M, Milewich L: Inhibition of carbamoyl phosphate synthetase-I by dietary dehydroepiandrosterone.

    J Steroid Biochem Mol Biol 1991, 38(5):599-609.

    [ http://view.ncbi.nlm.nih.gov/pubmed/1828177 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  31. Wydysh EA, Medghalchi SM, Vadlamudi A, Townsend CA: Design and synthesis of small molecule glycerol 3-phosphate acyltransferase inhibitors.

    J Med Chem 2009, 52(10):3317-27.

    [ http://view.ncbi.nlm.nih.gov/pubmed/19388675 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Mashima T, Oh-hara T, Sato S, Mochizuki M, Sugimoto Y, Yamazaki K, Hamada Ji, Tada M, Moriuchi T, Ishikawa Y, Kato Y, Tomoda H, Yamori T, Tsuruo T: p53-defective tumors with a functional apoptosome-mediated pathway: a new therapeutic target.

    J Natl Cancer Inst 2005, 97(10):765-77.

    [ http://view.ncbi.nlm.nih.gov/pubmed/15900046 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  33. Mitamura T, Palacpac NMQ: Lipid metabolism in Plasmodium falciparum-infected erythrocytes: possible new targets for malaria chemotherapy.

    Microbes Infect 2003, 5(6):545-52.

    [ http://view.ncbi.nlm.nih.gov/pubmed/12758284 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  34. Labaied M, Dagan A, Dellinger M, Gèze M, Egée S, Thomas SL, Wang C, Gatt S, Grellier P: Anti-Plasmodium activity of ceramide analogs.

    Malar J 2004, 3:49.

    [ http://view.ncbi.nlm.nih.gov/pubmed/15588325 webcite]

    PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  35. Keppler DO, Schulz-Holstege C, Fauler J, Reiffen KA, Schneider F: Uridylate trapping, induction of UTP deficiency, and stimulation of pyrimidine synthesis de novo by D-galactosone.

    Biochem J 1982, 206:139-46.

    [ http://view.ncbi.nlm.nih.gov/pubmed/7126188 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Kuhajda FP, Aja S, Tu Y, Han WF, Medghalchi SM, El Meskini R, Landree LE, Peterson JM, Daniels K, Wong K, Wydysh EA, Townsend CA, Ronnett GV: Pharmacological glycerol-3-phosphate acyltransferase inhibition decreases food intake and adiposity and increases insulin sensitivity in diet-induced obesity.

    Am J Physiol Regul Integr Comp Physiol 2011, 301:R116-30.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21490364 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Santiago TC, Zufferey R, Mehra RS, Coleman RA, Mamoun CB: The Plasmodium falciparum PfGatp is an endoplasmic reticulum membrane protein important for the initial step of malarial glycerolipid synthesis.

    J Biol Chem 2004, 279(10):9222-32.

    [ http://view.ncbi.nlm.nih.gov/pubmed/14668349 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  38. Palacpac NMQ, Hiramine Y, Mi-ichi F, Torii M, Kita K, Hiramatsu R, Horii T, Mitamura T: Developmental-stage-specific triacylglycerol biosynthesis, degradation and trafficking as lipid bodies in Plasmodium falciparum-infected erythrocytes.

    J Cell Sci 2004, 117(Pt 8):1469-80.

    [ http://view.ncbi.nlm.nih.gov/pubmed/15020675 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  39. Velanker SS, Ray SS, Gokhale RS, Suma S, Balaram H, Balaram P, Murthy MR: Triosephosphate isomerase from Plasmodium falciparum: the crystal structure provides insights into antimalarial drug design.

    Structure 1997, 5(6):751-61.

    [ http://view.ncbi.nlm.nih.gov/pubmed/9261072 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  40. Hübscher G, Mayer RJ, Hansen HJ: Glycolytic enzymes as a multi-enzyme system.

    J Bioenerg 1971, 2(2):115-8.

    [ http://view.ncbi.nlm.nih.gov/pubmed/4257578 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  41. Srinivasan P, Beatty WL, Diouf A, Herrera R, Ambroggio X, Moch JK, Tyler JS, Narum DL, Pierce SK, Boothroyd JC, Haynes JD, Miller LH: Binding of Plasmodium merozoite proteins RON2 and AMA1 triggers commitment to invasion.

    Proc Natl Acad Sci U S A 2011, 108(32):13275-80.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21788485 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. van Schalkwyk D A, Priebe W, Saliba KJ: The inhibitory effect of 2-halo derivatives of D-glucose on glycolysis and on the proliferation of the human malaria parasite Plasmodium falciparum.

    J Pharmacol Exp Ther 2008, 327(2):511-7.

    [ http://view.ncbi.nlm.nih.gov/pubmed/18713952 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  43. Santos de Macedo C, Gerold P, Jung N, Azzouz N, Kimmel J, Schwarz RT: Inhibition of glycosyl-phosphatidylinositol biosynthesis in Plasmodium falciparum by C-2 substituted mannose analogues.

    Eur J Biochem 2001, 268(23):6221-8.

    [ http://view.ncbi.nlm.nih.gov/pubmed/11733018 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  44. van Brummelen A C, Olszewski KL, Wilinski D, Llinás M, Louw AI, Birkholtz LM: Co-inhibition of Plasmodium falciparum S-adenosylmethionine decarboxylase/ornithine decarboxylase reveals perturbation-specific compensatory mechanisms by transcriptome, proteome, and metabolome analyses.

    J Biol Chem 2009, 284(7):4635-46.

    [ http://view.ncbi.nlm.nih.gov/pubmed/19073607 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  45. Zini M, Passariello CL, Gottardi D, Cetrullo S, Flamigni F, Pignatti C, Minarini A, Tumiatti V, Milelli A, Melchiorre C, Stefanelli C: Cytotoxicity of methoctramine and methoctramine-related polyamines.

    Chem Biol Interact 2009, 181(3):409-16.

    [ http://view.ncbi.nlm.nih.gov/pubmed/19576191 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  46. Déchamps S, Wengelnik K, Cerdan R, Vial HJ Gannoun-Zaki, Berry-Sterkers L: The Kennedy phospholipid biosynthesis pathways are refractory to genetic disruption in Plasmodium berghei and therefore appear essential in blood stages.

    Mol Biochem Parasitol 2010, 173(2):69-80.

    [ http://view.ncbi.nlm.nih.gov/pubmed/20478340 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  47. Saleh A, Friesen J, Baumeister S, Gross U, Bohne W: Growth inhibition of Toxoplasma gondii and Plasmodium falciparum by nanomolar concentrations of 1-hydroxy-2-dodecyl-4(1H)quinolone, a high-affinity inhibitor of alternative (type II) NADH dehydrogenases.

    Antimicrob Agents Chemother 2007, 51(4):1217-22.

    [ http://view.ncbi.nlm.nih.gov/pubmed/17242151 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Mita T, Tanabe K, Kita K: Spread and evolution of Plasmodium falciparum drug resistance.

    Parasitol Int 2009, 58(3):201-9.

    [ http://view.ncbi.nlm.nih.gov/pubmed/19393762 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  49. Tang W: Drug metabolite profiling and elucidation of drug-induced hepatotoxicity.

    Expert Opin Drug Metab Toxicol 2007, 3(3):407-20.

    [ http://view.ncbi.nlm.nih.gov/pubmed/17539747 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  50. Feist AM, Palsson BO: The biomass objective function.

    Curr Opin Microbiol 2010, 13(3):344-9.

    [ http://view.ncbi.nlm.nih.gov/pubmed/20430689 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Kolev NG, Tschudi C, Ullu E: RNA interference in protozoan parasites: achievements and challenges.

    Eukaryot Cell 2011, 10(9):1156-63.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21764910 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Bello AM, Poduch E, Liu Y, Wei L, Crandall I, Wang X, Dyanand C, Kain KC, Pai EF, Kotra LP: Structure-activity relationships of C6-uridine derivatives targeting plasmodia orotidine monophosphate decarboxylase.

    J Med Chem 2008, 51(3):439-48.

    [ http://view.ncbi.nlm.nih.gov/pubmed/18189347 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  53. Bello AM, Poduch E, Fujihashi M, Amani M, Li Y, Crandall I, Hui R, Lee PI, Kain KC, Pai EF, Kotra LP: A potent, covalent inhibitor of orotidine 5’-monophosphate decarboxylase with antimalarial activity.

    J Med Chem 2007, 50(5):915-21.

    [ http://view.ncbi.nlm.nih.gov/pubmed/17290979 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  54. Chait R, Craney A, Kishony R: Antibiotic interactions that select against resistance.

    Nature 2007, 446(7136):668-71.

    [ http://view.ncbi.nlm.nih.gov/pubmed/17410176 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  55. Olszewski KL, Morrisey JM, Wilinski D, Burns JM, Vaidya AB, Rabinowitz JD, Llins M: Host-parasite interactions revealed by Plasmodium falciparum metabolomics.

    Cell Host Microbe 2009, 5(2):191-9.

    [ http://view.ncbi.nlm.nih.gov/pubmed/19218089 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Azzarello J, Fung KM, Lin HK: Tissue distribution of human AKR1C3 and rat homolog in the adult genitourinary system.

    J Histochem Cytochem 2008, 56(9):853-61.

    [ http://view.ncbi.nlm.nih.gov/pubmed/18574251 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes.

    Nucleic Acids Res 2000, 28:27-30.

    [ http://view.ncbi.nlm.nih.gov/pubmed/10592173 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Scheer M, Grote A, Chang A, Schomburg I, Munaretto C, Rother M, Söhngen C, Stelzer M, Thiele J, Schomburg D: BRENDA, the enzyme information system in 2011.

    Nucleic Acids Res 2011, 39(Database issue):D670—6.

    [ http://view.ncbi.nlm.nih.gov/pubmed/21062828 webcite]

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Castro J, Cortés JP, Guzmán M: Properties of the mitochondrial membrane and carnitine palmitoyltransferase in the periportal and the perivenous zone of the liver. Effects of chronic ethanol feeding.

    Biochem Pharmacol 1991, 41(12):1987-95.

    [ http://view.ncbi.nlm.nih.gov/pubmed/2039548 webcite]

    PubMed Abstract | Publisher Full Text OpenURL

  60. Holzhütter HG: The principle of flux minimization and its application to estimate stationary fluxes in metabolic networks.

    Eur J Biochem 2004, 271(14):2905-22.

    [ http://view.ncbi.nlm.nih.gov/pubmed/15233787 webcite]

    PubMed Abstract | Publisher Full Text OpenURL