Exploration¶

biomassPrecursorCheck(model, checkCoupling, checkConservedQuantities, ATN)[source]¶

Checks if biomass precursors are able to be synthesized.

Usage

[missingMets, presentMets, coupledMets, missingCofs, presentCofs] = biomassPrecursorCheck(model, checkCoupling, checkConservedQuantities)

Input

model – COBRA model structure

Optional inputs

checkCoupling – Test, whether some compounds can only be produced if there is a sink for other biomass precursors (Default: false)

checkConservedQuantities – true to check whether the cofactor pairs containing conserved moieties (defined by the network structure) can be synthesized (e.g., ATP, NAD, NADPH, ACCOA, AA-charged tRNA, fatty acyl-ACP). They will otherwise be identified as missingMets (Default: false)

ATN – atom transition network outputed by buildAtomTransitionNetwork If provided, true internal conserved moieties will be identified and used for checking conserved quantities (default [])

Outputs

missingMets – List of biomass precursors that are not able to be synthesized

presentMets – List of biomass precursors that are able to be synthesized

coupledMets – List of metabolites which need an exchange reaction for at least one other biomass component because their production is coupled to it (returned only if checkCoupling = true)

missingCofs – List of cofactor pairs (defined by the network conserved moieties) that are not able to be synthesized (returned only if checkConservedQuantities = true)

presentCofs – List of cofactor pairs that are able to be synthesized (returned only if checkConservedQuantities = true)

essentialRxn4MultipleModels(allModels, objFun)[source]¶

Performs single reactions deletions across multiple models and identifies which reactions have variable essentiality across all models for the chosen objective function (eg. it identifies reactions whos deletion would result in a zero flux through ATP consumption reaction - ATPM).

Usage

[essentialRxn4Models, dataStruct] = essentialRxn4MultipleModels(allModels, objFun)

Inputs

allModels – directory of the structure with multiple COBRA model structures

objFun – string with objective function reaction (e.g. ‘ATPM’)

Outputs

essentialRxn4Models – Table with reaction fluxes (within the objective function reaction) after single deletion of model reaction (rows) across models (columns)

dataStruct – Structure with all models

Example

[essentialRxn4Models, dataStruct] = essentialRxn4MultipleModels(allModels, 'ATPM')

findBlockedReaction(model, method)[source]¶

Determines those reactions which cannot carry any flux in the given simulation conditions.

Usage

BlockedReaction = findBlockedReaction(model)

Input

model – COBRA model structure

Optional input

method – ‘FVA’ for flux variability analysis (default) ‘L2’ for 2-norm minimization

Output

blockedReactions – List of blocked reactions

findCarbonRxns(model, nCarbonThr)[source]¶

Returns the list of reactions that act of compounds which contain cabons greater than the thershhold set.

Usage

[hiCarbonRxns, zeroCarbonRxns, nCarbon] = findCarbonRxns(model, nCarbonThr)

Inputs

model – Structure containing all necessary variables to described a stoichiometric model

nCarbonThr – defines the min # of carbons that a metabolite, that is acted on in a reaction, can have in the final list of reactions

Outputs

hiCarbonRxns – The list of reactions that act on metabolites with greater than the threshold number of carbons

zeroCarbonRxns – Reactions with no carbon

nCarbon – The number of carbons in each metabolite in the model

findEpistaticInteractions(model, doubleDeletionFitness, lethalFlag, minEffect)[source]¶

Finds synthetic lethal and/or synthetic sick interactions based on double deletion analysis data

Usage

[interactions, epistaticEffect] = findEpistaticInteractions(model, doubleDeletionFitness, lethalFlag, minEffect)

Inputs

model – COBRA model structure

doubleDeletionFitness – A matrix of fitness (or growth rate) values for each of the double deletion strains. The diagonal of this matrix contains the single deletion fitness values.

Optional inputs

lethalFlag – Only consider SL interactions (Default = false)

minEffect – Minimum fitness effect considered to be significant (Default = 1e-2)

Outputs

interactions – A sparse binary matrix indicating a SL or SS interaction between two genes in the model

