Reconstruction


gapFind(model, findNCgaps, verbFlag)[source]

Identifies all blocked metabolites (anything downstream of a gap) in a model. MILP algorithm that finds gaps that may be missed by simple inspection of the S matrix. To find every gap in a model, change the rxn bounds on all exchange reactions to allow uptake of every metabolite.

USAGE

[allGaps, rootGaps, downstreamGaps] = gapFind (model, findNCgaps, verbFlag)

INPUT

model – a COBRA model

OPTIONAL INPUTS
  • findNCgaps – find no consupmption gaps as well as no production gaps (default false)

  • verbFlag – verbose flag (default false)

OUTPUTS
  • allGaps – all gaps found by GapFind

  • rootGaps – all root no production (and consumption) gaps

  • downstreamGaps – all downstream gaps

Based on Kumar, V. et al. BMC Bioinformatics. 2007 Jun 20;8:212.

testPathway(model, MetIn, MetOut, AdditionalMetsInorOut, ObjectiveOption)[source]

Allows the user to see if given one metabolite A, downstream metabolite B can be made made. Additional sinks can be added for co-factors if needed A –>–>–>–>–> B

USAGE

[Flux, FBAsolution, model] = testPathway (model, MetIn, MetOut, AdditionalMetsInorOut, ObjectiveOption)

INPUTS
  • model – COBRA model structure

  • MetIn – The input metabolite(s) (A)

  • MetOut – The output metabolite (B)

OPTIONAL INPUTS
  • AdditionalMetsInorOut – Additional metabolites for which sinks will be added

  • ObjectiveOption – Boolean:

    • 1 = objective will be production of B (default)

    • 0 = use objective in model

OUTPUTS
  • Flux – The rate of B production

  • FBAsolution – Solution to the FBA problem

  • model – COBRA model with sinks in it