TIde: a software for the systematic scanning of drug targets in kinetic network models
© Schulz et al; licensee BioMed Central Ltd. 2009
Received: 18 June 2009
Accepted: 19 October 2009
Published: 19 October 2009
During the stages of the development of a potent drug candidate compounds can fail for several reasons. One of them, the efficacy of a candidate, can be estimated in silico if an appropriate ordinary differential equation model of the affected pathway is available. With such a model at hand it is also possible to detect reactions having a large effect on a certain variable such as a substance concentration.
We show an algorithm that systematically tests the influence of activators and inhibitors of different type and strength acting at different positions in the network. The effect on a quantity to be selected (e.g. a steady state flux or concentration) is calculated. Moreover, combinations of two inhibitors or one inhibitor and one activator targeting different network positions are analysed. Furthermore, we present TIde (Target Identification), an open source, platform independent tool to investigate ordinary differential equation models in the common systems biology markup language format. It automatically assigns the respectively altered kinetics to the inhibited or activated reactions, performs the necessary calculations, and provides a graphical output of the analysis results. For illustration, TIde is used to detect optimal inhibitor positions in simple branched networks, a signalling pathway, and a well studied model of glycolysis in Trypanosoma brucei.
Using TIde, we show in the branched models under which conditions inhibitions in a certain pathway can affect a molecule concentrations in a different. In the signalling pathway we illuminate which inhibitions have an effect on the signalling characteristics of the last active kinase. Finally, we compare our set of best targets in the glycolysis model with a similar analysis showing the applicability of our tool.
In the current pharmaceutic development new drugs are often found by screening a library of small molecular entities (SME) against so-called 'blockbuster targets' which are supposed to play a relevant role in the onset of a certain disease. The development of drugs for new targets is in most cases less interesting for a pharmaceutical company due to the fact that the research is more expensive, they fail pre-clinical trials more often and are in most cases financially less successful [1, 2]. In order to increase the productivity of 'Research and Development' (R&D) when focussing on novel targets a possible way is to identify candidates which are likely to fail trials earlier in the drug development process .
One problem that drugs against novel targets can cause is their possible lack of efficacy. During the development, possible targets are validated via knock-out experiments which work in a totally different way than medication with competitive inhibitors against the corresponding enzymes. While in the first case the flux through a certain pathway can be completely shut down, in the second case it will only be partially decreased which the system can overcome, e.g. by substrate accumulation or feedback regulation. Therefore quantitative modelling should be incorporated into drug research.
A systematic approach to the identification of possible drug targets in a reaction network renders possible with the established tools and methods used in systems biology. Over the last years more and more mathematical models for chemical reaction networks have proven to be successful in predicting an microorganism's response to changes in its environment and to perturbations in its gene expression [4–6]. These models are being collected in steadily growing databases like BioModels  or JWS online . A promising approach to a systematic drug design is to simulate possible inhibitors to any reaction in a given network and to quantify their effects on a given observable. This observable can be defined as any inner variable of the system, e.g. the concentration of a substance or the flux through a certain reaction which is altered in the pathological state. For a more complex analysis this observable can also be defined as any function of these variables. From the time course of the observable, several characteristics can be extracted for later comparison. These include the steady state values, which are interesting for models of metabolic pathways, and several characteristics, which are relevant to signalling cascade models, namely the integrated concentration, the characteristic time, the signal duration, and the signal amplitude as described in [9, 10]. Given this information from several simulations using different inhibition targets, types, and inhibitor concentrations, one can select favourable modification scenarios. "Favourable" means here that in a certain scenario a given observable reaches a desired value while the system is perturbed by few inhibitors in small concentrations.
Other approaches to determine modified enzyme activities in order to achieve a certain change in a systems behaviour already exist (e.g.  and ), but none of them exhaustively searches the space of possible drug combinations, tries to minimise the amount of used inhibitor, or explicitly models different types of modification kinetics. A broad overview on similar methods can be found in .
In the following, these three steps of the tool's functionality will be explained.
Identification of reaction kinetics
Starting from an SBML model, TIde is supposed to simulate the effects of different modifications of single reactions and combinations thereof. Therefore, the first step of the tool is the identification of the present kinetic formulas for each reaction by numerically comparing them to known kinetic types collected in an internal database. The identity of two formulas is judged by testing all possible parameter matchings of those formulas, repeatedly replacing the parameters by random values, and checking whether the results are similar except for numerical inaccuracy. Using the information about the present kinetics, the database can be searched for known modifier formulas. The identification of the kinetic parameters makes it possible to automatically replace a kinetic by a formula representing a certain kind of modification mechanism acting on the corresponding enzyme.
Replacing reaction kinetics by corresponding modifications
Then, a set of new models is created in which systematically every reaction is replaced by a modification. During this step up to five new models per reaction can be created since four types of modifications are known to the database. These types are competitive, uncompetitive, and noncompetitive inhibition and nonessential activation. Additionally, a competitive inhibition for cofactors is available.
Unknown formulas (v i (S,p) with reaction velocity v i which is dependent on substance concentrations S and parameters p) will automatically be added to the internal database. Also a standard noncompetitive inhibition ( , with inhibitor concentration I and its dissociation constant K i ) and a standard nonessential activation ( ) with activator concentration A and its dissociation constant K a ) will be created. Because of that the number of new models created during this step will be at least twice as large as the number of reactions included in the original model.
In a step including user interaction, also different competitive and uncompetitive inhibitions can be included into the internal database. The user interaction involves the identification of K m values for substrates of reactions which are supposed to be inhibited. For the selected K m values new competitive inhibitions are created by multiplying the original K m by a factor . Uncompetitive inhibitions are created similarly by dividing the K m and V max values by the same factor. For SBML models including complex kinetics this step should be taken as the database can not cover all possible types of kinetic formulas.
TIde can also be used to check combinations of modifiers. These modifications do not need to be of the same type, but they have to affect different reactions.
Comparing results of different modifications
In order to examine the effects of different modifications the user has to define the objective function which is to be maximized or minimized. This function can be an arbitrary formula containing substance concentrations or fluxes or a combination of them if the steady values are of interest (e.g. ). In case that signalling characteristics should be evaluated the objective function can in the current implementation only be a single compound concentration.
The effects of different inhibitor concentrations can be evaluated by either testing several distinct effective concentrations ( or ) in separate simulations or one simulation of a continuous titration of the inhibitor. Optionally these titrations can also be performed until the objective function has reached a certain value.
After the simulations of the new models have been performed making use of either the SciPy  library, the SBML Ode Solver library , or Copasi , the absolute differences between the original value of the objective function and its new values in the modification scenarios are calculated and displayed.
Results and Discussion
We analysed four different models with the TIde tool (using Copasi for simulation) as discussed in the following. The first two models were simple metabolic pathways including a branch, the third was a simple signalling cascade, and the fourth one was a biologically relevant example of aerobic metabolism in Trypanosoma brucei .
The postbranch model (see Figure 2, middle) was a slight modification of the previous model in which the third pathway was not directed to carry biomass into the center of the model but away from it. Here the idea was to study whether this third pathway could be activated in order to compete for biomass from the pathway including species_7.
Simple signalling cascade
Relative changes of the signalling characteristics of "Zp" in the small cascade model.
Model of Trypanosoma glycolysis
Finally, in order to demonstrate the applicability of our tool, we investigated possible modification patterns in an updated model of the aerobic metabolism in Trypanosoma brucei  (BioModels BIOMD0000000211). For an older version of this model a similar analysis has already been performed manually . The updated model contains new and more detailed information on enzyme kinetics under more realistic conditions and therefore these new results could lead to more promising drug targets against sleeping sickness. We have chosen the steady state flux through the upper glycolysis as the variable to be minimized by different modifications.
Since the data produced by TIde was too extensive to display in paperform, only specific results could be shown in Figure 5 and Tables 2 and 3. The favourability of inhibition target and kinetic type over others depended on the degree to which the flux through the glycolysis should be reduced in the parasite (see Figure 5). Therefore, we focussed on inhibitor concentrations decreasing the flux by 50% which was supposed to be sufficient to kill the parasite. As proposed earlier , the most promising target was the glucose transporter (GT). The best targets inside the glycolysis to achieve this goal were the phosphoglycerate mutase (PGM), the glyceraldehyde-3 phosphate dehydrogenase (GAPDH), the glycorol-3-phosphate dehydrogenase (GPDH), the enolase (ENO), and the fructose bisphosphate aldolase (ALD) in this order. As depicted in Table 1, also in this group the quality of a target depended on the type of inhibition. In the case that only competitive inhibitors would be available the GAPDH made an equally potent target as the GT. Earlier results favoured the targets ALD, GAPDH, phosphoglycerate kinase (PGK), and GPDH in this order besides the GT. These findings do not differ that much from the current results except for the PGK. This difference can be explained by different enzyme kinetics, the fact that ENO was not included as a single reaction in the old model, and a slightly different analysis.
Effective inhibitor concentrations necessary to kill Trypanosoma brucei.
Strongest synergisms and antagonisms of combined modification in the Trypanosoma model.
We have introduced an algorithm to systematically test the effect of modifications or modifier combinations depending their target and type as well as TIde, an implementation working on ODE models in the SBML format. This process can be understood as an extension of metabolic control analysis to multiple and strong perturbations in an ODE model. Although TIde seems a bit difficult to use in the first place needing a reliable ODE model, being a command line tool, and having a limited kinetic database the tool still has many advantages. Since TIde is not platform dependend it can be run on any kind of server making it possible to analyse even large models, due to the extendible kinetic database it can be applied to any ODE model, and because of the flexible objective function it can be used for many different purposes. Furthermore, we have tested our tool with a new version of the glycolysis in Trypanosoma brucei. Our results are except for some explainable changes in good agreement with older findings proving the applicability of TIde.
Given a reliable ODE model and the information our tool provides, the process of determining possible drugs is more directed into efficacy of the candidates in pre-clinical studies. Inclusion of this analysis will increase the probability of a candidate to become a potent drug, and decrease the cost of the development of 'new target' drugs. This fact will hopefully draw the attention of the pharmaceutical industry towards the results produced by systems biology. Since this research area is still growing more and more ODE models, which are a necessary input to our method, will become available.
In the future we plan to extend our tool to handle constraints which have to be fulfilled during simulation. The idea behind this is that certain modifications and effective modifier concentrations could lead to the death of the simulated organism (e.g. if the ATP concentration drops below a certain level). Such lethal modifications should be considered in order to identify possible side effects of the potential drug.
Availability and requirements
Project name: TIde
Project home page: http://sysbio.molgen.mpg.de/tide/
Operating systems: Cross platform
Programming language: Python, C
Other requirements: Python 2.5 or higher, SciPy 0.5 or higher, NumPy 1.1 or higher, PyX 0.9 or higher, semanticSBML 1.0 or higher, libSBML 3.3
License: GNU GPL
Any restrictions to use by non-academics: None
Furthermore, we have attached a walkthrough for our tool (see Additional file 1), its source code (see Additional file 2), and a Windows installer for the required tools and libraries (see Additional file 3) to the article.
MS gratefully acknowledges funding and support from the International Max Planck Research School for Computational Biology and Scientific Computing, the German Research Foundation (CRC 618), and the Marie Curie EST project "Systems Biology" (MEST-CT-2004-514169). The work of BMB was funded by NWO Vernieuwingsimpuls and EK is supported by the EU and UNICELLSYS (HEALTH-2007-201142). The funders had no role in study design, data collection, analysis, and interpretation, and decision to submit the manuscript for publication.
- Ma P, Zemmel R: From the analyst's couch: Value of novelty? Nature Reviews Drug Discovery 2002, 1(8):571–572. 10.1038/nrd884View ArticlePubMedGoogle Scholar
- Kola I, Landis J: Can the pharmaceutical industry reduce attrition rates? Nat Rev Drug Discov 2004, 3(8):711–5. 10.1038/nrd1470View ArticlePubMedGoogle Scholar
- Goodall S, Ringel M, Tollman P: Rising to the Productivity Challenge: A Strategic Framework for Biopharma. Boston Consulting Group Report 2004. [http://www.bcg.com/publications/files/rising_to_the_productivity_challenge.pdf]Google Scholar
- Hynne F, Dano S, Sorenson P: Full-scale model of glycolysis in Saccharomyces cerevisiae. Biophysical Chemistry 2001, 94: 121–163. 10.1016/S0301-4622(01)00229-0View ArticlePubMedGoogle Scholar
- Hoefnagel M, Burgt A, Martens D, Hugenholtz J, Snoep J: Time Dependent Responses of Glycolytic Intermediates in a Detailed Glycolytic Model of Lactococcus Lactis During Glucose Run-Out Experiments. Molecular Biology Reports 2002, 29: 157–161. 10.1023/A:1020313409954View ArticlePubMedGoogle Scholar
- Helfert S, Estevez A, Bakker B, Michels P, Clayton C: Roles of triosephosphate isomerase and aerobic metabolism in Trypanosoma brucei. Biochem J 2001, 357(Pt 1):117–125. 10.1042/0264-6021:3570117PubMed CentralView ArticlePubMedGoogle Scholar
- Le Novère N, Bornstein B, Broicher A, Courtot M, Donizelli M, Dharuri H, Li L, Sauro H, Schilstra M, Shapiro B, Snoep J, Hucka M: BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems. Nucleic Acids Research 2006, 34: D689-D691. 10.1093/nar/gkj092PubMed CentralView ArticlePubMedGoogle Scholar
- Olivier B, Snoep J: Web-based kinetic modelling using JWS Online. Bioinformatics 2004, 20(13):2143–2144. 10.1093/bioinformatics/bth200View ArticlePubMedGoogle Scholar
- Heinrich R, Rapoport T: Mathematical analysis of multienzyme systems. II. Steady state and transient control. Biosystems 1975, 7: 130–6. 10.1016/0303-2647(75)90050-7View ArticlePubMedGoogle Scholar
- Heinrich R, Neel B, Rapoport T: Mathematical Models of Protein Kinase Signal Transduction. Molecular Cell 2002, 9(5):957–970. 10.1016/S1097-2765(02)00528-2View ArticlePubMedGoogle Scholar
- Gerber S, Aßmus H, Bakker B, Klipp E: Drug-efficacy depends on the inhibitor type and the target position in a metabolic network-A systematic study. Journal of Theoretical Biology 2008, 252(3):442–455. 10.1016/j.jtbi.2007.09.027View ArticlePubMedGoogle Scholar
- Hornberg J, Bruggeman F, Westerhoff H, Lankelma J: Cancer: A Systems Biology disease. BioSystems 2006, 83(2–3):81–90. 10.1016/j.biosystems.2005.05.014View ArticlePubMedGoogle Scholar
- Hornberg J, Bruggeman F, Binder B, Geest C, de Vaate A, Lankelma J, Heinrich R, Westerhoff H: ERK phosphorylation and kinase/phosphatase control. FEBS Journal 2005, 272: 244–258. 10.1111/j.1432-1033.2004.04404.xView ArticlePubMedGoogle Scholar
- Hornberg J, Westerhoff H: Oncogenes Are to Lose Control on Signaling Following Mutation: Should We Aim Off Target? Molecular Biotechnology 2006, 34(2):109. 10.1385/MB:34:2:109View ArticlePubMedGoogle Scholar
- Bakker B, Michels P, Opperdoes F, Westerhoff H: What controls glycolysis in bloodstream form Trypanosoma brucei. Journal of Biological Chemistry 1999, 274(21):14551–14559. 10.1074/jbc.274.21.14551View ArticlePubMedGoogle Scholar
- Hornberg J, Binder B, Bruggeman F, Schoeberl B, Heinrich R, Westerhoff H: Control of MAPK signalling: from complexity to what really matters. Oncogene 2005, 24: 5533–5542. 10.1038/sj.onc.1208817View ArticlePubMedGoogle Scholar
- Yang K, Bai H, Ouyang Q, Lai L, Tang C: Finding multiple target optimal intervention in disease-related molecular network. Molecular Systems Biology 2008., 4(228):Google Scholar
- Tveito A, Lines G: A note on a method for determining advantageous properties of an anti-arrhythmic drug based on a mathematical model of cardiac cells. Mathematical Biosciences 2009, 217(2):167–173. 10.1016/j.mbs.2008.12.001View ArticlePubMedGoogle Scholar
- Feala J, Cortes J, Duxbury P, Piermarocchi C, McCulloch A, Paternostro G: Systems approaches and algorithms for discovery of combinatorial therapies. Wiley Interdisciplinary Reviews: Systems Biology and Medicine 2009, in press.Google Scholar
- Hucka M, Finney A, Sauro H, Bolouri H, Doyle J, Kitano H, Arkin A, Bornstein B, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles E, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, Juty NS, Kasberger JL, Kremling A, Kummer U, Le Novèere N, Loew LM, Lucio D, Mendes P, Minch E, Mjolsness ED, Nakayama Y, Nelson MR, Nielsen PF, Sakurada T, Schaff JC, Shapiro BE, Shimizu TS, Spence HD, Stelling J, Takahashi K, Tomita M, Wagner J, Wang J: The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics 2003, 19(4):524–531. 10.1093/bioinformatics/btg015View ArticlePubMedGoogle Scholar
- Bornstein B, Keating S, Jouraku A, Hucka M: LibSBML: an API Library for SBML. Bioinformatics 2008, 24(6):880. 10.1093/bioinformatics/btn051PubMed CentralView ArticlePubMedGoogle Scholar
- Le Novère N: Model storage, exchange and integration. BMC Neuroscience 2006, 7(Suppl 1):S11. 10.1186/1471-2202-7-S1-S11PubMed CentralView ArticlePubMedGoogle Scholar
- Lu J, Muller S, Machne R, Flamm C: SBML Ode Solver library: Extensions for inverse analysis. Proceedings of the Fifth International Workshop on Computational Systems Biology, WCSB 2008.Google Scholar
- Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, Singhal M, Xu L, Mendes P, Kummer U: COPASI-a COmplex PAthway SImulator. Bioinformatics 2006, 22(24):3067. 10.1093/bioinformatics/btl485View ArticlePubMedGoogle Scholar
- Albert M, Haanstra J, Hannaert V, Van Roy J, Opperdoes F, Bakker B, Michels P: Experimental and in Silico Analyses of Glycolytic Flux Control in Bloodstream Form Trypanosoma brucei. Journal of Biological Chemistry 2005, 280(31):28306–28315. 10.1074/jbc.M502403200View ArticlePubMedGoogle Scholar
- Schulz M, Uhlendorf J, Klipp E, Liebermeister W: SBMLmerge, a System for Combining Biochemical Network Models. Genome Informatics 2006, 17: 62–71.PubMedGoogle Scholar
This article is published under license to 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.