epistaticEffect – Magnitude of the epistatic interaction defined as min(f1-f12, f2-f12) where f1 and f2 are the fitness values for the deletion strain of gene 1 and gene 2 respectively and f12 is the fitness value for the double deletion strain of genes 1 and 2

Note

The criteria for establishing a synthetic sick interaction are that the double deletion strain fitness must be at least minEffect lower than the fitness of either of the single deletion strains, i.e. f12 < f1 - minEffect and f12 < f2 - minEffect

The additional criterion for establishing a synthetic lethal interaction is that the double deletion fitness value is smaller than minEffect (i.e. essentially zero) f12 < minEffect

Note that the interactions matrix double counts all interactions

findExcRxns(model, inclObjFlag, irrevFlag)[source]¶

Finds exchange and uptake rxns

Usage

[selExc, selUpt] = findExcRxns(model, inclObjFlag, irrevFlag)

Input

model – COBRA model structure

Optional inputs

inclObjFlag – Include objective rxns in the exchange rxn set (1) or not (0) (Default = false)

irrevFlag – Model is in irreversible format (1) or not (Default = false)

Outputs

selExc – Boolean vector indicating whether each reaction in model is exchange or not

selUpt – Boolean vector indicating whether each reaction in model is nutrient uptake or not

Note

Exchange reactions only have one non-zero (+1 / -1) element in the corresponding column of the stoichiometric matrix. Uptake reactions are exchange reactions are exchange reactions with negative lower bounds.

findGPRFromRxns(model, rxnIDs)[source]¶

Get the Textual representations of the GPR rules for the indicated reactions. Usage

gprs = findGPRFromRxns(model,rxnIDs)

Inputs

model – The model to retrieve the GPR rules from

rxnIDs – The reaction IDs to obtain the GPR rules for

Output

gprs – The textual GPR Rules

Authors:

Thomas Pfau Dec 2017

findGeneIDs(model, geneList)[source]¶

Finds gene numbers in a model

Usage

geneID = findGeneIDs(model, geneList)

Inputs

model – COBRA model structure

geneList – List of genes

Output

geneID – List of gene IDs corresponding to geneList

findGenesFromRxns(model, reactions)[source]¶

Make a gene list of the genes that correspond to the selected reactions

Usage

[geneList] = findGenesFromRxns(model, reactions)

Inputs

model – COBRA model structure

reactions – Reactions to find the corresponding genes for

Output

geneList – List of genes corresponding to reactions

findHiCarbonRxns(model, nCarbonThr)[source]¶

Returns the list of reactions that act on compounds which contain cabons greater than the threshold set.

Usage

[hiCarbonRxns, nCarbon] = findHiCarbonRxns(model, nCarbonThr)

Inputs

model – Structure containing all necessary variables to describe a stoichiometric model

nCarbonThr – defines the min # of carbons that a metabolite, that is acted on in a reaction, can have in the final list of reactions

Outputs

hiCarbonRxns – The list of reactions that act on metabolites with greater than the threshold number of carbons

zeroCarbonRxns – Reactions with no carbon

nCarbon – The number of carbons in each metabolite in the model

findMatchingFieldEntries(field, searchString, isAnnotation, minSim)[source]¶

Find elements of the field that are similar to the provided identifier

Usage

[matchingIDs, positons, similarities] = findMatchingFieldEntries(field, searchString, isAnnotation, minSim)

Inputs

field – A Cell Array of strings

searchString – The item to search for

isAnnotation – Whether the strings in the field are annotation strings (multiple identifiers are concatenated with ; and will be separated)

minSim – minimum “similarity” between element and search item.

Outputs

matchingIDs – All IDs which are matching to the search string.

positions – The positions of these IDs in the supplied field.

similarities – Similarities of the matchingIDs with the search String

findMetFromCompartment(model, compartment)[source]¶

Finds all the metabolites and their identifiers in a compartment of interest.

Usage

[compartmentMetabolites] = findMetFromCompartment(model,Compartment)

Inputs

model – COBRA model strcture

compartment – compartment of interest (e.g.: ‘[m]’, ‘[n]’, ‘[e]’, etc.)

Output

compartmentMetabolites – List of metabolites in the compartment of interest

findMetIDs(model, metList)[source]¶

