Reconstruction and analysis of a genome-scale metabolic model for Scheffersomyces stipitis
© Balagurunathan et al; BioMed Central Ltd. 2012
Received: 21 December 2011
Accepted: 23 February 2012
Published: 23 February 2012
Fermentation of xylose, the major component in hemicellulose, is essential for economic conversion of lignocellulosic biomass to fuels and chemicals. The yeast Scheffersomyces stipitis (formerly known as Pichia stipitis) has the highest known native capacity for xylose fermentation and possesses several genes for lignocellulose bioconversion in its genome. Understanding the metabolism of this yeast at a global scale, by reconstructing the genome scale metabolic model, is essential for manipulating its metabolic capabilities and for successful transfer of its capabilities to other industrial microbes.
We present a genome-scale metabolic model for Scheffersomyces stipitis, a native xylose utilizing yeast. The model was reconstructed based on genome sequence annotation, detailed experimental investigation and known yeast physiology. Macromolecular composition of Scheffersomyces stipitis biomass was estimated experimentally and its ability to grow on different carbon, nitrogen, sulphur and phosphorus sources was determined by phenotype microarrays. The compartmentalized model, developed based on an iterative procedure, accounted for 814 genes, 1371 reactions, and 971 metabolites. In silico computed growth rates were compared with high-throughput phenotyping data and the model could predict the qualitative outcomes in 74% of substrates investigated. Model simulations were used to identify the biosynthetic requirements for anaerobic growth of Scheffersomyces stipitis on glucose and the results were validated with published literature. The bottlenecks in Scheffersomyces stipitis metabolic network for xylose uptake and nucleotide cofactor recycling were identified by in silico flux variability analysis. The scope of the model in enhancing the mechanistic understanding of microbial metabolism is demonstrated by identifying a mechanism for mitochondrial respiration and oxidative phosphorylation.
The genome-scale metabolic model developed for Scheffersomyces stipitis successfully predicted substrate utilization and anaerobic growth requirements. Useful insights were drawn on xylose metabolism, cofactor recycling and mechanism of mitochondrial respiration from model simulations. These insights can be applied for efficient xylose utilization and cofactor recycling in other industrial microorganisms. The developed model forms a basis for rational analysis and design of Scheffersomyces stipitis metabolic network for the production of fuels and chemicals from lignocellulosic biomass.
KeywordsGenome scale metabolic models Scheffersomyces stipitis Metabolic flux analysis Xylose utilization Anaerobic growth
Scheffersomyces stipitis (S. stipitis), formerly known as Pichia stipitis, is a hemiascomycetous yeast, closely related to several yeast endosymbionts of passalid beetles that inhabit and decay white-rotted hardwood [2, 3]. It has the highest native capacity for xylose fermentation of any known microbe [4, 5]. Fed batch cultures of S. stipitis produce around 47 g/l of ethanol with yields of 0.36 g/g xylose at 30°C . In addition to xylose, S. stipitis has the capability to ferment sugars from hydrolysates with yields equivalent to 80% of theoretical yield [6–8]. Auxotrophic strains have been created and methods for high efficiency transformation have been developed for S. stipitis[9, 10]. Genetic tools based on a loxP/Cre recombination system have been developed for functional genomics and metabolic engineering of this yeast . The availability of genetic tools and capability for fermentation of hydrolysates has made S. stipitis an attractive microorganism for bioconversion of lignocellulose to fuels and chemicals. S. stipitis has already been successfully engineered to produce lactic acid and xylitol [12, 13]. However, S. stipitis suffers from some drawbacks like lower fermentation rates, lower tolerance to ethanol and absence of anaerobic growth [5, 14, 15].
As a parallel approach, xylose utilization pathway from S. stipitis has been used to engineer xylose metabolism in Saccharomyces cerevisiae. Successive cycles of metabolic engineering have improved xylose utilization in recombinant S. cerevisiae[16, 17]. However, the ethanol productivity from xylose is still low. This has been attributed to: low substrate affinity of recombinant enzymes ; cofactor imbalance in the XR-XDH reactions [19, 20]; low xylose transport capacity [21, 22]; and failure to recognize xylose as a fermentable carbon source [23, 24]. The holistic analysis of metabolism in S. stipitis could provide useful insights to identify shortcomings in S. stipitis and S. cerevisiae metabolic networks.
The complete genome of S. stipitis has been sequenced . The functional annotation of the genome sequence showed numerous genes for lignocellulose bioconversion and systematic analysis of the genome sequence annotation is necessary to obtain useful insights. Genome scale metabolic models, which represent the link between the genotype and phenotype of the organism, can be reconstructed based on the genome sequence annotation and relevant biochemical and physiological information. These models have the ability to provide a holistic view of the metabolism of an organism. Once experimentally validated, these models can be used to characterize the metabolic resource allocation, generate experimentally testable predictions of cellular phenotypes, elucidate metabolic network evolution scenarios, design experiments that most effectively reveal the genotype-phenotype relationships, and design engineered microorganisms with desired properties like overproduction of commercially valuable chemicals [26–30]. Due to the genome wide-scale, these models enable systematic assessment of how perturbations in the metabolic network affect the organism as a whole which may not be possible by analyzing a set of enzymes or isolated pathways.
We have reported a framework for reconstruction of genome scale metabolic model of S. stipitis. In this study, a genome scale metabolic model has been developed for S. stipitis based on the proposed framework and a recently published protocol . Experimental procedure for the estimation of macromolecular composition of S. stipitis was standardized and used to obtain the biomass composition. Growth and non-growth associated maintenance energy requirements were also estimated from experimental data. The model was refined and validated based on the ability of S. stipitis to grow on different carbon, nitrogen, sulphur and phosphorus sources. In silico analysis of the model was used to identify biosynthetic requirements for anaerobic growth of S. stipitis in glucose and to analyze xylose utilization capability in S. stipitis. Model simulations were carried out to obtain insights on the recycling of nucleotide cofactors and mechanisms involved in mitochondrial respiration and oxidative phosphorylation.
Reconstruction of the genome scale metabolic model for Scheffersomyces stipitis
Macromolecular composition of S. stipitis Biomass
Biomass Composition of Scheffersomyces stipitis
Characteristics of the reconstructed network
Analysis of high-throughput substrate utilization
Biolog phenotyping results indicated that capric acid and caproic acid cannot by utilized by S. stipitis as a sole carbon source, but the model predicted growth. However, it has been observed that capric, caproic and other fatty acids were known to inhibit the growth of S. stipitis[35, 36] and other yeasts . Since the inhibition mechanisms are not incorporated in the metabolic model, the in silico computations predicted growth on these substrates. In the case of nitrogen source utilization, Biolog phenotyping results and experimental data reported by [38, 39] indicated growth on lysine as a sole nitrogen source. However, the in silico predictions did not predict growth as the pathway enzymes involved in metabolism of lysine has not been identified in S. stipitis genome. The in silico growth predictions for Biolog substrates evaluated as low confidence data were also compared with data available in literature for corresponding substrates. Experimental data was available from literature [38, 39] for 5 out of 40 low-confidence cases and the model could correctly predict the utilization of these substrates with 80% accuracy (4 cases)(Figure 4D). Incorrect prediction for glucuronic acid was due to the lack of homologs for pathway enzymes. The examples of lysine and glucuronic acid utilization illustrate the capability of the model to pinpoint potential gaps in the understanding of metabolism and to guide experimental design.
Metabolic requirements for anaerobic growth of Scheffersomyces stipitis
One drawback for using S. stipitis in industrial fermentation is its inability to grow under anaerobic conditions. The model developed was used to analyze the requirements for anaerobic growth. The model was simulated for anaerobic growth on a glucose based minimal media by reducing the oxygen uptake rate to zero, with an unconstrained uptake of sterols and unsaturated fatty acids (which have a known biosynthetic requirement for oxygen). No growth was predicted by the model simulations under these conditions. Growth under anaerobic conditions is a complex process. Several requirements need to be met for growth including biosynthetic requirements for oxygen, energy requirements, redox balance requirements and regulatory requirements.
List of reactions that lead to anaerobic growth on glucose identified by single reaction insertion analysis
h2o2[c] + h[c] + nadh[c] < = > 2 h2o[c] + nad[c]
nad[c] + 2 gthrd[c] < = > h[c] + nadh[c] + gthox[c]
h2o2[c] + h[c] + nadph[c] < = > 2 h2o[c] + nadp[c]
nadp[c] + 2 gthrd[c] < = > h[c] + nadph[c] + gthox[c]
o2[c] + pyr[c] + coa[c] < = > h2o2[c] + accoa[c] + co2[c]
o2[c] + lac-L[c] < = > h2o[c] + ac[c] + co2[c]
o2[c] + mal-L[c] < = > oaa[c] + h2o2[c]
o2[c] + glyclt[c] < = > glx[c] + h2o2[c]
asp-L[c] + o2[c] < = > h2o2[c] + iasp[c]
2 gthrd[c] < = > gthox[c]
h2o[c] + o2[c] + so3[c] < = > h2o2[c] + so4[c]
o2[c] + glyc3p[c] < = > h2o2[c] + dhap[c]
pyr[c] + pydam[c] < = > ala-L[c] + pydx[c]
oaa[c] + pydam[c] < = > asp-L[c] + pydx[c]
h2o[c] + o2[c] + hxan[c] < = > h2o2[c] + xan[c]
h2o[c] + cdpdag[c] < = > pa[c] + cmp[c]
pa[c] + ctp[c] < = > ppi[c] + cdpdag[c]
ser-L[c] + cdpdag[c] < = > cmp[c] + ps[c]
nadp[c] + dhor-S[c] < = > h[c] + nadph[c] + orot[c]
nad[c] + dhor-S[c] < = > h[c] + nadh[c] + orot[c]
akg[c] + o2[c] + duri[c] < = > co2[c] + succ[c] + uri[c]
atp[c] + pydxn[c] < = > adp[c] + pdx5p[c]
pi[c] + pydxn[c] < = > h2o[c] + pdx5p[c]
h2o[c] + o2[c] + xan[c] < = > h2o2[c] + urate[c]
amp[c] + gthox[c] + so3[c] < = > 2 gthrd[c] + aps[c]
chol[c] + cdpdag[c] < = > pc[c] + cmp[c]
o2[c] + h[c] + nadh[c] < = > h2o2[c] + nad[c]
o2[c] + h[c] + nadph[c] < = > h2o2[c] + nadp[c]
Model simulations were also carried out with xylose as the carbon source. The gene candidate identified for anaerobic growth on glucose resulted in a very low biomass flux. This flux was observed when xylose reductase activity was solely dependent on NADH and there was no flux when xylose reductase enzyme used NAPDH or a ratio of NADH and NADPH. S. stipitis strain with URA1 gene was not able to grow on xylose under anaerobic conditions . In addition to biosynthetic requirements, further analysis has to be performed on energy, redox balance and regulatory requirements to understand the limitations for anaerobic growth on xylose.
Xylose utilization by scheffersomyces stipitis
Reaction insertion analysis used for the identification of metabolic requirements for anaerobic growth was used to identify reactions that enable enhanced uptake of xylose under this condition. Several reactions were able to enhance xylose uptake rate and biomass flux [Data not shown]. A majority of these reactions were able to enhance the uptake rate of xylose by the effective production of the cofactors NADPH and NAD. This was evident when the reversible transdehydrogenase reaction which inter converts NADPH and NAD to NADH and NADP was introduced. This indicates that the metabolic network of S. stipitis lacks a sufficient NADPH forming transdehydrogenase reaction.
Flux variability analysis of the genome scale metabolic model
Flux variability analysis (FVA) was carried out for growth on glucose and xylose using the COBRA toolbox. Reactions known to result in loops within the major metabolic pathways were manually removed from the model before performing the flux variability analysis. The normal FVA calculates the minimum and maximum fluxes across various reactions in major metabolic pathways when maximizing the objective function (Biomass Flux). A variant of FVA called the sub-optimal FVA has been found to be more informative , wherein instead of fixing the objective value to an optimal value from the initial FBA, objective lower limit was chosen at 95% of the initial objective value. The normalized flux ranges (normalized with respect to substrate uptake rate) for the major reactions obtained using the suboptimal FVA is shown in supplementary Figures S1-S4 (Additional File 5). Experimental values for the metabolic flux distribution in S. stipitis are scarce. The only report on the metabolic flux profiling of S. stipitis compared the central carbon metabolism of this yeast with that of S. cerevisiae. This data cannot be directly compared with the model simulations as the substrate uptake rates and oxygen transfer rates were not reported.
NAD-dependent glutamate dehydrogenase and NADP-dependent glutamate dehydrogenase which causes a substrate shuffling between 2-oxoglutarate and glutamate which restores NADP from NADPH through the coupled conversion of NAD to NADH.
NADPH dehydrogenase which couples the oxidation of cytoplasmic NADPH to mitochondrial respiratory chain.
NAD-dependent alcohol dehydrogenase and NADP-dependent alcohol dehydrogenase which causes a substrate shuffling between ethanol and acetaldehyde which restores NADP from NADPH through the coupled conversion of NAD to NADH.
Comparison of the literature on the phosphoglucose isomerase mutant (pgi1-mutant) S. cerevisiae, Escherichia coli and Kluyveromyces lactis has indicated that the above mentioned pathways are either present in these organisms or when introduced has enhanced the growth of these organisms. K. lactis is reported to possess the mitochondrial NADPH dehydrogenase and a transdehydrogenase cycle involving the alcohol dehydrogenase [50, 51]. In E. coli and S. cerevisiae the introduction of a soluble transdehydrogenase gene was found to enhance the growth of pgi1-mutants on glucose [47, 52]. Further, in S. cerevisiae, over expression of the NAD-dependent glutamate dehydrogenase restored growth in these mutants . From the Biolog phenotype data presented earlier in this paper it can be observed that S. stipitis can grow on glutamate as the sole carbon source. The first step in glutamate utilization is NAD-dependent glutamate dehydrogenase and thus this pathway might be responsible for the growth of pgi1-mutant S. stipitis. However, all these pathways and the pgi1-mutant strains have to be evaluated experimentally to confirm their roles. Nevertheless, the analysis highlights the usefulness of the metabolic model developed in the designing microbial strains with desired properties.
Mechanism of mitochondrial respiration and oxidative phosphorylation
The role of the model in enhancing the understanding of cellular phenotypes and identifying requirements for metabolism were explained in the previous sections. In addition, the model can also be used to probe mechanisms in metabolism of a microorganism. In this section, the mechanisms involved in mitochondrial respiration and oxidative phosphorylation are analyzed using the model developed.
Effect of inhibition of various mitochondrial respiratory complexes on the growth of Scheffersomyces stipitis in glucose and xylose. (--) Complete Inhibition; (-) Partial Inhibition; (0) Negligible; (++) Enhanced; NA - Information not available
Effect on Growth
Effect on Growth Complex formation AOX and Complex III or IV (Predicted from model analysis)
( in Vivo /
( in Vivo /
( in Vivo /
( in Vivo /
Shi et al., 2002
Jeppsson et al., 1995
Jeppsson et al., 1995
Complex IV and
AOX (Sodium Azide)
Lighthelm et al., 1988
Complex IV(Cyanide) and AOX (SHAM)
Jeppsson et al., 1995
Complex I (Rotenone) + del AOX
Shi et al., 2002
Complex I (Rotenone) + del Complex IV
Shi et al., 2002
One of the well established ways by which microorganisms attain different phenotypic characteristics is by gene regulation. However, in the case of mitochondrial respiration, channeling the electron flow by the formation of super-complexes has been reported as a common mechanism [59–62]. Analysis of the inhibition data on glucose and xylose suggest that either the alternative oxidase or the alternate NADH dehydrogenase may be repressed when grown on glucose. However, studies on cytochrome-C mutants suggest that alternative oxidase may be expressed constitutively [56, 57]. Various hypothetical complexes were analyzed through model simulations and were then combined with regulation of gene expression (alternative NADH dehydrogenase) to result in a mechanism which explains all the experimental observations qualitatively. Based on the analysis, the complex formation between alternate oxidase and either Complex III or Complex IV is critical to explain the experimental observations. The proposed mechanism was also able to predict the observed increase in the growth yield in the presence of SHAM (an alternative oxidase inhibitor) during the growth of S. stipitis on xylose . This again highlights the utility of model in enhancing our understanding of metabolism of an organism.
An iterative procedure has been designed and used to develop the genome scale metabolic model of S. stipitis. The procedure begins with the reconstruction of genome scale metabolic network which is then converted to a fully functional in silico model by incorporating the experimentally determined macromolecular composition, maintenance coefficients and minimal medium requirements. Even though assumption of macromolecular composition based on other related microorganisms or compiling fragmented data from different sources is a common practice [46, 63–65], it may lead to improper predictions of essential genes and metabolic flux distribution. Determination of macromolecular composition is critical as it defines the minimal number of metabolites that has to be produced for growth and the relative contribution of these metabolites to growth. Further, it helps in refinement of the genome annotation as some pathways producing these metabolites might not be annotated properly. In the case of S. stipitis, higher content of chitin is observed as compared to other yeasts like S. cerevisiae (chitin is not a part of the biomass equation used for S. cerevisiae genome scale model ). Similarly, phosphatidyl-inositol was found in S. stipitis biomass but a pathway enzyme which converts myo-inositol to inositol was not annotated in S. stipitis genome. A homolog was identified for this enzyme (PICST_63214) and was incorporated into the model. The in silico model was validated for consistency against newly generated high-throughput substrate phenotyping data. The model was subjected to iterative refinements based on identified inconsistencies, leading to additional reactions incorporated into the network and other modifications to the model content. The net result is a biochemically and genetically detailed in silico model that consists of 1371 reactions that are catalyzed by 814 genes. To our knowledge this is the first genome scale metabolic model developed for S. stipitis and it qualitatively predicts the phenotypic behaviors for substrate utilization with 74% agreement (252 out of 339 cases).
The model developed has also generated several experimentally verifiable hypotheses that could provide insight into the metabolism of S. stipitis. Generation of new annotation for metabolic genes based on network-based gap analysis and high-throughput phenotyping is one such example. Another example is the insight obtained on the mechanism involved in mitochondrial respiratory chain of S. stipitis. The results of the inhibition experiments carried out for the various complexes in the mitochondrial respiration chain where explained based on a hypothetical complex formation using in silico metabolic flux analysis. The hypothetical complex formation has to be verified by designing appropriate experiments. Further, the formation of mitochondrial super complexes has already been reported in several yeasts like S. cerevisiae and Yarrowia lipolytica[61, 62].
The predictive capability of the genome scale metabolic network was demonstrated by identifying the biosynthetic requirements for anaerobic growth of S. stipitis. The gene insertion analysis performed on the metabolic model was able to identify the particular conversion essential for anaerobic growth on glucose and this was validated by comparing with literature. However, the identified reaction did not enable in silico growth on xylose and this was also observed in experiments . As mentioned in the results sections anaerobic growth is a complex phenotype with biosynthesis complemented by energy, redox balance and regulatory requirements. However, the ability to grow on glucose under anaerobic condition suggests that anaerobic growth on xylose can also be achieved if the other requirements are analyzed systematically. In silico analysis revealed that unconstrained supply of ATP could support anaerobic growth on xylose. Many pathways are known to generate ATP. For example, the bacterial acetate production pathway is commonly employed to improve ATP production . Since acetate has been reported to inhibit xylose fermentation in S. stipitis, other pathways for ATP generation needs to be investigated. Further redox balance requirements and regulatory requirements should also be considered when developing strategies to promote anaerobic growth.
Xylose uptake analysis carried out using the model indicated that S. stipitis metabolic network may not possess sufficient NADPH and/or NAD generation capability under lower oxygen uptake rates. Even though S. stipitis is considered to be an efficient xylose utilizing yeast, the metabolic model developed emphasizes the need for further characterization of this yeast to understand and improve its fermentation performance. High throughput techniques routinely used in the characterization of industrial yeasts like S. cerevisiae should also be developed for S. stipitis. Furthermore, analysis of pgi1-mutant shows that S. stipitis metabolic network possess many pathways to consume NADPH and recycle NADP. When grown on glucose the xylose utilization pathway itself can be considered as one such pathway. Thus by performing a detailed analysis on the various NADP recycling pathways and carefully designing suitable mutants the S. stipitis metabolic network can be optimized for simultaneous consumption of glucose and xylose. This trait has been a subject of active investigation and is of industrial significance .
In this study, we have reconstructed a genome scale metabolic model for S. stipitis by combining information from genome sequence annotation, pathway databases, literature and experimental data. The model was refined using high-throughput phenotyping data. The model predictions were in good agreement with experimental observations, thus allowing us to systematically investigate the physiological characteristics and metabolic capability of this yeast. In silico model analysis shows that S. stipitis possesses several pathways to recycle nucleotide cofactors and thus efficient xylose utilization. However, the flux through of these pathways needs experimental investigation. Analysis of mitochondrial respiration and identification of mitochondrial super-complexes demonstrate the novel applications of the model developed. Incorporation of thermodynamic constraints, enzyme kinetics information and high-throughput omics data can further improve the predictability of these models
Materials and methods
Metabolic network reconstruction
Figure 1 outlines the overall procedure for reconstruction and validation of the genome scale metabolic model for S. stipitis. The entire procedure is based on the recently published protocol for generation of genome scale metabolic models . The reconstruction process was initiated based on annotated genome of Scheffersomyces stipitis and ORF (open reading frame) information available on the National Centre for Biological Information http://www.ncbi.nlm.nih.gov/. The biochemical reactions corresponding to the ORFs were mainly compiled from KEGG Database [40, 69, 70] and MetaCyc Database , resulting in a draft reconstruction comprising of enzymes and metabolic reactions.
The next step is the refinement of the draft. For each reaction, gene-protein-reaction association, localization, cofactor specificity, and directionality were identified and assigned. The gene-protein-reaction associations include definitions for isoenzymes and enzyme complexes. The metabolic reactions in the model were organized into three compartments (Cytoplasm, Mitochondria and Extracellular) based on the localization of associated enzymes, which was obtained using the protein localization predictors [72, 73] (Lu et al., 2004; Claros and Vincens, 1996). The cofactor and substrate specificity and reaction directionality information were compiled from available literature information, completed genome scale reconstructions [33, 63, 74–76] and BRENDA database [77, 78]. The reactions were then organized into pathways/subsystems. For each metabolite, the charge, formula and identification information were compiled from KEGG database and BIGG database. In the next step, some spontaneous and non-gene associated reactions whose existence was supported by physiological or experimental data from the literature and databases were included. In addition to these reactions, exchange reactions and intracellular and extracellular transport reactions were added. The genes associated with the transport reactions were identified using the transport protein predictor . The final step is the incorporation of biomass equation and, growth and non-growth associated maintenance coefficient.
The final reconstruction was then loaded into Matlab and the in silico model obtained was evaluated for its capability to produce biomass precursors and known by-products from minimal media and gaps in the metabolic network were identified. The gaps were then filled based on pathway databases and published reconstructions. The dead-end reactions in the network were identified and appropriate exchange reactions were added when applicable. The model was then refined and validated iteratively using Biolog phenotyping data. The final model was used for the simulation studies.
Biomass macromolecular composition estimation
Media and cultivation conditions
Scheffersomyces stipitis CBS 6054 (CBS6054 = ATCC 58785 = NRRL Y-11545 = IFO 10063) was purchased from American Type Culture Collection. It was routinely cultured in yeast extract, peptone, dextrose (YPD) medium and samples were collected at the exponential growth phase for the estimation of biomass macromolecular composition. Experimental data for biomass composition estimation were collected in duplicates and average values obtained are reported. For growth rate and substrate uptake rate measurements, experiments were carried out in minimal media containing 0.17% yeast nitrogen base without amino acids, 0.5% ammonium sulphate and 2% glucose. Experiments were carried out in 1 liter Erlenmeyer flasks at 30°C and 200 rpm. The OD600 was measured periodically for several hours to calculate the growth rates and the supernatant was collected and analyzed using HPLC for calculating the substrate uptake rates. The macromolecular composition of S. stipitis biomass was estimated by measuring carbohydrate, protein, lipid and nucleic acid content and the components making up these macromolecules.
Biomass carbohydrate content
The cell wall polysaccharides (Glucan, Mannan and Chitin) and intracellular carbohydrates (Glycogen and Trehalose) together contribute to the total carbohydrate content. Cell wall isolation and quantification of polysaccharides in the cell wall was carried out as described by Francois, 2007 . Cells were harvested and then lyzed using glass beads in a homogenizer. The cell wall was separated by high speed centrifugation and dried. The supernatant was used for nucleic acid analysis. The dried cell wall was hydrolyzed with 72% sulphuric acid and the sugars liberated were analyzed using HPLC. Glucan and Mannan content was estimated from these liberated sugars. Laminarin and Mannan were used for calibration and galactose was used as the internal standard. For the analysis of chitin content, the cell wall was heated with 6% KOH to liberate chitin. The chitin was treated with chitinase to liberate glucosamine which was detected using Reissig's reagent. Glycogen content was determined as explained by Smolders et al., 1994 . Approximately 20 mg of lyophilised cells was resuspended in 10 ml of 0.6 M HCl and boiled on a heating block at 100°C for 1 h. Glucose liberated from glycogen hydrolysis was quantified using HPLC. For the measurement of trehalose content, cells were washed twice with cold water and resuspended in 3 ml of water for 15 min at 100°C. Sample was incubated overnight with trehalase in 60 ul of acetate buffer. Glucose liberated from trehalose hydrolysis was quantified using HPLC.
HPLC analysis of sugars
The sugars glucose and mannose was measured by HPLC using the Biorad Aminex HPX-87H column. 5 mM Sulphuric acid was used as mobile phase and sugars were detected using RID detector operating at 50°C.
Biomass protein and amino acid content
Total protein content of S. stipitis biomass is determined using nitric acid method . The harvested cells were washed in TE buffer and centrifuged. The cells were then solubilized in 70% nitric acid and incubated at 22°C for 24 h. The absorbance was measured at 358 nm to obtain the protein content. Calibration curves for protein estimation were obtained using BSA. The amino acid content of S. stipitis biomass was measured using standard protocols as explained in AOAC Official Method 994.12 and AOAC Official Method 985.28 .
Biomass nucleic acid content
The DNA and RNA content was determined using Orcinol reagent (0.1% Orcinol, 0.1% FeCl3.6H2O in concentrated HCl) . Freeze-dried cell lysate was resuspended and diluted in autoclaved MilliQ water for DNA and RNA analysis. Equal volumes of diluted sample and freshly prepared Orcinol reagent were mixed and incubated at 100°C for either 2 min (for DNA estimation) or 15 min (Total nucleic acid estimation). The mixtures were immediately cooled on ice. The mixture incubated for 2 min was further incubated at 37°C for 2 h after cooling on ice. The absorbance was measured at 600 nm to obtain the DNA and total nucleic acid content. The values obtained from DNA measurements were subtracted from the total nucleic acid values to obtain the RNA content of the cell lysate. Calibration curves were obtained using standard DNA and RNA.
Biomass lipid content
The total lipid content was determined as explained by Matyash et al., 2008 . Harvested cells were washed twice with 10 ml of 0.1% ammonium acetate solution and resuspended in 6 ml of ammonium acetate solution. Optical density was measured and 5 ml of cell suspension was added to 7.5 ml of methanol and vortexed vigorously. 25 ml of methyl-tert-butyl ether was added and the mixture was shaken at 250 rpm for 1 h. Cell debris was removed by filtration and 6.25 ml of water was added to the mixture. The organic phase was extracted twice with the solvent mixture and dried in a pre-weighed round bottom flask. Lipid content was then calculated based on the weight of the lipids extracted. The individual lipid classes were measured using TLC. Polar lipids were separated on a silica gel TLC plate with chloroform/methanol/water (65:25:4) as an eluent as described by Skipski et al., 1962 . Neutral lipids were separated on a silica gel TLC plate with hexane/Diethyl ether/formic acid (45:5:1) as an eluent as described by Low et al., 2009 . The separated lipids were quantified by densitometry using standard lipids. The fatty acid composition of S. stipitis was obtained using standard protocols as explained in AOAC Official Method 996.06 .
The composition of the various biomass macromolecules estimated were converted to a biomass synthesis equation to be incorporated into to model developed. The growth and substrates uptake data estimated in the present study, along with reported data from continuous cultivation experiments  were used to calculated the growth and non-growth associated maintenance coefficients (Refer to Additional File 1 for details).
In silico computations
The metabolic network was loaded into Matlab using functions available in the COBRA toolbox [89, 90]. The metabolic capabilities of S. stipitis network were calculated by using flux balance analysis and linear optimization. For growth simulations, biomass synthesis was selected as the objective to be maximized and the optimization was solved using the COBRA Toolbox. The external metabolites were allowed to freely cross the system boundary by having unconstrained exchange reaction. Units of all flux values are in mmol/gDCW/h. For the simulation of aerobic growth on minimal media, the following external metabolites were allowed to freely enter and leave the network: NH4, O2, H+, SO4, PO4, CO2 and H2O (exchange fluxs -1000-1000 mmol/gDCW/h). All other external metabolites, except the substrates tested were only allowed to leave the system (exchange fluxes 0-1000 mmol/gDCW/h). Growth on different substrates was simulated by allowing the corresponding external metabolite to enter the system with the definite exchange rate (10 mmol/gDCW/h). For the simulation of anaerobic growth, oxygen uptake was reduced to zero and unconstrained uptake of sterols and unsaturated fatty acids was allowed.
Phenotype microarray analysis of scheffersomyces stipitis
Biolog's Phenotype MicroArray™ technology  was used for the phenotypic analysis of Scheffersomyces stipitis. It permits assays of 190 carbon (PM1- and PM2A- microplates), 95 nitrogen (PM3BMicroplate), 59 phosphorus and 35 sulphur-source (PM4A-Microplate) utilizations at once. A defined medium containing 100 mM glucose, 5.0 mM NH4Cl, 2.0 mM NaH2PO4, 0.25 mM Na2SO4, 100 mM NaCl, 30 mM triethanolamine HCl (pH 7.1), 0.05 mM MgCl2, 1.0 mM KCl, 1.0 mM FeCl3, and 0.01% tetrazolium violet was used for the PM tests. The PM plates contained various carbon-, nitrogen-, phosphorus-, or sulphur-sources which are omitted from the defined medium. The microplates were incubated at 30°C and the dye reduction data were collected in 15-min intervals for 48 h. In addition to time profile dye reduction data, the absorbance was measured at 590 nm and 750 nm after 24 h and 48 h incubations. The colorimetric assay was considered positive when the absorbance corresponding to reduced dye was 1.2 times higher than the negative control. Similar threshold was applied for the absorbance at 750 nm. The threshold was set at 1.3 for sulphur source data as the negative control was observed to have a higher background. The confidence with which growth can be predicted based on these measurements was determined based on the number of positive reactions out of 4 absorbance measurements. High confidence growth (+) if all measurements were above threshold and high confidence no-growth (-) if below. The remaining cases were classified as low-confidence data.
We thank Dr. Dong-Yup Lee (Department of Chemical and Biomolecular Engineering, National University of Singapore) for valuable discussions on the reconstruction procedure. We also thank Dr Zhao Hua (Industrial Biotechnology Group, ICES, A*STAR) for suggestions on biomass composition estimation and Dr. Cletus P. Kurtzman (Microbial Genomics and Bioprocessing Research Unit, National Center for Agricultural Utilization Research, ARS, USDA) for sending Pichia stipitis strain to Biolog for phenotype microarray studies. We thank Dr. Michael Ziman (Biolog Inc.) for technical support during the design of phenotype microarray experiments. Financial support from ICES is gratefully acknowledged (project numbers ICES08-132002 and ICES/10-234A01).
- Kurtzman CP, Suzuki M: Phylogenetic analysis of the ascomycete yeasts that form coezyme Q-9 and the proposal of the new genera Babjeviella, Meyerozyma, Millerozyma, Priceomyces and Scheffersomyces. Mycoscience. 2010, 51: 2-14. 10.1007/s10267-009-0011-5.View Article
- Suh SO, Marshall CJ, McHugh JV, Blackwell M: Wood ingestion by passalid beetles in the presence of xylose-fermenting gut yeasts. Mol Ecol. 2003, 12: 3137-3145. 10.1046/j.1365-294X.2003.01973.x.View Article
- Nardi JB, Bee CM, Miller LA, Nyugen NH, Suh SO: Blackwell M Communities of microbes that inhabit the changing hindgut landscape of a subsocial beetle. Arthropod Struct Dev. 2006, 35: 57-68. 10.1016/j.asd.2005.06.003.View Article
- van Dijken JP, van den Bosch E, Hermans JJ, de Miranda LR, Scheffers WA: Alcoholic fermentation by 'non-fermentative' yeasts. Yeast. 1986, 2: 123-127. 10.1002/yea.320020208.View Article
- du Preez JC, van Dreissel B, Prior BA: Ethanol tolerance of Pichia stipitis and Candida shehatae strains in fed-batch cultures at controlled low dissolved oxygen levels. Appl Microbiol Biotechnol. 1989, 30: 53-58.View Article
- Ferrari MD, Neirotti E, Albornoz C, Saucedo E: Ethanol-production from eucalyptus wood hemicellulose hydrolysate by Pichia stipitis. Biotechnol Bioeng. 1992, 40: 753-759. 10.1002/bit.260400702.View Article
- Nigam JN: Development of xylose-fermenting yeast Pichia stipitis for ethanol production through adaptation on hardwood hemicellulose acid prehydrolysate. J Appl Microbiol. 2001, 90: 208-215. 10.1046/j.1365-2672.2001.01234.x.View Article
- Nigam JN: Ethanol production from hardwood spent sulphite liquor using an adapted strain of Pichia stipitis. J Ind Microbiol Biotechnol. 2001, 26: 145-150. 10.1038/sj.jim.7000098.View Article
- Yang VW, Marks JA, Davies BP, Jeffries TW: High-Efficiency transformation of Pichia stipitis based on its URA3 gene and a homologous autonomous replication sequence, ARS2. Appl Environ Microbiol. 1994, 60 (12): 4245-4254.
- Lu P, Davis BP, Hendrick J, Jeffries TW: Cloning and disruption of the β- isopropylmalate dehydrogenase gene (LEU2) of Pichia stipitis and with URA3 and recovery of the double auxotroph. Appl Microbiol Biotechnol. 1998, 49: 141-146. 10.1007/s002530051150.View Article
- Laplaza JM, Torres BR, Jin YS, Jeffries TW: Sh ble and Cre adapted for functional genomics and metabolic engineering of Pichia stipitis. Enzyme Microb Technol. 2006, 38: 741-747. 10.1016/j.enzmictec.2005.07.024.View Article
- Jin YS, Cruz J, Jeffries TW: Xylitol production by a Pichia stipitis D-xylulokinase mutant. Appl Microbiol Biotechnol. 2005, 68: 42-45. 10.1007/s00253-004-1854-5.View Article
- Ilmen M, Koivuranta K, Ruohonen L, Suominen P, Penttila M: Efficient production of l-lactic acid from xylose by Pichia stipitis. Appl Environ Microbiol. 2007, 73: 117-123. 10.1128/AEM.01311-06.View Article
- Grootjen DRJ, van der Lans GJM, Luyben KChAM: Effects of the aeration rate on the fermentation of glucose and xylose by Pichia stipitis CBS 5773. Enzyme Microb Technol. 1990, 12: 20-23. 10.1016/0141-0229(90)90174-O.View Article
- Shi XQ, Jeffries TW: Anaerobic growth and improved fermentation of Pichia stipitis bearing a URA1 gene from Saccharomyces cerevisiae. Appl Microbiol Biotechnol. 1998, 50: 339-345. 10.1007/s002530051301.View Article
- Johansson B, Hahn-Hagerdal B: The non-oxidative pentose phosphate pathway controls the fermentation rate of xylulose but not of xylose in Saccharomyces cerevisiae TMB3001. FEMS Yeast Res. 2002, 2: 277-282.
- Karhumaa K, Hahn-Hagerdahl B, Gorwa-Grauslund MF: Investigation of the limiting metabolic steps in the utilization of xylose by recombinant Saccharomyces cerevisiae using metabolic engineering. Yeast. 2005, 22: 359-368. 10.1002/yea.1216.View Article
- Kotter P, Ciriacy M: Xylose fermentation by Saccharomyces cerevisiae. Appl Microbiol Biotechnol. 1993, 38: 776-783. 10.1007/BF00167144.View Article
- Bruinenberg PM, Debot PH, van Dijken JP, Scheffers WA: The role of redox balances in the anaerobic fermentation of xylose by yeasts. Eur J Appl Microbiol Biotechnol. 1983, 18: 287-292. 10.1007/BF00500493.View Article
- Eliasson A, Christensson C, Wahlbom CF, Hahn-Hägerdal B: Anaerobic xylose fermentation by recombinant Saccharomyces cerevisiae carrying XYL, XYL, and XKS in mineral medium chemostat cultures. Appl Environ Microbiol. 2000, 66: 3381-3386. 10.1128/AEM.66.8.3381-3386.2000.View Article
- Saloheimo A, Rauta J, Stasyk OV, Sibirny AA, Penttila M, Ruohonen L: Xylose transport studies with xylose-utilizing Saccharomyces cerevisiae strains expressing heterologous and homologous permeases. Appl Microbiol Biotechnol. 2007, 74: 1041-1052. 10.1007/s00253-006-0747-1.View Article
- Runquist D, Fonseca C, Radstrom P, Spencer-Martins I, Hahn-Hagerdal B: Expression of the Gxf1 transporter from Candida intermedia improves fermentation performance in recombinant xylose-utilizing Saccharomyces cerevisiae. Appl Microbiol Biotechnol. 2009, 82: 123-130. 10.1007/s00253-008-1773-y.View Article
- Jin YS, Laplaza JM, Jeffries TW: Saccharomyces cerevisiae engineered for xylose metabolism exhibits a respiratory response. Appl Environ Microbiol. 2004, 70: 6816-6825. 10.1128/AEM.70.11.6816-6825.2004.View Article
- Souto-Maior AM, Runquist D, Hahn-Hagerdal B: Crabtree-negative characteristics of recombinant xylose-utilizing Saccharomyces cerevisiae. J Biotechnol. 2009, 143: 119-123. 10.1016/j.jbiotec.2009.06.022.View Article
- Jeffries TW, Grigoriev I, Grimwood J, Laplaza JM, Aerts A, Salamov A, Schumtz J, Linquist E, Dehal P, Shapiro H, Jin YS, Passoth V, Richardson PM: Genome sequence of the lignocellulose-bioconverting and xylose-fermenting yeast Pichia stipitis. Nat Biotechnol. 2007, 25 (3): 319-326. 10.1038/nbt1290.View Article
- Patil KR, Akesson M, Nielsen J: Use of genome-scale microbial models for metabolic engineering. Curr Opin Biotechnol. 2004, 15: 64-69. 10.1016/j.copbio.2003.11.003.View Article
- Jamshidi N, Palsson BO: Formulating genome-scale kinetic models in the post-genome era. Mol Syst Biol. 2008, 4: 171-View Article
- Oberhardt MA, Palsson BO, Papin JA: Applications of genome-scale metabolic reconstructions. Mol Syst Biol. 2009, 5: 320-View Article
- Park JM, Kim TY, Lee SY: Constraints-based genome-scale metabolic simulation for systems metabolic engineering. Biotechnol Adv. 2009, 27 (6): 979-988. 10.1016/j.biotechadv.2009.05.019.View Article
- Liu L, Agren R, Bordel S, Nielsen J: Use of genome-scale metabolic models for understanding microbial physiology. FEBS Lett. 2010, 584 (12): 2556-2564. 10.1016/j.febslet.2010.04.052.View Article
- Balagurunathan B, Jonnalagadda S, Lee DY, Srinivasan R: Reconstruction of genome scale metabolic model for naturally xylose fermenting yeast Pichia stipitis [abstract]. J Biotechnol. 2008, 136S: S49-View Article
- Thiele I, Palsson BO: A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010, 5 (1): 93-121. 10.1038/nprot.2009.203.View Article
- Forster J, Famili I, Fu P, Palsson BO, Nielsen J: Genome scale reconstruction of the Saccharomyces cerevisiae metabolic network. Genome Res. 2003, 13: 244-253. 10.1101/gr.234503.View Article
- Bochner BR, Gadzinski P, Panomitros E: Phenotype microarrays for high-throughput phenotypic testing and assay of gene function. Genome Res. 2001, 11: 1246-1255. 10.1101/gr.186501.View Article
- Tran AV, Chambers RP: Red oak wood derived inhibitors in the ethanol fermentation of xylose by Pichia stipitis CBS 5776. Biotechnol Lett. 1985, 7 (11): 841-846. 10.1007/BF01025567.View Article
- Tran AV, Chambers RP: Ethanol fermentation of red oak acid prehydrolysate by the yeast Pichia stipitis CBS 5776. Enzyme Microb Technol. 1986, 8: 439-444. 10.1016/0141-0229(86)90154-7.View Article
- Hunkova Z, Fencl Z: Toxic effects of fatty acids on yeast cells: dependence of inhibitory effects on fatty acid concentration. Biotechnol Bioeng. 1977, 19: 1623-1641. 10.1002/bit.260191103.View Article
- Kurtzman CP, Fell JW: The Yeasts: A Taxonomic Study. 1998, Amsterdam: Elsevier, 4
- Barnett JA, Payne RW, Yarrow D: Yeasts: Characteristics and Identification. 2000, Cambridge: Cambridge University Press, 3
- Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.View Article
- Jeffries TW: Engineering yeasts for xylose metabolism. Curr Opin Biotechnol. 2006, 17: 320-326. 10.1016/j.copbio.2006.05.008.View Article
- Hahn-Hagerdal B, Karhumaa K, Jeppsson M, Gorwa-Grauslund MF: Metabolic engineering for pentose utilization in Saccharomyces cerevisiae. Adv Biochem Eng Biotechnol. 2007, 108: 147-177.
- Verduyn C, van Kleef R, Frank J, Schreuder H, van Dijken JP, Scheffers WA: Properties of the NAD(P)H-dependent xylose reductase from xylose-fermenting yeast Pichia stipitis. Biochem J. 1985, 226: 669-677.View Article
- Skoog K, Hahn-Hagerdahl B: Effect of oxygenation on xylose fermentation by Pichia stipitis. Appl Environ Microbiol. 1990, 56 (11): 3389-3394.
- Agbogbo FK, Coward-Kelly G, Torry-Smith N, Wenger K, Jeffries TW: The effect of initial cell concentration on xylose fermentation by Pichia stipitis. Appl Biochem Biotechnol. 2007, 136-140: 653-662.
- Puchalka J, Oberhardt JA, Godinho M, Bielecka A, Regenhardt D, Timmis KN, Papin JA, Martins dos Santos VAP: Genome-scale reconstruction and analysis of the Pseudomonas putida KT2440 metabolic network facilitates applications in biotechnology. Plos Comput Biol. 2008, 4 (10): e1000210-10.1371/journal.pcbi.1000210.View Article
- Fiaux J, Cakar ZP, Sonderegger M, Wuthrich K, Szyperski T, Sauer U: Metabolic-flux profiling of the yeasts Saccharomyces cerevisiae and Pichia stipitis. Eukaryot cell. 2003, 2 (1): 170-180. 10.1128/EC.2.1.170-180.2003.View Article
- Aguilera A: Deletion of the phosphoglucose isomerase structural gene makes growth and sporulation glucose dependent in Saccharomyces cerevisiae. Mol Gen Genet. 1986, 204: 310-316. 10.1007/BF00425515.View Article
- Boles E, Lehnert W, Zimmermann FK: The role of the NAD-dependent glutamate dehydrogenase in restoring growth on glucose of a Saccharomyces cerevisiae phosphoglucose isomerase mutant. Eur J Biochem. 1993, 217: 469-477. 10.1111/j.1432-1033.1993.tb18266.x.View Article
- Overkamp KM, Bakker BM, Steensma HY, van Dijken JP, Pronk JT: Two mechanisms for the oxidation of cytosolic NADPH by Kluyveromyces lactis mitochondria. Yeast. 2002, 19: 813-824. 10.1002/yea.878.View Article
- Tarrio N, Becerra M, Cerdan ME, Gonzalez Siso MI: Reoxidation of cytosolic NAPDH in Kluyveromyces lactis. FEMS Yeast Res. 2006, 6: 371-380. 10.1111/j.1567-1364.2005.00021.x.View Article
- Sauer W, Canonaco F, Heri S, Perrenoud A, Fischer E: The soluble and membrane bound transdehydrogenase UdhA and PntAB have divergent functions in NAPDH Metabolism of Escherichia coli. J Biol Chem. 2004, 279 (8): 6613-6619.View Article
- Jeffries TW: Engineering the Pichia stipitis genome for fermentation of hemicelluloses hydrolysates. Bioenergy. Edited by: Wall JD, Harwood CS, Demain A. 2008, Washington: ASM press, 37-47.View Article
- Joseph-Horne T, Hollomon DW, Wood PM: Fungal respiration: a fusion of standard and alternative components. Biochim Biophys Acta. 2001, 1504: 179-195. 10.1016/S0005-2728(00)00251-6.View Article
- Lighthelm ME, Prior BA, du Preez JC: The effect of respiratory inhibitors on the fermentative ability of Pichia stipitis, Pachysolen tannopilus and Saccharomyces cerevisiae under various conditions of aerobiosis. Appl Microbiol Biotechnol. 1988, 29: 67-71. 10.1007/BF00258353.View Article
- Jeppsson H, Alexander NJ, Hahn-Hagerdahl B: Existence of cyanide-insensitive respiration in the yeast Pichia stipitis and its possible influence on product formation during xylose utilization. Appl Environ Microbiol. 1995, 61 (7): 2596-2600.
- Shi NQ, Davis B, Sherman B, Cruz J, Jeffries TW: Disruption of the Cytochrome c Gene in xylose utilizing yeast Pichia stipitis leads to higher ethanol production. Yeast. 1999, 15: 1021-1030. 10.1002/(SICI)1097-0061(199908)15:11<1021::AID-YEA429>3.0.CO;2-V.View Article
- Shi NQ, Cruz J, Sherman B, Jeffries TW: SHAM-sensitive alternative respiration in the xylose-metabolizing yeast Pichia stipitis. Yeast. 2002, 19: 1203-1220. 10.1002/yea.915.View Article
- Wittig I, Carrozzo R, Santorelli FM, Schagger H: Supercomplexes and supcomplexes of mitochondrial oxidative phosphorylation. Biochim Biophys Acta. 2006, 1757: 1066-1072. 10.1016/j.bbabio.2006.05.006.View Article
- Boekema EJ, Braun HP: Supramolecular structure of the mitochondrial oxidative phosphorylation system. J Biol Chem. 2007, 282: 1-4.View Article
- Heinemeyer J, Braun HP, Boekema EJ, Kouil R: A structural model of the cytochrome c reductase/oxidase supercomplex from yeast mitochondria. J Biol Chem. 2007, 282: 12240-12248.View Article
- Guerrero-Castillo S, Vazquez-Acevedo M, Gonzalez-Halphen D, Uribe-Carvajal S: In Yarrowia lipolytica mitochondria, the alternative NADH dehydrogenase interacts specifically with the cytochrome complexes of the classic respiratory pathway. Biochim Biophys Acta. 2009, 1787: 75-85. 10.1016/j.bbabio.2008.10.008.View Article
- Oh YK, Palsson BO, Park SM, Schilling CH, Mahadevan R: Genome scale reconstruction of metabolic network in Bacillus subtilis based on high-throughput phenotyping and gene essentiality data. J Biol Chem. 2007, 282: 28791-28799. 10.1074/jbc.M703759200.View Article
- Suthers PF, Dasika MS, Kumar VS, Denisov G, Glass JI, Maranas CD: A Genome-scale metabolic reconstruction of Mycoplasma genitalium iPS189. Plos Comput Biol. 2009, 5 (2): e1000285-10.1371/journal.pcbi.1000285.View Article
- Chung BK, Selvarasu S, Andrea C, Ryu J, Lee H, Ahn J, Lee H, Lee DY: Genome-scale metabolic reconstruction and in silico analysis of methylotrophic yeast Pichia pastoris for strain improvement. Microb Cell Fact. 2010, 9: 50-10.1186/1475-2859-9-50.View Article
- Sohn SB, Kim TY, Park JM, Lee SY: In silico genome-scale metabolic analysis of Pseudomonas putida KT2440 for polyhydroxyalkanoate synthesis, degradation of aromatics and anaerobic survival. Biotechnol J. 2010, 5 (7): 739-750. 10.1002/biot.201000124.View Article
- Nigam JN: Ethanol production from wheat straw hemicellulose hydrolysate by Pichia stipitis. J Biotechnol. 2001, 87: 17-27. 10.1016/S0168-1656(00)00385-0.View Article
- Kim JH, Block DE, Mills DA: Simultaneous consumption of pentose and hexose sugars: an optimal microbial phenotype for efficient fermentation of lignocellulose biomass. Appl Microbiol Biotechnol. 2010, 88: 1077-1085. 10.1007/s00253-010-2839-1.View Article
- Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 2006, 34: D354-D357. 10.1093/nar/gkj102.View Article
- Aoki-Kinoshita KF, Kanehisa M: Gene annotation and pathway mapping in KEGG. Methods Mol Biol. 2007, 396: 71-91. 10.1007/978-1-59745-515-2_6.View Article
- Caspi R, Altman T, Dale JM, Dreher K, Fulcher CA, Gilham F, Kaipa P, Karthikeyan AS, Kothari A, Krummenacker M, Latendresse M, Mueller LA, Paley S, Popescu L, Pujar A, Shearer A, Zhang P, Karp PD: The MetaCyc Database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic Acids Res. 2010, 38 (1): D473-D479.View Article
- Lu Z, Szafron D, Greiner R, Lu P, Wishart DS, Poulin B, Anvik J, Macdonell C, Eisner R: Predicting subcellular localization of proteins using machine-learned classifiers. Bioinformatics. 2004, 20: 547-556. 10.1093/bioinformatics/btg447.View Article
- Claros MG, Vincens P: Computational method to predict mitochondrially imported proteins and their targeting sequences. Eur J Biochem. 1996, 241: 779-786. 10.1111/j.1432-1033.1996.00779.x.View Article
- Feist AM, Henry CS, Reed JL, Krummenacker M, Joyce AR, Karp PD, Broadbelt LJ, Hatzimanikatis V, Palsson BO: A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol Syst Biol. 2007, 3: 121-View Article
- Schellenberger J, Park JO, Conrad TC, Palsson BO: BiGG: a Biochemical Genetic and Genomic knowledgebase of large scale metabolic reconstructions. BMC Bioinformatics. 2010, 11: 213-10.1186/1471-2105-11-213.View Article
- Cvijovic M, Olivares-Hernández R, Agren R, Dahr N, Vongsangnak W, Nookaew I, Patil KR, Nielsen J: BioMet Toolbox: genome-wide analysis of metabolism. Nucleic Acids Res. 2010, 38: W144-W149. 10.1093/nar/gkq404.View Article
- Schomburg I, Hofmann O, Baensch C, Chang A, Schomburg D: Enzyme data and metabolic information: BRENDA, a resource for research in biology, biochemistry, and medicine. Gene Funct Dis. 2000, 3-4: 109-118.View Article
- Schomburg I, Chang A, Hofmann O, Ebeling C, Ehrentreich F, Schomburg D: BRENDA: a resource for enzyme data and metabolic information. Trends Biochem Sci. 2002, 27 (1): 54-56. 10.1016/S0968-0004(01)02027-8.View Article
- Li H, Benedito VA, Udvardi MK, Zhao PX: TransportTP: a two-phase classification approach for membrane transport prediction and characterization. BMC Bioinformatics. 2009, 10: 418-10.1186/1471-2105-10-418.View Article
- Francois JM: A simple method for quantitative determination of polysaccharides in fungal cell walls. Nat Protoc. 2007, 1 (6): 2995-3000. 10.1038/nprot.2006.457.View Article
- Smolders GJ, van der Meij J, van Loosdrecht MC, Heijnen JJ: Stoichiometric model of the aerobic metabolism of the biological phosphorus removal process. Biotechnol Bioeng. 1994, 44 (7): 837-848. 10.1002/bit.260440709.View Article
- Bible KC, Boerner SA, Kaufmann SH: A one-step method for protein estimation in biological samples: nitration of tyrosine in nitric acid. Anal Biochem. 1999, 267: 217-221. 10.1006/abio.1998.2980.View Article
- Horwitz W: Official methods of analysis of the AOAC International, 2nd rev. 2003, Gaithersburg (MD): AOAC International, 17
- Almog R, Shirey TL: A modified Orcinol test for the specific determination of RNA. Anal Biochem. 1978, 91 (1): 130-137. 10.1016/0003-2697(78)90823-0.View Article
- Matyash V, Leibisch G, Kurzchalia TV, Shevchenko A, Schwudke D: Lipid extraction by methyl-tert-butyl ether for high-throughput lipidomics. J Lipid Res. 2008, 49 (5): 1137-1146. 10.1194/jlr.D700041-JLR200.View Article
- Skipski VP, Peterson RF, Barclay M: Separation of phosphatidylethanolamine, phosphatidylserine and other phospholipids by thin-layer chromatography. J Lipid Res. 1962, 3: 467-470.
- Low KL, Rao PS, Shui G, Bendt AK, Pethe K, Dick T, Wenk MR: Triacylglycerol utilization is required for regrowth of in vitro hypoxic nonreplicating Mycobacterium bovis Bacillus calmette-Guerin. J Bacteriol. 2009, 191 (16): 5037-5043. 10.1128/JB.00530-09.View Article
- van Urk H, Voll WSL, Scheffers WA, van Dijken JP: Transient-state analysis of metabolic fluxes in crabtree-positive and crabtree-negative yeasts. Appl Microbiol Biotechnol. 1990, 56 (1): 281-287.
- Becker SA, Feist AM, Mo ML, Hannum G, Palsson BO, Herrgard MJ: Quantitative prediction of cellular metabolism with constraint-based models: The COBRA Toolbox. Nat Protoc. 2007, 2: 727-738. 10.1038/nprot.2007.99.View Article
- Schellenberger J, Que R, Fleming RMT, Thiele I, Orth JD, Feist AM, Zielinski DC, Bordbar A, Lewis NE, Rahmanian S, Kang J, Hyduke DR, Palsson BO: Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011, 6 (9): 1290-307. 10.1038/nprot.2011.308.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.