Document Detail

Metabolic reconstruction identifies strain-specific regulation of virulence in Toxoplasma gondii.
Jump to Full Text
MedLine Citation:
PMID:  24247825     Owner:  NLM     Status:  MEDLINE    
Increasingly, metabolic potential is proving to be a critical determinant governing a pathogen's virulence as well as its capacity to expand its host range. To understand the potential contribution of metabolism to strain-specific infectivity differences, we present a constraint-based metabolic model of the opportunistic parasite, Toxoplasma gondii. Dominated by three clonal strains (Type I, II, and III demonstrating distinct virulence profiles), T. gondii exhibits a remarkably broad host range. Integrating functional genomic data, our model (which we term as iCS382) reveals that observed strain-specific differences in growth rates are driven by altered capacities for energy production. We further predict strain-specific differences in drug susceptibilities and validate one of these predictions in a drug-based assay, with a Type I strain demonstrating resistance to inhibitors that are effective against a Type II strain. We propose that these observed differences reflect an evolutionary strategy that allows the parasite to extend its host range, as well as result in a subsequent partitioning into discrete strains that display altered virulence profiles across different hosts, different organs, and even cell types.
Carl Song; Melissa A Chiasson; Nirvana Nursimulu; Stacy S Hung; James Wasmuth; Michael E Grigg; John Parkinson
Related Documents :
24307645 - The majority of a collection of us endocarditis enterococcus faecalis isolates from 197...
23114775 - 18β-glycyrrhetinic acid inhibits methicillin-resistant staphylococcus aureus survival a...
23468005 - A new enrichment method for isolation of bacillus thuringiensis from diverse sample types.
23395905 - Post-translational changes to prp alter transmissible spongiform encephalopathy strain ...
500565 - Chemotaxis of bdellovibrio bacteriovorus toward pure compounds.
8841725 - Video microscopy to quantitate the inhomogeneous equilibrium strain within articular ca...
Publication Detail:
Type:  Journal Article; Research Support, N.I.H., Intramural; Research Support, Non-U.S. Gov't     Date:  2013-11-19
Journal Detail:
Title:  Molecular systems biology     Volume:  9     ISSN:  1744-4292     ISO Abbreviation:  Mol. Syst. Biol.     Publication Date:  2013  
Date Detail:
Created Date:  2013-11-19     Completed Date:  2014-03-03     Revised Date:  2014-06-03    
Medline Journal Info:
Nlm Unique ID:  101235389     Medline TA:  Mol Syst Biol     Country:  England    
Other Details:
Languages:  eng     Pagination:  708     Citation Subset:  IM    
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Antiprotozoal Agents / pharmacology
Diphosphonates / pharmacology
Drug Resistance / drug effects
Fibroblasts / cytology,  parasitology*
Gene Expression Regulation*
Host Specificity
Host-Parasite Interactions
Metabolic Engineering
Metabolic Networks and Pathways*
Models, Genetic
Quinolines / pharmacology
Species Specificity
Toxoplasma / drug effects,  genetics,  metabolism*,  pathogenicity*
Grant Support
MOP #84556//Canadian Institutes of Health Research
Reg. No./Substance:
0/Antiprotozoal Agents; 0/Diphosphonates; 0/Quinolines

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Mol Syst Biol
Journal ID (iso-abbrev): Mol. Syst. Biol
ISSN: 1744-4292
Publisher: European Molecular Biology Organization
Article Information
Download PDF
Copyright © 2013, EMBO and Macmillan Publishers Limited
Received Day: 10 Month: 05 Year: 2013
Accepted Day: 10 Month: 10 Year: 2013
collection publication date: Year: 2013
Electronic publication date: Day: 19 Month: 11 Year: 2013
pmc-release publication date: Day: 19 Month: 11 Year: 2013
Volume: 9First Page: 708 Last Page: 708
PubMed Id: 24247825
ID: 4039375
Publisher Item Identifier: msb201362
DOI: 10.1038/msb.2013.62

Metabolic reconstruction identifies strain-specific regulation of virulence in Toxoplasma gondii
Carl Song12
Melissa A Chiasson3
Nirvana Nursimulu14
Stacy S Hung12
James Wasmuth16
Michael E Grigg3
John Parkinsona125
1Program in Molecular Structure and Function, The Hospital for Sick Children, Toronto, Ontario, Canada
2Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada
3Molecular Parasitology Section, Laboratory of Parasitic Diseases, NIAID, National Institutes of Health, Bethesda, MD, USA
4Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
5Department of Biochemistry, University of Toronto, Toronto, Ontario, Canada
6Current address: Department of Ecosystem and Public Health, Faculty of Veterinary Medicine, University of Calgary, 3280 Hospital Drive NW, Calgary, Alberta, Canada T2N 4Z6
aProgram in Molecular Structure and Function, The Hospital for Sick Children, 21.9709 Peter Gilgan Center for Research and Learning, 686 Bay Street, Toronto, Ontario, Canada M5G 0A4. Tel.:+1 416 813 5746; Fax:+1 416 813 5022;


Toxoplasma gondii is an opportunistic single-celled parasite with the capacity to infect any warm blooded animal. Thought to infect one in three people worldwide, infection by T. gondii results in toxoplasmosis, typically associated with flu-like symptoms in adults that resolve into a life-long chronic illness. More significantly, T. gondii can result in serious ocular disease in healthy adults and may be life threatening for pregnant women, transplant patients, and the immunocompromised (e.g., those living with HIV/AIDS) (Jeannel et al, 1988; Luft and Remington, 1992; Wong and Remington, 1994; Belanger et al, 1999). Despite its significance, few treatments are available and those that do exist do not promote sterile cure. The situation is further exacerbated with the emergence of new strains of parasites resistant to or tolerant of available prophylactics (Aspinall et al, 2002; Djaman et al, 2007). Indeed, it is increasingly evident that parasite strain ‘Type' is a predictor of virulent disease (Grigg et al, 2001b; McLeod et al, 2012). Among extant T. gondii lines, three strains (referred to as Types I, II, and III) dominate human infections in Europe and North America (Boothroyd and Grigg, 2002). Compared to Type II and Type III, Type I strains display relatively high growth rates and are acutely virulent in mice (LD100=1 parasite) (Howe and Sibley, 1995). Recent work has identified that murine virulence is highly dependent on the expression level of virulence factors, such as ROP18, GRA15, and SRS29C (Melo et al, 2011; Wasmuth et al, 2012), proteins that target host immune signalling pathways. At the same time, due to its importance in providing energy and the basic building blocks required for growth, metabolic potential is increasingly being viewed as a critical element governing a pathogen's virulence potential, as well as its ability to survive in infected hosts (McKinney et al, 2000; Olszewski et al, 2009; Willger et al, 2009; Ensminger et al, 2012). Through modulating metabolic capacity, parasites are able to tune growth in response to changes in host environment, offering a potential route to a broad host range. Genome comparisons reveal identical sets of genes encoding enzymes with the same predicted functional roles across the three strains. However, what is not known is how the differential expression of these genes across different Toxoplasma strains may influence their growth potential and hence virulence.

Genome-scale metabolic reconstruction has emerged as an effective strategy for systems-based investigations of an organism's metabolic potential, serving to crystallize current knowledge of an organism's metabolism as well as providing a framework for in silico investigation (Becker et al, 2007; Oberhardt et al, 2009; Thiele and Palsson, 2010). Metabolic reconstruction is an iterative process, beginning with the initial creation of a draft metabolic network based on the available enzyme annotation data. Subsequent rounds of simulation and refinement help resolve errors and fill gaps in otherwise incomplete networks (Green and Karp, 2004). With the increasing availability of high quality metabolic reconstructions, a variety of modeling procedures have been developed to analyse how these reconstructions are organized and operate. Arguably the most established method is flux balance analysis (FBA) (Kauffman et al, 2003; Lee et al, 2006), which solves for a steady-state distribution of reaction fluxes while satisfying a priori constraints (Orth et al, 2010). During an FBA simulation, the algorithm identifies sets of metabolic fluxes that optimize a specified function, for example, maximizing growth potential. This is achieved through the derivation of a ‘biomass equation', which details the proportions of all metabolites required for growth (e.g., DNA, RNA, protein, lipids, and cofactors). In the absence of constraints within the system, FBA can yield a range of optimal pathways to achieve maximal growth. To reduce the number of pathways, constraints can be placed on individual fluxes. Ideally these are obtained through systematic surveys of enzyme activities. However, in the absence of such data, mRNA expression data have been found to be an effective substitute (Colijn et al, 2009; Huthmacher et al, 2010). FBA has been successfully applied to a number of pathogens including Mycobacterium tuberculosis, Leishmania major, and Plasmodium falciparum (Raman et al, 2005; Chavali et al, 2008; Plata et al, 2010) to predict enzymes critical for growth and virulence.

Motivated by the need for a clearer understanding of the relationship between strain-specific metabolic capacity and a pathogen's ability to replicate and cause disease across a broad range of intermediate hosts, we present the first high quality metabolic reconstruction and constraint-based model of T. gondii. Our reconstruction, termed as iCS382, consists of 382 gene annotations, 282 enzymes, 384 unique metabolites, and 571 reactions, and provides a valuable reference resource detailing current state of knowledge of T. gondii metabolism. Integrating mRNA expression data, we use this model to identify strain-specific differences in metabolic potential that correlate with their observed growth rates. Drug inhibition assays are then applied to three enzymes to validate findings from our model.

A metabolic reconstruction crystallizes our current knowledge of T. gondii metabolism