Finds metabolite numbers in a model

Usage

metID = findMetIds(model, metList)

Inputs

model – COBRA model structure

metList – List of metabolites

Output

metID – List of metabolite IDs corresponding to metList

findMetabolicJunctions(model, nRxnsJnc)[source]¶

Finds metabolic branchpoints with different numbers of branches

Usage

validJunctionMets = findMetabolicJunctions(model, nRxnsJnc)

Inputs

model – COBRA model structure

nRxnJnc – Number of reactions to be considered a junction

Output

validJunctionMets – List of junction metabolites

findMetsFromRxns(model, reactions)[source]¶

Get all Metabolites from a set of reactions.

Usage

[metList, stoichiometries] = findMetsFromRxns(model, reactions) [metList] = findMetsFromRxns(model, reactions)

Inputs

model – COBRA model structure

reactions – Reaction IDs (Cell Array) or positions (Double array) to find the corresponding metabolites for

Outputs

metList – If only one output is requested, returns an ordered set of all metabolites involved in the provided reactions. Otherwise, this is a Cell Array of cell arrays containing the metabolites involved in each of the provided reactions.

stoichiometries – this is a Cell array of double arrays of the stoichiometric coefficients corresponding to the reactions in the order of provided reaction ids. If reactions not in the model are provided, those will be represented by empty arrays.

findNeighborRxns(model, rxns, asSingleArray, order, commonMets, withComp)[source]¶

Identifies the reactions and the corresponding genes that are adjacent (having a common metabolite) to a reaction of interest. Useful for characterizing the network around an orphan reaction.

Usage

[neighborRxns, neighborGenes, mets] = findNeighborRxns(model, rxns, asSingleArray, order, commonMets, withComp)

Inputs

model – COBRA model structure

rxns – the target reaction as a string or multiple reactions as cell array

Optional inputs

asSingleArray – If false, then return cell array of cell arrays with neighbor reactions for one particular connecting metabolite and input reaction combination. Else just return all neighbors in one cell array. (Default = false)

order – maximal order of neighbors to be returned (default = 1) order >= 2 only works with asSingleArray = true Neighborhoods of order >= 2 will usually also return the input reactions.

commonMets – Cell array of common metabolites, that should not count as edges between reactions. Use {‘’} if no such metabolite should be included (default = {‘atp’, ‘adp’, ‘h’, ‘h2o’, ‘pi’, ‘ppi’}).

withComp – if commonMets already have a compartment identifier, e.g. ‘atp[m]’, then true (default=false)

Outputs

neighborRxns – the neighboring rxns in the network, (having common metabolites)

neighborGenes – the genes associated with the neighbor rxns

mets – the metabolites in the target reaction

findOrphanRxns(model)[source]¶

Finds all orphan reactions in model (reactions with no associated genes), not including exchange rxns

Usage

orphans = findOrphanRxns(model)

Input

model – a COBRA model with GPRs

Output

orphans – all orphan reactions in the model

findRootNPmets(model, findNCmets)[source]¶

Finds the root no production (and no consumption) metabolites in a model, used by gapFind

Usage

gaps = findRootNPmets(model, findNCmets)

Input

model a COBRA model

Optional input

findNCmets – find no consumption mets as well as no production (default = false)

Output

gaps – all root no production metabolites

findRxnFromCompartment(model, compartment)[source]¶

Finds all the reactions and their identifiers in a compartment of interest.

Usage

[compartmentReactions] = findRxnFromCompartment(model,Compartment)

Inputs

model – COBRA model strcture

compartment – compartment of interest (e.g.: ‘[m]’, ‘[n]’, ‘[e]’, etc.)

Output

compartmentMetabolites – List of reactions in the compartment of interest

findRxnIDs(model, rxnList)[source]¶

Finds reaction numbers in a model

Usage

rxnID = findRxnIDs(model, rxnList)

Inputs

model – COBRA model strcture

rxnList – List of reactions

Output

rxnID – IDs for reactions corresponding to rxnList

findRxnsActiveWithGenes(model, genes)[source]¶

Finds all reactions for which the provided genes show sufficient evidence (i.e. are usfficient according to the GPRs)

Usage

