BACKGROUND
salmonellae are gram-negative bacterial pathogens with a wide host range.
the ability of salmonella enterica to cause disease is closely correlated with its ability to survive inside its host cells, which include epithelial cells and macrophages  <cit> .
infection of mice with salmonella enterica serovar typhimurium  serves as model for the systemic disease typhoid fever.
salmonella strains that are deficient for intra-macrophage survival cannot elicit disease in the mouse model.
many genomic and proteomic datasets exist for salmonella as a pathogen  <cit> .
most experiments have focused on identifying 'virulence' factors that might drive the disease process, such as surface structures or secreted products.
differential expression of metabolic gene-products has largely been left unexplored in terms of their role in pathogenesis.
however, it has now gained more attention that the metabolism of intracellular bacterial pathogens such as salmonella may hold important clues to understanding their intracellular life-style as well as to understanding host defense mechanisms  <cit> .
genome-scale metabolic in silico models that can integrate and exploit genomic and proteomic data may provide a first step towards a quantitative analysis of the metabolism of intracellular bacteria.
constraint-based reconstruction and analysis  of cellular metabolism has been employed successfully to develop organism-specific genome-scale in silico models across all the three major domains of the tree of life  <cit> .
the cobra approach seeks to differentiate between those network states that are achievable by the organism from those that are not, based on a detailed reconstruction of metabolism and incorporation of physiological parameters that are consistent with known experimental information  <cit> .
the solutions this analysis provides are based on the successive imposition of governing physicochemical constraints on genome-scale reconstructions.
cobra integrates and represents our current knowledge of network components and interactions and has been applied to several prokaryotic and eukaryotic organisms, including escherichia coli  <cit> , helicobacter pylori  <cit> , haemophilus influenzae  <cit> , saccharomyces cerevisiae  <cit>  and geobacter sulfurreducens  <cit> .
by systematizing and providing context for transcriptomic, proteomic, and metabolomic data, these models allow for the imposition of constraints, which together define the possible phenotypic behavior allowed by these in silico organisms.
the use of constrained based modeling has significant potential to identify new antibacterial targets  <cit> .
modeling allows simulation of single gene knockouts and, more importantly, the effect of combinatorial deletions, which might allow the rational development of combination drug therapy.
the use of such a modeling approach has been demonstrated for h. influenza  <cit>  and for mycobacterium tuberculosis  <cit> .
as the number of pathogens with such models grows  <cit> , one might expect that more basic principles of bacterial metabolism needed for pathogenesis can be derived.
here we describe the reconstruction of a detailed and validated genome-scale metabolic network for s. typhimurium, that reconciles established genomic and high-throughput phenotyping data and that is augmented by experimental data to facilitate iterative model building efforts.
furthermore, we demonstrate how the model can be used to reconcile existing virulence data and to resolve gaps in our understanding of pathogen metabolism and host-cell nutrient environment.
we also use gene expression data obtained for this study in combination with established microarray and proteomic data in the context of the model to understand pathogen metabolism in macrophage cell lines.
RESULTS
characteristics of the reconstructed metabolic network
the  <dig> -kilobase  chromosome of salmonella enterica serovar typhimurium lt <dig>  accounts for  <dig>  coding sequences  including  <dig> pseudogenes and a  <dig> kb plasmid encoding  <dig> orfs  <cit> .
the in silico metabolic reconstruction of the bacterium referred to as irr <dig> contains only annotated chromosomal genes.
the characteristics of the network are detailed in figure 1a and 1b .
it currently accounts for  <dig>  genes that encode proteins catalyzing  <dig>  unique metabolic and transport reactions in the bacterium.
the reconstruction includes gene to protein to reaction association relationships, which distinguish between isozymes, multimeric proteins and protein complexes  <cit> .
the reconstruction contains all the biosynthesis pathways to make biomass components that are required for growth.
a biomass production reaction was formulated, which is based on the measured biomass composition of the bacterium and which details the required metabolites and their relative amounts.
the reconstruction accounts for central metabolism, carbohydrate utilization, amino acid biosynthesis and degradation, nucleotide biosynthesis and salvage, and cofactor biosynthesis .
also included were anaerobic and aerobic electron transport chain reactions, transport reactions for various minerals and small molecules, and biosynthesis reactions for cell wall constituents and precursors, such as phospholipids and peptidoglycan.
lipid biosynthesis in bacteria involves  fatty acid biosynthesis,  attachment of the complete fatty acids to sn-glycerol-3-phosphate, and  addition and modification of head groups to yield the major cellular phospholipids.
fatty acids are important constituents of lipid a and the cell wall of bacteria.
in contrast to published models of other microbes that used compositional data of e. coli, we have directly measured amino acid, fatty acid and lipid profiles for salmonella lt <dig>  and incorporated them into our reconstruction in irr <dig>  since fatty acid synthesis has also become a recent target of new antimicrobial drugs, it is critical to represent fatty acids and lipids accurately for different pathogens.
in addition to accounting for these major aspects of cell metabolism and biosynthesis, the reconstruction also includes reactions required for drug efflux, proton pumps, and inhibitory effects of antibiotics on bacterial metabolism.
we have included the acrab-tolc multi-component drug efflux pump which has been characterized in s. typhimurium.
this efflux pump has been shown to be involved in the resistance to β-lactams , fluoroquinolone , cefotaxime, and chloramphenicol  <cit> .
an increase in the expression levels of this pump has been recognized as a significant determinant of antibiotic resistance in bacterial pathogens.
the inclusion of this efflux pump allows modeling of the reaction of salmonella to antibiotics, by varying efflux pump activity and/or limiting metabolic flux through the respective pathways.
for example, β-lactam antibiotics result in incomplete peptidoglycan synthesis , which can be modeled by restricting flux through peptidoglycan synthesis.
similarly, dihydropteroate synthase inhibition by sulfonamides and dihydrofolate reductase inhibition by trimethoprim can be simulated by constraining fluxes through their corresponding reactions.
another unique aspect of our reconstruction are the reactions involved in free radical quenching that are generally used by salmonella to counteract the host cell response to infection  <cit> .
in addition to the classically known mechanisms mediated by superoxide dismutases, catalases and the enzymes that make and recycle glutathione and thioredoxin, our reconstruction includes mechanisms for reactive nitrogen species  and reactive oxygen species  resistance and removal.
some of those include: a soluble flavohemoglobin  with a no dioxygenase activity that is induced solely by rns; and an alkyl hydroperoxide reductase  belonging to the peroxiredoxin family that is involved in dissipating both ros and rns  <cit> .
the inclusion of reactions for ros  and rns  species make this model a suitable platform to analyze their role in infection and pathogenesis.
growth predictions
flux balance analysis  was used to predict growth rates on six different carbon sources  using their measured substrate uptake rates.
these in silico predicted growth rates were then compared to experimental results of aerobic batch cultivation of s. typhimurium lt <dig> in m <dig> minimal media.
the calculated maximum specific growth rate  of salmonella for a given carbon source uptake rate was close to the experimentally determined values for s. typhimurium under similar conditions .
quantitative differences between calculated and observed growth rates could be attributed to differences in atp maintenance parameters  or secretion of by-products, which were not predicted in any of the simulations since carbon sources were assumed to be the limiting nutrient.
s. typhimurium is a mixed acid fermenter that metabolizes glucose via the embden-meyerhoff-parnas pathway  <cit> .
many of the predicted non-zero fluxes  corresponded to genes whose mrna transcript expression levels were measurable during growth of salmonella  on glucose minimal media.
most of the transcripts detected in our study concur with previous reports in literature.
during growth on glucose the irr <dig> model predicts that oxaloacetate can be produced by phosphoenolpyruvate carboxylase , which is known to be the only enzyme serving the anaplerotic role of replenishing oxaloacetate that is withdrawn for amino acid biosynthesis.
in agreement with our flux calculations, our experimental results also indicate a high abundance of the pepc mrna.
the data also indicate that glucose is converted to biomass and some part of the carbon is diverted into stoichiometric formation of formate, acetate, ethanol and succinate under anaerobic conditions  <cit> .
growth of s. typhimurium on  <dig> pentoses  was also tested experimentally.
growth was only observed on d-xylose and d-ribose, which is consistent with literature and model predictions  <cit> .
xylose is catabolized, as most other aldopentoses, by transport into the cell, isomerization and phosphorylation of ketopentose .
the phosphorylated derivative is then converted to xylulose-5-phosphate by an epimerase and channeled into the pentose phosphate pathway.
d-ribose is transported into the cell, concomitant with the disruption of a high-energy phosphate bond and enters the pentose phosphate pathway following epimerization.
these results highlight the utility of using modeling to predict growth on substrates a priori to experimentation and can help understand the different pathways that are operational during growth on specific carbon and nitrogen sources.
qualitative growth phenotypes
fba was also used to predict which of the  <dig> carbon containing compounds and  <dig> nitrogen-containing compounds, with transport reactions  that are included in the model, could be used for growth as sole carbon or nitrogen sources.
the in silico predictions were then compared to  existing high-throughput phenotyping data from biolog's phenotyping microarray™ technology  and  reports from experimental literature  <cit> ; .
comparisons between model predictions  and experimental data  taken from biolog and published literature
the model predicted that  <dig> out of  <dig> potential carbon sources could support growth as the sole carbon source, and  <dig> out of  <dig> potential nitrogen sources could support growth as the sole nitrogen source.
comparisons with experimental data could be made for  <dig> of the potential carbon sources and for  <dig> of the potential nitrogen sources.
overall, the model predicted carbon and nitrogen sources with 82% accuracy for carbon sources and 81% accuracy for nitrogen sources.
the relationship between model predictions and experimental growth phenotypes is statistically significant .
most of the inaccurately predicted carbon sources are compounds that can also serve as sole nitrogen sources experimentally .
this suggests that regulation or enzyme kinetic limitations likely prevent these compounds from being used by s. typhimurium lt <dig> as sole carbon sources.
s. typhimurium and e. coli are closely related microorganisms.
comparison of their physiological data and metabolic network structures suggest differences in carbon and nitrogen source utilization.
for example, e. coli can utilize d-alanine, l-asparagine, l-aspartate, l-glutamine, and d-lactate as sole carbon sources experimentally, while s. typhimurium can not, even though s. typhimurium has the necessary enzymes and transporters.
these phenotypic differences are likely to be caused by differences in how the organisms control metabolism through regulation and kinetics rather than by differences in metabolic network structure.
the same explanation can be made for carbon sources that only s. typhimurium can use even though both organisms can catalyze the required reactions .s.
typhimurium is also predicted to be able to utilize additional carbon sources due to the presence of unique transporters and/or enzymes .
analysis of gene deletions, essentiality and lethal genes
the impact of single gene deletions on growth of s. typhimurium lt <dig> was analyzed by simulating its growth on rich medium in the absence of the reactions catalyzed by the corresponding gene.
the results were compared for consistency with a compiled list of  <dig> essential genes from a large-scale gene essentiality screen.
of the  <dig> essential genes, only  <dig> were included in irr <dig> and could be compared with experimental results.
the total accuracy of essentiality predictions was around 40% .
the accuracy of predicting in vivo essential genes for growth in rich media is within the range reported for models of pseudomonas aeroginosa   <cit> , h. pylori  <cit> , and bacillus subtilis  <cit> .
discrepancies between model predictions and experimental results can be used to better identify what metabolites are present in the rich media and lead to model improvements.
for example, inclusion of charged trnas rather than free amino acids into the biomass reaction would lead to the correct prediction of amino acid trna synthetases as essential genes.
to further explore the metabolic differences between e. coli and s. typhimurium we compared gene essentiality predictions for in silico aerobic growth in glucose minimal media using the most recent e. coli metabolic reconstruction, iaf <dig>  <cit>  and irr <dig> .
the number of predicted essential genes was fairly similar for the two organisms, with e. coli having  <dig>  essential metabolic genes and s. typhimurium having  <dig>  essential metabolic genes.
the two models have  <dig> overlapping orthologs, meaning there are  <dig> genes included in irr <dig> which do not have an ortholog included in iaf <dig> .
for  <dig> of the  <dig> overlapping orthologs both models predict the genes to be non-essential for in silico growth in a glucose minimal medium under aerobic conditions.
the models also agree on  <dig> of the  <dig> orthologs being essential.
comparing the predictions to experimental e. coli knock-out mutant phenotypes .
the discrepancies between essentiality predictions made by the two models can be useful in identifying model differences  and metabolic network differences  for the two microorganisms.
in  <dig> of the  <dig> overlapping orthologs, the e. coli model predicted genes as essential and the s. typhimurium model predicted those genes as non-essential; in  <dig> cases the s. typhimurium model predicted genes as essential and the e. coli model predicted those as non-essential.
most of these discrepancies  can be attributed to differences in the biomass formulations that were used for the two models.
the s. typhimurium model either has additional compounds  that cannot be produced in a gene's absence , or the e. coli model has additional compounds  that cannot be produced in a gene's absence .
reaction reversibility differences are also responsible for some of the discrepancies.
in the s. typhimurium model threonine aldolase is reversible so that threonine can be produced from glycine and acetaldehyde; however, this reaction was changed to irreversible in the updated e. coli model, iaf <dig>  as a result the biosynthetic genes for threonine production encoding homoserine kinase  and threonine synthase  are predicted to be essential by the e. coli model.
in contrast, the s. typhimurium model predicts that in their absence threonine can be synthesized from glycine and acetaldehyde.
since both thrb and thrc are essential for e. coli growth in minimal media conditions  <cit>  it is likely that threonine aldolase does not produce threonine.
glycine hydroxymethyltransferase also differs in its reversibility between the two models, where it is reversible in iaf <dig> and irreversible in irr <dig>  the reversible reaction provides another pathway in the e. coli model for serine production in the absence of phosphoserine phosphatase  or 3-phosphoglycerate dehydrogenase , and explains why these two enzymes are predicted to be essential for s. typhimurium growth and non-essential for e. coli growth.
since both sera and serb are essential experimentally in e. coli  <cit> , it is likely that neither organism uses glycine hydroxymethyltransferase for serine production, which is consistent with the s. typhimurium model predictions.
gene essentiality discrepanices caused by metabolic network differences, such as additional reactions and/or isozymes, provide interesting insights.
for example, s. typhimurium has genes encoding for two separate 3-isopropylmalate dehydratases  while e. coli only has genes for one .
as a result, when stm <dig> or stm <dig> are deleted the irr <dig> model predicts that stm <dig> and stm <dig> can be used instead to produce leucine, which is not possible in e. coli since it lacks the second isozyme.
s. typhimurium also carries out reactions that e. coli can not.
for example, an nadph dependent sulfite reductase  is predicted to be used for hydrogen sulfide production  when the nadh dependent sulfite reductase  is absent.
there are also cases where e. coli has additional isozymes or metabolic reactions, that explain why genes are predicted to be non-essential in e. coli but essential in s. typhimurium.
for example, there are multiple carbonic anhydrases and dihydrofolate reductases in e. coli while s. typhimurium has only one copy of each of these essential enzymes.
e. coli also contains a gene encoding ribose  <dig> -bisphosphokinase, for which there is no ortholog in s. typhimurium.
this e. coli specific enzyme is predicted by the iaf <dig> model to compensate for ribose-phosphate pyrophosphokinase in 5-phospho-alpha-d-ribose 1-diphosphate production, which is needed for histidine and nucleotide biosynthesis in the absence of ribose-phosphate pyrophosphokinase  <cit> .
as shown here, comparisons between gene essentiality predictions can be useful in identifying model discrepancies which can lead to updated and improved models, as well as identification of unique metabolic features, which allow for discovery of possible phenotypic differences between two organisms.
metabolic gene essentiality for virulence
the identification of genes essential for survival is also important for understanding the minimal requirements for pathogen survival in the host cell and to subsequently identify targets for new antibiotics.
a list of genes predicted to be essential for in silico growth of irr <dig> was compiled and compared with a set of in vivo essential and non-essential genes from the literature .
genes essential in vivo are those that are reported to be required for growth.
hence, salmonella mutants defective in these essential genes are avirulent in macrophage or animal models, while mutants defective in non-essential genes are still virulent and can grow.
an attenuated mutant strain is capable of growth in vivo at a lower rate than the wildtype.
this is consistent with the physiological definition of virulent, avirulent and attenuated strains in literature  <cit> .
since the availability of nutrients in the macrophage or animal host has not been completely delineated we have used an approximate representation based on that indirectly inferred from existing literature .
model predicted growth rates are highly sensitive to nutrient uptake rates which are unknown in the host celll.
the model was not used to distinguish between attenuated and virulent phenotypes as this would require detailed nutrient uptake rate measurements.
instead mutants were only predicted to be virulent  or avirulent .
the irr <dig> reconstruction as described earlier only contains metabolic genes, so there are additional non-metabolic essential genes in vivo that were not included in our virulence prediction study.
the in silico virulence predictions matched in vivo virulence characterization for 92%  of the in vivo virulent mutants and 56%  of the in vivo avirulent mutants.
most attenuated strains were also correctly predicted by the model to still be able to grow .
again a statistically significant relationship between model predictions and experimental observations was found .
the model predicted incorrectly that mutants defective in pmrf and udg were avirulent  <cit> .
both of these genes are involved in the biosynthesis of the aminoarabinose component of lipopolysaccharide .
however, this aminoarabinose component of the lps is not always present.
if alternate lps structures were also included in the model then the model would accurately predict the virulence of pmrf and udg mutant strains.
this would increase the final virulence prediction accuracy to 100%.
the consistencies and discrepancies between in silico and in vivo mutant phenotypes deserve further discussion and may generate potential hypotheses that can subsequently be tested.
some genes that are involved in major central metabolic pathways and are non-essential for intracellular survival  are discussed below.
acea encodes isocitrate lyase, which catalyzes the first step in the glyoxylate shunt, an anaplerotic carbon assimilatory pathway that synthesizes c <dig> dicarboxylates from c <dig> compounds like acetate or short chain fatty acids.
acea is required only for chronic salmonella infections and not for acute infections  <cit> .
this enzyme is generally used when substrates such as glucose or pyruvate are absent and acetate or fatty acids are present.
acka and pta are two enzymes in salmonella that lead to acetyl phosphate formation and the double deletion strain is virulent.
acetyl-phosphate is thought to activate expression of invasion genes, which can also be activated via a sensory kinase bara.
however, virulence is decreased only when pta, acka, and bara are deleted  <cit> .
phagocytic cells are able to generate superoxide, hydrogen peroxide and other reactive oxygen species  and pathogens possess multiple defenses to such oxidative stresses.
the model also accurately predicted that catalase  mutants  <cit>  and super oxide dismutase mutants   <cit>  are virulent.
mutants lacking a third superoxide dismutase  are attenuated and the model incorrectly predicts that a sodci mutant would be virulent; however, if the host-cell environment included superoxide and none of the other superoxide dismutases were expressed or active  then the bacterium would have no way of dissipating superoxide and the model would predict an avirulent phenotype.
one of the incorrect predictions for the experimentally determined avirulent strains was for the naga mutant.
naga is a n-acetylglucosamine-6-phosphate deacetylase, and it's substrate acts as an inducer of the nag operon, which relieves repression by nagc.
it has been proposed that naga mutants indirectly lead to reduction in aminosugar biosynthesis, which is needed in membrane structures  <cit> .
the model's inability to identify this mutant correctly as being avirulent could be due to the absence of any regulatory mechanisms not accounted for in the model.
inclusion of simple regulatory rules as has been done for the e. coli model, would allow us to address the role of regulation of metabolism during salmonella-host interaction.
of the experimentally determined attenuated mutants, the model predicted that mutants defective in aromatic amino acid biosynthesis are avirulent .
aro salmonella strains have been used as oral vaccines  <cit> .
the attenuation of this strain, which is accomplished by interrupting pathways for biosynthesis of aromatic amino acids, is consistent with our in silico analysis.
strains of s. typhimurium with complete non-reverting blocks in aromatic biosynthesis are avirulent  <cit>  as predicted by our model.
a large portion of the experimentally determined avirulent mutants which were incorrectly predicted by the model to be virulent are associated with nucleotide biosynthesis: guaab, purab, and thya .
these mutants were all predicted to be virulent due to the inclusion of thymidine, guanosine, adenosine, and deoxyadenosine into simulations of the host-cell environment.
if these metabolites are not available, then the model correctly predicts that these strains are avirulent.
other mutant strains defective in purine biosynthesis genes  are attenuated rather than avirulent  <cit> .
the experimental results indicate that intermediates in the purine biosynthesis pathway  are available in the host-cell environment and that guanosine and adenosine are likely produced from these pathway intermediates rather than provided by the host .
better accounting for nutrient availability in vivo is then able to resolve  <dig> % of the mutants that were incorrectly predicted in silico to be virulent.
this approach can thus also serve as an indirect probe of the host cell intracellular environment in which the pathogen finds itself.
pseudogenes and their role in infection among serovars of salmonella
sixty-five genes that are included in irr <dig>  and thus present in salmonella serovar typhimurium, are missing or are pseudogenes in other salmonella serovars   <cit> .
this implies that these genes could be dispensable for virulence.
based on the in silico analyses of the effect of single gene deletions on growth of s. typhimurium in both rich media and host-cell environments, all  <dig> genes are predicted to be dispensable.
only one gene out of these  <dig>  pand, was predicted to be required for growth in minimal media.
pand is required for panthothenate synthesis, an intermediate in the biosynthesis of coa.
in e. coli this enzyme is not required for growth in lb but is required for growth in minimal media  <cit> .
presumably then, serovars lacking pand would require panthothenate from its host.
compounds that affect pantothenate transport could potentially be used as antibiotics.
salmonella metabolism in infection
the transcriptional status of in vivo derived bacteria may correlate with mechanisms of pathogen survival and proliferation within host cells.
such genome-scale expression profiles are indicative of the host cell environment as well as the pathogen's adaptation to the intracellular niche.
based on a complete transcriptional profile of intracellular salmonella following short-term  infection of macrophage cell lines, only 21% of the total genes in the genome are differentially expressed .
further, only about 6% of these are characterized metabolic genes of known function.
since it is known that some metabolic genes are essential for intracellular replication  <cit>  and that salmonella's niche in the host has limited nutrient potential, the pathogen appears to have efficiently adapted to minimize its metabolic capacity in order to survive and grow in the phagosome.
we sought to identify genes  and select pathways that may be crucial for salmonella's infection of and growth inside macrophage cell lines by quantitatively monitoring the expression of  <dig> selected metabolic genes  that are involved in different key pathways.
since only a small percentage of metabolic genes were hypothesized to change mrna abundance levels during infection, we use the genomelab expression profiler  for estimating gene expression quantitatively.
this technique combines rt-pcr and capillary gel electrophoresis in a high-throughput format to obtain accurate, consistent mrna abundance levels for expressed transcripts of chosen genes rapidly.
consistent with metabolic gene essentiality hypothesis suggested in literature, only 24% of the genes chosen in this study were up-regulated  <dig> hrs post infection.
our data indicates that isocitrate lyase , and fatty acid degradation genes fada, fadd/fade are not expressed.
although these genes have been implicated for virulence in chronic salmonella infection, they are known to be unchanged during short-term infection experiments  <cit> .
our results suggest the presence of glucose or pyruvate and not acetate or short chain fatty acids in vivo during the early stages of infection as suggested in literature  <cit> .
this also explains why genes from the oxidative branch of the tca cycle  are over expressed.
the simultaneous up regulation of succinate dehydrogenase and electron transport chain genes cyoe and cyoa is suggestive of an aerobic environment in vivo, at least during the early stages of infection  <cit> .
increased levels of f1f <dig> atp synthase subunits measured within our data set also support this.
atp synthase is primarily responsible for the oxidative phosphorylation of adp to atp.
it also plays a role in maintaining the proton gradient to maintain intracellular ph and drive efflux pumps.
the high levels of mrna transcripts for succinate dehydrogenase and protoheme farnesyl transferase genes also suggest heme biosynthesis occurs in s. typhimurium during infection.
these are in concurrence with other existing data sets  <cit> .
genes belonging to the histidine biosynthesis cluster  were up regulated and those involved in aromatic amino acid synthesis pathways  were down regulated.
these findings are consistent with the literature  <cit> .
genes involved in nucleotide  and nad biosynthesis  were down regulated or not detected.
the down regulated genes also included those involved in glycolysis and carbon metabolite repression, certain amino acid decarboxylases and some transporter genes including drug efflux genes.
these results suggest that salmonella acquires a significant amount of essential nutrients directly from the host cell, thus shutting down several biosynthetic pathways when inside the host.
the genes that cannot be detected include some genes involved in glycolysis/gluconeogenesis  and anaplerotic reactions .
other genes are detectable, but do not change in expression level over the course of infection, including those in the pentose phosphate pathway , entner duoduroff pathway  and glutathione cycle .
fatty acid biosynthesis and degradation genes are also not detected during the monitored early stages of infection.
metabolic capacity of s. typhimurium in the host cell
in silico prediction of metabolic pathways utilized during infection allows us to identify redundant metabolism that the pathogen could exploit in order to survive and replicate.
during conditions of infection it is reasonable to assume that salmonella cells are under an energy limitation.
hence, atp generated is fully coupled to biosynthesis and growth  <cit> .
given the possible nutrients provided in the host-cell environment, we used flux variability analysis  to identify which reactions can be utilized  <cit>  and found that  <dig> reactions span a large range of flux values  when biomass production is optimal.
these  <dig> reactions  encompass all the possible reactions  that could be used by salmonella during infection.
when the biomass production is allowed to be suboptimal , an additional  <dig> reactions  could be used by salmonella in the host cell environment .
the constraining conditions of the simulation , were derived by literature mining to identify a media composition that could closely mimic nutrient availability for s. typhimurium in the host cell.
we also used fva to identify blocked reactions and thus pathways that will never be used by the pathogen under the described conditions in vivo during infection.
these reactions are typically those catalyzed by proteins that cannot be detected in the proteome or those genes that are down regulated under similar conditions.
our model simulation results of these blocked reactions  along with our gene expression data identify correctly fatty acid degradation metabolic genes  to be non-operational during the above conditions.
thus, there is some correlation between blocked reactions found by fva and unexpressed genes during infection.
we compiled the list of all proteins that comprise the in silico reactome and compared that to previously determined proteome of s. typhimurium extracted from macrophage cells after infection .
of the  <dig> proteins identified with a high level of confidence in s. typhimurium  <cit> ,  <dig> have an annotated metabolic function  and in our model irr <dig>  eighty of these proteins were associated with reactions present in the reactome calculated by fva analysis during optimal growth, while only  <dig> were identified in the host-cell reactome when growth was sub-optimal.
in another experiment, where the proteome of s. typhimurium was probed under different in vitro environmental conditions, a total of  <dig> proteins were identified  <cit> , of which  <dig> were accounted for in our metabolic reconstruction.
out of these  <dig> proteins,  <dig> and  <dig> were associated with the optimal and suboptimal host cell reactome, respectively, as calculated using fva.
most of the proteins that were found in both proteomic data sets and associated with blocked reactions, as calculated by fva, were trna synthetases .
incorporating charged trnas instead of free amino acids in the biomass reaction would change the trna synthetase reactions from blocked to part of the optimal reactome and would result in higher agreement between the optimal reactome and identified proteins.
using a similar strategy we compared micro-array gene expression data obtained from s. typhimurium extracted from macrophage cells  <cit> .
of the  <dig> genes that were differentially expressed,  <dig> genes had a metabolic function, thus comprising the metabolic transcriptome.
when we compared those with the gene set represented by our fva analysis reactome, we identified 80% of the genes in the metabolic transcriptome were associated with reactions in the reactome.
of the remaining  <dig> that were not present in our in silico reactome, 65% of the genes were down-regulated and formed a subset of the blocked reactions from our fva analysis .
although the accuracy of prediction is very good, it could be increased by improving our knowledge of existing nutrients in the macrophage or incorporating regulation into the model.
by identifying preferred metabolic pathways we can better understand how the pathogen responds to its microenvironment during invasion and infection.
discussion
this study describes a genome-scale metabolic model, irr <dig>  of the gram-negative intracellular pathogen s. typhimurium.
the model has been reconstructed by integration of existing genomic, proteomic, and phenotypic experimental data.
comparison between model predictions and experimental data demonstrated an accuracy of more than 80% for growth and virulence phenotypes.
the model allowed us to use flux balance analysis to analyze growth and virulence genes given the environment the intracellular pathogen encounters.
further, gene expression experiments and computational analysis allowed us to identify metabolic enzymes and pathways that are potentially active during infection.
model-development is an iterative process that moves from a genome annotation and biochemistry of the organism to the development of a genome-scale metabolic network that can be used as a fully functional in silico organism to compute physiological phenotypes.
our model irr <dig> serves also as a framework to integrate further gene expression and proteomic data and explain the expression level of specific gene products during infection.
computing cellular functions from reconstructed networks allowed us to understand the relationship between environmental factors  and genetic factors  and their integrated function to orchestrate cell phenotypes .
knowledge of the functional status of metabolic genes under defined constraints allows us to interpret the role of critical nodes in a cellular metabolic network.
the function of redundant pathways under internally perturbed conditions, such as in the host cell during infection, provides important insights into metabolic network structure.
although this internal robustness has been experimentally demonstrated, building a complete catalog of genes essential for growth under a myriad of conditions for any organism is time-consuming and laborious.
a computational approach that predicts gene essentiality with high accuracy is thus of great value.
model predictions of virulence for different mutants improved from 80% to 91% upon altering host-cell environment components  and this shows how the model can be used as a tool for providing insights into the role of metabolism in infection and pathogenesis.
to date, constraint-based models cannot predict effects of toxic intermediates and methods are needed that allow for prediction of concentrations within constraint-based models.
a biomass reaction  is needed to accurately predict the effects of gene deletions on growth and viability.
although predictions have been made in the past with approximations of the biomass function to that of e. coli, there are some major differences in fatty acid and amino acid composition that require experimental measurements to determine the drain of different metabolites and the cost of their production on the function of the organism.
our study has addressed this issue by directly determining the components of salmonella biomass, thus increasing the validity and predictive accuracy of model irr <dig>
s. typhimurium is a pathogen that is the leading cause of human gastroenteritis and is used as a mouse model for human typhoid fever.
although more than  <dig> serovars exist, the lt <dig> strain, whose genome sequence has been used here for building model irr <dig>  is one of the principal strains for cellular and molecular biology of salmonella.
attenuated mutants have been used not only as live oral vaccines against salmonella infection but also to express antigens from other pathogens and to deliver proteins to solid tumors.
one reason for attenuation or avirulence is the deletion or functional silencing of a genetic element in the bacterium resulting in the impairment of the intracellular survival capacity of s. typhimurium.
our in silico analysis shows the ability to predict such effects with good accuracy.
at this point the host cell metabolic environment has been inferred indirectly from gene-expression data.
improved accuracy in directly determining what metabolites are available in the host cell environment and incorporation of transcriptional regulatory rules into irr <dig> will allow a more precise analysis of metabolic selective pressures a pathogen encounters during infection.
in addition to pathogenicity and virulence analysis, the s. typhimurium metabolic reconstruction can be a valuable tool to interpret several 'omics' data types.
although there is not a direct correlation between gene or protein expression levels and metabolic fluxes, the genome-scale model serves as a platform to visualize these large data sets and interpret them in the context of operational, redundant or incomplete metabolic pathways.
this provides a better understanding of cellular phenotypes and offers insights into how system level properties affect survival or sensitivity of pathogens to factors in the environment provided by the host cell.
for example, the model prediction of fad genes not being operational during early stages of infection  and operational during sub-optimal growth can be explained by the fact that fatty acids are not utilized in the initial stages of infection by s. typhimurium  <cit> .
this observation along with the calculated host cell reactome from irr <dig> could suggest a switch from optimal to sub-optimal growth or a change in nutrients available for salmonella inside host cell tissue.
this switch is concomitant with a marked decrease in growth during persistence or chronic infection.
growth rates are presumably higher  during the initial stage of infection.
the entire reactome of salmonella as predicted by the model provides many hypotheses that need to be further investigated to identify the implications of metabolic pathways and functions in pathogenesis.
metabolomic data have become more readily available and easier to measure.
they can be used to apply additional constraints to the model  <cit> , thus improving model predictions of cell function and global phenotypes.
the genome-scale model of s. typhimurium we have developed here also presents us with the opportunity to interrogate several other serovars or species that belong to the salmonella genus.
for example, pand is the only gene that is non-functional in salmonella typhi, and is predicted to be required for virulence in specific conditions.
the presence of pseudogenes in other serovars can be analyzed in the context of the genome-scale model to present us with systemic similarities and differences that may or may not play a role in the pathogenic potential of all the four species.
serovars lacking pand would require panthothenate from its host.
thus, compounds that affect pantothenate transport could potentially be used as antibiotics.
such analyses across different serovars require genome-scale models and might accelerate the identification of rational drug targets.
although many of the reactions in the irr <dig> model for s. typhimurium are conserved in models for e. coli mg <dig>  there are unique metabolic and genetic differences that affect growth phenotypes in response to different environmental or genetic perturbations.
the two organisms are each capable of utilizing unique carbon sources, which are related to differences in metabolic networks  or differences in regulation or kinetics that affect flow through shared reactions common to both organisms .
we also found that there are differences in the number of isozymes or alternative reactions/pathways that enable one organism to be more robust against gene deletions in a glucose minimal media condition.
we anticipate that comparisons of model-predicted essential genes in more environmental conditions will lead to the identification of more metabolic differences between the two organisms, as the analysis was only performed for aerobic growth in glucose minimal media.
these metabolic differences might shed light onto the evolutionary differences between the two organisms and the nutrients they are or were exposed to.
the genome-scale metabolic network reconstruction presented here supports the notion of a robust minimal set of pathways that salmonella requires for intracellular survival during infection  <cit> .
our model allows in silico simulation of different nutritional requirements and constraints as they may be present in different host cell types or tissues that are difficult to directly address experimentally.
our model can generate hypotheses that can be validated by using mutant strains, high-throughput data, and selective analysis of pathways.
such an endeavor has the potential to lead to an improved understanding of how intracellular pathogens utilize the existing microenvironment of the host.
this in turn might lead to the identification of conserved metabolic pathways that can be targeted by novel interfering compounds.
CONCLUSIONS
we have presented here extensive datasets within a reconstruction of a genome-scale metabolic model for s. typhimurium .
this model has been validated by experimental data and by legacy data.
it demonstrates an accuracy of 80% for growth and virulence phenotypes.
thus the model irr <dig> provides a suitable platform for computing cellular phenotypes for salmonella and for further integration of high-throughput data.
biomass composition of salmonella was determined experimentally, which increased the accuracy of our model building and also showed distinctive difference to that of e. coli.
constraint-based analysis using the metabolic reconstructions irr <dig> for s. typhimurium and mg <dig> for e. coli identified unique metabolic differences between these two bacteria.
expression profiling of selective key metabolic enzymes in salmonella during infection of macrophages identified operational and non-operational metabolic pathways.
flux variabity analysis of metabolic pathways identified a minimal set of reactions required for salmonella intracellular survival.
this study thus provides another high-resolution reconstruction of microbial metabolism of an intracellular pathogen next to that of m. tuberculosis and demonstrates its suitability for analyzing bacterial metabolic networks during infection.
from further validation and testing of hypotheses generated from these reconstruction efforts we expect to arrive at a conserved set of metabolic reactions that interdigitate with those of the host cell and might thus provide novel targets for therapeutic intervention.

 salmonellae are gram-negative bacterial pathogens with a wide host range. the ability of salmonella enterica to cause disease is closely correlated with its ability to survive inside its host cells, which include epithelial cells and macrophages. infection of mice with salmonella enterica serovar typhimurium serves as model for the systemic disease typhoid fever. salmonella strains that are deficient for intra-macrophage survival cannot elicit disease in the mouse model. many genomic and proteomic datasets exist for salmonella as a pathogen. most experiments have focused on identifying'virulence' factors that might drive the disease process, such as surface structures or secreted products. differential expression of metabolic gene-products has largely been left unexplored in terms of their role in pathogenesis. however, it has now gained more attention that the metabolism of intracellular bacterial pathogens such as salmonella may hold important clues to understanding their intracellular life-style as well as to understanding host defense mechanisms. genome-scale metabolic in silico models that can integrate and exploit genomic and proteomic data may provide a first step towards a quantitative analysis of the metabolism of intracellular bacteria. constraint-based reconstruction and analysis of cellular metabolism has been employed successfully to develop organism-specific genome-scale in silico models across all the three major domains of the tree of life. the cobra approach seeks to differentiate between those network states that are achievable by the organism from those that are not, based on a detailed reconstruction of metabolism and incorporation of physiological parameters that are consistent with known experimental information. the solutions this analysis provides are based on the successive imposition of governing physicochemical constraints on genome-scale reconstructions. cobra integrates and represents our current knowledge of network components and interactions and has been applied to several prokaryotic and eukaryotic organisms, including escherichia coli, helicobacter pylori, haemophilus influenzae, saccharomyces cerevisiae and geobacter sulfurreducens. by systematizing and providing context for transcriptomic, proteomic, and metabolomic data, these models allow for the imposition of constraints, which together define the possible phenotypic behavior allowed by these in silico organisms. the use of constrained based modeling has significant potential to identify new antibacterial targets. modeling allows simulation of single gene knockouts and, more importantly, the effect of combinatorial deletions, which might allow the rational development of combination drug therapy. the use of such a modeling approach has been demonstrated for h. influenza and for mycobacterium tuberculosis. as the number of pathogens with such models grows, one might expect that more basic principles of bacterial metabolism needed for pathogenesis can be derived. here we describe the reconstruction of a detailed and validated genome-scale metabolic network for s. typhimurium, that reconciles established genomic and high-throughput phenotyping data and that is augmented by experimental data to facilitate iterative model building efforts. furthermore, we demonstrate how the model can be used to reconcile existing virulence data and to resolve gaps in our understanding of pathogen metabolism and host-cell nutrient environment. we also use gene expression data obtained for this study in combination with established microarray and proteomic data in the context of the model to understand pathogen metabolism in macrophage cell lines.  characteristics of the reconstructed metabolic network the -kilobase chromosome of salmonella enterica serovar typhimurium lt accounts for coding sequences including pseudogenes and a kb plasmid encoding orfs. the in silico metabolic reconstruction of the bacterium referred to as irr contains only annotated chromosomal genes. the characteristics of the network are detailed in figure 1a and 1b. it currently accounts for genes that encode proteins catalyzing unique metabolic and transport reactions in the bacterium. the reconstruction includes gene to protein to reaction association relationships, which distinguish between isozymes, multimeric proteins and protein complexes. the reconstruction contains all the biosynthesis pathways to make biomass components that are required for growth. a biomass production reaction was formulated, which is based on the measured biomass composition of the bacterium and which details the required metabolites and their relative amounts. the reconstruction accounts for central metabolism, carbohydrate utilization, amino acid biosynthesis and degradation, nucleotide biosynthesis and salvage, and cofactor biosynthesis. also included were anaerobic and aerobic electron transport chain reactions, transport reactions for various minerals and small molecules, and biosynthesis reactions for cell wall constituents and precursors, such as phospholipids and peptidoglycan. lipid biosynthesis in bacteria involves fatty acid biosynthesis, attachment of the complete fatty acids to sn-glycerol-3-phosphate, and addition and modification of head groups to yield the major cellular phospholipids. fatty acids are important constituents of lipid a and the cell wall of bacteria. in contrast to published models of other microbes that used compositional data of e. coli, we have directly measured amino acid, fatty acid and lipid profiles for salmonella lt and incorporated them into our reconstruction in irr since fatty acid synthesis has also become a recent target of new antimicrobial drugs, it is critical to represent fatty acids and lipids accurately for different pathogens. in addition to accounting for these major aspects of cell metabolism and biosynthesis, the reconstruction also includes reactions required for drug efflux, proton pumps, and inhibitory effects of antibiotics on bacterial metabolism. we have included the acrab-tolc multi-component drug efflux pump which has been characterized in s. typhimurium. this efflux pump has been shown to be involved in the resistance to β-lactams, fluoroquinolone, cefotaxime, and chloramphenicol. an increase in the expression levels of this pump has been recognized as a significant determinant of antibiotic resistance in bacterial pathogens. the inclusion of this efflux pump allows modeling of the reaction of salmonella to antibiotics, by varying efflux pump activity and/or limiting metabolic flux through the respective pathways. for example, β-lactam antibiotics result in incomplete peptidoglycan synthesis, which can be modeled by restricting flux through peptidoglycan synthesis. similarly, dihydropteroate synthase inhibition by sulfonamides and dihydrofolate reductase inhibition by trimethoprim can be simulated by constraining fluxes through their corresponding reactions. another unique aspect of our reconstruction are the reactions involved in free radical quenching that are generally used by salmonella to counteract the host cell response to infection. in addition to the classically known mechanisms mediated by superoxide dismutases, catalases and the enzymes that make and recycle glutathione and thioredoxin, our reconstruction includes mechanisms for reactive nitrogen species and reactive oxygen species resistance and removal. some of those include: a soluble flavohemoglobin with a no dioxygenase activity that is induced solely by rns; and an alkyl hydroperoxide reductase belonging to the peroxiredoxin family that is involved in dissipating both ros and rns. the inclusion of reactions for ros and rns species make this model a suitable platform to analyze their role in infection and pathogenesis. growth predictions flux balance analysis was used to predict growth rates on six different carbon sources using their measured substrate uptake rates. these in silico predicted growth rates were then compared to experimental results of aerobic batch cultivation of s. typhimurium lt in m minimal media. the calculated maximum specific growth rate of salmonella for a given carbon source uptake rate was close to the experimentally determined values for s. typhimurium under similar conditions. quantitative differences between calculated and observed growth rates could be attributed to differences in atp maintenance parameters or secretion of by-products, which were not predicted in any of the simulations since carbon sources were assumed to be the limiting nutrient. s. typhimurium is a mixed acid fermenter that metabolizes glucose via the embden-meyerhoff-parnas pathway. many of the predicted non-zero fluxes corresponded to genes whose mrna transcript expression levels were measurable during growth of salmonella on glucose minimal media. most of the transcripts detected in our study concur with previous reports in literature. during growth on glucose the irr model predicts that oxaloacetate can be produced by phosphoenolpyruvate carboxylase, which is known to be the only enzyme serving the anaplerotic role of replenishing oxaloacetate that is withdrawn for amino acid biosynthesis. in agreement with our flux calculations, our experimental results also indicate a high abundance of the pepc mrna. the data also indicate that glucose is converted to biomass and some part of the carbon is diverted into stoichiometric formation of formate, acetate, ethanol and succinate under anaerobic conditions. growth of s. typhimurium on pentoses was also tested experimentally. growth was only observed on d-xylose and d-ribose, which is consistent with literature and model predictions. xylose is catabolized, as most other aldopentoses, by transport into the cell, isomerization and phosphorylation of ketopentose. the phosphorylated derivative is then converted to xylulose-5-phosphate by an epimerase and channeled into the pentose phosphate pathway. d-ribose is transported into the cell, concomitant with the disruption of a high-energy phosphate bond and enters the pentose phosphate pathway following epimerization. these results highlight the utility of using modeling to predict growth on substrates a priori to experimentation and can help understand the different pathways that are operational during growth on specific carbon and nitrogen sources. qualitative growth phenotypes fba was also used to predict which of the carbon containing compounds and nitrogen-containing compounds, with transport reactions that are included in the model, could be used for growth as sole carbon or nitrogen sources. the in silico predictions were then compared to existing high-throughput phenotyping data from biolog's phenotyping microarray™ technology and reports from experimental literature;. comparisons between model predictions and experimental data taken from biolog and published literature the model predicted that out of potential carbon sources could support growth as the sole carbon source, and out of potential nitrogen sources could support growth as the sole nitrogen source. comparisons with experimental data could be made for of the potential carbon sources and for of the potential nitrogen sources. overall, the model predicted carbon and nitrogen sources with 82% accuracy for carbon sources and 81% accuracy for nitrogen sources. the relationship between model predictions and experimental growth phenotypes is statistically significant. most of the inaccurately predicted carbon sources are compounds that can also serve as sole nitrogen sources experimentally. this suggests that regulation or enzyme kinetic limitations likely prevent these compounds from being used by s. typhimurium lt as sole carbon sources. s. typhimurium and e. coli are closely related microorganisms. comparison of their physiological data and metabolic network structures suggest differences in carbon and nitrogen source utilization. for example, e. coli can utilize d-alanine, l-asparagine, l-aspartate, l-glutamine, and d-lactate as sole carbon sources experimentally, while s. typhimurium can not, even though s. typhimurium has the necessary enzymes and transporters. these phenotypic differences are likely to be caused by differences in how the organisms control metabolism through regulation and kinetics rather than by differences in metabolic network structure. the same explanation can be made for carbon sources that only s. typhimurium can use even though both organisms can catalyze the required reactions.s. typhimurium is also predicted to be able to utilize additional carbon sources due to the presence of unique transporters and/or enzymes. analysis of gene deletions, essentiality and lethal genes the impact of single gene deletions on growth of s. typhimurium lt was analyzed by simulating its growth on rich medium in the absence of the reactions catalyzed by the corresponding gene. the results were compared for consistency with a compiled list of essential genes from a large-scale gene essentiality screen. of the essential genes, only were included in irr and could be compared with experimental results. the total accuracy of essentiality predictions was around 40%. the accuracy of predicting in vivo essential genes for growth in rich media is within the range reported for models of pseudomonas aeroginosa, h. pylori, and bacillus subtilis. discrepancies between model predictions and experimental results can be used to better identify what metabolites are present in the rich media and lead to model improvements. for example, inclusion of charged trnas rather than free amino acids into the biomass reaction would lead to the correct prediction of amino acid trna synthetases as essential genes. to further explore the metabolic differences between e. coli and s. typhimurium we compared gene essentiality predictions for in silico aerobic growth in glucose minimal media using the most recent e. coli metabolic reconstruction, iaf and irr. the number of predicted essential genes was fairly similar for the two organisms, with e. coli having essential metabolic genes and s. typhimurium having essential metabolic genes. the two models have overlapping orthologs, meaning there are genes included in irr which do not have an ortholog included in iaf. for of the overlapping orthologs both models predict the genes to be non-essential for in silico growth in a glucose minimal medium under aerobic conditions. the models also agree on of the orthologs being essential. comparing the predictions to experimental e. coli knock-out mutant phenotypes. the discrepancies between essentiality predictions made by the two models can be useful in identifying model differences and metabolic network differences for the two microorganisms. in of the overlapping orthologs, the e. coli model predicted genes as essential and the s. typhimurium model predicted those genes as non-essential; in cases the s. typhimurium model predicted genes as essential and the e. coli model predicted those as non-essential. most of these discrepancies can be attributed to differences in the biomass formulations that were used for the two models. the s. typhimurium model either has additional compounds that cannot be produced in a gene's absence, or the e. coli model has additional compounds that cannot be produced in a gene's absence. reaction reversibility differences are also responsible for some of the discrepancies. in the s. typhimurium model threonine aldolase is reversible so that threonine can be produced from glycine and acetaldehyde; however, this reaction was changed to irreversible in the updated e. coli model, iaf as a result the biosynthetic genes for threonine production encoding homoserine kinase and threonine synthase are predicted to be essential by the e. coli model. in contrast, the s. typhimurium model predicts that in their absence threonine can be synthesized from glycine and acetaldehyde. since both thrb and thrc are essential for e. coli growth in minimal media conditions it is likely that threonine aldolase does not produce threonine. glycine hydroxymethyltransferase also differs in its reversibility between the two models, where it is reversible in iaf and irreversible in irr the reversible reaction provides another pathway in the e. coli model for serine production in the absence of phosphoserine phosphatase or 3-phosphoglycerate dehydrogenase, and explains why these two enzymes are predicted to be essential for s. typhimurium growth and non-essential for e. coli growth. since both sera and serb are essential experimentally in e. coli, it is likely that neither organism uses glycine hydroxymethyltransferase for serine production, which is consistent with the s. typhimurium model predictions. gene essentiality discrepanices caused by metabolic network differences, such as additional reactions and/or isozymes, provide interesting insights. for example, s. typhimurium has genes encoding for two separate 3-isopropylmalate dehydratases while e. coli only has genes for one. as a result, when stm or stm are deleted the irr model predicts that stm and stm can be used instead to produce leucine, which is not possible in e. coli since it lacks the second isozyme. s. typhimurium also carries out reactions that e. coli can not. for example, an nadph dependent sulfite reductase is predicted to be used for hydrogen sulfide production when the nadh dependent sulfite reductase is absent. there are also cases where e. coli has additional isozymes or metabolic reactions, that explain why genes are predicted to be non-essential in e. coli but essential in s. typhimurium. for example, there are multiple carbonic anhydrases and dihydrofolate reductases in e. coli while s. typhimurium has only one copy of each of these essential enzymes. e. coli also contains a gene encoding ribose -bisphosphokinase, for which there is no ortholog in s. typhimurium. this e. coli specific enzyme is predicted by the iaf model to compensate for ribose-phosphate pyrophosphokinase in 5-phospho-alpha-d-ribose 1-diphosphate production, which is needed for histidine and nucleotide biosynthesis in the absence of ribose-phosphate pyrophosphokinase. as shown here, comparisons between gene essentiality predictions can be useful in identifying model discrepancies which can lead to updated and improved models, as well as identification of unique metabolic features, which allow for discovery of possible phenotypic differences between two organisms. metabolic gene essentiality for virulence the identification of genes essential for survival is also important for understanding the minimal requirements for pathogen survival in the host cell and to subsequently identify targets for new antibiotics. a list of genes predicted to be essential for in silico growth of irr was compiled and compared with a set of in vivo essential and non-essential genes from the literature. genes essential in vivo are those that are reported to be required for growth. hence, salmonella mutants defective in these essential genes are avirulent in macrophage or animal models, while mutants defective in non-essential genes are still virulent and can grow. an attenuated mutant strain is capable of growth in vivo at a lower rate than the wildtype. this is consistent with the physiological definition of virulent, avirulent and attenuated strains in literature. since the availability of nutrients in the macrophage or animal host has not been completely delineated we have used an approximate representation based on that indirectly inferred from existing literature. model predicted growth rates are highly sensitive to nutrient uptake rates which are unknown in the host celll. the model was not used to distinguish between attenuated and virulent phenotypes as this would require detailed nutrient uptake rate measurements. instead mutants were only predicted to be virulent or avirulent. the irr reconstruction as described earlier only contains metabolic genes, so there are additional non-metabolic essential genes in vivo that were not included in our virulence prediction study. the in silico virulence predictions matched in vivo virulence characterization for 92% of the in vivo virulent mutants and 56% of the in vivo avirulent mutants. most attenuated strains were also correctly predicted by the model to still be able to grow. again a statistically significant relationship between model predictions and experimental observations was found. the model predicted incorrectly that mutants defective in pmrf and udg were avirulent. both of these genes are involved in the biosynthesis of the aminoarabinose component of lipopolysaccharide. however, this aminoarabinose component of the lps is not always present. if alternate lps structures were also included in the model then the model would accurately predict the virulence of pmrf and udg mutant strains. this would increase the final virulence prediction accuracy to 100%. the consistencies and discrepancies between in silico and in vivo mutant phenotypes deserve further discussion and may generate potential hypotheses that can subsequently be tested. some genes that are involved in major central metabolic pathways and are non-essential for intracellular survival are discussed below. acea encodes isocitrate lyase, which catalyzes the first step in the glyoxylate shunt, an anaplerotic carbon assimilatory pathway that synthesizes c dicarboxylates from c compounds like acetate or short chain fatty acids. acea is required only for chronic salmonella infections and not for acute infections. this enzyme is generally used when substrates such as glucose or pyruvate are absent and acetate or fatty acids are present. acka and pta are two enzymes in salmonella that lead to acetyl phosphate formation and the double deletion strain is virulent. acetyl-phosphate is thought to activate expression of invasion genes, which can also be activated via a sensory kinase bara. however, virulence is decreased only when pta, acka, and bara are deleted. phagocytic cells are able to generate superoxide, hydrogen peroxide and other reactive oxygen species and pathogens possess multiple defenses to such oxidative stresses. the model also accurately predicted that catalase mutants and super oxide dismutase mutants are virulent. mutants lacking a third superoxide dismutase are attenuated and the model incorrectly predicts that a sodci mutant would be virulent; however, if the host-cell environment included superoxide and none of the other superoxide dismutases were expressed or active then the bacterium would have no way of dissipating superoxide and the model would predict an avirulent phenotype. one of the incorrect predictions for the experimentally determined avirulent strains was for the naga mutant. naga is a n-acetylglucosamine-6-phosphate deacetylase, and it's substrate acts as an inducer of the nag operon, which relieves repression by nagc. it has been proposed that naga mutants indirectly lead to reduction in aminosugar biosynthesis, which is needed in membrane structures. the model's inability to identify this mutant correctly as being avirulent could be due to the absence of any regulatory mechanisms not accounted for in the model. inclusion of simple regulatory rules as has been done for the e. coli model, would allow us to address the role of regulation of metabolism during salmonella-host interaction. of the experimentally determined attenuated mutants, the model predicted that mutants defective in aromatic amino acid biosynthesis are avirulent. aro salmonella strains have been used as oral vaccines. the attenuation of this strain, which is accomplished by interrupting pathways for biosynthesis of aromatic amino acids, is consistent with our in silico analysis. strains of s. typhimurium with complete non-reverting blocks in aromatic biosynthesis are avirulent as predicted by our model. a large portion of the experimentally determined avirulent mutants which were incorrectly predicted by the model to be virulent are associated with nucleotide biosynthesis: guaab, purab, and thya. these mutants were all predicted to be virulent due to the inclusion of thymidine, guanosine, adenosine, and deoxyadenosine into simulations of the host-cell environment. if these metabolites are not available, then the model correctly predicts that these strains are avirulent. other mutant strains defective in purine biosynthesis genes are attenuated rather than avirulent. the experimental results indicate that intermediates in the purine biosynthesis pathway are available in the host-cell environment and that guanosine and adenosine are likely produced from these pathway intermediates rather than provided by the host. better accounting for nutrient availability in vivo is then able to resolve % of the mutants that were incorrectly predicted in silico to be virulent. this approach can thus also serve as an indirect probe of the host cell intracellular environment in which the pathogen finds itself. pseudogenes and their role in infection among serovars of salmonella sixty-five genes that are included in irr and thus present in salmonella serovar typhimurium, are missing or are pseudogenes in other salmonella serovars. this implies that these genes could be dispensable for virulence. based on the in silico analyses of the effect of single gene deletions on growth of s. typhimurium in both rich media and host-cell environments, all genes are predicted to be dispensable. only one gene out of these pand, was predicted to be required for growth in minimal media. pand is required for panthothenate synthesis, an intermediate in the biosynthesis of coa. in e. coli this enzyme is not required for growth in lb but is required for growth in minimal media. presumably then, serovars lacking pand would require panthothenate from its host. compounds that affect pantothenate transport could potentially be used as antibiotics. salmonella metabolism in infection the transcriptional status of in vivo derived bacteria may correlate with mechanisms of pathogen survival and proliferation within host cells. such genome-scale expression profiles are indicative of the host cell environment as well as the pathogen's adaptation to the intracellular niche. based on a complete transcriptional profile of intracellular salmonella following short-term infection of macrophage cell lines, only 21% of the total genes in the genome are differentially expressed. further, only about 6% of these are characterized metabolic genes of known function. since it is known that some metabolic genes are essential for intracellular replication and that salmonella's niche in the host has limited nutrient potential, the pathogen appears to have efficiently adapted to minimize its metabolic capacity in order to survive and grow in the phagosome. we sought to identify genes and select pathways that may be crucial for salmonella's infection of and growth inside macrophage cell lines by quantitatively monitoring the expression of selected metabolic genes that are involved in different key pathways. since only a small percentage of metabolic genes were hypothesized to change mrna abundance levels during infection, we use the genomelab expression profiler for estimating gene expression quantitatively. this technique combines rt-pcr and capillary gel electrophoresis in a high-throughput format to obtain accurate, consistent mrna abundance levels for expressed transcripts of chosen genes rapidly. consistent with metabolic gene essentiality hypothesis suggested in literature, only 24% of the genes chosen in this study were up-regulated hrs post infection. our data indicates that isocitrate lyase, and fatty acid degradation genes fada, fadd/fade are not expressed. although these genes have been implicated for virulence in chronic salmonella infection, they are known to be unchanged during short-term infection experiments. our results suggest the presence of glucose or pyruvate and not acetate or short chain fatty acids in vivo during the early stages of infection as suggested in literature. this also explains why genes from the oxidative branch of the tca cycle are over expressed. the simultaneous up regulation of succinate dehydrogenase and electron transport chain genes cyoe and cyoa is suggestive of an aerobic environment in vivo, at least during the early stages of infection. increased levels of f1f atp synthase subunits measured within our data set also support this. atp synthase is primarily responsible for the oxidative phosphorylation of adp to atp. it also plays a role in maintaining the proton gradient to maintain intracellular ph and drive efflux pumps. the high levels of mrna transcripts for succinate dehydrogenase and protoheme farnesyl transferase genes also suggest heme biosynthesis occurs in s. typhimurium during infection. these are in concurrence with other existing data sets. genes belonging to the histidine biosynthesis cluster were up regulated and those involved in aromatic amino acid synthesis pathways were down regulated. these findings are consistent with the literature. genes involved in nucleotide and nad biosynthesis were down regulated or not detected. the down regulated genes also included those involved in glycolysis and carbon metabolite repression, certain amino acid decarboxylases and some transporter genes including drug efflux genes. these results suggest that salmonella acquires a significant amount of essential nutrients directly from the host cell, thus shutting down several biosynthetic pathways when inside the host. the genes that cannot be detected include some genes involved in glycolysis/gluconeogenesis and anaplerotic reactions. other genes are detectable, but do not change in expression level over the course of infection, including those in the pentose phosphate pathway, entner duoduroff pathway and glutathione cycle. fatty acid biosynthesis and degradation genes are also not detected during the monitored early stages of infection. metabolic capacity of s. typhimurium in the host cell in silico prediction of metabolic pathways utilized during infection allows us to identify redundant metabolism that the pathogen could exploit in order to survive and replicate. during conditions of infection it is reasonable to assume that salmonella cells are under an energy limitation. hence, atp generated is fully coupled to biosynthesis and growth. given the possible nutrients provided in the host-cell environment, we used flux variability analysis to identify which reactions can be utilized and found that reactions span a large range of flux values when biomass production is optimal. these reactions encompass all the possible reactions that could be used by salmonella during infection. when the biomass production is allowed to be suboptimal, an additional reactions could be used by salmonella in the host cell environment. the constraining conditions of the simulation, were derived by literature mining to identify a media composition that could closely mimic nutrient availability for s. typhimurium in the host cell. we also used fva to identify blocked reactions and thus pathways that will never be used by the pathogen under the described conditions in vivo during infection. these reactions are typically those catalyzed by proteins that cannot be detected in the proteome or those genes that are down regulated under similar conditions. our model simulation results of these blocked reactions along with our gene expression data identify correctly fatty acid degradation metabolic genes to be non-operational during the above conditions. thus, there is some correlation between blocked reactions found by fva and unexpressed genes during infection. we compiled the list of all proteins that comprise the in silico reactome and compared that to previously determined proteome of s. typhimurium extracted from macrophage cells after infection. of the proteins identified with a high level of confidence in s. typhimurium, have an annotated metabolic function and in our model irr eighty of these proteins were associated with reactions present in the reactome calculated by fva analysis during optimal growth, while only were identified in the host-cell reactome when growth was sub-optimal. in another experiment, where the proteome of s. typhimurium was probed under different in vitro environmental conditions, a total of proteins were identified, of which were accounted for in our metabolic reconstruction. out of these proteins, and were associated with the optimal and suboptimal host cell reactome, respectively, as calculated using fva. most of the proteins that were found in both proteomic data sets and associated with blocked reactions, as calculated by fva, were trna synthetases. incorporating charged trnas instead of free amino acids in the biomass reaction would change the trna synthetase reactions from blocked to part of the optimal reactome and would result in higher agreement between the optimal reactome and identified proteins. using a similar strategy we compared micro-array gene expression data obtained from s. typhimurium extracted from macrophage cells. of the genes that were differentially expressed, genes had a metabolic function, thus comprising the metabolic transcriptome. when we compared those with the gene set represented by our fva analysis reactome, we identified 80% of the genes in the metabolic transcriptome were associated with reactions in the reactome. of the remaining that were not present in our in silico reactome, 65% of the genes were down-regulated and formed a subset of the blocked reactions from our fva analysis. although the accuracy of prediction is very good, it could be increased by improving our knowledge of existing nutrients in the macrophage or incorporating regulation into the model. by identifying preferred metabolic pathways we can better understand how the pathogen responds to its microenvironment during invasion and infection. discussion this study describes a genome-scale metabolic model, irr of the gram-negative intracellular pathogen s. typhimurium. the model has been reconstructed by integration of existing genomic, proteomic, and phenotypic experimental data. comparison between model predictions and experimental data demonstrated an accuracy of more than 80% for growth and virulence phenotypes. the model allowed us to use flux balance analysis to analyze growth and virulence genes given the environment the intracellular pathogen encounters. further, gene expression experiments and computational analysis allowed us to identify metabolic enzymes and pathways that are potentially active during infection. model-development is an iterative process that moves from a genome annotation and biochemistry of the organism to the development of a genome-scale metabolic network that can be used as a fully functional in silico organism to compute physiological phenotypes. our model irr serves also as a framework to integrate further gene expression and proteomic data and explain the expression level of specific gene products during infection. computing cellular functions from reconstructed networks allowed us to understand the relationship between environmental factors and genetic factors and their integrated function to orchestrate cell phenotypes. knowledge of the functional status of metabolic genes under defined constraints allows us to interpret the role of critical nodes in a cellular metabolic network. the function of redundant pathways under internally perturbed conditions, such as in the host cell during infection, provides important insights into metabolic network structure. although this internal robustness has been experimentally demonstrated, building a complete catalog of genes essential for growth under a myriad of conditions for any organism is time-consuming and laborious. a computational approach that predicts gene essentiality with high accuracy is thus of great value. model predictions of virulence for different mutants improved from 80% to 91% upon altering host-cell environment components and this shows how the model can be used as a tool for providing insights into the role of metabolism in infection and pathogenesis. to date, constraint-based models cannot predict effects of toxic intermediates and methods are needed that allow for prediction of concentrations within constraint-based models. a biomass reaction is needed to accurately predict the effects of gene deletions on growth and viability. although predictions have been made in the past with approximations of the biomass function to that of e. coli, there are some major differences in fatty acid and amino acid composition that require experimental measurements to determine the drain of different metabolites and the cost of their production on the function of the organism. our study has addressed this issue by directly determining the components of salmonella biomass, thus increasing the validity and predictive accuracy of model irr s. typhimurium is a pathogen that is the leading cause of human gastroenteritis and is used as a mouse model for human typhoid fever. although more than serovars exist, the lt strain, whose genome sequence has been used here for building model irr is one of the principal strains for cellular and molecular biology of salmonella. attenuated mutants have been used not only as live oral vaccines against salmonella infection but also to express antigens from other pathogens and to deliver proteins to solid tumors. one reason for attenuation or avirulence is the deletion or functional silencing of a genetic element in the bacterium resulting in the impairment of the intracellular survival capacity of s. typhimurium. our in silico analysis shows the ability to predict such effects with good accuracy. at this point the host cell metabolic environment has been inferred indirectly from gene-expression data. improved accuracy in directly determining what metabolites are available in the host cell environment and incorporation of transcriptional regulatory rules into irr will allow a more precise analysis of metabolic selective pressures a pathogen encounters during infection. in addition to pathogenicity and virulence analysis, the s. typhimurium metabolic reconstruction can be a valuable tool to interpret several'omics' data types. although there is not a direct correlation between gene or protein expression levels and metabolic fluxes, the genome-scale model serves as a platform to visualize these large data sets and interpret them in the context of operational, redundant or incomplete metabolic pathways. this provides a better understanding of cellular phenotypes and offers insights into how system level properties affect survival or sensitivity of pathogens to factors in the environment provided by the host cell. for example, the model prediction of fad genes not being operational during early stages of infection and operational during sub-optimal growth can be explained by the fact that fatty acids are not utilized in the initial stages of infection by s. typhimurium. this observation along with the calculated host cell reactome from irr could suggest a switch from optimal to sub-optimal growth or a change in nutrients available for salmonella inside host cell tissue. this switch is concomitant with a marked decrease in growth during persistence or chronic infection. growth rates are presumably higher during the initial stage of infection. the entire reactome of salmonella as predicted by the model provides many hypotheses that need to be further investigated to identify the implications of metabolic pathways and functions in pathogenesis. metabolomic data have become more readily available and easier to measure. they can be used to apply additional constraints to the model, thus improving model predictions of cell function and global phenotypes. the genome-scale model of s. typhimurium we have developed here also presents us with the opportunity to interrogate several other serovars or species that belong to the salmonella genus. for example, pand is the only gene that is non-functional in salmonella typhi, and is predicted to be required for virulence in specific conditions. the presence of pseudogenes in other serovars can be analyzed in the context of the genome-scale model to present us with systemic similarities and differences that may or may not play a role in the pathogenic potential of all the four species. serovars lacking pand would require panthothenate from its host. thus, compounds that affect pantothenate transport could potentially be used as antibiotics. such analyses across different serovars require genome-scale models and might accelerate the identification of rational drug targets. although many of the reactions in the irr model for s. typhimurium are conserved in models for e. coli mg there are unique metabolic and genetic differences that affect growth phenotypes in response to different environmental or genetic perturbations. the two organisms are each capable of utilizing unique carbon sources, which are related to differences in metabolic networks or differences in regulation or kinetics that affect flow through shared reactions common to both organisms. we also found that there are differences in the number of isozymes or alternative reactions/pathways that enable one organism to be more robust against gene deletions in a glucose minimal media condition. we anticipate that comparisons of model-predicted essential genes in more environmental conditions will lead to the identification of more metabolic differences between the two organisms, as the analysis was only performed for aerobic growth in glucose minimal media. these metabolic differences might shed light onto the evolutionary differences between the two organisms and the nutrients they are or were exposed to. the genome-scale metabolic network reconstruction presented here supports the notion of a robust minimal set of pathways that salmonella requires for intracellular survival during infection. our model allows in silico simulation of different nutritional requirements and constraints as they may be present in different host cell types or tissues that are difficult to directly address experimentally. our model can generate hypotheses that can be validated by using mutant strains, high-throughput data, and selective analysis of pathways. such an endeavor has the potential to lead to an improved understanding of how intracellular pathogens utilize the existing microenvironment of the host. this in turn might lead to the identification of conserved metabolic pathways that can be targeted by novel interfering compounds.  we have presented here extensive datasets within a reconstruction of a genome-scale metabolic model for s. typhimurium. this model has been validated by experimental data and by legacy data. it demonstrates an accuracy of 80% for growth and virulence phenotypes. thus the model irr provides a suitable platform for computing cellular phenotypes for salmonella and for further integration of high-throughput data. biomass composition of salmonella was determined experimentally, which increased the accuracy of our model building and also showed distinctive difference to that of e. coli. constraint-based analysis using the metabolic reconstructions irr for s. typhimurium and mg for e. coli identified unique metabolic differences between these two bacteria. expression profiling of selective key metabolic enzymes in salmonella during infection of macrophages identified operational and non-operational metabolic pathways. flux variabity analysis of metabolic pathways identified a minimal set of reactions required for salmonella intracellular survival. this study thus provides another high-resolution reconstruction of microbial metabolism of an intracellular pathogen next to that of m. tuberculosis and demonstrates its suitability for analyzing bacterial metabolic networks during infection. from further validation and testing of hypotheses generated from these reconstruction efforts we expect to arrive at a conserved set of metabolic reactions that interdigitate with those of the host cell and might thus provide novel targets for therapeutic intervention. 
