Metabolic modeling and analysis of the metabolic switch in Streptomyces coelicolor
© Alam et al; licensee BioMed Central Ltd. 2010
Received: 22 October 2009
Accepted: 26 March 2010
Published: 26 March 2010
The transition from exponential to stationary phase in Streptomyces coelicolor is accompanied by a major metabolic switch and results in a strong activation of secondary metabolism. Here we have explored the underlying reorganization of the metabolome by combining computational predictions based on constraint-based modeling and detailed transcriptomics time course observations.
We reconstructed the stoichiometric matrix of S. coelicolor, including the major antibiotic biosynthesis pathways, and performed flux balance analysis to predict flux changes that occur when the cell switches from biomass to antibiotic production. We defined the model input based on observed fermenter culture data and used a dynamically varying objective function to represent the metabolic switch. The predicted fluxes of many genes show highly significant correlation to the time series of the corresponding gene expression data. Individual mispredictions identify novel links between antibiotic production and primary metabolism.
Our results show the usefulness of constraint-based modeling for providing a detailed interpretation of time course gene expression data.
The transition from exponential growth to stationary phase is a major event in microbial physiology . During the exponential phase of growth, bacterial cells produce metabolites necessary for growth and grow rapidly. Once essential nutrients have been depleted, cells switch to stationary phase, stop growing, reorganize their energy metabolism and often start producing a new set of secondary metabolites, including antibiotics .
In this study, we have explored the metabolic switch in Streptomyces coelicolor, the model organism of the antibiotics producing genus Streptomyces. The genome of this soil bacterium has been sequenced and contains about 7825 genes, one of the largest numbers for any bacterium . More than 20 clusters coding for the 4 known and several predicted antibiotics or related compounds have been identified in the genome . To optimize the production of valuable secondary metabolites, understanding the shift from primary to secondary metabolism during the transition phase will play a key role.
We constructed a constraints-based genome-scale stoichiometric model of S. coelicolor metabolism, based on earlier similar models [5, 6], and integrated the model predictions with a large gene expression dataset . The constraints-based approach, in particular flux balance analysis, has been shown to be highly predictive of growth phenotypes in many microbial systems [8, 9] and can be used to construct large scale metabolic models based on genome sequences in the absence of kinetic information, making it particularly attractive for less well-studied organisms like S. coelicolor.
Predictions from constraint-based models usually hold for steady-state assumptions [10, 11]. To enable the incorporation of experimental information from time-series measurements, we extend the approach by applying a dynamically changing input function (specifying constrains on nutrient uptake) and objective function (specifying the shift of cellular resources from cellular growth to antibiotics production). The predicted flux profiles are then compared to the gene expression profiles of the corresponding enzyme-coding genes to validate the model.
We observe a surprisingly good correlation between predicted fluxes and measured gene expression, indicating both the power of the constraint-based modeling approach and the tight regulation of gene expression in S. coelicolor. A small number of incorrectly predicted fluxes indicate the need for including additional gene regulatory constraints to the model [12, 13], but also allows the sensitive identification of misannotations and putative novel reactions involved in secondary metabolite biosynthesis.
Results and Discussion
We have reconstructed a genome-scale model of Streptomyces coelicolor metabolism with recent updated annotations as discussed in the Methods section. Our aim was to study the metabolic switch between the primary phase and secondary phase of growth.
Initial model validation
Comparison of experimentally observed dilution rates from chemostat data  and predicted specific growth rates
Actinorhodin (μ g/g.h)
Observed dilution rate D
Predicted specific growth rate μ
Global metabolic switching from primary phase to secondary phase of growth
For a more detailed understanding of the metabolic transition phase, we then modeled flux changes happening during fermentation culture on phosphate limiting medium. For this growth condition we had earlier collected a detailed gene expression time series.
A large set of genes does not show correlation (64% of genes; -0.5 < r < 0.5). These are mostly genes that do not change expression (nor predicted flux) along the time course. In these cases of constant expression no correlation information is present in the data, leading to correlation coefficient close to zero. Of course, there will also be cases where gene expression levels and flux levels do not correlate for other reasons, for instance due to post-transcriptional and post-translational regulation mechanisms.
List of anticorrelated genes
fatty acid biosynthesis
fatty acid biosynthesis
asparagine synthetase (glutamine-hydrolysing)
pentadecanoyl-[acyl-carrier protein] synthesis
fatty acid biosynthesis
pentadecanoyl-[acyl-carrier protein] synthesis
fatty acid biosynthesis
pentadecanoyl-[acyl-carrier protein] synthesis
fatty acid biosynthesis
pentadecanoyl-[acyl-carrier protein] synthesis
fatty acid biosynthesis
GTP cyclohydrolase II
GTP cyclohydrolase II
One large group of anticorrelated genes is the set of 10 genes located the middle of the calcium dependent antibiotics (CDA) biosynthesis gene cluster (SCO3210-SCO3249) . SCO3210 and SCO3221 are annotated as 2-dehydro-3-deoxyheptonate aldolase and prephenate dehydrogenase respectively, part of the shikimate pathway (tryptophan biosynthesis). Tryptophan is a precursor for CDA, and there are four anticorrelated genes (SCO3211-3214) which encode for enzymes TrpC2, TrpD2, TrpG, and TrpE2. It seems obvious that these genes are involved in the biosynthesis of tryptophan for CDA biosynthesis and not in the production of tryptophan for general primary metabolism. Indeed it has been shown that these genes do not complement a deficiency in central tryptophan biosynthesis . SCO3249 encodes an ACP homolog, and the adjacent genes SCO3246 and SCO3248 along with SCO3228 are proposed to be involved in the biosynthesis of the N-terminal epoxyhexanoyl fatty acid side chain . While the direct involvement in CDA biosynthesis has not yet been established for all of these genes, the non-complementation as well as the clear anticorrelation in our model analysis point to the existence of strong regulatory constraints on the expression of these genes. Such regulatory constraints are not routinely included in flux balance analysis, but can substantially enhance its predictive accuracy [12, 13]. Our result shows that a lack of regulatory information can be efficiently compensated by the integration of transcriptomics information, which quite specifically highlights this group of genes for further study.
Another group of anticorrelating genes is seen in the middle of the undecylprodigiosin biosynthesis gene cluster [4, 18]. Three genes (SCO5886, SCO5887 and SCO5888) in this cluster were automatically annotated in our model as fatty acid biosynthesis genes on the basis of sequence similarity with fatty acid genes (3-oxoacyl-[acyl-carrier-protein] synthase II, acyl carrier protein, 3-oxoacyl-[acyl carrier protein] synthase III). However these three genes are well known to be involved in undecylprodigiosin production under the gene names redR (SCO5886), redQ (SCO5887) and redP (SCO5888). This is a clear example of a misannotation that is revealed by the correlation analysis and can easily be fixed in the model.
A third example of strongly anticorrelated genes listed in Table 2 are three alkaline phosphatases - SCO2286 (phoA), SCO0828 (phoC) and SCO2068 (phoD) - which are assigned in the KEGG database (and consequently in our model) to the folate biosynthesis pathway. Their expression pattern, which shows strong induction upon phosphate depletion, is consistent with earlier reports on their control by PhoR/PhoP  and a potential role in secondary metabolism, but is less easy to reconcile with a putative function in folate biosynthesis, which is based only on sequence homology.
In all three of these cases, the integration of gene expression and model flux predictions highlighted groups of genes involved in antibiotics production. A small set of additional anticorrelated genes (Table 2) are widely scattered through out the genome (Figure 5). Each of them is a potential candidate from model correction and for the identification of new secondary metabolite biosynthesis genes with specifically constrained gene expression patterns.
Our biological understanding of S. coelicolor metabolism is further enhanced by a more detailed analysis of the reactions for which the flux balance analysis predicted zero flux. When clustering the measured gene expression profiles for the genes encoding the enzymes of these zero-flux reactions, a substantial number of genes showed consistent changes in gene expression along the time course, suggesting that the corresponding reactions are in fact active (Additional file 1). Striking examples include a large number of genes for vitamin B12 (cobalamin) biosynthesis, a group of ten genes involved in calcium-dependent antibiotic (CDA) biosynthesis, and three genes involved in ectoine biosynthesis (Additional files 2, 3, 4 and 5). Each of these cases provides important insights: the first one shows that vitamin B12 is likely to be produced by S. coelicolor under the growth conditions of our experiment, even if it is not essential due to the availability of cobalamine-independent enzymes . The second one highlights that CDA biosynthesis genes are coherently induced in expression during the metabolic switch, similar to undecylprodigiosin and actinorhodin and concordant with the results of the correlation analysis discussed above. This could indicate that this additional antibiotic compound is potentially also produced in phosphate starvation conditions, contrary to previous expectations . Finally, the case of ectoine biosynthesis genes suggests that this novel osmoprotectant metabolite is produced by S. coelicolor. This has in fact been experimentally confirmed recently . In each of these cases, the activity of the pathway was not predicted, based on the biological evidence incorporated in the stoichiometric model and the expected biomass composition, and the comparison of flux balance predictions and gene expression data indicated relevant modifications of our metabolic model. A complete list of genes that have zero predicted flux but show gene expression is included in the supplementary material (Additional files 2, 3, 4 and 5).
Conversely, our model can be used to identify those genes that are predicted to be essential for growth (non-zero flux under all conditions), but show no or very low gene expression. There are 159 predicted essential genes in our model, which have a median log gene expression level of 7.47, compared to 6.83 for the non-essential genes and 4.66 for the negative controls. This indicates that on average the essential genes have a 60% higher expression than the non-essential genes. There is only one predicted essential gene with a detected median expression level below 5.0, compared to 23 non-essential genes with such low expression levels. This non-expressed essential gene is panB (SCO2256), a 3-methyl-2-oxobutanoate hydroxymethyltransferase of pantothenate and coenzyme A biosynthesis, which has a maximum log expression signal of only 5.53. Its apparent non-expression can be due to insufficient hybridization of the gene-specific probes on the microarray, but it could also indicate the existence of another isoenzyme or additional metabolic pathways that would make this reaction redundant. In both of these cases, this gene might warrant further detailed study.
The observed good correlation between gene expression and predicted metabolic flux is not necessarily expected; expression levels can show little correlation to protein levels, enzyme activity and metabolic flux for many reasons . It could be that the relationship between expression and flux is tighter in prokaryotes like S. coelicolor, than in multicellular eukaryotic model organisms [15, 16]. However, we cannot exclude that the group of non-correlated genes contains not only reactions with constant flux, but also reactions with dynamic flux little correlation between gene expression and protein activity or metabolic flux. In a next step, it will be interesting to directly incorporate the gene expression information in the model, providing additional constraints on the maximum flux [23, 24].
Our study demonstrates the ability of flux balance analysis to not only study classical steady-state conditions but also to predict microbial behaviour in dynamic growth conditions provided that sufficiently detailed measurements of the changing growth conditions (nutrient uptake) and cellular objective (antibiotic production rate) are available. In combination with detailed gene expression information, these dynamic model predictions can help identifying potential new players in the metabolic switch, including putative new genes for antibiotic synthesis.
The gene expression dataset used in this study has been described in detail in . Briefly, S. coelicolor was cultivated in a phosphate limiting defined medium containing glucose as a carbon source and glutamate as a nitrogen as well as carbon source. Samples for transcriptomics and off-line analysis were taken every hour from 20 to 44 hours after inoculation (25 sample points), and subsequently every second hour from 46 to 60 hours after inoculation. Cell dry weight was measured on samples collected every third hour between 20 and 40 hours. The last sample, collected at the end of the fermentation (68 hours after inoculation), was used for analysis of remaining nutrients and total production levels of red and blue pigments. Only one sample was collected at each time point and no re-samplings were performed. Gene expression was measured on custom-made Affymetrix gene chips as described in . Expression data have been deposited in the GEO database under accession number GSE18489. Measurements for all known or predicted enzyme-coding genes were extracted and matched to the corresponding reactions in the constraint-based model.
Constraints-based genome-scale metabolic model reconstruction
A genome-scale stoichiometric metabolic model of Streptomyces coelicolor was reconstructed from different sources of data, including KEGG pathways, ScoCyc pathways, biochemistry textbooks, an extensive literature survey and available genome-scale models of other organisms. The initial stoichiometric matrix was generated based on KEGG and ScoCyc and manually curated to refine the S. coelicolor-specific parts of the metabolic network (e.g., antibiotic biosynthesis), to specify the correct reversibility constraints of reactions, and to add missing essential reactions. Missing essential reactions were identified iteratively; a minimum set of hypothetical reactions was added to the model if an essential metabolite could not be produced otherwise. Reversibility and essentiality of reactions were also compared to other published genome-scale models of S. coelicolor and other organisms [6, 25–27]. The resulting model is very similar to the model of Borodina et al. [5, 6], and differs mainly in the more comprehensive inclusion of antibiotic pathways.
In the final curated model, one lumped reaction is added to produce the biomass of the cell. Information about biomass composition and growth and non-growth associated ATP maintenance were taken from Borodina et al.  and Ingraham et al.  and complemented with literature information [29, 30]. Some of the biomass precursor biosynthesis reactions are also lumped reactions, e.g. protein translation, and were specified according to the literature and published genome-scale models . The full model in SBML format is available in the supplementary material (Additional file 2).
Analysis of the model was based on standard flux balance analysis (FBA) to predict optimal in silico growth and metabolic flux distribution using the COBRA tool . Uptake fluxes for metabolites not available in the medium were set to zero, while metabolic by-products were always allowed to leave the metabolic system. Observed nutrient uptake rates from the fermenter culture used for the transcriptome analysis were used to define the constraints of nutrients uptake for the model (input function). The objective function was defined as maximizing the growth rate. Beginning at 34 hours, we dynamically varied the biomass composition by adding increasing amounts of antibiotics, based on the observed antibiotics production rate.
Comparing transcriptome data and predicted flux
Our computational model contains 643 metabolites and 1015 reactions: 747 reactions for metabolite biosynthesis and degradation, 152 transport reactions, and 116 additional input and output constraints of the system. 666 reactions were annotated as enzyme-catalyzed reactions and could be matched to an enzyme-coding gene. Some reactions were annotated as potentially catalyzed by more than one gene and some genes catalyze more than one reaction. If one gene catalyzes multiple reactions, we matched its expression profile to the reaction with the maximum predicted flux, hypothesizing that this reaction will dominate the expression behavior. In total, 789 genes are assigned to 666 enzymatic reactions. Of these, 558 genes are predicted to have non-zero flux (the remaining 231 genes are not used for biomass production according to the model). Out of these 558 genes, 9 genes were involved in cell maintenance with constant flux and zero standard deviation; these were excluded from the further analysis. In total we therefore considered 549 enzyme-coding genes with non-zero predicted flux. For each of these genes, we compared the predicted flux profile and the observed gene expression levels using Pearson's correlation, testing whether gene expression was indeed upregulated when a much higher flux through a particular reaction was required at a certain growth phase.
We are very grateful to the STREAM consortium (funded by ERA-NET SySMO [GEN2006-27745-E/SYS]: (P-UK-01-11-3i) and the Research Council of Norway [project no. 181840/I30]) for providing the expression data prior to publication. The STREAM consortium is an international project funded by a SysMO grant of the ERASysBio framework http://www.erasysbio.net studying "Global metabolic switching in Streptomyces coelicolor". Members of the consortium include Florian Battke, Sven Even Borgos, Per Bruheim, Nigel Burroughs, Gregory Challis, Lubbert Dijkhuizen, Trond E Ellingsen, William Gaze, Alexander Herbig, Sunniva Hoel, Øyvind Jakobsen, Ritsert C. Jansen, Brent Kiernan, Preben Krabben, Roxane Legaie, Juan Francisco Martín, Jay Moore, Kay Nieselt, Walid Omara, David Rand, Jens Reuther, Antonio Rodríguez-García, Håvard Sletta, Maggie Smith, Louise Thomas, Alexander Wentzel, David Wild, Wolfgang Wohlleben, Arouna Woukeu, Anders Øverby, and the authors https://www.wsbc.warwick.ac.uk/groups/sysmopublic/. MEM was funded by a 4 × 4 Ubbo Emmius scholarship and ET by a Rosalind Franklin Fellowship, both from the University of Groningen. RB is supported by an NWO-Vidi fellowship.
- Kolter R, Siegele DA, Tormo A: The stationary phase of the bacterial life cycle. Annu Rev Microbiol. 1993, 47: 855-874. 10.1146/annurev.mi.47.100193.004231.PubMedView ArticleGoogle Scholar
- Roszak DB, Colwell RR: Survival strategies of bacteria in the natural environment. Microbiol Rev. 1987, 51 (3): 365-379.PubMed CentralPubMedGoogle Scholar
- Bentley SD, Chater KF, Cerdeno-Tarraga AM, Challis GL, Thomson NR, James KD, Harris DE, Quail MA, Kieser H, Harper D: Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2). Nature. 2002, 417 (6885): 141-147. 10.1038/417141a.PubMedView ArticleGoogle Scholar
- Challis GL, Hopwood DA: Synergy and contingency as driving forces for the evolution of multiple secondary metabolite production by Streptomyces species. Proc Natl Acad Sci USA. 2003, 100 (Suppl 2): 14555-14561. 10.1073/pnas.1934677100.PubMed CentralPubMedView ArticleGoogle Scholar
- Borodina I, Siebring J, Zhang J, Smith CP, van Keulen G, Dijkhuizen L, Nielsen J: Antibiotic overproduction in Streptomyces coelicolor A3 2 mediated by phosphofructokinase deletion. J Biol Chem. 2008, 283 (37): 25186-25199. 10.1074/jbc.M803105200.PubMedView ArticleGoogle Scholar
- Borodina I, Krabben P, Nielsen J: Genome-scale analysis of Streptomyces coelicolor A3(2) metabolism. Genome Res. 2005, 15 (6): 820-829. 10.1101/gr.3364705.PubMed CentralPubMedView ArticleGoogle Scholar
- Nieselt K, Battke F, Herbig A, Bruheim P, Wentzel A, Jakobsen OM, Sletta H, Alam MT, Merlo ME, Moore J: The dynamic architecture of the metabolic switch in Streptomyces coelicolor. BMC Genomics. 2010, 11: 10-10.1186/1471-2164-11-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Price ND, Papin JA, Schilling CH, Palsson BO: Genome-scale microbial in silico models: the constraints-based approach. Trends Biotechnol. 2003, 21 (4): 162-169. 10.1016/S0167-7799(03)00030-1.PubMedView ArticleGoogle Scholar
- Price ND, Reed JL, Palsson BO: Genome-scale models of microbial cells: evaluating the consequences of constraints. Nat Rev Microbiol. 2004, 2 (11): 886-897. 10.1038/nrmicro1023.PubMedView ArticleGoogle Scholar
- Durot M, Bourguignon PY, Schachter V: Genome-scale models of bacterial metabolism: reconstruction and applications. FEMS Microbiol Rev. 2009, 33 (1): 164-190. 10.1111/j.1574-6976.2008.00146.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee JM, Gianchandani EP, Papin JA: Flux balance analysis in the era of metabolomics. Brief Bioinform. 2006, 7 (2): 140-150. 10.1093/bib/bbl007.PubMedView ArticleGoogle Scholar
- Covert MW, Schilling CH, Palsson B: Regulation of gene expression in flux balance models of metabolism. J Theor Biol. 2001, 213 (1): 73-88. 10.1006/jtbi.2001.2405.PubMedView ArticleGoogle Scholar
- Lee JM, Gianchandani EP, Eddy JA, Papin JA: Dynamic analysis of integrated signaling, metabolic, and regulatory networks. PLoS Comput Biol. 2008, 4 (5): e1000086-10.1371/journal.pcbi.1000086.PubMedView ArticleGoogle Scholar
- Melzoch K, de Mattos MJ, Neijssel OM: Production of actinorhodin by Streptomyces coelicolor A3(2) grown in chemostat culture. Biotechnol Bioeng. 1997, 54 (6): 577-582. 10.1002/(SICI)1097-0290(19970620)54:6<577::AID-BIT8>3.0.CO;2-G.PubMedView ArticleGoogle Scholar
- Wagner A: Energy costs constrain the evolution of gene expression. J Exp Zool B Mol Dev Evol. 2007, 308 (3): 322-324. 10.1002/jez.b.21152.PubMedView ArticleGoogle Scholar
- Wagner A: Energy constraints on the evolution of gene expression. Mol Biol Evol. 2005, 22 (6): 1365-1374. 10.1093/molbev/msi126.PubMedView ArticleGoogle Scholar
- Hojati Z, Milne C, Harvey B, Gordon L, Borg M, Flett F, Wilkinson B, Sidebottom PJ, Rudd BA, Hayes MA: Structure, biosynthetic origin, and engineered biosynthesis of calcium-dependent antibiotics from Streptomyces coelicolor. Chem Biol. 2002, 9 (11): 1175-1187. 10.1016/S1074-5521(02)00252-1.PubMedView ArticleGoogle Scholar
- Cerdeno AM, Bibb MJ, Challis GL: Analysis of the prodiginine biosynthesis gene cluster of Streptomyces coelicolor A3(2): new mechanisms for chain initiation and termination in modular multienzymes. Chem Biol. 2001, 8 (8): 817-829. 10.1016/S1074-5521(01)00054-0.PubMedView ArticleGoogle Scholar
- Apel AK, Sola-Landa A, Rodriguez-Garcia A, Martin JF: Phosphate control of phoA, phoC and phoD gene expression in Streptomyces coelicolor reveals significant differences in binding of PhoP to their promoter regions. Microbiology. 2007, 153 (Pt 10): 3527-3537. 10.1099/mic.0.2007/007070-0.PubMedView ArticleGoogle Scholar
- Martens JH, Barg H, Warren MJ, Jahn D: Microbial production of vitamin B12. Appl Microbiol Biotechnol. 2002, 58 (3): 275-285. 10.1007/s00253-001-0902-7.PubMedView ArticleGoogle Scholar
- Kim HB, Smith CP, Micklefield J, Mavituna F: Metabolic flux analysis for calcium dependent antibiotic (CDA) production in Streptomyces coelicolor. Metab Eng. 2004, 6 (4): 313-325. 10.1016/j.ymben.2004.04.001.PubMedView ArticleGoogle Scholar
- Kol S, Merlo ME, Scheltema RA, de Vries M, Vonk RJ, Kikkert NA, Dijkhuizen L, Breitling R, Takano E: Metabolomic characterization of the salt stress response in Streptomyces coelicolor. Appl Environ Microbiol.
- Akesson M, Forster J, Nielsen J: Integration of gene expression data into genome-scale metabolic models. Metab Eng. 2004, 6 (4): 285-293. 10.1016/j.ymben.2003.12.002.PubMedView ArticleGoogle Scholar
- Colijn C, Brandes A, Zucker J, Lun DS, Weiner B, Farhat MR, Cheng TY, Moody DB, Murray M, Galagan JE: Interpreting expression data with metabolic flux models: predicting Mycobacterium tuberculosis mycolic acid production. PLoS Comput Biol. 2009, 5 (8): e1000489-10.1371/journal.pcbi.1000489.PubMed CentralPubMedView ArticleGoogle Scholar
- Jamshidi N, Palsson BO: Investigating the metabolic capabilities of Mycobacterium tuberculosis H37Rv using the in silico strain iNJ661 and proposing alternative drug targets. BMC Syst Biol. 2007, 1: 26-10.1186/1752-0509-1-26.PubMed CentralPubMedView ArticleGoogle Scholar
- Reed JL, Vo TD, Schilling CH, Palsson BO: An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). Genome Biol. 2003, 4 (9): R54-10.1186/gb-2003-4-9-r54.PubMed CentralPubMedView ArticleGoogle Scholar
- 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 (39): 28791-28799. 10.1074/jbc.M703759200.PubMedView ArticleGoogle Scholar
- Ingraham JL, Maaløe O, Neidhardt FC: Growth of the bacterial cell. 1983, Sinauer, Sunderland, MAGoogle Scholar
- Shahab N, Flett F, Oliver SG, Butler PR: Growth rate control of protein and nucleic acid content in Streptomyces coelicolor A3(2) and Escherichia coli B/r. Microbiology. 1996, 142 (Pt 8): 1927-1935. 10.1099/13500872-142-8-1927.PubMedView ArticleGoogle Scholar
- Zuneda MC, Guillenea JJ, Dominguez JB, Prado A, Goni FM: Lipid composition and protoplast-forming capacity of Streptomyces antibioticus. Lipids. 1984, 19 (3): 223-228. 10.1007/BF02534801.PubMedView ArticleGoogle Scholar
- 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 (3): 727-738. 10.1038/nprot.2007.99.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.