[Reaclist] = findRxnsActiveWithGenes(model, genes)

Inputs

model – COBRA model structure

genes – A list of gene identifiers

Output

Reaclist – Cell array of reactions which are supported by the provided genes

Note

Only reactions which do have a GPR are considered for this function. Reactions without GPRs are ignored, as we don’t have evidence.

findRxnsFromGenes(model, genes, numericFlag, ListResultsFlag)[source]¶

Print every reaction associated with a gene of interest

Usage

[results, ListResults] = findRxnsFromGenes(model, genes, numericFlag, ListResultsFlag)

Inputs

model – COBRA model structure

genes – string of single gene or cell array of multiple genes for which rxns are desired.

Optional inputs

numericFlag – 1 if using Human Recon (Default = 0)

ListResultsFlag – 1 if you want to output ListResults (Default = 0)

Outputs

results – structure containing cell arrays for each gene. Each cell array has one column of rxn abbreviations and one column containing the reaction formula

ListResults – same as above, but in a cell array

findRxnsFromMets(model, metList, varargin)[source]¶

Returns a list of reactions in which at least one metabolite listed in metList participates.

Usage

[rxnList, rxnFormulaList] = findRxnsFromMets(model, metList, printFlag)

Inputs

model – COBRA model structure

metList – Metabolite list

Optional inputs

printFlag – Print reaction formulas to screen (Default = false)

Property/Value – Allowed Properties are: * containsAll - If true only reactions containing all metabolites in metList are returned (Default: false) * printFlag - as above, will overwrite a printFlag set individually (Default: false) * producersOnly - Only return reactions which produce any of the given metabolites. (Default: false) * consumersOnly - Only return reactions which consume any of the given metabolites. (Default: false) * exclusive - Only return those reactions, which do not contain any metabolites but those given (Default: false)

Outputs

rxnList – List of reactions

rxnFormulaList – Reaction formulas coresponding to rxnList

findRxnsFromSubSystem(model, subSystem)[source]¶

Find all reactions which are part of a given subsystem.

Usage

[reactionNames,rxnPos] = findRxnsFromSubSystem(model,subSystem)

Inputs

model – A COBRA model struct with at least rxns and subSystems fields

subSystem – A String identifying a subsystem

Outputs

reactionNames – A Cell array of reaction names

rxnPos – A double array of positions of the reactions in reactionNames in the model (same order).

Usage

%Obtain all reactions with Glycolysis in their respective subSystems

field.

[reactionNames,rxnPos] = findRxnsFromSubSystem(model,’Glycolysis’)

findRxnsInActiveWithGenes(model, genes)[source]¶

Finds all reactions for which the provided genes show sufficient evidence of their absence. (i.e. make the corresponding GPRs always be zero)

Usage

[Reaclist] = findRxnsActiveWithGenes(model, genes)

Inputs

model – COBRA model structure

genes – A list of gene identifiers

Output

Reaclist – Cell array of reactions which are supported by the provided genes

Note

Only reactions which do have a GPR are considered for this function. Reactions without GPRs are ignored, as we don’t have evidence.

findSubSysGen(model)[source]¶

Lists the subsystem that a reaction occurs in encoded by a gene. Returns list of subsystems. If multiple reactions are associated with gene, subsystem of first occurance will be listed.

Usage

[GenSubSystem] = findSubSysGen(model)

Input

model – COBRA model structure

Output

GenSubSystem – array listing genes and subsystmes

findSubSystemsFromGenes(model, genes, varargin)[source]¶

Returns the subsystems associated with the provided genes.

Usage

[geneSubSystems,singleList] = findSubSystemsFromGenes(model,genes,...)

Input

model – COBRA model structure

Optional inputs

genes – The genes to find subSystems for (Default: model.genes)

varargin – Optional arguments as parameter/value pairs, or parameter struct. * structResult - return individual Gene Associations as a struct instead of a Cell Array with each gene represented by a field. Gene names might be adapted to fit to field names and all genes will be prefixed with gene_. * onlyOneSub - Only return one subSystem for each gene.

Output

geneSubSystems – All subsystems associated with the provided genes.

Optional output

