Fluxconsistency

findFluxConsistentSubset(model, param, printLevel)[source]

Finds the subset of S that is flux consistent using various algorithms, but fastcc from fastcore by default

USAGE:

[fluxConsistentMetBool, fluxConsistentRxnBool, fluxInConsistentMetBool, fluxInConsistentRxnBool, model] = findFluxConsistentSubset (model, param, printLevel)

INPUTS:

model – structure with field:

  • .S - m x n stoichiometric matrix

OPTIONAL INPUTS:

  • param – can contain: * param.LPsolver - the LP solver to be used * param.epsilon - minimum nonzero flux, default feasTol*10

    Note that fastcc is very sensitive to the value of parm.epsilon

    • param.modeFlag - {(0),1} 1 = return flux modes

    • param.method - {‘swiftcc’, (‘fastcc’), ‘dc’,’fastB’}

    • param.reduce - {(0),1} 1 = return fluxConsistModel

    • param.signedCyclicFluxConsistency -

    {(0),1} 1 = identify reactions that are forward flux consistent for a forward cycle, or reverse flux consistent for a reverse cycle

  • printLevel – verbose level

  • model.rev - the 0-1 vector with 1’s corresponding to the reversible reactions (if using swiftcc)

OUTPUTS:

  • fluxConsistentMetBool – m x 1 boolean vector indicating flux consistent mets

  • fluxConsistentRxnBool – n x 1 boolean vector indicating flux consistent rxns

  • fluxInConsistentMetBool – m x 1 boolean vector indicating flux inconsistent mets

  • fluxInConsistentRxnBool – n x 1 boolean vector indicating flux inconsistent rxns

  • model – structure with fields duplicating the single output arguments:

    • .fluxConsistentMetBool

    • .fluxConsistentRxnBool

    • .fluxInConsistentMetBool

    • .fluxInConsistentRxnBool

findSignedCyclicFluxConsistentSubset(model, param)[source]

Find the subset of reactions that is signed, cyclically flux consistent Reaction j is said to be signed, cyclically flux consistent if

  1. there exists a v, such that N*v = 0

and either (2.1) or (2.2) hold, where

(2.1) v(j) >= param.epsilon (2.2) v(j) <= -param.epsilon

and either (3.1) or (3.2) hold, where

(3.1) v(k) >= 0, for all k in 1..n, except k = j (3.2) v(k) <= 0, for all k in 1..n, except k = j

where N is an (internal) stoichiometric matrix

USAGE:

[signedCyclicFluxConsistentMetBool, signedCyclicFluxConsistentRxnBool, model, fluxConsistModel] = findSignedCyclicFluxConsistentSubset (model, param)

INPUTS:

model – structure with field:

  • .S - m x n stoichiometric matrix

OPTIONAL INPUTS:

param – can contain: * param.LPsolver - the LP solver to be used * param.epsilon - minimum nonzero flux, default feasTol*10

Note that fastcc is very sensitive to the value of parm.epsilon

  • param.modeFlag - {(0),1} 1 = return flux modes

  • printLevel - verbose level

OUTPUTS:

  • signedCyclicFluxConsistentMetBool – m x 1 boolean vector indicating flux consistent mets

  • signedCyclicFluxConsistentRxnBool – n x 1 boolean vector indicating flux consistent rxns

  • model – structure with fields duplicating the single output arguments: * .signedCyclicFluxConsistentMetBool * .signedCyclicFluxConsistentRxnBool