Here, we present a systematic high quality reconstruction of the metabolic capabilities of T. gondii. We began by constructing an initial set of 258 enzymatic reactions (unique EC identifiers) that captures: (1) expert curation provided through the established T. gondii knowledgebase—ToxoDB (Gajria et al, 2008); (2) the curated Braunschweig Enzyme Database (BRENDA) (Barthelmes et al, 2007); (3) 43 previously published studies in T. gondii metabolism supporting 89 enzymatic reactions; and (4) 132 high confidence predictions of enzymes using the DETECT pipeline (Hung et al, 2010) (Figure 1A; Supplementary Table S1). Among these predictions were two genes that appear misannotated in ToxoDB (Supplementary Figure S1). TGME49_088450 is annotated as an aldehyde dehydrogenase (EC, but is predicted by DETECT to be 1-pyrroline-5-carboxylate dehydrogenase (EC, completing the pathway for autotrophic L-proline biosynthesis from L-glutamate. TGME49_109730 is annotated as a glutathione-disulfide reductase (EC, but is predicted by DETECT to be a thioredoxin-disulfide reductase (EC, which allows for the regeneration of thioredoxin and consequently the reduction in RNA nucleotides to their DNA counterparts. Note that reactions EC and EC are associated with other genes and are therefore also included in our reconstruction (Supplementary Table S2). During initial FBA simulations (see below), we found it is necessary to include an additional 24 so-called ‘gap-filling' enzymes to produce the full complement of biomass components required for parasite growth. Literature surveys were used to assign 282 enzymes (258 curated+24 gap filling) to five defined subcellular compartments: apicoplast; cytosol; endoplasmic reticulum; mitochondrion and mitochondrial intermembrane space. Some reactions were assigned to more than one compartment; thus 282 unique enzymatic activities (i.e., EC identifiers) were assigned to 400 reactions (Supplementary Table S1). Of these, 352 reactions are predicted through ToxoDB and DETECT to be encoded by 382 genes. Of the remaining 48 reactions, 19 represent 11 unique enzyme activities for which there is biochemical support only (see below), while the remaining 29 reactions are provided by the 24 ‘gap-filling' enzymes.

In addition to the 400 enzymatic reactions, we included 107 transport reactions comprising 85 organellar transport reactions that shuttle metabolites between different subcellular compartments, as well as 22 extracellular transport reactions that allow metabolites to enter and leave the system. For extracellular transport reactions, 21 represent known transporter proteins or auxotrophies supported by literature evidence, while one (tyrosine transport) represents an auxotrophy predicted by our model. In addition we include seven reactions that allow the diffusion of small molecules previously shown not to require a catalytic transporter (Boyle and Radke, 2009), as well as six sink reactions: metabolite exchanges with no support of prior knowledge but required for the production of defined biomass components. The 85 organellar transport reactions include both those catalyzed by active transporters as well as reactions that permit the passive diffusion of metabolites between different cell locations (Chavali et al, 2008; Thiele and Palsson, 2010).

In addition to enzymatic and transport reactions, our reconstruction includes three spontaneous molecular interconversions, three reactions required for recycling of currency metabolites through undefined cofactors, one reaction accounting for oxidative stress and nine reactions representing artificial conversions of biomass components to cell growth required by the model. Finally, our model includes an additional 33 so-called ‘dead-end reactions', which involve metabolites that are neither produced nor consumed by other metabolic reactions in the network and additionally lack evidence supporting the import or the export of the metabolites implying that the flux for the reaction must be zero. Enzymatic reactions, as defined by EC identifiers, were cross-referenced with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Kanehisa et al, 2006), to obtain reaction details including substrates and product metabolites, as well as reaction stoichiometry and direction. From an overall set of 571 reactions, the metabolic reconstruction is represented as a network incorporating 492 (384 unique) metabolites (Figure 1B); details of the reconstruction are provided in Table I.

To provide a guide for future hypothesis generation and model refinement, we provide two independent categories of confidence associated with the annotation of each reaction: Genomic confidence indicates that the reaction is supported by a gene association; Biochemical confidence indicates that the reaction is supported by experimental evidence (Figure 1B). Genomic confidence is assigned on the basis of either DETECT predictions (194 reactions; 132 unique EC IDs) or annotations provided by ToxoDB (348 reactions; 244 unique EC IDs). Similarly, Biochemical confidence is assigned on the basis of entries in the BRENDA database (63 reactions; 35 unique EC IDs) or the general literature (126 reactions; 89 unique EC IDs). During the reconstruction, we identified 19 reactions, involving 11 enzymes, for which there is biochemical evidence but for which a gene has yet to be assigned (Table II). These include 2-methylisocitrate dehydratase (EC, the only member of the 2-methylcitrate cycle (2-MCC), which generates acetyl CoA from the degradation of branched chain amino acids, for which no gene has currently been assigned (Seeber et al, 2008). Diamine N-acetyltransferase (EC is a member of the pathway required for the production of spermidine and putrescine from uptake of host spermine. Intriguingly we note two genes, TGME49_219760 and TGME49_289900, both of which are annotated by ToxoDB as N-acetyltransferase family proteins in ToxoDB with the EC identifier corresponding to peptide alpha-N-acetyltransferase ( We propose that since peptide alpha-N-acetyltransferase activity is not required by our model, these genes are potential candidates for diamine N-acetyltransferase (EC activity. Finally, we note that pyridoxine 5-phosphate synthase (EC consists of a heterodimer composed of two subunits, pdx1 and pdx2, which have been assigned to TGME49_237140 and TGME49_281490, respectively (Muller and Kappes, 2007). Both genes are annotated in ToxoDB as involved in pyridoxal phosphate (Vitamin B6) metabolism but without associated EC identifiers. These examples serve to illustrate the use of metabolic reconstructions to refine future annotation efforts.

Among the six classes of enzymes, transferases represent the most abundant followed by oxidoreductases (Figure 2A). The 492 metabolites involved in these reactions are predicted to be localized in five different subcellular compartments (Figure 2B). The large number of reactions in the cytosol is likely an overestimate as it is the default compartment for reactions without detailed localization data. Figure 2C shows the breakdown of reactions by general metabolic activity. Contrary to a previous metabolic reconstruction of an unrelated unicellular parasite, L. major, which highlighted the importance of lipid and amino-acid metabolism (Chavali et al, 2008), T. gondii, like P. falciparum (Huthmacher et al, 2010), possesses a large number of reactions involved in nucleic-acid metabolism. This reflects the impact of purine auxotrophy on apicomplexan parasites, resulting in the retention of multiple purine salvage pathways that can utilize a variety of purine precursors obtained from the host environment (Chaudhary et al, 2006). The final model is named ‘iCS382' consistent with metabolic reconstruction naming conventions (Reed et al, 2003); further details are provided in Supplementary Table S1. In the next section, we illustrate the use of this model through the application of constraint-based modelling to predict strain-specific metabolic behaviour.

FBA applied to iCS382 reveals strain-specific differences in metabolic potential

Previous studies have shown that T. gondii can be classified into three dominant clonal lineages, with each possessing different virulence profiles (Howe and Sibley, 1995). Given the relationship between parasite growth and virulence potential, we examined whether strain-specific differences in metabolic capabilities could account for changes in growth rate and hence virulence. Since genome sequencing reveals identical enzyme complements across strains, any metabolic differences will occur either at the level of gene/protein expression or through sequence variation. Here, we explore the former by integrating gene expression data into a series of FBA simulations to model the metabolic capabilities of the proliferative, tachyzoite form of the parasite across different strains. To our knowledge, this is the first time FBA has been applied to examine the potential impact of enzyme expression across different strains of pathogens.

To determine flux distributions for each reaction, the FBA framework requires a so-called objective function, defined here as maximizing the production of biomass components (i.e., growth; Table III). In the absence of a complete set of biochemical data detailing the kinetics of each T. gondii enzyme, flux constraints were assigned on the basis of gene expression profile data (see Materials and methods). These constraints define upper and lower bounds for each reaction. While gene-expression profiles provide only an approximation of enzyme activity, recent studies demonstrate that they significantly enhance model predictions through defining maximum flux constraints through reactions (Colijn et al, 2009; Plata et al, 2010). Our model relies on the assumption that flux associated with an enzyme is dependent on its mRNA expression, for example, strains expressing higher levels of enzyme transcripts will possess a higher flux for that reaction. Consequently, we integrated mRNA expression data previously generated for four strains of T. gondii (the Type I strains: RH and GT1, and the two Type II strains: Me49 and Prugniaud) as flux constraints in our FBA model (see Materials and methods). Flux constraints based on the expression data could be added for 209 reactions. These were further scaled to yield a predicted doubling time of 11.8 h for strain Me49 (see Materials and methods). On the basis of the scaling applied to Me49 and consistent with previous experimental studies (Radke et al, 2001; Saeij et al, 2005), our model correctly predicts that Type I strains of T. gondii have much higher growth rates (3.0 and 4.5 h for RH and GT, respectively) than Type II strains (11.8 and 14.0 h for Me49 and Prugniaud, respectively; Figure 3A). On the basis of the expression data, we find that among the pathways upregulated in RH relative to Me49 are pyrimidine biosynthesis, the TCA cycle and the pentose phosphate shunt (Supplementary Figure S2A). These differences suggest that Type I strains produce higher levels of ATP (driving an increased production of biomass) relative to Type II strains. In subsequent simulations, we explore these links in more detail (see below).

It should be appreciated that FBA rarely results in a single solution of optimal reaction fluxes. Nevertheless, by performing flux variability analysis (FVA) (Mahadevan and Schilling, 2003), it is possible to identify reactions that operate only at their maximum constraint for all possible solutions. Such reactions may be deemed ‘bottleneck' reactions, since the expectation is that by relaxing the constraint on these reactions, biomass production will similarly increase. We therefore applied FVA to identify potential bottleneck reactions in our four strains that appear to be constrained in different parts of the network (Supplementary Figure S2B). Consistent with our observations on growth rate, four of the eight bottlenecks for the production of biomass for all four strains involve the production of energy. Both ubiquinol-cytochrome-c reductase (EC, an essential enzyme in the production of ATP via oxidative phosphorylation, and phosphoglycerate kinase (EC and b) are found to be limiting across three strains. Strain RH is predicted to be only constrained by acetyl-CoA carboxylase (EC:, involved in one of the first steps of fatty-acid biosynthesis in the apicoplast. While bottleneck analysis identified the various chokepoints in the model that give rise to the optimal growth rates specific to each strain, it should be appreciated that these findings are reliant on the steady-state assumptions implied in the model, and may therefore not represent the true bottlenecks in a dynamic biological environment. Furthermore, removal of a bottleneck (e.g., by enzyme overexpression) may demonstrate little impact on organism growth in vitro/in vivo, as this may simply result in the identification of a new bottleneck reaction restricting growth to a similar (albeit lesser) degree. On the other hand, since these bottlenecks dictate the growth rate of the model, they allow reactions to be prioritized for further in silico and in vitro experiments that seek to refine associated flux constraints and consequently increase the accuracy of the model.

FBA of iCS382 identifies 39 new potential drug targets

An advantage of in silico modelling is the ability to rapidly investigate the impact of gene knockouts. On the basis of the Me49 reconstruction, just less than half of the enzymatic, extracellular, and organellar transport reactions (242 of 507 total) are predicted to be essential for biomass production and may therefore represent suitable targets for therapeutic intervention (Figure 3B and C; Supplementary Table S3). An additional 27 (5.3%) are predicted to have a ‘major' impact on growth rate (<80% optimal growth). Of the essential reactions, 171 are encoded by 129 metabolic enzymes. The model correctly predicts the essentiality of dihydrofolate reductase (EC in the folate biosynthetic pathway (targeted by the current anti-Toxoplasma therapeutic, pyrimethamine). However, dihydropteroate synthase (EC, the target of sulfadiazine, was not predicted to be essential, presumably due to the presence of thymidylate synthase activity (EC offering an alternative route for the production of dihydrofolate. Further model predictions of essentiality consistent with previous studies include two enzymes involved in fatty-acid synthesis (acetyl-CoA carboxylase, EC and enoyl reductase, EC (Zuther et al, 1999; Tipparaju et al, 2010). Due to a reliance on purine salvage, adenosine kinase (EC has been suggested as a putative target due to its high level of expression (Rodriguez and Szajnman, 2012). However, the presence of alternative routes for the production of AMP leads the model to predict a non-essential role for this enzyme, with its knockout predicted to have only a modest impact on growth rate (99.6% optimal growth). Finally, contrary to previous experimental studies (Fox and Bzik, 2003), our model does not predict the knockout of carbamoyl phosphate synthetase II (CPSII, EC to impact Toxoplasma growth. CPSII catalyses the production of carbamoyl phosphate, a substrate that is acted on by both ornithine carbamoyltransferase (EC and carbamate kinase (EC However, both reactions were set as reversible in the model, allowing them to provide alternate routes to the production of carbamoyl phosphate. But when both reactions were set as irreversible, our model subsequently predicted CPSII to be essential, confirming the experimental data and suggesting that both ornithine carbamoyltransferase and carbamate kinase reactions are irreversible under physiological conditions. In addition, our simulations also correctly predict tryptophan, cysteine and arginine auxotrophy (Pfefferkorn et al, 1986; Fox et al, 2004). Contrary to previous findings (Chaudhary and Roos, 2005), our model predicts Toxoplasma to be an auxotroph for tyrosine. Recent work has identified a putative phenylalanine hydroxylase (EC (Gaskell et al, 2009); inclusion of this reaction results in tyrosine import no longer being essential. These results highlight the iterative process of refining metabolic reconstructions through integration of additional experimental data to increase the accuracy of simulations.

To further prioritize this list, we identified 39 enzymes (unique EC identifiers), catalyzing 67 essential reactions, which lack homologs in humans (Supplementary Table S3). For example, we identified three major biosynthetic pathways essential to the parasite, with enzymes lacking homology in human: isoprenoid biosynthesis, pantothenate biosynthesis, and lysine biosynthesis. Previous studies have already implicated isoprenoid and pantothenate pathways for therapeutic intervention (Muller and Kappes, 2007; Moreno and Li, 2008), however, for lysine biosynthesis, we note mRNA expression levels of genes encoding five enzymes in the pathway are low, suggesting a reliance on host uptake.

Moving beyond single reaction knockouts, the availability of a robust metabolic model enables the rapid investigation of epistatic interactions within the network (i.e., pairs of reactions with redundant activities). The identification of such reaction pairs offers additional opportunities for therapeutic intervention through targeted drug combination strategies. Ignoring trivial reaction combinations in which one of the reactions is predicted to be essential, we systematically explored the impact of 26 106 pairs of reaction knockouts based on the Me49 reconstruction (Supplementary Table S4). We derived a genetic interaction score, calculated as the ratio of double knockout growth ratio to the product of single knockout growth ratios (see Materials and methods). Of all enzyme pairs simulated for double knockout, 1721 (6.6%) demonstrated any degree of epistasis, of which only 322 (1.2%) were predicted to be lethal (Figure 4). Our results predict that the TCA cycle and the pentose phosphate pathway are functionally redundant for a vital metabolic process, namely the regeneration of NAD(P)H and the subsequent production of ATP via oxidative phosphorylation. Other pathway-centric epistatic interactions of note are the glycolysis pathway in the apicoplast with both the TCA cycle and the pentose phosphate shunt, as well as reactions involved in the inter-conversion of nucleic-acid bases.

Comparisons of flux balance models identify both species- and strain-specific metabolic dependencies

In previous work, we showed that core metabolic functions encoded by different apicomplexans are provided both by sets of highly conserved enzymes, together with those that are lineage specific (Hung and Parkinson, 2011). The former group of enzymes is of particular interest as they represent putative pan-apicomplexan therapeutic targets. With the recent availability of a metabolic reconstruction for Plasmodium falciparum (Huthmacher et al, 2010), as well as a non-exhaustive list of anti-malarial enzyme drug targets determined through previous biochemical studies, we were therefore interested in examining how the different strategies adopted by the parasites to perform similar core metabolic activities might impact reaction essentiality. Figure 3D shows the overlap in predicted essential enzymes (see also Supplementary Figure S2C). Interestingly, while some (64 unique EC identifiers) of the 346 enzymes conserved between T. gondii and P. falciparum are either predicted (Hung and Parkinson, 2011) or confirmed to be essential in both species, we also note species-specific dependencies. For example, 33 of the conserved enzymes are predicted to be essential only in T. gondii, while 31 are either predicted or confirmed to be essential only in P. falciparum. Pathways predicted to be essential only in T. gondii include lysine biosynthesis and beta-alanine biosynthesis. Of the shared enzymes predicted to be essential in both species are components of glycolysis, pantothenate metabolism, heme biosynthesis, isoprenoid biosynthesis, pyrimidine metabolism, and fatty-acid metabolism—the latter four validated in P. falciparum. Due to the conserved nature of predictions, these pathways may be considered as most likely to yield targets for pan-apicomplexan therapeutics.

Moving beyond obvious lifestyle differences at the species level, we were interested in examining how changes in metabolism between different strains of the same species might impact metabolic dependencies. While Type I and Type II strains of T. gondii share identical enzyme complements, they possess markedly different growth rates. Here, we hypothesize that differences in enzyme expression result in significant changes in enzyme essentiality. Strains Me49 and RH are predicted to share 228 essential reactions (enzymatic, organellar, extracellular transport and excluding sink reactions), of which 162 are catalyzed by 121 metabolic enzymes (Supplementary Table S3). A further 14 reactions are predicted to be essential for strain Me49 but not RH. Ten of these involve metabolic enzymes or extracellular transporters, including six encoded by metabolic enzymes associated with glycolysis (Figure 3C; Table IV; Supplementary Table S3). These include lactate dehydrogenase (EC and diphosphate-fructose-6-phosphate 1-phosphotransferase (EC, both predicted to reduce growth to 42% of optimal for strain RH. In addition, while not essential to Me49, knockout of glucose-6-phosphate isomerase (EC, also involved in glycolysis, was predicted to reduce Me49 growth to 0.2% of optimal compared to 60% of optimal for strain RH. Conversely, our model does not predict any enzymes to be essential for strain RH. Model predictions also identified 18 enzymatic reactions, in which knockout was predicted to result in a significantly greater effect on growth rate (difference of >10% in optimal growth rate between the two strains) for strain Me49 compared to strain RH (Figure 3C; Table IV; Supplementary Table S3). These include enzymes in the pentose phosphate, glycolysis, and TCA cycle pathways. In contrast, 18 enzymatic and 2 extracellular transport reactions were predicted to result in a significantly greater effect on growth rate (difference of >10% in optimal growth rate between the two strains) for strain RH compared to strain Me49 (Figure 3C; Table IV; Supplementary Table S3). These include enzymes in the leucine degradation, oxidative phosphorylation, and TCA cycle pathways. For example, knockout of 3-methyl-2-oxobutanoate dehydrogenase (EC, hydroxymethylglutaryl-CoA lyase (EC, methylglutaconyl-CoA hydratase (EC, and methylcrotonoyl-CoA carboxylase (EC in the leucine degradation pathway were predicted to reduce growth of strain RH to 83% of optimal, with no predicted effect on Me49. Similarly, knockout of citrate synthase, EC, was predicted to result in 32% optimal growth for strain RH, compared to 94% optimal growth in strain Me49. Together, these findings indicate how the strains differ in their reliance on energy production pathways, with strain Me49 predicted to be more susceptible to the knockout of enzymes involved in glycolysis compared to RH, while both strains appear to be more reliant on different subsets of enzymes in the pentose phosphate pathway and TCA cycle.

To investigate the robustness of our conclusions and their dependence on the constraints assigned to individual reactions, we performed a systematic set of sensitivity analyses (Supplementary Figure S3). We first investigated the impact of increasing the maximum constraints for each individual metabolic reaction and examined their effect on the single reaction knockout predictions for strain Me49. Of the 400 reactions investigated, changing the constraints of only four reactions had any significant impact on the growth rates predicted for the single reaction knockouts (Supplementary Figure S3A). We also explored the dependence of the observed strain differences in growth rate, on the assigned constraints for eight reactions, mostly involved in energy production pathways (Supplementary Figure S3B). Results from these analyses suggest that our model predictions are robust to changes in the constraints assigned to individual reactions; observed strain differences appear to be driven by relative differences in global enzyme expression across the entire set of energy production pathways.

In the next section, we apply drug intervention screens to validate the observed strain-specific differences for two enzymes involved in glycolysis.

Drug sensitivity assays validate predicted strain-specific metabolic behavior

The strain-specific differences predicted by the model raise important questions regarding the impact of metabolic regulation on parasite virulence. In an attempt to validate these predictions, we selected three enzyme targets for drug inhibition assays: coproporphyrinogen oxidase (EC, fructose-bisphosphate aldolase (EC, and diphosphate-fructose-6-phosphate 1-phosphotransferase (EC Coproporphyrinogen oxidase is predicted to be essential for both strains Me49 and RH, while knockout of fructose-bisphosphate aldolase or diphosphate-fructose-6-phosphate 1-phosphotransferase is predicted to be lethal for strain Me49, but to only reduce growth to 42% for strain RH (Table IV). Five drug candidates from a list of 400 possibilities were prioritized based on previous work showing inhibition of growth in other microorganisms (Bai et al, 1982; Camadro et al, 1986; Bruchhaus et al, 1996; Rukseree et al, 2008). The five compounds used were 5,5-dithiobis(2-nitrobenzoic acid), alendronate, clodronate, 8-hydroxyquinoline, and 5-chloro-8-quinolinol. Two assays were conducted to measure growth of parasites under drug treatment: a 24-h doubling kinetic assay using microscopy and a 7-day growth assay using fluorescence. Assays at two different time points were necessary; it has been shown that many drugs, including atovaquone and clindamycin, do not act immediately to inhibit growth (Camps et al, 2002). Thus, while an inhibitory effect may not be observable at 24 h, at 7 days such an effect should be detectable. To assay for the predicted strain-specific differences in growth, drugs were tested in both RH GFP (Type I) and Me49 RFP (Type II) strains of the parasite.

5,5-dithiobis(2-nitrobenzoic acid) is predicted to inhibit coproporphyrinogen oxidase in the heme pathway (Bogard et al, 1989). Under treatment with 5,5-dithiobis(2-nitrobenzoic acid) for 24 h, no inhibitory effect was seen in RH GFP or Me49 RFP. However, after 7 days of treatment, both strains displayed a similar 50% inhibition of growth at a concentration of 900 μM of 5,5-dithiobis(2-nitrobenzoic acid) (Supplementary Figure S4). We next targeted the two enzymes in the glycolytic pathway. First, we applied the bisphosphonate, alendronate, to inhibit the activity of diphosphate-fructose-6-phosphate 1-phosphotransferase (Bruchhaus et al, 1996), which our model predicts should have a greater impact on growth of Me49 compared to RH. Consistent with these predictions we found that by day 7, while the growth of RH was not inhibited by alendronate up to a concentration of 100 μM, growth of Me49 was completely inhibited at the same concentration (Figure 5A). The 7-day assay confirmed that alendronate inhibits parasite growth in a strain-specific manner. Previous studies have shown that alendronate can also target the enzyme farnesyl pyrophosphate synthase in the mevalonate pathway (Montalvetti et al, 2001). To confirm that the observed strain-specific inhibition was the result of targeting diphosphate-fructose-6-phosphate 1-phosphotransferase in the fructose/mannose pathway, and not a result of off-target effects, we also examined the impact of clodronate (also a known inhibitor of diphosphate-fructose-6-phosphate 1-phosphotransferase; Bruchhaus et al, 1996) on parasite growth. Consistent with our findings for alendronate, we observed no effect on RH growth at 7 days up to 300 μM, whereas growth of Me49 was completely inhibited at this concentration (Figure 5A). Finally to determine the importance of the glycolytic pathway on parasite growth, we tested the impact of two additional drugs (8-hydroxyquinoline and 5-chloro-8-quinolinol) that target the enzyme fructose-bisphosphate aldolase immediately downstream of diphosphate-fructose-6-phosphate 1-phosphotransferase (Bai et al, 1982; Rukseree et al, 2008). After 7 days, treatment with 8-hydroxyquinoline at 25 μM and higher concentrations inhibited growth completely in both RH and Me49 strains, and treatment with 5-chloro-8-quinolinol at 12.5 μM and higher concentrations inhibited growth completely in both RH and Me49 strains (Figure 5B and C). Note parasite growth inhibition is not simply a consequence of the compounds impacting host cell viability (Supplementary Figure S5).

To confirm the relative expression levels of the targeted enzymes, quantitative PCR (qPCR) assays were performed. Four genes (TGME49_226960, 240890, 281390, and 281410) are annotated with diphosphate-fructose-6-phosphate 1-phosphotransferase activity, and are predicted to be targeted by alendronate and clodronate. qPCR revealed that only one of these genes was expressed at an appreciable level, confirming a strain-specific difference in its expression, with RH having significantly higher expression than Me49 (Figure 5D). Two genes (TGME49_236040 and 236050) are associated with fructose-bisphosphate aldolase activity that 8-hydroxyquinoline and 5-chloro-8-quinolinol are predicted to inhibit; both genes were expressed in a strain-specific manner (although expression of TGME_236040 was greater than TGME_236050) with RH significantly higher than Me49 (Figure 5D). Together, these results demonstrate the usefulness of our metabolic reconstruction, through validating the predicted essentiality of all three enzymes targeted here. In addition, we validated the predicted strain-specific sensitivity of Me49 to inhibition of diphosphate-fructose-6-phosphate 1-phosphotransferase (EC by both alendronate and clondronate. However, we were unable to validate a similar prediction for fructose-bisphosphate aldolase (EC As we note in Discussion, this finding may be related to the translocation of this enzyme during host cell egress (Starnes et al, 2006).


There is an increasing recognition that metabolic potential has the capacity to impact pathogen virulence and its ability to expand or restrict its host range (McKinney et al, 2000; Olszewski et al, 2009; Willger et al, 2009). For example, the accumulation of novel mutations in the lysine biosynthetic pathway of the opportunistic pathogen, Legionella pneumophila, was recently shown to confer a selective advantage in the mammalian host, at the expense of displayed reduced fitness in its natural host, amoeba (Ensminger et al, 2012). In this current study, we have demonstrated that changes in the expression of enzymes in different strains of T. gondii can likewise impact parasite virulence, which may impact its ability to expand its host range.

T. gondii can be classified into three distinct clonal lineages, designated as Types I, II, and III (McLeod et al, 2012). Each type exhibits distinct characteristics in terms of virulence during murine infection. For example, for strain RH; a single parasite is capable of lethal infection in mice, while for strain Me49, the LD50 for mice is 2000 parasites (Howe and Sibley, 1995). Furthermore, different strains exhibit partitioning in terms of their animal host, organs in which they persist and even the host cell type that they target (Wendte et al, 2011; Pan et al, 2012). Although it is known that different strains differentially target the host immune system (Melo et al, 2011), it is still not clear how the biochemical potential of different strains impacts their rate of growth and potential for clonal expansion and success.

To examine variation in metabolic behaviour across the different strains of T. gondii, we constructed a manually curated, high quality metabolic reconstruction (iCS382) that captures current biochemical knowledge of the parasite's metabolism. Applying a constraint-based framework, we systematically explored strain-specific models based on the differential expression of mRNA data. Although this method does not take into account any post-translational modifications or the specific activity of enzymes, mRNA expression has previously been shown to be useful in determining relative metabolic flux capacity (Colijn et al, 2009; Plata et al, 2010). Integrating these data we were able to demonstrate a faster growth phenotype for the highly virulent strain RH relative to Me49 consistent with experimental observations (Radke et al, 2001; Saeij et al, 2005). Our bottleneck and knockout analyses predicted the increased growth observed for strain RH, which appears due to greater production of ATP through upregulation of enzymes in the glycolytic, pentose phosphate, and TCA cycle pathways. Using the drugs alendronate and clodronate, known inhibitors of pyrophosphate-fructose-6-phosphate 1-phosphotransferase (EC, we show a greater inhibition of growth in strain Me49 relative to RH. Interestingly, drug inhibition studies focused on fructose-bisphosphate aldolase (EC did not confirm predictions from our model. This was all the more surprising given that: (1) EC was predicted to be essential for Me49; (2) drug treatment inhibited parasite growth; and (3) qPCR confirmed differential expression of the enzyme between the two strains. To explain this contradiction, we cannot exclude the possibility that Toxoplasma possess alternative metabolic routes, not captured in iCS382, that render the function of EC redundant. It is also possible that both 8-hydroxyquinoline and 5-chloro-8-quinolinol operate on an alternative target (the two chemical moieties share similar structures). However, a further explanation may concern the translocation of fructose-bisphosphate aldolase (Starnes et al, 2006). Specifically, during host cell egress, fructose-bisphosphate aldolase translocates to the parasite's pellicle as a route for optimizing ATP delivery to processes critical for extracellular survival (Pomel et al, 2008). Interference with the aldolase enzyme activity has been shown to impact parasite invasion, a necessary stage of the parasite's lifecycle. We therefore speculate that our inability to validate the strain-specific differences predicted by the glycolytic role of the enzyme is masked by the drug's interference post-translocation, impacting the parasite's ability to invade. Nonetheless, the strain-specific behaviour observed in the inhibition of EC demonstrates that while both strains of T. gondii possess the same complement of enzymes, changes in their regulation impact not only growth rate, but also the parasite's reliance on critical pathways involved in the production of energy.

These findings may have significant implications for the remarkably broad host range associated with the parasite. T. gondii is thought to be capable of infecting any nucleated cell from any warm blooded animal. During this stage of its life cycle, the parasite switches between a proliferative tachyzoite stage and a slow growing tissue cyst bradyzoite stage capable of transmissible infection to other hosts. In order to facilitate its development into the infective bradyzoite form, the parasite must strike a balance between: (a) growing too rapidly such that it kills its host before it has an opportunity to differentiate into the bradyzoite form; and (b) growing too slowly such that the host's immune system sterilizes the infection, again before it has an opportunity to differentiate into the bradyzoite form. Consequently, to be able to successfully colonize a wide range of hosts, the parasite must be able to rapidly modify its growth potential to account for the variety of immune responses and nutrient availability it encounters.

Previous studies suggest that regulatory changes in core components are a major evolutionary force driving the generation of viable selectable phenotypic variation that may operate at a faster rate than the innovation of new genes, which typically occurs through gene duplication and subsequent sub- or neo-functionalization (Carroll, 2000; Wray et al, 2003; He and Zhang, 2005; Gerhart and Kirschner, 2007; Han et al, 2009). On the basis of our current findings, we propose that changes in the regulation of enzymes required for the production of energy, as opposed to changes in the enzyme coding sequence itself as observed for L. pneumophila (Ensminger et al, 2012), offer a rapid evolutionary vehicle that allows T. gondii to modify its rate of growth in response to changes in host environment. Hence, observed differences in enzyme expression present the parasite with a fundamental route to adapt to changes in host environment through optimizing exploitation of nutrient availability. Recent studies have shown that sexual crosses of T. gondii in the definitive feline host result in the generation of a vast genetic repertoire of progeny (Grigg et al, 2001a). As subtle differences in the regulation of metabolic enzyme expression arise, the parasite is able to explore a vast landscape of metabolic potential, each with the possibility of increasing the efficient use of host nutrient availability. At the same time, in areas of restricted host heterogeneity, it might be expected that the parasite's metabolism becomes increasingly optimized for a limited host range, resulting in a subsequent partitioning into discrete strains that display markedly altered virulence profiles across different hosts, different organs, and even cell types (Wendte et al, 2011; Pan et al, 2012). As additional strains of T. gondii are sequenced, we will be able to explore these questions in more detail to examine both the diversity and the extent of metabolic potential that may underlie this drive toward optimizing host exploitation.

Beyond providing insights into the potential for changes in metabolic capacity to drive strain partitioning, iCS382 serves as a valuable platform to drive the design of anti-parasitic therapeutics. Studies of T. gondii metabolism, capitalizing on genomic data sets, are beginning to identify pathways that may represent additional routes for therapeutic investigations, including the biosynthesis of isoprenoids (Mol and Oudega, 1996), vitamins (Muller and Kappes, 2007), and polyamines (Cook et al, 2007). On the other hand, several biosynthetic processes previously thought to be critical for parasite survival have subsequently been found to be complemented by salvage pathways (i.e., scavenged from the host) (Massimine et al, 2005; Crawford et al, 2006). In a systematic set of single reaction deletion simulations, iCS382 predicts 242 of the 507 reactions to be essential for parasite growth including components of the beta-alanine, pantothenate, and glycolytic pathways. Integration of human homology data allows the prioritization of candidates suitable for therapeutic intervention. With predicted differences in essentiality and growth, it is clear that strain-specific differences in chemotherapeutic sensitivity must be considered during the development of any novel therapeutic. This is consistent with a previous study of Trypanosoma brucei rhodesiense infections in a murine model, which revealed significant strain-specific responses to drug treatment, with only six of eight clinical isolates susceptible to MDL73811, a potent inhibitor of S-adenosylmethionine decarboxylase (Bacchi et al, 1992).

In addition to strain-specific considerations, a valuable goal in these types of investigations is to identify candidates that represent new broad spectrum therapeutic targets. Comparison of our model predictions with those from a previous constraint-based modelling study of P. falciparum metabolism (Huthmacher et al, 2010) revealed a reasonably high level of overlap in essential reactions (57 of 116 common EC identifiers). Remarkably, despite their relatively close evolutionary relationship, comparisons of enzyme complements between the two reconstructions find only 193 enzymes (of the 609 enzymes included in both reconstructions) in common. Despite this lack of overlap, both models demonstrate the ability to support growth, lending support to the idea that each parasite lineage has evolved distinct metabolic strategies to exploit differences in nutrient availability provided by their unique host profiles (Hung and Parkinson, 2011). For example, while T. gondii has the biosynthetic capacity to produce phospholipids as well as several amino acids including serine, alanine, and glycine, P. falciparum relies on the breakdown of host heme and scavenging of host lipids. Moving beyond single drug therapies to reduce the rate of development in drug resistance, often a combination of enzyme inhibitors is administered to disrupt pathogen metabolism synergistically. This approach has widely been used for antimalarial and HIV treatments, and has also been shown to overcome drug resistance in Leishmania (Perez-Victoria et al, 2006). Our model identified 322 combinations that might serve as a useful basis for developing combination therapies in T. gondii.

In this study, we present a metabolic reconstruction of T. gondii and demonstrate how it might be usefully exploited not only to gain insights into parasite metabolism, but also to inform the choice of new targets for therapeutic intervention. In addition to validating predictions of strain-specific dependencies, this work unveils a host of new testable hypotheses concerning the reliance on specific pathways by different strains as well as a large number of new predicted drug targets. We acknowledge that the more streamlined nature of our T. gondii model may reflect the minimalist approach used in the reconstruction process and that there is much potential for future model refinement. These refinements include the incorporation of additional host-scavenging functions, which through the provision of key metabolites has the potential to render otherwise auxotrophic reactions non-essential. In addition, due to the incompleteness of the initial network reconstructed from current knowledge, a number of gap-filling reactions were added to the model to complete essential metabolic pathways. Furthermore, we have identified 33 so-called dead-end reactions that produce metabolites which are neither transported nor used as part of the biomass equation. As additional biochemical investigations are performed, we may expect these reactions to be confirmed and extended or alternative routes to be discovered. Future iterations of the model are also expected to incorporate additional compartments and pathways. For example, we do not currently include either the nucleus or the extracellular milieu as separate compartments. In terms of pathways, while we include a framework for glycan metabolism, the lack of detailed data concerning the role of glycans in the production of biomass for Toxoplasma (e.g., abundance and composition of GPI-anchors) has resulted in these pathways being largely incomplete. Similarly, despite evidence of 15-lipoxygenase activity in Toxoplasma, lack of additional evidence for supporting reactions precluded the incorporation of a complete arachidonic acid pathway, instead the production of arachidonic acid is captured by a single ‘black box' reaction. Finally, it should be appreciated that expression data provide only a crude approximation for reaction constraints used in the model. Detailed kinetic data, preferably obtained from in vivo experiments, would significantly contribute to the accuracy of the model. These caveats aside, by crystallizing current knowledge of T. gondii metabolism, we have shown how this resource may be leveraged as a valuable platform for integrating and organizing additional meta-datasets such as mRNA and protein expression data, SNPs and other comparative data sets. To facilitate access to iCS382, we make the model reconstruction freely available in the accepted reporting standard, systems biology markup language (SBML) level 2 format, available for download from our project website ( Our vision is that subsequent investigations will serve to both refine the model itself and shed light on the organization of T. gondii metabolism and enzyme dependencies throughout the entire life cycle of the parasite.

Materials and methods
Metabolic network reconstruction

An initial list of enzymes was compiled from the following resources: (1) ToxoDB gene annotations (Gajria et al, 2008); (2) BRENDA (Barthelmes et al, 2007) entries with T. gondii annotations; and (3) DETECT predictions (Hung et al, 2010) from T. gondii gene models (version 2008-07-23). Segregated reactions that do not link to other parts of the network and are responsible for the production of non-essential metabolites or whose substrates are not considered metabolites (e.g., polypeptides or other polymers) were removed. Subsequent manual analysis identified pathways with missing reactions by comparing the network to known pathway schemes in KEGG. A search in the literature was conducted to provide evidence for the presence of the pathway in T. gondii. Missing reactions were then added to the pathway if: (1) the pathway was reported or predicted with high confidence in the literature; (2) intermediate metabolites associated with the missing reactions were reported in the literature; or (3) a sizeable fraction of other reactions associated with the pathway were present. Lack of evidence resulted in the pathway being removed. For each enzyme, additional details such as substrate and product metabolites, direction of reaction, subcellular localization, and pathway annotation were retrieved from KEGG (Kanehisa et al, 2006), BRENDA (Barthelmes et al, 2007), and other literature sources. Intracellular transport reactions were added to move currency metabolites (H2O, ATP, etc.) between compartments. These include both previously known transporter proteins catalyzing intracellular transport reactions for specific metabolites and reactions involving passive diffusion. Negatively charged metabolites translocating to a relatively basic environment (i.e., from cytosol to mitochondrion) were accompanied by hydrolysis of ATP. The reconstruction is maintained as a spreadsheet (Supplementary Table S1) in a standard format (Thiele and Palsson, 2010). Gene to enzyme mappings, metabolite mappings and full citations are listed in Supplementary Tables S2, S5, and S6 respectively. Note here we represent enzymes performing the same reaction in different compartments as different reactions. The metabolic network was visualized as a bipartite graph using Cytoscape (Shannon et al, 2003), with nodes representing reactions and metabolites.

Biomass estimations

The mass of a single T. gondii cell together with its protein and lipid composition was obtained from a previous study of membrane fluidity (Gallois et al, 1988). The DNA fraction was calculated from the genome size, and the RNA fraction was estimated using a ratio of 1:7 DNA:RNA observed in other organisms (Beste et al, 2005). The amino-acid composition of the protein fraction and the NTP composition of RNA were estimated from the codon usage table available on ToxoDB (, while the dNTP composition of DNA was tabulated from the genomic sequence. For the remaining 3.76% of unaccounted total cell mass was unaccounted for by amino acids, lipids, DNA or RNA and was simply evenly distributed to cofactors and macromolecule subunits essential for parasite growth. Amounts of each biomass component were converted to mmol/gDW (Thiele and Palsson, 2010). We also included two additional biomass terms: growth-associated maintenance (GAM) and non-GAM (NGAM), which define the energy cost in terms of ATP hydrolysis for the organism to grow and sustain, respectively. GAM is estimated from the number of peptide and nucleotide bonds needed to polymerize the DNA, RNA, and protein fractions from their respective subunits. For NGAM, we define an invariable flux of 5 mmol/gDW·h (Thiele and Palsson, 2010).

Flux balance analysis

FBA was performed using the COBRA Toolbox (version 1.3.4) in MATLAB (Becker et al, 2007). Each reaction in the model is supplied upper and lower constraints for its flux. By default, irreversible reactions have their fluxes constrained to 0–1000, mmol/gDW·h, and reversible reactions −1000 to 1000, mmol/gDW·h. For reactions with single-gene associations, mRNA expression values derived from the Gene Expression Omnibus (GEO accession: GSE22315) were used to derive flux constraints. Each reaction receiving a flux constraint based on its associated gene expression relative to the highest gene expression value in the data set (which received an initial constraint of 1000, mmol/gDW·h). For irreversible reactions, the lower bound was set to 0 whereas for reversible reactions, the lower bound was set at the negative of the upper bound. For example, a reversible reaction encoded by a gene expressed at 20% of the level of the most highly expressed gene received constraints of −200 to 200 mmol/gDW·h. On the other hand, an irreversible reaction encoded by a gene expressed at 20% of the level of the most highly expressed gene received constraints of 0 to 200 mmol/gDW·h. The entire set of constraints was then scaled linearly (with the exception of NGAM which was set at a constant flux of 5 mmol/gDW·h consistent with previous studies; Thiele and Palsson, 2010) so that the predicted doubling time for Me49 matched the in vivo observation of 11.8 h (Blader et al, 2001). The calculation of scaling factors was performed only once, for the Me49 expression data. The same scaling of constraints based on the Me49 doubling time was then applied to expression data from the other three strains (RH, GT1, and Prugniaud).

Bottleneck analysis

For a given FBA solution, the optimization of the objective function is constrained by one or more bottleneck reactions in the network. Such bottlenecks were identified as those reactions with fluxes which equal to their maximum constraint (or minimum, in the case of reversible reactions) for all possible solutions leading to the optimized objective function. The solution space was computed using FVA, part of the COBRA Toolbox package. Reactions in the model were flagged as bottlenecks if the minimum flux of the reaction was equal to its maximum constraint.

Gene deletion predictions

Single knockouts were simulated for each reaction in the model by setting the constraints of the reaction to 0. Knockout effects were assessed by computing a growth ratio, which is the biomass production rate of the knockout divided by that of the wild type. Non-lethal single knockouts (having non-zero growth ratios) were permutated in pairs for double knockout simulations. Each double knockout simulation is in essence two single knockouts performed in tandem. Double knockouts were assessed by computing a genetic interaction score that is defined as follows:

[Formula ID: eq1]

where gr is the growth ratio of the respective knockouts.

Cell and parasites

Two strain-specific transgenic lines of Toxoplasma gondii were engineered. eGFP driven by the GRA1 promoter was stably integrated into the Type I RHΔhxgprt strain whereas DSRed 2.0 (RFP), also driven by the GRA1 promoter, was stably integrated into the Type II Me49 strain by Restriction Enzyme Mediated Insertion (REMI). Parasites with stable integration of the GFP or RFP cassette were selected by fluorescence. The point of insertion is unknown, but the parasites are identical in phenotype to wild type in terms of attachment, invasion, and replication. Parasites were routinely passaged in human foreskin fibroblast (HFF) cells.

Drug preparations

All drugs were purchased from Sigma (St Louis, MO). Both alendronate and clodronate were dissolved in RPMI media free of phenol red to form 5.4 mM stocks, and 8-quinoline and 5-chloro-8-quinolinol were dissolved in dimethyl sulfoxide (DMSO) to form 0.54 M stocks. Media prepared with mycophenolic acid (50 mg/ml) and xanthine (50 mg/ml) that inhibits the growth of Δhxgprt Toxoplasma strains was used as a no growth control for assays performed using the RHΔhxgprt GFP strain. All drug stocks were stored at 4°C.

Fluorescence assays

Black 96-well tissue culture-treated plates with opaque bottoms were purchased from BD Falcon (Bridgeport, NJ). Each well was seeded with cells in a volume of 200 μl. Before infection, plates were washed with RPMI to remove all traces of phenol red, to reduce interference and to increase the sensitivity of fluorescence detection. Freshly lysed parasites were plated in parasite culture medium without phenol red at a concentration of 104 parasites per well. Parasites were allowed to invade for 4 h before drug treatment was applied. Drug stocks were diluted three-fold through six wells, from 2.7 to 11 μM for alendronate, clodronate, and 5,5-dithiobis, and from 100 to 3.125 μM for 8-hydroxyquinoline and 5-chloro-8-quinolinol. Plates were kept in a humidified incubator at 37°C with 5% CO2 and read in a Perkin-Elmer Wallac fluorescent plate reader. The following excitation (485 nm for both GFP and RFP) and emission (535 nm for GFP; 595 nm for RFP) values were used. Fluorescence readings were measured daily, the assay proceeded until day 7 post infection.

Microscopy assays

Clear 96-well, Costar tissue culture-treated plates were purchased from Corning (Sigma-Aldrich, St Louis, MO). Cells were grown and infected as described above for the fluorescence assays. After 24 h of drug treatment, the average number of parasites present in ∼100 vacuoles across 5–7 fields of view (1, 2, 4, 8, 16, and, in rare cases 32 parasites per vacuole) was calculated using a fluorescence microscope at × 32 magnification. Photos were also taken to show growth or inhibition visually.

Quantitative PCR

Total RNA (2 μg) isolated from tachyzoites by the RNeasy mini kit (Qiagen) was reverse transcribed using random primers and SuperScript II (Invitrogen). Gene expression was measured by Taqman qPCR using an Applied Biosystems 7900HT Real-Time PCR System. The cycling program included 2 min at 50°C, 10 min incubation at 95°C followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Toxoplasma 18S rRNA and LDH1 were used as reference genes to normalize the quantity of transcripts (Livak and Schmittgen, 2001). Transcript levels were represented as 2−ΔCT to show absolute levels of transcript relative to every gene examined.

Host cell viability assays

Host cell viability assays under drug treatment were performed using the CellTiter-Glo Luminescent Cell Viability Assay (Promega). The kit's reagent was prepared according to the protocol, then added to 96-well plates in which HFF cells had been left under drug treatment for either 24 h or 7 days. Plates were read using a Perkin-Elmer Wallac fluorescent plate reader. Data analysis was performed using Prism.

Supplementary Material Supplementary Information

Supplementary Figures S1–5

Supplementary Table S1

Summary of iCS382

Click here for additional data file (msb201362-s2.xlsx)

Supplementary Table S2

List of enzymes with predicted gene mappings.

Click here for additional data file (msb201362-s3.xlsx)

Supplementary Table S3

Predicted impact of single enzyme knockouts on parasite growth.

Click here for additional data file (msb201362-s4.xlsx)

Supplementary Table S4

Predicted impact of double enzyme knockouts on parasite growth.

Click here for additional data file (msb201362-s5.xlsx)

Supplementary Table S5

List of metabolites that map to KEGG chemical IDs listed in Supplementary Table S1.

Click here for additional data file (msb201362-s6.xlsx)

Supplementary Table S6

List of citations that map to PubMed IDs listed in Supplementary Table S1

Click here for additional data file (msb201362-s7.xlsx)

Model file in SBML format
Click here for additional data file (msb201362-s8.xml)

MatLab files for Model
Click here for additional data file (

Review Process File


The authors declare that they have no conflict of interest.

This study was funded by the Canadian Institutes for Health Research (CIHR—MOP #84556 to JP and MEG) and the Natural Sciences and Engineering Research Council (NSERC #188266-04 to JP). JP also acknowledges support from the Ontario Ministry of Research and Innovation. MEG and MAC also acknowledge support from the Intramural Research Program of the National Institutes of Health (NIH) and National Institute of Allergy and Infectious Diseases (NIAID). JP is a member of the Center for the analysis of genome evolution and function (CAGEF). MEG is a Scholar of the Canadian Institute for Advanced Research (CIFAR) Program for Integrated Microbial Biodiversity. Computing resources were provided by the SciNet HPC Consortium.

Author contributions: JP and MEG conceived and designed the study. CS, JW, NN, and SSH performed reconstruction. CS and NN performed simulations. MAC performed drug assays. CS, JP, MAC, MEG, NN, and SSH analysed results. CS, JP, JW, MEG, MAC, NN, and SSH drafted manuscript.

Aspinall TV,, Joynson DH,, Guy E,, Hyde JE,, Sims PF, (Year: 2002) The molecular basis of sulfonamide resistance in Toxoplasma gondii and implications for the clinical management of toxoplasmosis. J Infect Dis185: 1637–164312023770
Bacchi CJ,, Nathan HC,, Yarlett N,, Goldberg B,, McCann PP,, Bitonti AJ,, Sjoerdsma A, (Year: 1992) Cure of murine Trypanosoma brucei rhodesiense infections with an S-adenosylmethionine decarboxylase inhibitor. Antimicrob Agents Chemother36: 2736–27401482141
Bai NJ,, Pai MR,, Murthy PS,, Venkitasubramanian TA, (Year: 1982) Fructose-bisphosphate aldolases from mycobacteria. Methods Enzymol90(Pt E): 241–2506818421
Barthelmes J,, Ebeling C,, Chang A,, Schomburg I,, Schomburg D, (Year: 2007) BRENDA, AMENDA and FRENDA: the enzyme information system in 2007. Nucleic Acids Res35: D511–D51417202167
Becker SA,, Feist AM,, Mo ML,, Hannum G,, Palsson BO,, Herrgard MJ, (Year: 2007) Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox. Nat Protoc2: 727–73817406635
Belanger F,, Derouin F,, Grangeot-Keros L,, Meyer L, (Year: 1999) Incidence and risk factors of toxoplasmosis in a cohort of human immunodeficiency virus-infected patients: 1988-1995. HEMOCO and SEROCO Study Groups. Clin Infect Dis28: 575–58110194081
Beste DJ,, Peters J,, Hooper T,, Avignone-Rossa C,, Bushell ME,, McFadden J, (Year: 2005) Compiling a molecular inventory for Mycobacterium bovis BCG at two growth rates: evidence for growth rate-mediated regulation of ribosome biosynthesis and lipid metabolism. J Bacteriol187: 1677–168415716438
Blader IJ,, Manger ID,, Boothroyd JC, (Year: 2001) Microarray analysis reveals previously unknown changes in Toxoplasma gondii-infected human cells. J Biol Chem276: 24223–2423111294868
Bogard M,, Camadro JM,, Nordmann Y,, Labbe P, (Year: 1989) Purification and properties of mouse liver coproporphyrinogen oxidase. Eur J Biochem181: 417–4212540974
Boothroyd JC,, Grigg ME, (Year: 2002) Population biology of Toxoplasma gondii and its relevance to human infection: do different strains cause different disease?Curr Opin Microbiol5: 438–44212160866
Boyle JP,, Radke JR, (Year: 2009) A history of studies that examine the interactions of Toxoplasma with its host cell: Emphasis on in vitro models. Int J Parasitol39: 903–91419630139
Bruchhaus I,, Jacobs T,, Denart M,, Tannich E, (Year: 1996) Pyrophosphate-dependent phosphofructokinase of Entamoeba histolytica: molecular cloning, recombinant expression and inhibition by pyrophosphate analogues. Biochem J316(Pt 1): 57–638645233
Camadro JM,, Chambon H,, Jolles J,, Labbe P, (Year: 1986) Purification and properties of coproporphyrinogen oxidase from the yeast Saccharomyces cerevisiae. Eur J Biochem156: 579–5873516695
Camps M,, Arrizabalaga G,, Boothroyd J, (Year: 2002) An rRNA mutation identifies the apicoplast as the target for clindamycin in Toxoplasma gondii. Mol Microbiol43: 1309–131811918815
Carroll S, (Year: 2000) Endless forms: the evolution of gene regulation and morphological diversity. Cell101: 57710892643
Chaudhary K,, Darling JA,, Fohl LM,, Sullivan WJ Jr.,, Donald RG,, Pfefferkorn ER,, Ullman B,, Roos DS. (Year: 2004) Purine salvage pathways in the apicomplexan parasite Toxoplasma gondii. J Biol Chem279: 31221–3122715140885
Chaudhary K,, Roos DS, (Year: 2005) Protozoan genomics for drug discovery. Nat Biotechnol23: 1089–109116151400
Chaudhary K,, Ting LM,, Kim K,, Roos DS, (Year: 2006) Toxoplasma gondii purine nucleoside phosphorylase biochemical characterization, inhibitor profiles, and comparison with the Plasmodium falciparum ortholog. J Biol Chem281: 25652–2565816829527
Chavali AK,, Whittemore JD,, Eddy JA,, Williams KT,, Papin JA, (Year: 2008) Systems analysis of metabolism in the pathogenic trypanosomatid Leishmania major. Mol Syst Biol4: 17718364711
Colijn C,, Brandes A,, Zucker J,, Lun DS,, Weiner B,, Farhat MR,, Cheng TY,, Moody DB,, Murray M,, Galagan JE, (Year: 2009) Interpreting expression data with metabolic flux models: predicting Mycobacterium tuberculosis mycolic acid production. PLoS Comput Biol5: e100048919714220
Cook T,, Roos D,, Morada M,, Zhu G,, Keithly JS,, Feagin JE,, Wu G,, Yarlett N. (Year: 2007) Divergent polyamine metabolism in the Apicomplexa. Microbiology153: 1123–113017379721
Crawford MJ,, Thomsen-Zieger N,, Ray M,, Schachtner J,, Roos DS,, Seeber F, (Year: 2006) Toxoplasma gondii scavenges host-derived lipoic acid despite its de novo synthesis in the apicoplast. EMBO J25: 3214–322216778769
Dautu G,, Ueno A,, Munyaka B,, Carmen G,, Makino S,, Kobayashi Y,, Igarashi M, (Year: 2008) Molecular and biochemical characterization of Toxoplasma gondii beta-hydroxyacyl-acyl carrier protein dehydratase (FABZ). Parasitol Res102: 1301–130918274777
Djaman JA,, Mazabraud A,, Basco L, (Year: 2007) Sulfadoxine-pyrimethamine susceptibilities and analysis of the dihydrofolate reductase and dihydropteroate synthase of Plasmodium falciparum isolates from Cote d'Ivoire. Ann Trop Med Parasitol101: 103–11217316496
Ensminger AW,, Yassin Y,, Miron A,, Isberg RR, (Year: 2012) Experimental evolution of Legionella pneumophila in mouse macrophages leads to strains with altered determinants of environmental survival. PLoS Pathog8: e100273122693450
Fox BA,, Bzik DJ, (Year: 2003) Organisation and sequence determination of glutamine-dependent carbamoyl phosphate synthetase II in Toxoplasma gondii. Int J Parasitol33: 89–9612547350
Fox BA,, Gigley JP,, Bzik DJ, (Year: 2004) Toxoplasma gondii lacks the enzymes required for de novo arginine biosynthesis and arginine starvation triggers cyst formation. Int J Parasitol34: 323–33115003493
Gajria B,, Bahl A,, Brestelli J,, Dommer J,, Fischer S,, Gao X,, Heiges M,, Iodice J,, Kissinger JC,, Mackey AJ,, Pinney DF,, Roos DS,, Stoeckert CJ Jr.,, Wang H,, Brunk BP, (Year: 2008) ToxoDB: an integrated Toxoplasma gondii database resource. Nucleic Acids Res36: D553–D55618003657
Gallois Y,, Foussard F,, Girault A,, Hodbert J,, Tricaud A,, Mauras G,, Motta C, (Year: 1988) Membrane fluidity of Toxoplasma gondii: a fluorescence polarization study. Biol Cell62: 11–153365515
Gaskell EA,, Smith JE,, Pinney JW,, Westhead DR,, McConkey GA, (Year: 2009) A unique dual activity amino acid hydroxylase in Toxoplasma gondii. PLoS One4: e480119277211
Gerhart J,, Kirschner M, (Year: 2007) The theory of facilitated variation. Proc Natl Acad Sci USA104(Suppl 1): 8582–858917494755
Green ML,, Karp PD, (Year: 2004) A Bayesian method for identifying missing enzymes in predicted metabolic pathway databases. BMC Bioinformatics5: 7615189570
Grigg ME,, Bonnefoy S,, Hehl AB,, Suzuki Y,, Boothroyd JC, (Year: 2001a) Success and virulence in Toxoplasma as the result of sexual recombination between two distinct ancestries. Science294: 161–16511588262
Grigg ME,, Ganatra J,, Boothroyd JC,, Margolis TP, (Year: 2001b) Unusual abundance of atypical strains associated with human ocular toxoplasmosis. J Infect Dis184: 633–63911474426
Han MV,, Demuth JP,, McGrath CL,, Casola C,, Hahn MW, (Year: 2009) Adaptive evolution of young gene duplicates in mammals. Genome Res19: 859–86719411603
He X,, Zhang J, (Year: 2005) Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution. Genetics169: 1157–116415654095
Howe DK,, Sibley LD, (Year: 1995) Toxoplasma gondii comprises three clonal lineages: correlation of parasite genotype with human disease. J Infect Dis172: 1561–15667594717
Hung SS,, Parkinson J, (Year: 2011) Post-genomics resources and tools for studying apicomplexan metabolism. Trends Parasitol27: 131–14021145790
Hung SS,, Wasmuth J,, Sanford C,, Parkinson J, (Year: 2010) DETECT--a density estimation tool for enzyme classification and its application to Plasmodium falciparum. Bioinformatics26: 1690–169820513663
Huthmacher C,, Hoppe A,, Bulik S,, Holzhutter HG, (Year: 2010) Antimalarial drug targets in Plasmodium falciparum predicted by stage-specific metabolic network analysis. BMC Syst Biol4: 12020807400
Jeannel D,, Niel G,, Costagliola D,, Danis M,, Traore BM,, Gentilini M, (Year: 1988) Epidemiology of toxoplasmosis among pregnant women in the Paris area. Int J Epidemiol17: 595–6023264821
Kanehisa M,, Goto S,, Hattori M,, Aoki-Kinoshita KF,, Itoh M,, Kawashima S,, Katayama T,, Araki M,, Hirakawa M, (Year: 2006) From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res34: D354–D35716381885
Kauffman KJ,, Prakash P,, Edwards JS, (Year: 2003) Advances in flux balance analysis. Curr Opin Biotechnol14: 491–49614580578
Lee JM,, Gianchandani EP,, Papin JA, (Year: 2006) Flux balance analysis in the era of metabolomics. Brief Bioinform7: 140–15016772264
Livak KJ,, Schmittgen TD, (Year: 2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods25: 402–40811846609
Luft BJ,, Remington JS, (Year: 1992) Toxoplasmic encephalitis in AIDS. Clin Infect Dis15: 211–2221520757
Mahadevan R,, Schilling CH, (Year: 2003) The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab Eng5: 264–27614642354
Massimine KM,, Doan LT,, Atreya CA,, Stedman TT,, Anderson KS,, Joiner KA,, Coppens I, (Year: 2005) Toxoplasma gondii is capable of exogenous folate transport. A likely expansion of the BT1 family of transmembrane proteins. Mol Biochem Parasitol144: 44–5416159678
McKinney JD,, Honer zu Bentrup K,, Munoz-Elias EJ,, Miczak A,, Chen B,, Chan WT,, Swenson D,, Sacchettini JC,, Jacobs WR Jr.. , Russell DG, (Year: 2000) Persistence of Mycobacterium tuberculosis in macrophages and mice requires the glyoxylate shunt enzyme isocitrate lyase. Nature406: 735–73810963599
McLeod R,, Boyer KM,, Lee D,, Mui E,, Wroblewski K,, Karrison T,, Noble AG,, Withers S,, Swisher CN,, Heydemann PT,, Sautter M,, Babiarz J,, Rabiah P,, Meier P,, Grigg ME, (Year: 2012) Prematurity and severity are associated with Toxoplasma gondii alleles (NCCCTS, 1981-2009). Clin Infect Dis54: 1595–160522499837
Melo MB,, Jensen KD,, Saeij JP, (Year: 2011) Toxoplasma gondii effectors are master regulators of the inflammatory response. Trends Parasitol27: 487–49521893432
Mol O,, Oudega B, (Year: 1996) Molecular and structural aspects of fimbriae biosynthesis and assembly in Escherichia coli. FEMS Microbiol Rev19: 25–528916554
Montalvetti A,, Bailey BN,, Martin MB,, Severin GW,, Oldfield E,, Docampo R, (Year: 2001) Bisphosphonates are potent inhibitors of Trypanosoma cruzi farnesyl pyrophosphate synthase. J Biol Chem276: 33930–3393711435429
Moreno SN,, Li ZH, (Year: 2008) Anti-infectives targeting the isoprenoid pathway of Toxoplasma gondii. Expert Opin Ther Targets12: 253–26318269336
Muller S,, Kappes B, (Year: 2007) Vitamin and cofactor biosynthesis pathways in Plasmodium and other apicomplexan parasites. Trends Parasitol23: 112–12117276140
Oberhardt MA,, Palsson BO,, Papin JA, (Year: 2009) Applications of genome-scale metabolic reconstructions. Mol Syst Biol5: 32019888215
Olszewski KL,, Morrisey JM,, Wilinski D,, Burns JM,, Vaidya AB,, Rabinowitz JD,, Llinas M, (Year: 2009) Host-parasite interactions revealed by Plasmodium falciparum metabolomics. Cell Host Microbe5: 191–19919218089
Orth JD,, Thiele I,, Palsson BO, (Year: 2010) What is flux balance analysis?Nat Biotechnol28: 245–24820212490
Pan S,, Thompson RC,, Grigg ME,, Sundar N,, Smith A,, Lymbery AJ, (Year: 2012) Western Australian marsupials are multiply infected with genetically diverse strains of Toxoplasma gondii. PLoS One7: e4514723028812
Perez-Victoria JM,, Cortes-Selva F,, Parodi-Talice A,, Bavchvarov BI,, Perez-Victoria FJ,, Munoz-Martinez F,, Maitrejean M,, Costi MP,, Barron D,, Di Pietro A. , Castanys S,, Gamarro F, (Year: 2006) Combination of suboptimal doses of inhibitors targeting different domains of LtrMDR1 efficiently overcomes resistance of Leishmania spp. to Miltefosine by inhibiting drug efflux. Antimicrob Agents Chemother50: 3102–311016940108
Pfefferkorn ER,, Eckel M,, Rebhun S, (Year: 1986) Interferon-gamma suppresses the growth of Toxoplasma gondii in human fibroblasts through starvation for tryptophan. Mol Biochem Parasitol20: 215–2243093859
Plata G,, Hsiao TL,, Olszewski KL,, Llinas M,, Vitkup D, (Year: 2010) Reconstruction and flux-balance analysis of the Plasmodium falciparum metabolic network. Mol Syst Biol6: 40820823846
Pomel S,, Luk FC,, Beckers CJ, (Year: 2008) Host cell egress and invasion induce marked relocations of glycolytic enzymes in Toxoplasma gondii tachyzoites. PLoS Pathog4: e100018818949028
Radke JR,, Striepen B,, Guerini MN,, Jerome ME,, Roos DS,, White MW, (Year: 2001) Defining the cell cycle for the tachyzoite stage of Toxoplasma gondii. Mol Biochem Parasitol115: 165–17511420103
Raman K,, Rajagopalan P,, Chandra N, (Year: 2005) Flux balance analysis of mycolic Acid pathway: targets for anti-tubercular drugs. PLoS Comput Biol1: e4616261191
Reed JL,, Vo TD,, Schilling CH,, Palsson BO, (Year: 2003) An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). Genome Biol4: R5412952533
Rodriguez JB,, Szajnman SH, (Year: 2012) New antibacterials for the treatment of toxoplasmosis; a patent review. Expert Opin Ther Pat22: 311–33322404108
Rukseree K,, Thammarongtham C,, Palittapongarnpim P, (Year: 2008) One-step purification and characterization of a fully active histidine-tagged Class II fructose-1,6-bisphosphate aldolase from Mycobacterium tuberculosis. Enzyme Microbial Technol43: 500–506
Saeij JP,, Boyle JP,, Boothroyd JC, (Year: 2005) Differences among the three major strains of Toxoplasma gondii and their specific interactions with the infected host. Trends Parasitol21: 476–48116098810
Seeber F,, Limenitakis J,, Soldati-Favre D, (Year: 2008) Apicomplexan mitochondrial metabolism: a story of gains, losses and retentions. Trends Parasitol24: 468–47818775675
Shannon P,, Markiel A,, Ozier O,, Baliga NS,, Wang JT,, Ramage D,, Amin N,, Schwikowski B,, Ideker T, (Year: 2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res13: 2498–250414597658
Smith TK,, Kimmel J,, Azzouz N,, Shams-Eldin H,, Schwarz RT, (Year: 2007) The role of inositol acylation and inositol deacylation in the Toxoplasma gondii glycosylphosphatidylinositol biosynthetic pathway. J Biol Chem282: 32032–3204217804418
Starnes GL,, Jewett TJ,, Carruthers VB,, Sibley LD, (Year: 2006) Two separate, conserved acidic amino acid domains within the Toxoplasma gondii MIC2 cytoplasmic tail are required for parasite survival. J Biol Chem281: 30745–3075416923803
Thiele I,, Palsson BO, (Year: 2010) A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc5: 93–12120057383
Tipparaju SK,, Muench SP,, Mui EJ,, Ruzheinikov SN,, Lu JZ,, Hutson SL,, Kirisits MJ,, Prigge ST,, Roberts CW,, Henriquez FL,, Kozikowski AP,, Rice DW,, McLeod RL, (Year: 2010) Identification and development of novel inhibitors of Toxoplasma gondii enoyl reductase. J Med Chem53: 6287–630020698542
Wasmuth JD,, Pszenny V,, Haile S,, Jansen EM,, Gast AT,, Sher A,, Boyle JP,, Boulanger MJ,, Parkinson J,, Grigg ME, (Year: 2012) Integrated bioinformatic and targeted deletion analyses of the SRS gene superfamily identify SRS29C as a negative regulator of Toxoplasma virulence. MBio3:
Wendte JM,, Gibson AK,, Grigg ME, (Year: 2011) Population genetics of Toxoplasma gondii: new perspectives from parasite genotypes in wildlife. Vet Parasitol182: 96–11121824730
Willger SD,, Grahl N,, Cramer RA Jr., (Year: 2009) Aspergillus fumigatus metabolism: clues to mechanisms of in vivo fungal growth and virulence. Med Mycol47(Suppl 1): S72–S7919253141
Wong SY,, Remington JS, (Year: 1994) Toxoplasmosis in pregnancy. Clin Infect Dis18: 853–861 quiz 8628086543
Wray GA,, Hahn MW,, Abouheif E,, Balhoff JP,, Pizer M,, Rockman MV,, Romano LA, (Year: 2003) The evolution of transcriptional regulation in Eukaryotes. Mol Biol Evol20: 1377–141912777501
Zuther E,, Johnson JJ,, Haselkorn R,, McLeod R,, Gornicki P, (Year: 1999) Growth of Toxoplasma gondii is inhibited by aryloxyphenoxypropionate herbicides targeting acetyl-CoA carboxylase. Proc Natl Acad Sci USA96: 13387–1339210557330

Article Categories:
  • Article

Keywords: flux balance analysis, metabolic reconstruction, strain differences, Toxoplasma gondii.

Previous Document:  Renal function in normal and disordered pregnancy.
Next Document:  LIM kinase 1 is required for insulin?dependent cell growth of osteosarcoma cell lines.