singleList – A Cell array of size(nGenes,2), with the first column representing the gene and the second column containing the subSystems associated with the gene.

findTrspRxnFromMet(model, metList, compFlag)[source]¶

Find transport reactions for defined metabolites. Option available to define the compartment of the transport

Usage

[TrspRxns] = findTrspRxnFromMet(model, metList, compFlag)

Inputs

model – COBRA model structure

metList – metabolites list

Optional input

compFlag – compartment of the transport (e.g. ‘c’, cytosol)

Output

TrspRxns – List of transporters reactions

getMIRIAMAnnotations(model, type, varargin)[source]¶

Get the MIRIAM Annotations a struct array of annotation structs for the given type

Usage

annotations = getMIRIAMAnnotations(model,field,...)

Inputs

model – The COBRA model structure.

type – the basic field to get annotations for (e.g. rxn, met, or model

Optional input

varagin – Additional arguments as parameter/value pairs. * bioQualifiers - A Cell array of BioQualifiers to look for if not provided, or empty, all bioQualifiers defined in getBioQualifiers() will be used * ids - A specific ID or IDs to get the annotation data for. Cannot be combined with type model.(Default: model.([type ‘s’])).

Output

annotations – A struct array with the following structure: * .annotation -> A Struct array with one element per element of the given type. * .annotation.cvterms -> A Struct array of controlled vocabulary terms with one element per qualifier used * .annotation.cvterms.qualifier -> the bioQualifier for this all ressourcesof this cvterm * .annotation.cvterms.qualifierType -> the Qualifier type (modelQualifier or bioQualifier) for all ressources of this cvterm * .annotation.cvterms.ressources -> struct with the following fields: * .annotation.cvterms.ressources.database -> the database for this ressource. * .annotation.cvterms.ressources.id-> The ID in the database for this ressource.

getModelSizes(model)[source]¶

Get the sizes of the basic fields of the model structure

Usage

[nMets, nRxns, nCtrs, nVars, nGenes, nComps] = getModelSizes(model)

Input

model – A COBRA model structure

Outputs

nMets – The number of metabolites in the model

nRxns – The number of reactions in the model

nCtrs – The number of constraints in the model

nVars – The number of variables in the model

nGenes – The number of genes in the model

nComps – The number of compartments in the model

getModelSubSystems(model)[source]¶

Get all subSystems present in a model

Usage

[subSystems] = getModelSubSystems(model)

Input

model – A COBRA model struct with at least the subSystems fields

Output

subSystems – A Cell Array of strings containing all subSystems in the model

Usage

%Get all subSystems present in the model. [subSystems] = getModelSubSystems(model)

getObjectiveSense(model)[source]¶

Get the objective sense of the model (both the osenseStr (‘max’ or ‘min’) and the correspdoning osense value (-1 or 1)

Usage

[osenseStr,osense] = getObjectiveSense(model)

Input

model – The model to obtain the sense for. If the model has a osenseStr field, it has to be either ‘min’ or ‘max’ (case insensitive) otherwise this function will error.

OUTPUS:: osenseStr: ‘min’ or ‘max’ osense: 1 or -1;

identifyFastBlockedRxns(model, rxnList, printLevel)[source]¶

This function evaluates the presence of blocked reactions in a metabolic model

Usage

[BlockedRxns] = identifyFastBlockedRxns(model,rxnList)

Inputs

organisms – model in COBRA model structure format

rxnList – nx1 cell array with reactions to test

printLevel – Verbose level (default: printLevel = 1)

Output

BlockedRxns – nx1 cell array containing blocked reactions

isReactionInSubSystem(model, reactions, subSystem)[source]¶

Determine whether a reaction is in a given subSystem.

Usage

[subSystems] = getModelSubSystems(model)

Inputs

model – A COBRA model struct with at least rxns and subSystems fields

reactions – Either a string identifying a reaction, or a cell array of strings identifying multiple reactions, or a double vector identifying the positions.

subSystem – A String identifying a subsystem.

Output

present – a boolean vector for each provided reaction.

Usage

%Get all subSystems present in the model. [subSystems] = getModelSubSystems(model)

plotEssentialRxns(essentialRxn4Models, essentialityRange, numModelsPresent)[source]¶

Identifies and plots (heatmap) a set of essential reactions (reactions of interest) based on conditional inputs.

Usage

rxnInterest4Models = plotEssentialRxns(essentialRxn4Models, essentialityRange, numModelsPresent)

Inputs

essentialRxn4Models – Table with reaction fluxes (within the objective function reaction) after single deletion of model reaction (rows) across models (columns) This input can be obtained from essentialRxn4MultipleModels.m

essentialityRange – Range of fluxes (e.g. [-100,100])

numModelsPresent – Minimum number of models where a reaction is essential in order to be ploted.

Output

rxnInterest4Models – Table with reaction fluxes (within the objective function reaction) for reactions of interest (rows) across models (columns)

Example

rxnInterest4Models = plotEssentialRxns(essentialRxn4Models, essentialityRange, numModelsPresent)

printConstraints(model, minInf, maxInf)[source]¶

Print all network constraints that are not -Inf (minInf) or +Inf (maxInf)

Usage

printConstraints(model, minInf, maxInf)

Inputs

model – COBRA model structure

minInf – value that is considered as -Inf (or desired minimum cutoff value)

maxInf – value that is considered as +Inf (or desired maximum cutoff value)

printFluxBounds(model, rxns)[source]¶

Prints the reactionID and upper/lower flux bounds.

Usage

printFluxBounds(model, rxns)

Input

model – The model to print

Optional input

rxns – a string array of reaction ids for which flux bounds need to be printed

printFluxVector(model, fluxData, nonZeroFlag, excFlag, sortCol, fileName, headerRow, formulaFlag)[source]¶

Prints a flux vector with reaction labels

Usage

printFluxVector(model, fluxData, nonZeroFlag, excFlag, sortCol, fileName, headerRow, formulaFlag)

Inputs

model – COBRA model structure

fluxData – Data matrix/vector (for example, solution.x)

Optional inputs

nonZeroFlag – Only print nonzero rows (Default = false)

excFlag – Only print exchange fluxes (Default = false)

sortCol – Column used for sorting (-1, none; 0, labels; >0, data columns;) (Default = -1)

fileName – Name of output file (Default = [])

headerRow – Header (Default = [])

formulaFlag – Print reaction formulas (Default = false)

printGPRForRxns(model, rxnIDs, printLevel)[source]¶

Print the GPRs (in textual format) for the given RxnIDs Usage

gprs = printGPRForRxns(model,rxnID)

Inputs

model – The model to retrieve the GPR rules from

rxnIDs – The reaction IDs to obtain the GPR rules for

Optional input

printLevel – Whether to print out the GPRs. If printLevel is 0, the function will only return the Strings in a cell array.

Output

gprs – The textual GPR Rules

Authors:

Thomas Pfau Dec 2017

printLabeledData(labels, data, nonzeroFlag, sortCol, fileName, headerRow, sortMode)[source]¶

Print a matrix of data with labels

Usage

printLabeledData(labels, data, nonzeroFlag, sortCol, fileName, headerRow, sortMode)

Inputs

labels – Row labels

data – Data matrix/vector

nonzeroFlag – Only print nonzero rows (opt)

sortCol – Column used for sorting (-1, none; 0, labels; >0, data columns; opt)

fileName – Name of output file (opt)

headerRow – Header (opt)

sortMode – Sort mode, ‘ascend’ or ‘descend’ (opt, default ‘ascend’)

printMatrix(A, format, file)[source]¶

Prints matrix into a file or screen

Usage

printMatrix(A, format, file)

Inputs

A – Matrix

format – Format string (opt, default ‘%6.4ft’)

file – File name (opt)

printRxnFormula(model, varargin)[source]¶

Prints the reaction formulas for a list of reactions

Reactions that have an upperbound <= 0 and lowerbound < 0 will have its directionality reversed unless directionFlag = false.

Usage

formulas = printRxnFormula(model, varargin)

Input

model – COBRA model structure

Optional input

varargin – Optional Inputs provided as ‘ParameterName’, Value

pairs. the following parameternames are available:

rxnAbbrList: Cell array of rxnIDs to be printed (Default = print all reactions)

printFlag: Print formulas or just return them (Default = true)

lineChangeFlag: Append a line change at the end of each line

(Default = true)

metNameFlag: Print full met names instead of abbreviations

(Default = false)

fid: Optional file identifier for printing in files

(default 1, i.e. stdout)

directionFlag: Checks directionality of reaction. See Note.

(Default = false)

gprFlag: Print gene protein reaction association

(Default = false)

proteinFlag: Print the protein names associated with the genes in the

GPRs associated with the reactions. (Default = false)

printBounds: Print the upper and lower Bounds of the reaction (Default = false)

Output

formulas – Cell array containing formulas of specified reactions

Example

1) print only ATPM and TKT1:
   printRxnFormula(model,'rxnAbbrList',{'ATPM','TKT1'});
2) print the reactions of the model with the metabolite Names instead of ids.
   printRxnFormula(model, 'metNameFlag',true);
