Dynamic FBA¶
- dynamicFBA(model, substrateRxns, initConcentrations, initBiomass, timeStep, nSteps, plotRxns, exclUptakeRxns)¶
Performs dynamic FBA simulation using the static optimization approach
USAGE:
[concentrationMatrix, excRxnNames, timeVec, biomassVec] = dynamicFBA(model, substrateRxns, initConcentrations, initBiomass, timeStep, nSteps, plotRxns, exclUptakeRxns)
- INPUTS:
model: COBRA model structure substrateRxns: List of exchange reaction names for substrates
initially in the media that may change (e.g. not h2o or co2)
- initConcentrations: Initial concentrations of substrates (in the same
structure as substrateRxns)
initBiomass: Initial biomass (must be non zero) timeStep: Time step size nSteps: Maximum number of time steps
- OPTIONAL INPUTS:
plotRxns: Reactions to be plotted (Default = {‘EX_glc(e)’, ‘EX_ac(e)’, ‘EX_for(e)’}) exclUptakeRxns: List of uptake reactions whose substrate concentrations do not change
(Default = {‘EX_co2(e)’, ‘EX_o2(e)’, ‘EX_h2o(e)’, ‘EX_h(e)’})
- OUTPUTS:
concentrationMatrix: Matrix of extracellular metabolite concentrations excRxnNames: Names of exchange reactions for the EC metabolites timeVec: Vector of time points biomassVec: Vector of biomass values
If no initial concentration is given for a substrate that has an open uptake in the model (i.e. model.lb < 0) the concentration is assumed to be high enough to not be limiting. If the uptake rate for a nutrient is calculated to exceed the maximum uptake rate for that nutrient specified in the model and the max uptake rate specified is > 0, the maximum uptake rate specified in the model is used instead of the calculated uptake rate.
Note
The dynamic FBA method implemented in this function is essentially the same as the method described in [Varma, A., and B. O. Palsson. Appl. Environ. Microbiol. 60:3724 (1994)]. This function does not implement the dynamic FBA using dynamic optimization approach described in [Mahadevan, R. et al. Biophys J, 83:1331-1340 (2003)].