sybilDynFBA (version 0.0.2)

dynamicFBA: dynamic flux balance analysis

Description

Calculate concentrations of metabolites of exchange reactions at defined time points given the initial concentrations. To accomplish this task this function calls optimizeProb function to get the fluxes then update the concentrations and the reaction boundaries ..etc.

Usage

dynamicFBA(model, substrateRxns, initConcentrations, initBiomass, timeStep, nSteps, exclUptakeRxns, fld = FALSE, verboseMode = 2, ...)

Arguments

model

An object of class modelorg.

substrateRxns

List of exchange reaction names for substrates initially in the media that may change (e.g. not h2o or co2)

initConcentrations

The given start concentrations of substrates

initBiomass

The start value of biomass (must be nonzero)

timeStep

Define the points of time to evaluate the problem at.

nSteps

The maximum number of steps, the procedure may stop before completing this number when the substrate run out.

exclUptakeRxns

List of uptake reactions whose substrate concentrations do not change (Default ={'EX_co2(e)','EX_o2(e)','EX_h2o(e)','EX_h(e)'})

fld

Boolean. Save the resulting flux distribution. Default: FALSE

verboseMode

An integer value indicating the amount of output to stdout: 0: nothing, 1: status messages, 2: like 1 plus a progress indicator, 3: a table containing the reaction id's and the corresponding min max values. Default: 2.

...

Further arguments passed to sysBiolAlg. Argument solverParm is a good candidate.

Value

• returns optsol_dynamicFBA

encoding

utf8

References

Varma, A. and Palsson, B.O. 1994. Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110. Appl Environ Microbiol 60: 3724-3731. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox. Nat Protoc 2, 727--738.

See Also

modelorg, optsol_dynamicFBA, optimizeProb, sysBiolAlg, SYBIL_SETTINGS

Examples

## The examples here require the package glpkAPI to be
	## installed. If that package is not available, you have to set
	## the argument 'solver' (the default is: solver = "glpk").

	## load the example data set
	data(Ec_core)
        lowbnd(Ec_core)[react_id(Ec_core)=='EX_glc(e)']=-10;
        lowbnd(Ec_core)[react_id(Ec_core)=='EX_o2(e)']=-18;
	## run dynamicFBA(), Ec_df will be an object of class \code{\link{optsol_dynamicFBA}}
	Ec_df <- dynamicFBA(Ec_core,substrateRxns={'EX_glc(e)'},initConcentrations=10,initBiomass=.035,timeStep=.25,nSteps=20,verbose=3)

	## plot biomass and reactions
	plot(Ec_df,plotRxns=c('EX_glc(e)','EX_ac(e)'));