3) print all reactions with the metabolite names to a given fileID
   printRxnFormula(model, 'metNameFlag',true, 'fid', fileID);

printUptakeBound(model)[source]¶

Returns the indices to the substrate uptake reactions and prints substrate uptake bounds

Usage

uptakeInd = printUptakeBound(model)

Input

model – CORBRA model structure

Output

upInd – Vector containing indecies of uptake reactions

searchModel(model, searchTerm, varargin)[source]¶

Search for the specified term in the given model. The performed search is fuzzy if similarity is lower than 1.

Usage

results = searchModel(model,searchTerm,...)

Inputs

model – The model to search in

searchTerm – The term to search for

varargin – Additional parameters as parameter/value pairs or value struct. Available parameters are:

printLevel - Printout switch, 0 silent, 1 print results. (Default = 1) .

similarity - Minimum similarity (as provided in calcSim) for matches.

Output

results – The results struct build as follows:

.field - a struct array containing of the basic fields (e.g. rxns, mets etc) of the model in which matches were found with each field indicating the matches as detailed below.

.field.id - The id of the matching element

.field.matches - information about the matches of the respective ID.

.field.matches.value - the value that matched.

.field.matches.source - the field that contained the matching value.

surfNet(model, metrxn, metNameFlag, flux, nonzeroFluxFlag, showMets, field2print, nCharBreak, iterOptions)[source]¶