gene expression profiling of s. typhimurium isolated from macrophage cell lines was used to constrain the model to predict metabolic pathways that are likely to be operational during infection. we have developed and validated a genome-scale metabolic network of salmonella typhimurium lt. this model accounts for genes that encode proteins catalyzing unique metabolic and transport reactions in the bacterium. this model also serves as platform for the integration of high-throughput data. replication of salmonella typhimurium inside its host cell is a model system for studying the pathogenesis of intracellular bacterial infections. we employed flux balance analysis and in silico gene essentiality analysis to investigate growth under a wide range of conditions that mimic in vitro and host cell environments. our analysis suggests that there is a robust minimal set of metabolic pathways that is required for successful replication of salmonella inside the host cell. its computational power allows identification of networked metabolic pathways and generation of hypotheses about metabolism during infection, which might be used for the rational design of novel antibiotics or vaccine strains. genome-scale modeling of bacterial metabolic networks provides a powerful tool to identify and analyze pathways required for successful intracellular replication during host-pathogen interaction.
 infections with salmonella cause significant morbidity and mortality worldwide. replication of salmonella typhimurium inside its host cell is a model system for studying the pathogenesis of intracellular bacterial infections. genome-scale modeling of bacterial metabolic networks provides a powerful tool to identify and analyze pathways required for successful intracellular replication during host-pathogen interaction.  we have developed and validated a genome-scale metabolic network of salmonella typhimurium lt. this model accounts for genes that encode proteins catalyzing unique metabolic and transport reactions in the bacterium. we employed flux balance analysis and in silico gene essentiality analysis to investigate growth under a wide range of conditions that mimic in vitro and host cell environments. gene expression profiling of s. typhimurium isolated from macrophage cell lines was used to constrain the model to predict metabolic pathways that are likely to be operational during infection.  our analysis suggests that there is a robust minimal set of metabolic pathways that is required for successful replication of salmonella inside the host cell. this model also serves as platform for the integration of high-throughput data. its computational power allows identification of networked metabolic pathways and generation of hypotheses about metabolism during infection, which might be used for the rational design of novel antibiotics or vaccine strains. 