Predicting and characterizing selective multiple drug treatments for metabolicdiseases and cancer
© Facchetti et al.; licensee BioMed Central Ltd. 2012
Received: 18 May 2012
Accepted: 13 August 2012
Published: 29 August 2012
In the field of drug discovery, assessing the potential of multidrug therapies isa difficult task because of the combinatorial complexity (both theoretical andexperimental) and because of the requirements on the selectivity of the therapy.To cope with this problem, we have developed a novel method for the systematic insilico investigation of synergistic effects of currently available drugs ongenome-scale metabolic networks.
The algorithm finds the optimal combination of drugs which guarantees theinhibition of an objective function, while minimizing the side effect on the othercellular processes. Two different applications are considered: finding drugsynergisms for human metabolic diseases (like diabetes, obesity and hypertension)and finding antitumoral drug combinations with minimal side effect on the normalhuman cell. The results we obtain are consistent with some of the availabletherapeutic indications and predict new multiple drug treatments. A clusteranalysis on all possible interactions among the currently available drugsindicates a limited variety on the metabolic targets for the approved drugs.
The in silico prediction of drug synergisms can represent an important tool forthe repurposing of drugs in a realistic perspective which considers also theselectivity of the therapy. Moreover, for a more profitable exploitation ofdrug-drug interactions, we have shown that also experimental drugs which have adifferent mechanism of action can be reconsider as potential ingredients of newmulticompound therapeutic indications. Needless to say the clues provided by acomputational study like ours need in any case to be thoroughly evaluatedexperimentally.
KeywordsMetabolic network Drug synergism Flux balance analysis Metabolic diseases Cancer
In spite of the advances in molecular and computational biology, the discovery of newdrugs still remains a very challenging task which requires a very long period ofresearch and development before any new compound can be commercialized. A possiblealternative to the search of new active compounds is to make use of the unexploitedproperties of already available drugs, since a wide knowledge about both theirtherapeutic and toxicity effects has already been gathered during the study for theirapproval. In this perspective, a natural approach to broaden the range of applicationsof the existing drugs is to try to combine them in multiple drug therapies [1–3]. However, even though both the financial burden of conducting trials as wellas the risk of adverse events in trial populations is expected to be sensibly lower foralready approved drugs, so far the experimental investigations of multicomponenttherapies have been quite limited [4, 5]. Major obstacles to this approach are the high number of possiblecombinations but also our limited understanding of the complex mode of action of amultidrug treatment. Indeed, multiple perturbations can show three types of interaction,which have been classified as synergistic, antagonistic and additive  (alternatively called, aggravating, buffering and non-epistatic ). We have focused our attention on the first type, where the use of drugcombinations represents an enhancement with respect to the superposition of the singleperturbations.
In order to identify synergistic effects, Ref.  investigated all pairs of a set of known drugs at different doses, obtaininga functional classification of the interactions by looking at their inhibitory effect.Various computational approaches, based on reconstructed genome-scale metabolicnetworks, have been also developed in Refs. [7–10], in order to identify the synergistic effects triggered by multiple drugs ormultiple genetic perturbations (for example the so-called “syntheticlethality”). Unfortunately, a systematic evaluation of the effects of all possiblecombinations of drugs is unfeasible, because their number scales exponentially with thenumber of chemicals taken into account (for an exhaustive search over 40 drugs, morethan a trillion possible combinations should be tested). Moreover, a drug profile isgiven by both its therapeutic effect and by its side effect, the latter being related toits selectivity. For example, drugs such as anticancer agents have to selectively actonly on tumor cells . Similarly, metabolic diseases are induced by the imbalance of key metabolicpathways which, if modulated without affecting other vital functions, can rescue fromthe pathology [12, 13]. Referring to these specific needs, we have developed an algorithm based onthe metabolic network of humans and on the comparison between the metabolic networks ofhuman and cancer cells (as reconstructed recently in ) aimed at expanding the spectrum of applications of the existing drugs to newselective treatments against metabolic diseases and tumors.
In spite of the limitations of restricting our analysis only to metabolic networks, theresults presented in the following provide a promising example of how in silicomodels can be used as practical tools for exploring the functional interactions betweendrugs and of the (little explored) potential offered by synergistic drugcombinations.
Results and discussion
the drug combination leads to a new inhibition since no singledrug can stop the objective reaction;
the objective reaction can be stopped also by a single drug but the drugcombination is more selective (has a minor side effect);
the objective reaction can be stopped also by a single drug and themultiple drug solution is less selective (this solution is not interestingbecause it triggers a larger side effect).
List of all drug synergisms
Rosiglitazone (# 7) - Quinacrine (# 36) -Cerulenin (# 62) - Tyloxapol (# 85)
Rosiglitazone (# 7) - Quinacrine (# 36) - Orlistat(# 65) - Tyloxapol (# 85)
Rosiglitazone (# 7) - Indomethacin (# 22) - Cerulenin (# 62)- Tyloxapol (# 85)
Rosiglitazone (# 7) - Diclofenac (# 35) -Cerulenin (# 62) - Tyloxapol (# 85)
Rosiglitazone (# 7) - Indomethacin (# 22) - Orlistat(# 65) - Tyloxapol (# 85)
Rosiglitazone (# 7) - Diclofenac (# 35) - Orlistat(# 65) - Tyloxapol (# 85)
Rosiglitazone (# 7) - Cerulenin (# 62)
Rosiglitazone (# 7) - Orlistat (# 65) -
Indomethacin (# 22) - Fomepizole (# 75)
Naftifine (# 43) - Acetylsalicylic acid (# 55)
Acetylsalicylic acid (# 55) - Tioconazole (# 60)
Simvastatin/Pravastatin (# 4) - Acetylsalicylicacid (# 55)
Rosiglitazone (# 7) - Tioconazole (# 60)
Rosiglitazone (# 7) - Naftifine (# 43)
Simvastatin/Pravastatin (# 4) - Rosiglitazone(# 7)
Carbidopa (# 6) - Droxidopa (# 24)
Droxidopa (# 24) - Selegiline (# 45)
Droxidopa (# 24) - Minaprine (# 49)
Droxidopa (# 24) - Zonisamide (# 54)
Mycophenolic acid (# 42) - Mercaptopurine (# 58)
Ribavirin (# 51) - Mercaptopurine (# 58)
Udenafil (# 10) - Mycophenolic acid (# 42) -Mercaptopurine (# 58) -
Mycophenolic acid (# 42) - Dipyridamole (# 57) -Mercaptopurine (# 58)
Udenafil (# 10) - Ribavirin (# 51) - Mercaptopurine(# 58)
Ribavirin (# 51) - Dipyridamole (# 57) -Mercaptopurine (# 58)
Theophylline (# 18) - Mycophenolic acid (# 42) -Mercaptopurine (# 58)
Mycophenolic acid (# 42) - Pentoxifylline (# 50) -Mercaptopurine (# 58)
Theophylline (# 18) - Ribavirin (# 51) -Mercaptopurine (# 58)
Pentoxifylline (# 50) - Ribavirin (# 51) -Mercaptopurine (# 58)
Pentoxifylline (# 50) - Arsenic trioxide (# 72)
Cladribirne (# 16) - Pentoxifylline (# 50)
Gemcitabine (# 30) - Pentoxifylline (# 50)
The complete list of the objective reactions, with relative pathways and synergisticinhibitions, is reported in Additional file 1: Table S2:this list contains many inhibitions in the fatty acid, cholesterol and carnitinetransport pathways, which may represent solutions for obesity. In particular,concerning the case of hyperlipidemia diseases, the algorithm finds the combinationof Rosiglitazone and Cerulenin (class A) as inhibitor of many reactions in thecarnitine transferase and fatty acid desaturase pathways. This synergism has beenreported in the literature for being active versus the biosynthesis of fatty acids inprostate tumors : the mean IC50 are 45μ M and 32μ Mfor Rosiglitazone and Cerulenin alone, whereas it reduces to 5μ M whenthey are combined: the authors claimed that this effect comes from the reducedproduction of fatty acids preventing the growth and differentiation of prostatecells. It is worth noting that we predict this combination also in three anticancersolutions (together with an additional target on palmitate conversion, see nextsubsection). Moreover, this pair is part of other six synergisms on the humanmetabolism, all concerning the same pathways (i.e. a very similar mechanism ofinteraction).
Another significant example is represented by the inhibition of Dihydroceramidedesaturase. Indeed, ceramide is the hydrophobic membrane anchor of sphingolipids andis involved as a bioactive molecule in cell growth regulation, apoptosis, senescence,and diverse cell responses, particularly those linked to stress situations [28, 29]; moreover, recent studies have shown the role for ceramide biosynthesis inbody weight regulation, energy expenditure, hence in the metabolic obesity syndrome . For these reasons Dihydroceramide desaturase has been proposed as apromising potential target for metabolic diseases. Currently some specific inhibitorsof Dihydroceramide desaturase are under investigation and development (for exampleGT11, XM462 and analogous [31, 32]) although there is no approved drug yet. Indeed in our model no singledrug can stop this reaction. Our algorithm finds some possible multidrug treatmentsthat block this reaction: among them, there are the synergistic pairs ofRosiglitazone plus Simvastatin (Pravastatin), and Acetylsalicyclic acid plusAtorvastatin (Pravastatin). Concerning the first synergism, clinical experiments haveshown that combining these two drugs a significant reduction (about 30% less) of theintracellular accumulation of lipid is achieved . Moreover, the authors of  investigate the adverse effect of single and combined therapies(hypoglycemia, body weight increase) and claim that adverse events are generallysimilar (the safety profile of Rosiglitazone was not adversely affected by theaddition of Atorvastatin). Also our results predict a limited worsening of theadverse effect: indeed, after the combination with Atorvastatin the side effect ofRosiglitazone passes from 252.7 (Additional file 1: TableS1) to 288.8 (Table 2), i.e. it increases of about 14% only.For the same therapeutic purpose, the pair of Acetylsalicyclic acid and Atorvastatinhas been also studied. Clinical trials are currently ongoing  and some of them have already shown promising results . The rationale for this approach is based on the restoration of plateletsensitivity by reduction of the cholesterol levels. Moreover, as mentioned above, itis known that ceramide is involved in apoptosis. Indeed, this combination has beentested for the treatment of prostate cancer: the results have shows a linearsynergism between these two drugs .
Human vs Cancer
A similar approach can be used to improve the selectivity andspecificity of the treatment when dealing simultaneously with more thanone type of cells. With a small adjustment, our procedure can force the solution topreserve the metabolism of one cell while inhibiting an objective reaction of another(see Additional file 1 for a detailed formulation). Indeed,a drug interaction can explore the differences in the topologies of the metabolicnetworks and, in this way, bypass the restrictions caused for instance by targetshomology. This is crucial in case of anticancer therapy since tumoral and normalcells share the same genes.
Solutions for the cancer vs human problem
Floxuridine (# 20)
Mycophenolic acid (# 42)
Trimethoprim (# 29)
Methotrexate (# 11)
Atovaquone (# 69)
Tyloxapol (# 85)
Ezetimibe (# 56)
Pemetrexed (# 41)
Ribavirin (# 51)
Quinacrine (# 36)
Myo-Inositol (# 82)
Tioconazole (# 60)
Naftifine (# 43)
Simvastatin (# 4)
Leflunomide (# 66)
Auranofin (# 59) - Fomepizole (# 75)
Indomethacin (# 22)
Diclofenac (# 35)
Hydroxyurea (# 16)
Arsenic trioxide (# 72)
Gemcitabine (# 30)
In order to increase the range of putative synergisms, we have tried to use thealgorithm to set up a search for targets potentially interacting with the currentlyavailable drugs. Instead of searching for the reactions which are syntheticallylethal (as partially done, for instance, in ) and proposing them as potential new targets, we look for the reactionswhose single inhibition may give a lethal synergism with any combination of theapproved drugs. The problem is formally equivalent to the one we have alreadydescribed: we search again for the optimal drug combination after having deleted areaction in the cancer network (the same reaction is removed also from the humannetwork; this reaction is called “additional target” since it will be thetarget of an additional new drug). We systematically consider each reaction of thecancer network as an additional target; the results of this screening are reported inAdditional file 3: Table S3. After a search in literatureof possible inhibitors of the additional targets of the results, we identify someinteresting solutions. For instance, the inhibition of methenyltetrahydrofolatecyclohydrolase (combined with the use of Mimosine) can be induced by the experimentaldrug 5,6,7,8- tetrahydro-N5,N10- carbonylfolic acid ; therefore, its combination with Mimosine could represent a potentialantitumor therapy. Also for other four additional targets there exist experimentalinhibitors whose activity is reported in the literature [51–53] (details are shown in Additional file 3: TableS3). Concerning the drug solutions identified by the algorithm, the pairRosiglitazone plus Cerulenin is proposed in combination with three possibleadditional targets. One is the palmitate fatty acid conversion: it is worth notingthat the validated antitumoral activity of this pairs versus prostate cancer cells(as mentioned above) is due to the reduction of the synthesis of fatty acids . In our prediction, indeed, among the metabolites which are no longeravailable because of this inhibition, there are cholesterol ester, mono-, di- andtri-acylglycerol. The other two additional targets are related tophosphatidylcholine. Part of the phosphatidyilcholine pathway has been alreadyidentified as synthetic lethal , but without mentioning any possible exploitation. Our results suggest theuse of a combined drug therapy as possible way to take advantage of this syntheticlethality.
The field of drug combinatorics is largely unexplored experimentally and the potentialof combined drug therapies is difficult to assess, mostly for lack of suitablesystematic methodologies. To try to fill this gap, we have developed an algorithm whichis capable of exploring efficiently the optimal synergisms among all possible drugcombinations and of characterizing them in terms of side effect and selectivity. Indeed,the success of a drug discovery process depends on multiple aspects, not least on thefulfillment of requirements regarding selectivity and toxicity: for metabolic diseasesthe modulation of the key pathways without affecting the other vital functions can beinstrumental for rescuing from the pathology. Similarly, anticancer compounds shouldonly kill cancer cells without affecting normal ones. These requirements are rarelytaken into account in standard computational approaches. The results we obtained byapplying our algorithm to the human and tumor vs human metabolic networks showthe possibility to take advantage of drug synergisms in proposing new therapies: thepotentialities lay in the possibility to intervene with a different mechanism of actionwith respect to those that are currently available. In this enlarged repertoire ofpossibilities, we have identified examples of drug repurposing (some of them werepreviously demonstrated experimentally), a procedure which is becoming more and moreattractive thanks to the reduced costs on the preclinical and clinical steps.
One of the main features of FBA-based knockout studies is that metabolic networks appearto be robust [1, 17], meaning that there seem to be an high degree of redundancy of the pathwaysinside a network (property alternatively reported as “nonessentiality” ofthe gene in Refs. [8, 19]). In the context of drug synergism, this property reflects into the presenceof optimal solutions consisting of many drugs (in our case up to four, or even more ifwe consider the number of inhibition targets of each drug). For the same reason, theresults show the necessity to extend the search to all possible drug combinationswithout limiting to those of low cardinality. This fact becomes significant especiallywhen the drugs to combine present a high similarity in terms of inhibited targets;indeed, the characterization of the synergisms we have found shows a limited variety ofpossible interactions between the available drugs (only six classes were identified, seeFigure 2). However, beside increasing the cardinality of thesolutions, a very strong robustness may also reduce the total number of solutionsbecause it makes more difficult to induce the simultaneous inhibition of all theredundant pathways; indeed, if the screening we have run on the human metabolism isapplied also to the less robust cancer network (seen as stand-alone network), the numberof possible inhibitions is much higher (see Additional file 1:Figure S2 in comparison with Figure 3). Moreover, when we want tostop the biomass reaction of the cancer, mainly single drug solutions are found: this isagain an index of the low redundancy of the cancer network (see Table 3) and of the limited variety in the metabolic targets for the availabledrugs.
Nevertheless, adding the possibility to inhibit an extra target, we could identifiedsome experimental compounds (other than the drugs from DrugBank) which may be used asanticancer in a combination with the approved drugs (see Additional file 3: Table S3). These examples show that predictive tools like themethod we are proposing become more important if one considers also the possibility ofcombining active compounds which are not yet approved but for which a minimalcharacterization of the mechanism of action is available. In this perspective,experimental compounds which inhibit additional targets that are different from thoseaffected by the approved drugs may represent a good chance for improving, throughsynergism, the spectrum of the whole set of the currently available drugs. Moreover, theapplication of our method can be extended to situations where multiple networks arecompared and contrasted. It is expected that problems like this will become important assoon as tissue-specific networks of human metabolism and cancer-specific networks willbecome available in the near future.
In a broader perspective, if instead of confining our study only to the human network weconsider also the metabolism of microorganisms, the exploitation of drug synergismsobtained with our algorithm can be useful in investigating a wide range of situations:(i) when specific enzyme inhibitors are not currently available, multiple drug solutionscould represent an example of the reprofiling of existing drugs for newtherapeutic indications; (ii) when the target enzyme has undergone a mutation renderingineffective the original therapy, a synergistic solution may bypass the resistanceacting on other enzymes and therefore help in fighting resistance; (iii) whenthe optimal synergism has no lethal impact while the single drug solution has. Thischange in lethality can be important for instance in cases of human-hostedbacteria producing toxic by-products: in order to save the useful symbiosis with thesecommensal bacteria, a selective (but not lethal) inhibition of the toxic processes mustbe pursued [54, 55].
Apart from the specific cases we have studied, the main objective of this paper is topropose an efficient method to single out drug combinations with potentially interestingtherapeutic effect. Given the exponential character of the combinatorics involved,unguided “fishing expeditions” are intrinsically ill posed. Hence we expectthat methods like the one proposed in this paper should be of help in dealing with sucha complex problem. Needless to say, as the predictions are based on in silicomodels of metabolic pathways, the validity of the results must be assessedexperimentally, but this is beyond the scope of the work.
Then, the optimal synergism problem can be stated as follows:
a metabolic network, which means a stoichiometric matrix and a vector of fluxes v with upper-bounds U , both laying in ;
an objective reaction flux v obj which has to bestopped;
the set of drugs together with their inhibition targets;
we want to find the subset of drugs such that D blocks v obj causing the minimal side effect, i.e. a minimum perturbation on the overallreaction fluxes.
Of course, we are not interested in procedures which perform an exhaustive search in thespace of all drugs combinations.
This reduces the space of feasible steady state fluxes to a subsetH(D) ⊂ H (notice thatH(∅) ≡ H). We can now introduce adefinition which quantifies the side effect, σ(D), of adrug combination D on the reaction fluxes.
Clearly, there are many possible definitions. If we were dealing with microorganisms wecould adapt the analogous MOMA (Minimization of Metabolic Adjustment) or ROOM(Regulatory On/Off Minimization) [20, 21] criteria modeling the perturbation effect (induced by knockout effect, in theliterature) on the fluxes with respect to the “wild type” fluxes. Thisrequires however to know the unperturbed fluxes (the wild type reference) which for amicroorganism corresponds to using the biomass production as cost function of the FBAproblem . Unfortunately, for human metabolic network such a commonly accepted FBAcriterion is unavailable, hence it is not possible to determine the metabolic fluxes forthe unperturbed network (“wild type”). Consequently, also the calculation ofthe perturbed fluxes and the quantification of the side effect are more ambiguous. Ourchoice in this paper is to quantify the side effect as number of stopped reaction (moreprecisely the number of reactions that cannot take place because of the perturbation).This bypasses the lack of the reference fluxes and allows to quantify the number ofcellular functions which are no longer available, regardless to the type oftissue to which the human cell belongs.
where the numerical coefficient captures the spread of the perturbation across the non-metabolicsystems. In analogy with the metabolic network, is equal to the mean number of metabolic reactions that are stoppedwhen a single metabolic target is inhibited. Referring to the human metabolic networkand averaging over all drugs we selected, we obtain .
The optimal solution can be described as follows:
The resulting bilevel optimization is a min-max integer linear programming problem . The inner problem adjusts the fluxes in order to achieve themaximum flow for the objective reaction when all fluxes are subjected to theinhibitions (drugs) imposed by the outer problem and to the stoichiometric constraints.The outer problem selects the combination of drugs which minimizes the sideeffect, restricting to those solutions of the inner problem which guarantee no flow forthe objective reaction.
where the parameter b ≪ 1(b = 0. 001 in this paper) is introduced in theobjective function of the outer problem in order to exclude the combinations containingredundant inhibitions and therefore avoiding an “over-selection” ofdrugs.
the final optimization problem becomes:
where is the upper bound for the dual variable δ i (chosen arbitrarily big in the implementation).
The codes for the algorithm have been developed in MATLAB (MathWorks R2010b) and areavailable in Additional files 4, 5,6, 7, 8and 9. All Mixed Integer Linear Optimizations have beenperformed using the ILOG-IBM CPLEX 12.1, under free academic license.
Selection procedure for the drugs
the whole database contains 6708 drugs;
only approved drugs have been picked out, restricting the search to 1570 drugs;
we filter for drugs which act on human enzymes (identified by the EC number): the set reduces to 473 drugs;
we select only drugs for which an inhibitory effect on at least one enzyme of the human metabolism has been experimentally proven. The set reduces to 267 drugs (with the EC numbers of the inhibited enzymes). For these drugs, we count the number of non-metabolic targets, including also cases of agonism, antagonism and activation, since they all represent a perturbation to the regular functioning of the target (activation is not considered for metabolic targets because it does not affect the number of stopped reaction in our FBA formulation).
the reactions of the metabolic network directly inhibited by each drug are identified through the available correspondence between EC numbers of the inhibited enzymes and the gene codes first, and then through the correspondence between gene codes and metabolic reactions. During this step, it may happen that many genes, and hence many reactions, are associated to the same EC number. Therefore, although in the original database a drug inhibits only a single or a few targets, the number of metabolic reactions affected by the drug can be high. For instance, this is the case of Rosiglitazone which inhibits a single target, the long-chain-fatty-acid-CoA ligase an enzyme responsible for the binding of the acyl-CoA group to a long fatty acid chain. Since the substrate of this enzyme can be any carbon chain, regardless of the unsaturation (presence or not of double bonds C=C) and of the exact length (it just requires a chain longer than 12 carbons atoms), we end up with a drug which inhibits up to 60 metabolic targets.
drugs which have exactly the same metabolic targets are grouped; the final 85 groups are listed in Additional file 1: Table S1. For each group only a representative is reported, namely the drug which has the minimal number of targets outside the metabolism. The full list of drugs in each group is reported in the Additional file 2.
Flux Balance Analysis
Minimization of Metabolic Adjustment.
We would like to thank the anonymous Reviewers for their valuable comments andconstructive suggestions. We would like to thank also Eytan Ruppin, Tomer Shlomi andtheir collaborators for discussions and for providing us the cancer network. Thiswork is supported by a grant from MIUR (Ministero dell’Istruzione,Università e Ricerca)
- Lehar J, Krueger A, Zimmermann G, Borisy A: High-order combination effect and biological robustness. Mol Syst Biol. 2008, 4: 215.-10.1038/msb.2008.51.View ArticleGoogle Scholar
- Chou T: Theoretical basis, experimental design, and computerized simulation of synergismand antagonism in drug combination studies. Pharmacol Rev. 2006, 58 (3): 621-681. 10.1124/pr.58.3.10.View ArticleGoogle Scholar
- Csermely P, Agoston V, Pongor S: The efficiency of multi-target drugs: the network approach might help drugdesign. Trends Pharmacol Sci. 2005, 26: 178-182. 10.1016/j.tips.2005.02.007.View ArticleGoogle Scholar
- Yeh P, Tschumi A, Kishony R: Functional classification of drugs by properties of their pairwise interaction. Nature Gen. 2006, 38 (4): 489-494. 10.1038/ng1755.View ArticleGoogle Scholar
- Keith C, Borisy A, Stockwell B: Multicomponent therapeutics for networked systems. Nat Rev Drug Discov. 2005, 4: 1-8.View ArticleGoogle Scholar
- Yeh P, Kishony R: Networks from drug-drug surface. Mol Syst Biol. 2007, 3: 85-87.View ArticleGoogle Scholar
- Segré D, De Luna A, Church G, Kishony R: Modular epistasis in yeast metabolism. Nat Genet. 2005, 37 (1): 77-83.Google Scholar
- Suthers P, Zomorrodi A, Maranas C: Genome-scale gene/reaction essentiality and synthetic lethality analysis. Mol Syst Biol. 2005, 5: 301-Google Scholar
- Burgard A, Pharkya P, Maranas C: OptKnock: A bilevel programming framework for identifying gene knockout strategiesfor microbial strain optimization. Biotechnol Bioeng. 2003, 84 (6): 647-657. 10.1002/bit.10803.View ArticleGoogle Scholar
- Sridhar P, Song B, Kahveci T, Ranka S: Mining metabolic networks for optimal drug targets. Pac Symp Biocomput. 2008, 13: 291-302.Google Scholar
- Kaelin W: The concept of synthetic lethality in the context of anticancer therapy. Nat Rev Cancer. 2005, 5: 689-698. 10.1038/nrc1691.View ArticleGoogle Scholar
- Goh K, Cusick M, Valle D, Childs B, Vidal M, Barabási A: The human disease network. Proc Natl Acad Sci USA. 2007, 104 (21): 8685-8690. 10.1073/pnas.0701361104.View ArticleGoogle Scholar
- Ferrannini E: Physiological and metabolic consequences of obesity. Metabolism. 1995, 44 (3): 15-17.View ArticleGoogle Scholar
- 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.-10.1038/msb.2011.35.View ArticleGoogle Scholar
- Palsson B: Systems Biology: Properties of Reconstructed Networks. 2006, New York: Cambridge University PressView ArticleGoogle Scholar
- Raman K, Chandra N: Flux balance analysis of biological systems: applications and challenges. Brief Bioinform. 2009, 10 (4): 435-449. 10.1093/bib/bbp011.View ArticleGoogle Scholar
- Deutscher D, Meilijson I, Kupiec M, Ruppin E: Multiple knockout analysis of genetic robustness in the yeast metabolicnetwork. Nat Gen. 2006, 38 (9): 993-998. 10.1038/ng1856.View ArticleGoogle Scholar
- Vazquez A: Optimal drug combinations and minimal hitting sets. BMC Syst Biol. 2009, 3: 81-87. 10.1186/1752-0509-3-81.View ArticleGoogle Scholar
- Imielinski M, Belta B: Deep epistasis in human metabolism. Chaos. 2010, 20 (2): 026104-10.1063/1.3456056.View ArticleGoogle Scholar
- Segré D, Vitkup D, Church G: Analysis of optimality in natural and perturbed metabolic networks. Proc Natl Acad Sci USA. 2002, 99 (23): 15112-15117. 10.1073/pnas.232349399.View ArticleGoogle Scholar
- Shlomi T, Berkman O, Ruppin E: Regulatory on/off minimization of metabolic flux changes after geneticperturbations. Proc Natl Acad Sci USA. 2005, 102 (21): 7695-7700. 10.1073/pnas.0406346102.View ArticleGoogle Scholar
- Snitkin E, Mackay T, Segrè D: Epistatic interaction maps relative to multiple metabolic phenotypes. PLoS Genet. 2011, 7 (2): e1001294-10.1371/journal.pgen.1001294.View ArticleGoogle Scholar
- Salter A: Overview. Trimethoprim-sulfamethoxazole: an assessment of more than 12 years ofuse. Rev Infect Dis. 1982, 4 (2): 196-236. 10.1093/clinids/4.2.196.View ArticleGoogle Scholar
- Weinreb O, Amit T, Bar-Am O, Youdim M: Rasagiline: a novel anti-Parkinsonian monoamine oxidase-B inhibitor withneuroprotective activity. Prog Neurobiol. 2010, 92 (3): 330-344. 10.1016/j.pneurobio.2010.06.008.View ArticleGoogle Scholar
- Duarte N, Becker S, Jamshidi N, Thiele I, et al.: Global reconstruction of the human metabolic network based on genomic andbibliomic data. PNAS. 2007, 104 (6): 1777-1782. 10.1073/pnas.0610772104.View ArticleGoogle Scholar
- Knox C, Law V, Jewison T, Liu P, Ly S, Frolkis A, Pon A, Banco K, Mak C, et al.: DrugBank 3.0: a comprehensive resource for ‘omics’ research ondrugs. Nucleic Acids Res. 2011, 39 (Suppl 1): D1035-D1041.View ArticleGoogle Scholar
- Mansour M, Schwarts D, Judd R, Akingbemi B, Braden T, Morrison E, Dennis J, Bartol F, Hazi A, Napier I, Abdel-Mageed A: Thiazolidinediones/PPARγ agonists and fatty acid synthase inhibitors as anexperimental combination therapy for prostate cancer. Int J Oncol. 2010, 38 (2): 537-546.. 103892/ijo2010877. 103892/ijo2010877 103892/ijo2010877Google Scholar
- Luberto C, Hannun Y: Ceramide in the eukaryotic stress response. Trends Cell Biol. 2000, 10 (2): 73-80. 10.1016/S0962-8924(99)01694-3.View ArticleGoogle Scholar
- Obeid L, Hannun Y: The Ceramide-centric universe of lipid-mediated cell regulation: stress encountersof the lipid kind. J Biol Chem. 2002, 277: 25847-25850. 10.1074/jbc.R200008200.View ArticleGoogle Scholar
- Yang G, Badeanlou L, Bielawski J, Roberts A, Hannun Y, Samad F: Central role of ceramide biosynthesis in body weight regulation, energymetabolism, and the metabolic syndrome. Am J Physiol Endocrinol Metab. 2009, 297: E211—E224-[http://ajpendo.physiology.org/content/297/1/E211.abstract],View ArticleGoogle Scholar
- Triola G, Fabrias G, Dragusin M, Niederhausen L, Broere R, et al.: Specificity of the Dihydroceramide Desaturase InhibitorN-[(1R,2S)-2-Hydroxy-1-hydroxymethyl-2-(2-tridecyl-1-cyclopropenyl)ethyl]octanamide (GT11) in Primary Cultured Cerebellar Neurons. Mol Pharmacol. 2004, 66 (6): 1671-1678. 10.1124/mol.104.003681.View ArticleGoogle Scholar
- Munoz-Olaya J, Matabosch X, Bedia C, Egido-Gabs M, Casas J, et al.: Synthesis and biological activity of a novel inhibitor of DihydroceramideDesaturase. Chem Med Chem. 2009, 3: 946-953.View ArticleGoogle Scholar
- Llaverias G, Lacasa D, Viñals M, Vázquez-Carrera M, Sánchez R, Laguna J, Alegret M: Reduction of intracellular cholesterol accumulation in THP-1 macrophages by acombination of rosiglitazone and atorvastatin. Biochem pharmacol. 2004, 68: 155-163. 10.1016/j.bcp.2004.03.009.View ArticleGoogle Scholar
- Freed M, Ratner R, Marcovina S, Kreider M, Biswas N, Cohen B, Brunzell J: Effects of Rosiglitazone alone and in combination with Atorvastatin on themetabolic abnormalities in type 2 diabetes mellitus. Ame J Cardiology. 2002, 90 (9): 947-952. 10.1016/S0002-9149(02)02659-0.View ArticleGoogle Scholar
- De Berardis G, Sacco M, Evangelista V, Filippi A, Giorda C, Tognoni G, Valentini U, Nicolucci A, et al.: Aspirin and Simvastatin Combination for Cardiovascular Events Prevention Trial inDiabetes (ACCEPT-D): design of a randomized study of the efficacy of low-doseaspirin in the prevention of cardiovascular events in subjects with diabetesmellitus treated with statins. Trials. 2007, 8: 21-10.1186/1745-6215-8-21.View ArticleGoogle Scholar
- Armstrong P, Lorell B, Nissen S, Borer J: 97th Meeting of the Cardiovascular and Renal Drugs Advisory Committee July 18thand 19th, 2002. Circulation. 2002, 106 (6): e9011-e9012.Google Scholar
- Murtola T, Pennanen P, Syvälä H, Bläuer M, Ylikomi T, Tammela T: Effects of simvastatin, acetylsalicylic acid, and rosiglitazone on proliferationof normal and cancerous prostate epithelial cells at therapeuticconcentrations. Prostate. 2009, 69 (9): 1017-1023. 10.1002/pros.20951.View ArticleGoogle Scholar
- Weiss G, Garnick M, Osteen R, Steele G, Wilson R, Schade D, Kaplan W, Boxt L, Kandarpa K, Mayer R: Long-term hepatic arterial infusion of 5-fluorodeoxyuridine for liver metastasesusing an implantable infusion pump. J Clin Oncology. 1983, 1 (5): 337-344.Google Scholar
- Carter S, Franklin T, Jones D, Leonard B, Mills S, Turner R, Turner W: Mycophenolic acid: an anti-cancer compound with unusual properties. Nature. 1969, 223: 848-850. 10.1038/223848a0.View ArticleGoogle Scholar
- Bleich H, Boro E, Frei III E, Jaffe N, Tattersall M, Pitman S, Parker L: New approaches to cancer chemotherapy with methotrexate. New England J Med. 1975, 292 (16): 846-851. 10.1056/NEJM197504172921607.View ArticleGoogle Scholar
- Hazarika M, White R, Johnson J, Pazdur R: FDA drug approval summaries: Pemetrexed (Alimta®). The oncologist. 2004, 9 (5): 482-488. 10.1634/theoncologist.9-5-482.View ArticleGoogle Scholar
- Borden K, Culjkovic-Kraljacic B: Ribavirin as an anti-cancer therapy: acute myeloid leukemia and beyond?. Leukemia Lymphoma. 2010, 51 (10): 1805-10.3109/10428194.2010.496506.View ArticleGoogle Scholar
- Lam S, McWilliams A, MacAulay C, Wattenberg L, Szabo E, et al.: A phase I study of myo-inositol for lung cancer chemoprevention. Cancer Epidemiol Biomarkers Prev. 2006, 15 (8): 1526-1531. 10.1158/1055-9965.EPI-06-0128.View ArticleGoogle Scholar
- Cho S, Kim J, Kim J, Lee J, Jung H, Song I: Simvastatin induces apoptosis in human colon cancer cells and in tumor xenografts,and attenuates colitis-associated colon cancer in mice. Int J Cancer. 2008, 123 (4): 951-957. 10.1002/ijc.23593.View ArticleGoogle Scholar
- Baumann P, Mandl-Weber S, Völkl A, Adam C, Bumeder I, Oduncu F, Schmidmaier R: Dihydroorotate dehydrogenase inhibitor A771726 (leflunomide) induces apoptosis anddiminishes proliferation of multiple myeloma cells. Mol cancer ther. 2009, 8 (2): 366-10.1158/1535-7163.MCT-08-0664.View ArticleGoogle Scholar
- Zhang G, Tu C, Zhang G, Zhou G, Zheng W: Indomethacin induces apoptosis and inhibits proliferation in chronic myeloidleukemia cells. Leukemia res. 2000, 24 (5): 385-392. 10.1016/S0145-2126(99)00198-8.View ArticleGoogle Scholar
- Schrell U, Rittig M, Anders M, Kiesewetter F, Marschalek R, Koch U, Fahlbusch R: Hydroxyurea for treatment of unresectable and recurrent meningiomas. I. Inhibitionof primary human meningioma cells in culture and in meningioma transplants byinduction of the apoptotic pathway. J Neurosurg. 1997, 86 (5): 845-852. 10.3171/jns.1997.86.5.0845.View ArticleGoogle Scholar
- Antman K: Introduction: the history of arsenic trioxide in cancer therapy. The oncologist. 2001, 6 (Suppl 2): 1-2.View ArticleGoogle Scholar
- von der Maase H, Hansen S, Roberts J, Dogliotti L, Oliver T, Moore M, Bodrogi I, Albers P, Knuth A, Lippert C, et al.: Gemcitabine and cisplatin versus methotrexate, vinblastine, doxorubicin, andcisplatin in advanced or metastatic bladder cancer: results of a large,randomized, multinational, multicenter, phase III study. J Clin Oncology. 2000, 18 (17): 3068-3077.Google Scholar
- Temple Jr C, Bennett Jr L, Rose J, Elliott R, Montgomery J: Synthesis of pseudocofactor analogs as potential inhibitors of the folateenzymes. J Med Chem. 1982, 25 (2): 161-166. 10.1021/jm00344a014.View ArticleGoogle Scholar
- Mahfouz M, Johnson S, Holman R: Inhibition of desaturation of palmitic, linoleic and eicosa-8, 11, 14-trienoicacids in vitro by isomeric cis-octadecenoic acids. Biochimica et Biophysica Acta (BBA)-Lipids and Lipid Metabolism. 1981, 663: 58-68. 10.1016/0005-2760(81)90194-6.View ArticleGoogle Scholar
- Borges T, Glauert H, Robertson L: Perfluorodecanoic acid noncompetitively inhibits the peroxisomal enzymes enoyl-CoAhydratase and 3-hydroxyacyl-CoA dehydrogenase. Toxicol Appl Pharmacol. 1993, 118: 8-15. 10.1006/taap.1993.1003.View ArticleGoogle Scholar
- Tserng K, Jin S, Hoppel C: Spiropentaneacetic acid as a specific inhibitor of medium-chain acyl-CoAdehydrogenase. Biochemistry. 1991, 30 (44): 10755-10760. 10.1021/bi00108a021.View ArticleGoogle Scholar
- Mazumdar V, Snitkin E, Amar S, Segré D: Metabolic network model of a human oral pathogen. J Bacteriol. 2009, 191 (1): 74-90. 10.1128/JB.01123-08.View ArticleGoogle Scholar
- Hooper L, Gordon J: Commensal host–bacterial relationships in the gut. Science. 2001, 292: 1115-1118. 10.1126/science.1058709.View ArticleGoogle Scholar
- Palsson B, Varma A: Metabolic capabilities of Escherichia coli II: optimal growth pattern. J Theor Biol. 1993, 165: 503-522. 10.1006/jtbi.1993.1203.View ArticleGoogle Scholar
- Dempe S: Foundations of bilevel programming. 2010, New York: Kluwer Academic PublisherGoogle Scholar
- Matoušek J Gärtner B: Understanding and Using Linear Programming. 2000, Berlin: SpringerGoogle Scholar
- Schrijver A: Theory of Linear and Integer Programming. 1986, New York: John Wiley SonsGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative CommonsAttribution License (http://creativecommons.org/licenses/by/2.0), whichpermits unrestricted use, distribution, and reproduction in any medium, provided theoriginal work is properly cited.