A simple and convenient tool to navigate through the metabolic network interactively, possibly with given flux vectors, in the Matlab command window using mouse once the starting point is called by command.

Usage

surfNet(model, metrxn, metNameFlag, flux, NonzeroFluxFlag, showMets,field2print, nCharBreak)

Input

model – COBRA model

Optional inputs

metrxn – mets or rxns names, can be a cell array of multiple rxns and mets, or simply a string (default: objective reactions or a random reaction)

metNameFlag – print model.metNames in reaction formulas. (default: false)

flux – flux vectors for the model, a #reactions x #columns matrix. (default: []) If given, the producing and consuming reactions displayed are in accordance with the flux direction.

nonzeroFluxFlag – show only reactions with nonzero fluxes if flux is given (default: true)

showMets – show metabolites in a list when printing reactions (default: true)

field2print – print fields associated with mets or rxns, e.g. metNames, rxnNames, grRules. Either: (i) a cell array of two cells, 1st being a character array for met fields and 2nd for rxn fields or (ii) a character array of field names recognizable from the field names or the sizes (default: {{‘metNames’,’metFormulas’}, {‘rxnNames’,’lb’,’ub’}})

nCharBreak – max. no. of characters per line for reaction formulas (default 0, the current width of the command window)

Example

surfNet(model, 'glc-D[c]')  % for starting the navigation at D-glucose
surfNet(model, 'EX_glc-D(e)')  % for starting the navigation at the glucose exchange reaction
surfNet(model, 'glc-D[c]', 1)  % for printing model.metNames in reaction formulas
surfNet(model, 'glc-D[c]', 0, flux, 1)  % for starting at glucose and show reactions with non-zero flux in 'flux' only
surfNet(model, 'EX_glc-D(e)', 0, [], [], 0)  % to show the reaction formulas only but not the details of metabolites during navigation
surfNet(model, {'glc-D[c]'; 'fru[c]'})  % to view several mets and rxns