The silicon trypanosome

Edinburgh Research Explorer The silicon trypanosome Citation for published version: Bakker, BM, Krauth-Siegel, RL, Clayton, C, Matthews, K, Girolami, M, Westerhoff, HV, Michels, PAM, Breitling, R & Barrett, MP 2010, 'The silicon trypanosome', Parasitology, vol. 137, no. 9, pp. 1333-1341. https://doi.org/10.1017/S0031182010000466 Digital Object Identifier (DOI): 10.1017/S0031182010000466 Link: Link to publication record in Edinburgh Research Explorer Document Version: Publisher's PDF, also known as Version of record Published In: Parasitology Publisher Rights Statement: RoMEO green General rights Copyright for the publications made accessible via the Edinburgh Research Explorer is retained by the author(s) and / or other copyright owners and it is a condition of accessing these publications that users recognise and abide by the legal requirements associated with these rights. Take down policy The University of Edinburgh has made every reasonable effort to ensure that Edinburgh Research Explorer content complies with UK legislation. If you believe that the public display of this file breaches copyright please contact [email protected] providing details, and we will remove access to the work immediately and investigate your claim. Download date: 25. Sep. 2019 1333 The silicon trypanosome BARBARA M. BAKKER 1*, R. LUISE KRAUTH-SIEGEL 2, CHRISTINE CLAYTON 3, KEITH MATTHEWS 4, MARK GIROLAMI 5, HANS V. WESTERHOFF 6, PAUL A. M. MICHELS 7, RAINER BREITLING 8 and MICHAEL P. BARRETT 9 1 Department of Pediatrics, Center for Liver, Digestive and Metabolic Diseases, University Medical Center Groningen, University of Groningen, Hanzeplein 1, 9713 GZ Groningen, The Netherlands 2 Biochemie-Zentrum der Universität Heidelberg, 69120 Heidelberg, Germany 3 Zentrum für Molekulare Biologie der Universität Heidelberg, ZMBH-DKFZ Alliance, D69120 Heidelberg, Germany 4 School of Biological Sciences, University of Edinburgh, Ashworth Laboratories, Edinburgh EH9 3JT, United Kingdom 5 University of Glasgow, Department of Computing Science & Department of Statistics, Glasgow, G12 8QQ, United Kingdom 6 Department of Molecular Cell Physiology, VU University Amsterdam, 1081 HV Amsterdam, The Netherlands ; and Manchester Centre for Integrative Systems Biology, Manchester Interdisciplinary BioCentre, The University of Manchester, Manchester M1 7ND, United Kingdom 7 Research Unit for Tropical Diseases, de Duve Institute and Laboratory of Biochemistry, Université catholique de Louvain, Brussels, Belgium and Faculty of Biomolecular and Life Sciences, University of Glasgow, Glasgow, G12 8QQ, UK 8 Groningen Bioinformatics Centre, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, 9751 NN Haren, The Netherlands 9 Faculty of Biomolecular and Life Sciences and Wellcome Centre of Molecular Parasitology, University of Glasgow, Glasgow Biomedical Research Centre, Glasgow G12 8TA, United Kingdom (Received 8 December 2009; revised 17 February 2010; accepted 18 February 2010; first published online 6 May 2010) SUMMARY African trypanosomes have emerged as promising unicellular model organisms for the next generation of systems biology. They offer unique advantages, due to their relative simplicity, the availability of all standard genomics techniques and a long history of quantitative research. Reproducible cultivation methods exist for morphologically and physiologically distinct life-cycle stages. The genome has been sequenced, and microarrays, RNA-interference and high-accuracy metabolomics are available. Furthermore, the availability of extensive kinetic data on all glycolytic enzymes has led to the early development of a complete, experiment-based dynamic model of an important biochemical pathway. Here we describe the achievements of trypanosome systems biology so far and outline the necessary steps towards the ambitious aim of creating a ‘ Silicon Trypanosome ’, a comprehensive, experiment-based, multi-scale mathematical model of trypanosome physiology. We expect that, in the long run, the quantitative modelling enabled by the Silicon Trypanosome will play a key role in selecting the most suitable targets for developing new anti-parasite drugs. Key words: Trypanosoma brucei, metabolism, gene expression, differentiation, Silicon Cell. THE AMBITION OF SYSTEMS BIOLOGY Systems biology seeks to understand how functional properties of living systems, such as biological rhythms, cellular differentiation or the adaptation of organisms to changes in their environment, emerge from interactions between the components in the underlying molecular networks (Bruggeman and Westerhoff, 2007). In the case of parasites with multiple hosts, differentiation and adaptation to drugs may be particularly relevant. Current systems biology is to a large extent (but not exclusively) focused on single-cell systems. These are more amenable to global molecular analysis than multicellular organisms. This is partly because high-throughput post-genomic technologies (transcriptomics, proteomics and metabolomics) make it relatively easy to * Corresponding author : B. M. Bakker, Tel : +31 (0) 50 361 1542. E-mail : [email protected] measure many components of a homogeneous population of cells simultaneously. Furthermore, dynamic measurements of the response of a cell population to a shared stimulus allow insight into the functional connectivity between components (Richard et al. 1996 ; Hynne et al. 2001 ; Nikerel et al. 2006, 2009 ; Schmitz et al. 2009). Mathematical methods that enable quantitative descriptions of the dynamic interplay between the molecules in living cells are being developed and, for the first time, it is possible to envisage a comprehensive molecular description of the functional circuitry of cellular systems. The Silicon Cell project (Snoep et al. 2006 ; Westerhoff et al. 2009) involves an international consortium of researchers aiming at a mathematical description of life at the cellular level on the basis of complete and quantitative genomic, transcriptomic, proteomic, metabolomic and phenotypic information. So far, the most ambitious wholecell modelling efforts have targeted the model Parasitology (2010), 137, 1333–1341. f Cambridge University Press 2010 doi:10.1017/S0031182010000466 1334 Barbara M. Bakker and others organisms Escherichia coli and Saccharomyces cerevisiae. The advanced state of understanding and the enormous amount of data relating to these organisms have made them obvious candidates for such a comprehensive description. Yet, the variety of organisms in the ‘ JWS Online ’ model repository (www.jjj.bio. vu.nl) demonstrates that the silicon cell initiative is not limited to these organisms. THE UNIQUE ADVANTAGES OF TRYPANOSOME SYSTEMS BIOLOGY The African trypanosome, Trypanosoma brucei, the causative agent of human African trypanosomiasis and Nagana cattle disease (Barrett et al. 2003), has emerged as a front runner in systems biology analysis. The relative simplicity of the energy metabolism of its bloodstream form and the early availability of a comprehensive and uniform set of kinetic data of the enzymes involved, were crucial factors for the successful construction of a detailed computer model of trypanosome glycolysis (Bakker et al. 1997). The obvious questions for this model were initially in the realm of drug-target selection ; the first studies analysed in depth how sensitive the pathway overall would be to varying extents of inhibition of each enzyme (Bakker et al. 1999 a). Another important factor stimulating further development of trypanosome systems biology was the possibility of reproducible in vitro cultivation, first of the procyclic insect stage, but later also of the long slender bloodstream form (Hirumi and Hirumi, 1989 ; Haanstra, 2009). Transitions between distinct life-cycle stages can be studied in a tractable and synchronous differentiation system (Fenn and Matthews, 2007). Moreover, the complete genome of T. brucei has been sequenced and annotated, and a metabolic pathway database has been developed (Berriman et al. 2005 ; Chukualim et al. 2008). An extension of the scope of trypanosome systems biology to include gene expression is facilitated by the absence of transcriptional regulation in trypanosomes. This implies that, unlike in most other organisms, the gene-expression cascade is regulated only post-transcriptionally. The genes of African trypanosomes – as well as those of the other kinetoplastids – are arranged in polycistronic transcription units which can be hundreds of kilobases long (Berriman et al. 2005 ; Siegel et al. 2009). All evidence so far indicates that RNA polymerase II transcribes constitutively, without intervention of regulatory factors (Lee et al. 2009 ; Palenchar and Bellofatto, 2006). Individual mRNAs are excised by a transsplicing complex which places identically capped 39 nt leaders at the 5k end of every mRNA (Liang et al. 2003); this splicing is co-ordinated with polyadenylation of the RNA located immediately upstream. Indeed, regulation of mRNA biogenesis may well be restricted to the processing steps (Lustig et al. 2007 ; Stern et al. 2009), while steady-state levels are further influenced by the rate of mRNA degradation (Clayton and Shapira, 2007). In fact, the majority of evidence concerning regulation of gene expression has implicated mRNA decay as the dominant factor (Clayton and Shapira, 2007 ; Haanstra et al. 2008 b), and this is the only step for which mechanistic details of the regulation are available. After the complete sequencing of the trypanosome genome (Berriman et al. 2005), mRNA microarray analyses of the differentiation from the bloodstream to the procyclic form have demonstrated that the expression of whole sets of mRNAs is coordinately regulated (Queiroz et al. 2009). When gene-expression is studied during synchronous differentiation, accurate time profiles of extremely homogenous cell populations can be obtained (Kabani et al. 2009). Most results so far suggest that such regulation is mediated by RNA-binding proteins that bind to specific sequences in the 3kuntranslated regions of mRNAs (Archer et al. 2009 ; Clayton and Shapira, 2007 ; Estévez, 2008). The rate of mRNA translation and protein turnover are other factors influencing the steady-state protein levels. Here too, key regulatory proteins have been identified (Paterou et al. 2006 ; Walrad et al. 2009). MILESTONES OF TRYPANOSOME SYSTEMS BIOLOGY A quantitative mathematical model of energy metabolism in the long slender form of the trypanosome (i.e. the form that replicates in the mammalian bloodstream) has been developed (Bakker et al. 1997) and iteratively updated after experimental testing (Bakker et al. 1999 a, b ; Albert et al. 2005 ; Haanstra et al. 2008 a). This model yields quantitative predictions of the flux through glycolysis, the concomitant ATP production flux, and the concentrations of glycolytic metabolites, at steady state as well as following a perturbation. Input data for the model are kinetic equations and parameters of enzymes and their concentrations. Through this model, the effects of drugs on the glycolytic pathway can be assessed quantitatively starting from their effects on the individual enzymes. Free-energy metabolism in the bloodstream form of T. brucei has been a logical starting point for the ‘ bottom-up ’ construction of a ‘ Silicon Trypanosome ’. Bloodstream forms of the parasite depend exclusively on substrate-level phosphorylation for ATP production through glycolysis, which proceeds as far as pyruvate (Flynn and Bowman, 1973) (Fig. 1). Pyruvate is the end product and is secreted from the cell. Many of the glycolytic enzymes differ in terms of allosteric regulation from their mammalian counterparts, and this probably relates to the fact that in T. brucei the first seven enzymes of the pathway reside within membrane-bounded, peroxisome-like The silicon trypanosome Fig. 1. The glycolytic pathway in Trypanosoma brucei. Reaction numbers indicate : 1. glucose transport ; 2. hexokinase ; 3. phosphoglucose isomerase ; 4. phosphofructokinase ; 5. aldolase; 6. triose-phosphate isomerase ; 7. glyceraldehyde-3-phosphate dehydrogenase ; 8. phosphoglycerate kinase ; 9. phosphoglycerate mutase; 10. enolase; 11. pyruvate kinase ; 12. pyruvate transport ; 13. glycerol-3-phosphate dehydrogenase ; 14. glycerol-3-phosphate oxidase (a combined process of mitochondrial glycerol-3phosphate dehydrogenase and trypanosome alterative oxidase ; 15. glycerol kinase ; 16. combined ATP utilization ; 17. glycosomal adenylate kinase ; 18. cytosolic adenylate kinase. Question marks indicate uncharacterized transport processes. Abbreviations of metabolite names : Glc-6-P: glucose 6-phosphate ; Fru-6-P : fructose 6-phosphate ; Fru-1,6-BP : fructose 1,6-bisphosphate ; DHAP : dihydroxyacetone phosphate ; Gly-3-P : glycerol 3-phosphate ; GA-3-P: glyceraldehyde 3-phosphate ; 1,3-BPGA : 1,3-bisphosphoglycerate ; 3-PGA : 3-phosphoglycerate ; 2-PGA : 2-phosphoglycerate ; PEP : phospho-enolpyruvate. organelles called glycosomes (Opperdoes and Borst, 1977 ; Parsons, 2004 ; Michels et al. 2006 ; Haanstra et al. 2008), which isolate most of the glycolytic pathway from the rest of the metabolic network involved in consumption of ATP and NAD(H). Even in growing and dividing trypanosomes, virtually all glucose is converted to pyruvate, as the amount required for biosynthesis is quantitatively negligible (Haanstra, 2009). This finding supports the initial choice to model the glycolytic pathway without any 1335 branches other than the one to glycerol. Glycerol production is crucial under anaerobic conditions (Fairlamb et al. 1977). Since the publication of the first version of the glycolysis model (Bakker et al. 1997), there have been two major updates (Helfert et al. 2001 ; Albert et al. 2005). Both of these involved updates and extensions of the enzyme kinetic dataset, e.g. the explicit inclusion of individual enzymes that were previously grouped into a net multi-step conversion. In the second update (Albert et al. 2005) the enzyme expression levels (Vmax) were adapted to reflect the concentrations observed in trypanosomes obtained from controlled state-of-the-art in vitro cultivation. Key missing pieces of information remain the mechanism and kinetics of the transport of glycolytic metabolites across the glycosomal membrane. The identification of semi-selective pores in peroxisomal membranes suggests that the smaller metabolites equilibrate across the glycosomal membrane, while bulkier molecules like ATP or NADH require specific transporters (Grunau et al. 2009 ; Rokka et al. 2009). This idea justifies, with hindsight, the choice to model the transport of a number of small intermediates as rapid-equilibrium processes. A number of basic and applied biological questions have been addressed using the glycolysis model. For example, it was predicted and then experimentally confirmed (Bakker et al. 1999 a, b) that the uptake of glucose across the plasma membrane was a major flux controlling step and therefore an interesting drug target. Enzymes that have been suggested to control glycolysis in mammalian cells, like hexokinase, phosphofructokinase and pyruvate kinase (Schuster and Holzhütter, 1995), exerted little control in trypanosomes, according to the model (Bakker et al. 1999 a ; Albert et al. 2005). Experiments, in which the expression of these enzymes was knocked down, confirmed this prediction qualitatively. However, the enormous overcapacity of some enzymes, which was predicted by the model, was shown to be exaggerated (Albert et al. 2005). This suggests that there are in vivo regulation mechanisms affecting these enzymes in a currently unknown fashion. Protein phosphorylation may contribute, since a glycosomal phosphatase has been identified in developmental signalling (Szoor and Matthews, unpublished data). The inhibition of anaerobic glycolysis by glycerol was also reproduced by the model, first qualitatively and then quantitatively (Bakker et al. 1997 ; Albert et al. 2005). An interesting biological feature that was revealed by the model was the relationship between compartmentation of glycolysis in glycosomes and the virtual absence of allosteric regulation of the glycolytic enzymes. Glycolysis models predict that glycolytic intermediates accumulate readily due to the investment of ATP at the beginning of the pathway (Teusink et al. 1998 ; Bakker et al. 2000). This risky ‘ turbo ’ effect can be avoided either by allosteric feedback 1336 Barbara M. Bakker and others Fig. 2. A. The positive feedback from the ATP produced by glycolysis to the initial kinase reactions can lead to toxic accumulation of hexose phosphates. In many organisms this is prevented by a negative feedback from the hexose phosphates to hexokinase. B. In trypanosomes, the compartmentation of glycolysis prevents the positive feedback. This renders the negative feedback unnecessary, and indeed there is no evidence for such feedback in trypanosomes. regulation of hexokinase or by compartmentation of the pathway in glycosomes (Fig. 2). Compartmentation prevents the accumulation of intermediates, because the net ATP production occurs outside the glycosome and this excess of ATP cannot activate the first enzymes of glycolysis. This model prediction was recently confirmed experimentally (Haanstra et al. 2008 a), providing a clear example of model-driven experimental design and hypothesis-driven systems biology. According to model predictions the glycolytic intermediates glucose 6phosphate, fructose 6-phosphate and fructose1,6-bisphosphate should accumulate on addition of glucose if the glycolytic enzymes are not properly located in the glycosome. Indeed, accumulation of glucose 6-phosphate could be measured in a PEX14RNAi mutant in which protein import into the glycosomes is disturbed. A similar phenotype was observed on glycerol addition, which led to accumulation of glycerol 3-phosphate, both in the model and in the PEX14-RNAi cells. Also in accordance with model predictions, a down-regulation of the expression of the genes encoding hexokinase and glycerol kinase rescues the PEX14-RNAi cells on glucose and glycerol, respectively (Kessler and Parsons, 2005 ; Haanstra et al. 2008 a). More recently, a model of the gene-expression cascade, based on quantitative knowledge of transcription, RNA precursor degradation, trans-splicing and mRNA degradation for phosphoglycerate kinase (PGK) has been generated (Haanstra et al. 2008 b). The model allowed a quantitative analysis of the control and regulation of the expression of the PGK isoenzymes. It was shown that regulation of mRNA degradation explains 80–90 % of the regulation of mature mRNA levels, while precursor degradation and trans-splicing make only minor contributions. In spite of the success of the model, it covers to date only a small part of trypanosome metabolism. This relates, for instance, to the fact that even the compartmentalised glycolysis does branch into other pathways, for example towards the biosynthesis of glycoconjugates and the pentose phosphate pathway. Although the fluxes into these branches may be small, they are vital for trypanosomes. Sufficient kinetic data have become available to enable extension of the model to include the pentose phosphate pathway which provides NADPH for reductive biosyntheses and also reducing equivalents to sustain cellular redox balance. Since redox balance is intimately related to the biosynthesis of trypanothione (from polyamine and glutathione precursors), a natural next step in a bottom-up systems biology approach to trypanosome metabolism would be the inclusion of the trypanothione–pentose phosphate pathway and related areas of redox metabolism (Fig. 3). GROWTH STAGES OF BUILDING A SILICON TRYPANOSOME Our current level of knowledge of trypanosome redox metabolism, as well as its biological importance (Krauth-Siegel and Comini, 2008), render it a natural choice for a next model extension (Fig. 3). The inclusion of redox metabolism is particularly The silicon trypanosome 1337 Fig. 3. The glycolytic and trypanothione pathways are linked through the oxidative pentose phosphate pathway. Metabolites are presented in abbreviated form within rectangles. Enzymes and transporters are circles. Kinetic data is available for those shaded grey. Met=methionine ; Arg=arginine ; Gly=glycine ; Glu=glutamate; Cys=cysteine ; ATP=adenosine triphosphate ; SAM=S-adenosylmethionine ; dcSAM=decarboxylated S-adenosylmethionine ; MTA=methylthioadenosine ; Orn=ornithine ; Put=putrescine ; Spd=spermidine ; c-GC=c-glutamylcysteine ; GSH=glutathione ; GSpd=glutathionylspermidine ; T(SH2)=reduced trypanothione ; TS2=oxidised trypanothione ; NADP=nicotinamide adenine dinucleotide phosphate ; NADPH=reduced nicotinamide adenine dinucleotide phosphate ; Glc=glucose ; G6P=glucose 6-phosphate ; Ru5P=ribulose 5-phosphate ; Ox=oxidised cellular metabolites ; Red=reduced cellular metabolites. 1.=methionine transport ; 2.=arginine transport ; 3.=glycine transport ; 4.=glutamate transport ; 5=cysteine transport ; 6.=methionine adenosyltransferase ; 7.=arginase (N.B., a robust arginase gene orthologue has not been annotated in the T. brucei genome project, raising the possibility that arginine does not serve as a source of ornithine in these cells) ; 8.=S-adenosylmethionine decarboxylase ; 9.=prozyme ; 10.=ornithine decarboxylase ; 11.=spermidine synthase ; 12.=c-glutamylcysteine synthetase ; 13.=glutathione synthetase ; 14.=c-glutathionylspermidine synthetase ; 15.=trypanothione synthetase/amidase (in T. brucei 14. & 15. are catalysed by a single protein) ; 16.=trypanothione reductase ; 17. oxidative pentose phosphate pathway (glucose 6-phosphate dehydrogenase, 6-phopshogluconolactonase & 6-phosphoglconate dehydrogenase) ; 18.=glucose transporter ; 19.=hexokinase (this enzyme links the redox pathway to glycolysis) ; 20.=pyruvate transporter ; 21.=methionine cycle enzymes ; 22. The pathway of electrons from reduced trypanothione for final acceptance on oxidised cellular metabolites or macromolecules is complex, involving transfers via other redox active intermediates including tryparedoxin (thioredoxin-like proteins) and peroxidoxin. interesting as trypanosome redox metabolism is sufficiently different from its human counterpart to offer perspectives for drug discovery. The unusual polyamine–glutathione conjugate trypanothione or bis(glutathionyl)spermidine (Fairlamb and Cerami, 1992) takes on the majority of roles served by glutathione in most other cell types. In addition, work in the last few years revealed that the enzymes involved in the synthesis and reduction of trypanothione are essential for the parasite (Krauth-Siegel and Comini, 2008). The trypanocidal drug eflornithine exerts its trypanocidal activity as an irreversible inhibitor of the enzyme ornithine decarboxylase (Bacchi et al. 1980), an enzyme involved in trypanothione biosynthesis (enzyme 10 in Fig. 3). A significant amount of information is available on kinetic parameters of that pathway, too. Preliminary attempts to model trypanothione metabolism have been made (Xu Gu, University of Glasgow PhD thesis, unpublished). Information available on the abundance of key metabolites measured in bloodstream form T. brucei grown in vitro (Fairlamb et al. 1987) and in vivo (Xiao et al. 2009), before and after exposure to eflornithine, was used to determine whether predicted behaviour under those perturbed conditions emulated the measured behaviour. The scarcity of kinetic data describing the whole pathway, however, has presented many challenges to constructing a model that captures observed behaviour. The acquisition of new kinetic data and the implementation of new mathematical tools to fill gaps in the data (Nikerel et al. 2006 ; Smallbone et al. 2007 ; Resendis-Antonio, 2009) should improve this. An extension of the glycolysis model to include the pentose phosphate pathway (Hanau et al. 1996 ; Barrett, 1997 ; Duffieux et al. 2000) and trypanothione metabolism should be a suitable next step in the modular approach that we envisage towards a complete Silicon Trypanosome. Initial efforts in this direction (not published) have indicated the importance of the compartmentation of the pentose phosphate pathway. Although most of the enzymes of the pathway have a peroxisome targeting 1338 Barbara M. Bakker and others sequence (PTS1), a significant fraction of their activity is often found in the cytosol (Michels et al. 2006 ; Heise and Opperdoes, 1999 ; Duffieux et al. 2000). A correct localisation of the enzymes as well as good estimates of the transport of intermediates across the glycosomal membrane will be key to a good model of the pentose phosphate pathway. CHALLENGES OF TRYPANOSOME SYSTEMS BIOLOGY The first initiatives in systems biology of trypanosomes as well as of other organisms dealt with enzymatic sub-systems, such as glycolysis. The models have depended largely on kinetic data for isolated enzymes. However, the abundance of these enzymes can, in principle, be controlled by the rates of transcription, RNA processing, translation, protein modification and turnover. These processes themselves may be regulated through complex signalling networks in response to both internal and external conditions (Westerhoff et al. 1990). To include gene expression in a Silicon Trypanosome requires a dramatic increase in model complexity – as well as the acquisition of new types of data on a large scale. Fortunately, the absence of transcriptional control of most individual open reading frames makes trypanosome gene expression simpler than that of yeast or even E. coli, rendering it uniquely amenable to mathematical modelling. It may well be possible to describe much of trypanosome mRNA metabolism using the following parameters : the rate constant of processing of the precursor RNA, i.e. of trans-splicing ; the rate constant of degradation of the precursor (which competes with its trans-splicing) ; the rate constant of polyadenylation; and the rate constant of mRNA degradation. The rates of degradation of the precursor and the mature mRNA can be measured by inhibiting splicing and transcription. To measure the rate of mRNA processing two approaches are possible. First, one can inhibit transcription alone, and assay precursor decay ; this approach is, however, compromised by practical constraints since splicing is very rapid. Second, the splicing rate can be calculated based on the steady-state abundance of the precursor mRNA, and the half-life and abundance of the products. This methodology has already been applied to the mRNA encoding PGK and it was demonstrated that splicing occurred within less than one minute ; mRNA decay was the primary determinant of mRNA abundance (Haanstra et al. 2008 b). Previous microarray studies with yeast have yielded estimates of the half-lives and polysomal loading of many mRNAs (e.g. Grigull et al. 2004). Deep sequencing technology – being more sensitive – should allow measurement of the abundances of all mRNAs and precursors on a genome-wide scale and to the accuracy required for the modelling ; from these data, it should be possible to derive the steadystate abundances and half-lives of all RNAs, revealing regulated or inefficient processing. This – combined with global polysome profiling – will provide quantitative data which allow quantifying the regulation of the processing, degradation and translation of each mRNA (Daran-Lapujade et al. 2007). The next challenge would then be to integrate such measurements with metabolic modelling in order to provide a complete model of pathways, from DNA to metabolic end-products. ANTICIPATED OUTCOMES FROM A SILICON TRYPANOSOME So far, the systems biology approach to trypanosomes has contributed to a fundamental understanding of cellular regulation (Bakker et al. 1999 a ; Haanstra et al. 2008 a), as well as to improvements in the drugtarget selection process (Bakker et al. 1999 a, b ; Albert et al. 2005 ; Hornberg et al. 2007). Since the initial systems biology analysis only addressed processes associated with less that 1 % of the organism’s genome, we would expect many more new insights to lie ahead. Drugs currently used against human African trypanosomiasis are unsatisfactory for a number of reasons, including their extreme adverse effects in the patient and the emergence of resistant parasites. New drugs are urgently needed and there is hope that a better understanding of the control points of the metabolic network can guide the selection of optimal drug targets. This has already been achieved for enzymes of the glycolytic pathway (Hornberg et al. 2007). This information can be used alongside enhanced chemoinformatics (Frearson et al. 2007) in order to determine those components of the trypanosome that are most amenable to drug targeting. As a consortium we have embarked on the construction of a Silicon Trypanosome. In this review we have discussed the current status and future directions of trypanosome systems biology that form the context of this endeavour. Our ambition is to achieve a comprehensive, quantitative description of the flow of information from gene, through transcript and protein, to metabolism and back. This will allow prediction of how the parasite responds to changes in its environment, with respect to nutrients, temperature and/or chemical inhibitors. It will also assist the deciphering of complex phenotypes generated by genetic perturbations in the laboratory or in the field. Thus, model predictions will improve our biological understanding of the differentiation and adaptation of the parasite as well as stimulate the discovery of inhibitors that attack processes which control trypanosome growth. The latter should contribute to the development of new optimised drugs for trypanocidal chemotherapy. Pioneering efforts The silicon trypanosome have focused on energy metabolism and recently started to include adaptations of the parasite via gene expression (Haanstra, 2009). The construction of a complete Silicon Trypanosome, which integrates metabolism, gene expression and signal transduction is an ambitious project. Clearly the route towards this objective will be long, and many challenges will emerge as the datasets required to build such a model are collected and analysed. However, the emergence of methods to allow collection of massive datasets, at every level, suggests that we may, in time, be able to generate a reasonably complete mathematical description of trypanosome cellular biology. Even if completion is not feasible, the evolving description will always represent the best conceivable dynamic representation of our knowledge of trypanosome biology. As a result, drug development programmes will have at their disposal a predictive model of the trypanosome to help identify those parts of metabolism most amenable to targeting by novel drugs and to controlling vital functions of the parasite. The project will be strengthened by parallel world-wide systems biology projects of human metabolism, in which some of us will be involved. After all, killing trypanosomes is easy. The difficulty is to kill the trypanosome without harming its host (Bakker et al. 2002). A careful comparison of the behaviour of our Silicon Trypanosome to quantitative knowledge of the control of human metabolism, will allow the identification of selective targets. ACKNOWLEDGEMENTS The work of BMB was funded by NWO Vernieuwingsimpuls and by a Rosalind Franklin Fellowship. RB was supported by an NWO-Vidi fellowship. HVW thanks BBSRC and EPSC for support through the MCISB grant (http://www.systembiology.net/ support/ ). MPB is grateful to the BBSRC for their support of the BBSRC-ANR ‘‘ Systryp ’’ consortium. The Silicon Trypanosome consortium is supported by a grant from SysMO2 (www.sysmo.net). REFERENCES Albert, M. A., Haanstra, J. R., Hannaert, V., Van Roy, J., Opperdoes, F. R., Bakker, B. M., and Michels, P. A. M. (2005). Experimental and in silico analyses of glycolytic flux control in bloodstream form Trypanosoma brucei. Journal of Biological Chemistry 280, 28306–28315. Archer, S. K., Luu, V.- D., de Queiroz, R., Brems, S. and Clayton, C. E. (2009). Trypanosoma brucei PUF9 regulates mRNAs for proteins involved in replicative processes over the cell cycle. PLoS Pathogens 5, e1000565. Bacchi, C. J., Nathan, H. C., Hutner, S. H., McCann, P. P. and Sjoerdsma, A. (1980). Polyamine metabolism : a potential therapeutic target in trypanosomes. Science 210, 332–334. Bakker, B. M., Aßmus, H. E., Bruggeman, F., Haanstra, J., Klipp, E. and Westerhoff, H. V. (2002). 1339 Network-based selectivity of antiparasitic inhibitors. Molecular Biology Reports 29, 1–5. Bakker, B. M., Mensonides, F. I. C., Teusink, B., Michels, P. A. M. and Westerhoff, H. V. (2000). Compartmentation protects trypanosomes from the dangerous design of glycolysis. Proceedings of the National Academy of Sciences, USA 97, 2087–2092. Bakker, B. M., Michels, P. A. M., Opperdoes, F. R. and Westerhoff, H. V. (1997). Glycolysis in bloodstream form Trypanosoma brucei can be understood in terms of the kinetics of the glycolytic enzymes. Journal of Biological Chemistry 272, 3207–3215. Bakker, B. M., Michels, P. A. M., Opperdoes, F. R. and Westerhoff, H. V. (1999 a). What controls glycolysis in bloodstream form Trypanosoma brucei? Journal of Biological Chemistry 274, 14551–14559. Bakker, B. M., Walsh, M. C., ter Kuile, B. H., Mensonides, F. I., Michels, P. A. M., Opperdoes, F. R. and Westerhoff, H. V. (1999 b). Contribution of glucose transport to the control of the glycolytic flux in Trypanosoma brucei. Proceedings of the National Academy of Sciences, USA 96, 10098–10103. Barrett, M. P. (1997). The pentose phosphate pathway and parasitic protozoa. Parasitology Today 13, 11–16. Barrett, M. P., Burchmore, R. J., Stich, A., Lazzari, J. O., Frasch, A. C., Cazzulo, J. J. and Krishna, S. (2003). The trypanosomiases. Lancet 362, 1469–1480. Berriman, M., Ghedin, E., Hertz-Fowler, C., Blandin, G., Renauld, H. et al. (2005). The genome of the African trypanosome Trypanosoma brucei. Science 309, 416–422. Bruggeman, F. J. and Westerhoff, H. V. (2007). The nature of systems biology. Trends in Microbiology 15, 45–50. Chukualim, B., Peters, N., Hertz-Fowler, C. and Berriman, M. (2008). TrypanoCyc – a metabolic pathway database for Trypanosoma brucei. BMC Bioinformatics 9 (Suppl 10), P5. Clayton, C. and Shapira, M. (2007). Post-transcriptional regulation of gene expression in trypanosomes and leishmanias. Molecular and Biochemical Parasitology 156, 93–101. Daran-Lapujade, P., Rossell, S., van Gulik, W. M., Luttik, M. A. H., de Groot, M. J. L., Slijper, M., Heck, A. J. R., Daran, J. M., de Winde, J. H., Westerhoff, H. V., Pronk, J. T. and Bakker, B. M. (2007). The fluxes through glycolytic enzymes in Saccharomyces cerevisiae are predominantly regulated at posttranscriptional levels. Proceedings of the National Academy of Sciences, USA 104, 15753–15758. Duffieux, F., Van Roy, J., Michels, P. A. M., Opperdoes, F. R. (2000). Molecular characterization of the first two enzymes of the pentose-phosphate pathway of Trypanosoma brucei. Glucose-6-phosphate dehydrogenase and 6-phosphogluconolactonase. Journal of Biological Chemistry 275, 27559–27565. Estévez, A. (2008). The RNA-binding protein TbDRBD3 regulates the stability of a specific subset of mRNAs in trypanosomes. Nucleic Acids Research 36, 4573–4586. Fairlamb, A. H. and Cerami, A. (1992). Metabolism and functions of trypanothione in the Kinetoplastida. Annual Review of Microbiology 46, 695–729. Barbara M. Bakker and others Fairlamb, A. J., Opperdoes, F. R. and Borst, P. (1977). New approach to screening drugs for activity against African trypanosomes. Nature 265, 270–271. Fairlamb, A. H., Henderson, G. B., Bacchi, C. J. and Cerami, A. (1987). In vivo effects of difluoromethylornithine on trypanothione and polyamine levels in bloodstream forms of Trypanosoma brucei. Molecular and Biochemical Parasitology 24, 185–191. Fenn, K. and Matthews, K. R. (2007). The cell biology of Trypanosoma brucei differentiation. Current Opinion in Microbiology 10, 539–546. Flynn, I. W. and Bowman, I. B. (1973). The metabolism of carbohydrate by pleomorphic African trypanosomes. Comparative Biochemistry and Physiology B 45, 25–42. Frearson, J. A., Wyatt, P. G., Gilbert, I. H. and Fairlamb, A. H. (2007). Target assessment for antiparasitic drug discovery. Trends in Parasitology 23, 589–595. Grigull, J., Mnaimneh, S., Pootoolal, J., Robinson, M. and Hughes, T. (2004). Genome-wide analysis of mRNA stability using transcription inhibitors and microarrays reveals posttranscriptional control of ribosome biogenesis factors. Molecular and Cellular Biology 24, 5534–5547. Grunau, S., Mindthoff, S., Rottensteiner, H., Sormunen, R. T., Hiltunen, J. K., Erdmann, R. and Antonenkov, V. D. (2009) Channel-forming activities of peroxisomal membrane proteins from the yeast Saccharomyces cerevisiae. FEBS Journal 276, 1698–1708. Haanstra, J. (2009). The power of network-based drug design and the interplay between metabolism and gene expression in Trypanosoma brucei. PhD thesis Vrije Universiteit Amsterdam, ISBN 978-90-8659331-6. Haanstra, J. R., Stewart, M., Luu, V. D., van Tuijl, A., Westerhoff, H. V., Clayton, C. and Bakker, B. M. (2008b). Control and regulation of gene expression : quantitative analysis of the expression of phosphoglycerate kinase in bloodstream form Trypanosoma brucei. Journal of Biological Chemistry 283, 2495–2507. Haanstra, J. R., van Tuijl, A., Kessler, P., Reijnders, W., Michels, P. A. M., Westerhoff, H. V., Parsons, M. and Bakker, B. M. (2008 a). Compartmentation prevents a lethal turbo-explosion of glycolysis in trypanosomes. Proceedings of the National Academy of Sciences, USA 105, 17718–17723. Hanau, S., Rippa, M., Bertelli, M., Dallocchio, F. and Barrett, M. P. (1996). 6-Phosphogluconate dehydrogenase from Trypanosoma brucei. Kinetic analysis and inhibition by trypanocidal drugs. European Journal of Biochemistry 240, 592–599. Heise, N. and Opperdoes, F. R. (1999). Purification, localisation and characterisation of glucose-6-phosphate dehydrogenase of Trypanosoma brucei. Molecular and Biochemical Parasitology 99, 21–32. Helfert, S., Bakker, B. M., Michels, P. A. M. and Clayton, C. (2001). An essential role of triosephosphate isomerase and aerobic metabolism in trypanosomes. Biochemical Journal 357, 117–125. Hirumi, H. and Hirumi, K. (1989). Continuous cultivation of Trypanosoma brucei blood stream forms in 1340 a medium containing a low concentration of serum protein without feeder cell layers. Journal of Parasitology 75, 985–989. Hornberg, J. J., Bruggeman, F. J., Bakker, B. M. and Westerhoff, H. V. (2007). Metabolic control analysis to identify optimal drug targets. Progress in Drug Research 64, 173–189. Hynne, F., Danø, S. and Sørensen, P. G. (2001). Full-scale model of glycolysis in Saccharomyces cerevisiae. Biophysical Chemistry 94, 121–163. Kabani, S., Fenn, K., Ross, A., Ivens, A., Smith, T. K., Ghazal, P. and Matthews, K. (2009). Genome-wide expression profiling of in vivo-derived bloodstream parasite stages and dynamic analysis of mRNA alterations during synchronous differentiation in Trypanosoma brucei. BMC Genomics 10, 427. Kessler, P. S. and Parsons, M. (2005). Probing the role of compartmentation of glycolysis in procyclic form Trypanosoma brucei : RNA interference studies of PEX14, hexokinase and phosphofructokinase. Journal of Biological Chemistry 280, 9030–9036. Krauth-Siegel, R. L. and Comini, M. A. (2008). Redox control in trypanosomatids, parasitic protozoa with trypanothione-based thiol metabolism. Biochimica et Biophysica Acta 1780, 1236–1248 Lee, J. H., Jung, H. S. and Gunzl, A. (2009). Transcriptionally active TFIIH of the early-diverged eukaryote Trypanosoma brucei harbors two novel core subunits but not a cyclin-activating kinase complex. Nucleic Acids Research 37, 3811–3820. Liang, X., Haritan, A., Uliel, S. and Michaeli, S. (2003). Trans and cis splicing in trypanosomatids : mechanism, factors, and regulation. Eukaryotic Cell 2, 830–840. Lustig, Y., Sheiner, L., Vagima, Y., Goldshmidt, H., Das, A., Bellofatto, V. and Michaeli, S. (2007). Spliced-leader RNA silencing : a novel stress-induced mechanism in Trypanosoma brucei. EMBO Reports 8, 408–413. Michels, P. A. M., Bringaud, F., Herman, M. and Hannaert, V. (2006). Metabolic functions of glycosomes in trypanosomatids. Biochimica et Biophysica Acta 1763, 1463–1477. Nikerel, I. E., van Winden, W. A., van Gulik, W. M. and Heijnen, J. J. (2006). A method for estimation of elasticities in metabolic networks using steady state and dynamic metabolomics data and linlog kinetics. BMC Bioinformatics 7, 540. Nikerel, I. E., van Winden, W. A., Verheijen, J. and Heijnen, J. J. (2009). Model reduction and a priori kinetic parameter identifiability analysis using metabolome time series for metabolic reaction networks with linlog kinetics. Metabolic Engineering 11, 20–30. Opperdoes, F. R. and Borst, P. (1977). Localization of nine glycolytic enzymes in a microbody-like organelle in Trypanosoma brucei : the glycosome. FEBS Letters 80, 360–364. Parsons, M. (2004). Glycosomes : parasites and the divergence of peroxisomal purpose. Molecular Microbiology 53, 717–724. Palenchar, J. B. and Bellofatto, V. (2006). Gene transcription in trypanosomes. Molecular and Biochemical Parasitology 146, 135–141. The silicon trypanosome Paterou, A., Walrad, P., Craddy, P., Fenn, K. and Matthews, K. (2006). Identification and stage-specific association with the translational apparatus of TbZFP3, a ccch protein that promotes trypanosome life cycle development. Journal of Biological Chemistry 281, 39002–39013. Queiroz, R., Benz, C., Fellenberg, K., Hoheisel, J. and Clayton, C. (2009). Transcriptome analysis of differentiating trypanosomes reveals the existence of multiple post-transcriptional regulons. BMC Genomics 10, 495. Resendis-Antonio, O. (2009). Filling kinetic gaps : dynamic modeling of metabolism where detailed kinetic information is lacking. PLoS One 4, e4967. Richard, P., Teusink, B., Hemker, M. B., van Dam, K. and Westerhoff, H. V. (1996). Sustained oscillations in free-energy state and hexose phosphates in yeast. Yeast 12, 731–740. Rokka, A., Antonenkov, V. D., Soininen, R., Immonen, H. L., Pirilä, P. L., Bergmann, U., Sormunen, R. T., Weckström, M., Benz, R. and Hiltunen, J. K. (2009). Pxmp2 is a channel-forming protein in Mammalian peroxisomal membrane. PLoS One 4, e5090. Schmitz, J. P., Van Riel, N. A., Nicolay, K., Hilbers, P. A. and Jeneson, J. A. (2009). Silencing of glycolysis in muscle : experimental observation and numerical analysis. Experimental Physiology 95, 380–397. Schuster, R. and Holzhütter, H. G. (1995). Use of mathematical models for predicting the metabolic effect of large-scale enzyme activity alterations. Application to enzyme deficiencies of red blood cells. European Journal of Biochemistry 229, 403–418. Siegel, T., Hekstra, D., Kemp, L., Figueiredo, L., Lowell, J., Fenyo, D., Wang, X., Dewell, S. and Cross, G. (2009). Four histone variants mark the boundaries of polycistronic transcription units in Trypanosoma brucei. Genes & Development 23, 1063–1076. 1341 Smallbone, K., Simeonidis, E., Broomhead, D. S. and Kell, D. B. (2007). Something from nothing : bridging the gap between constraint-based and kinetic modelling. FEBS Journal 274, 5576–5585. Snoep, J. L., Bruggeman, F., Olivier, B. G. and Westerhoff, H. V. (2006). Towards building the silicon cell : a modular approach. Biosystems 83, 207–216. Stern, M., Gupta, S., Salmon-Divon, M., Haham, T., Barda, O., Levi, S., Wachtel, C., Nilsen, T. and Michaeli, S. (2009). Multiple roles for polypyrimidine tract binding (PTB) proteins in trypanosome RNA metabolism. RNA 15, 648–665. Teusink, B., Walsh, M. C., Van Dam, K. and Westerhoff, H. V. (1998). The danger of metabolic pathways with turbo design. Trends in Biochemical Sciences 23, 162–169. Walrad, P., Paterou, A., Acosta-Serrano, A. and Matthews, K. R. (2009). Differential trypanosome surface coat regulation by a CCCH protein that coassociates with procyclin mRNA cis-elements. PLoS Pathogens 5, e1000317. Westerhoff, H. V., Kolodkin, A., Conradie, R., Wilkinson, S. J., Bruggeman, F. J., Krab, K., Van Schuppen, J. H., Hardin, H., Bakker, B. M., Moné, M. J., Rybakova, K. N., Eijken, M., Van Leeuwen, H. J. and Snoep, J. L. (2009). Sytems biology towards life in silico : mathematics of the control of living cells. Journal of Mathematical Biology 58, 7–34. Westerhoff, H. V., Koster, J. G., Van Workum, M. and Rudd, K. E. (1990). On the control of gene expression. In Control of Metabolic Processes (ed. Cornish-Bowden, A. and Cardenas, M.-L.), pp. 399–412. Plenum Press, New York. Xiao, Y., McCloskey, D. E. and Phillips, M. A. (2009). RNA interference-mediated silencing of ornithine decarboxylase and spermidine synthase genes in Trypanosoma brucei provides insight into regulation of polyamine biosynthesis. Eukaryotic Cell 8, 747–755.