Team:Alberta/Project/Modeling/FindEssentials

From 2009.igem.org

(Difference between revisions)
(New page: %start.m script changeCobraSolver('glpk'); model = readCbModel() %('Ec_iAF1260_flux1.xml') solution1 = optimizeCbModel(model, 'max', 'false', 'false'); essgenes=cell(1) fValue = so...)
 
(2 intermediate revisions not shown)
Line 1: Line 1:
-
  %start.m script
+
  %Any single gene knockout resulting in a growth rate less than CutOffCoeff*InitialGrowthRate will be put into a list
-
  changeCobraSolver('glpk');  
+
   
-
  model = readCbModel() %('Ec_iAF1260_flux1.xml')
+
function [Essentials,grRatio,grRateKO,grRateWT,hasEffect,fluxSolution] = FindEssentials(Model,CutOffCoeff)
-
  solution1 = optimizeCbModel(model, 'max', 'false', 'false');
+
   
-
  essgenes=cell(1)
+
  Essentials=cell(1);
-
  fValue = solution1.f
+
b=1;
-
  startSays = 'If comparing researched/suspected essential genes, dont forget to add s0001 to the essgenes list!'
+
  solution1 = optimizeCbModel(Model, 'max', 'false', 'false');
 +
  InitialGrowthRate = solution1.f
 +
CutOffGrowthRate=CutOffCoeff*InitialGrowthRate;
 +
[grRatio,grRateKO,grRateWT,hasEffect,fluxSolution]=singleGeneDeletion(Model);
 +
   
 +
%This part uses the results of singleGeneDeletion to make a list
 +
%containing all of the genes whose deletion caused a growth rate less than
 +
%CutOffGrowthRate
 +
 +
for a=(1:length(grRateKO))
 +
    if grRateKO(a)<CutOffGrowthRate
 +
        Essentials(b) = Model.genes(a);
 +
        b=b+1;
 +
    end
 +
end
 +
 +
  NumberOfEssentialGenes = length(Essentials)
 +
Essentials = Essentials';

Latest revision as of 00:45, 18 October 2009

%Any single gene knockout resulting in a growth rate less than CutOffCoeff*InitialGrowthRate will be put into a list

function [Essentials,grRatio,grRateKO,grRateWT,hasEffect,fluxSolution] = FindEssentials(Model,CutOffCoeff)

Essentials=cell(1);
b=1;
solution1 = optimizeCbModel(Model, 'max', 'false', 'false');
InitialGrowthRate = solution1.f
CutOffGrowthRate=CutOffCoeff*InitialGrowthRate;
[grRatio,grRateKO,grRateWT,hasEffect,fluxSolution]=singleGeneDeletion(Model);

%This part uses the results of singleGeneDeletion to make a list
%containing all of the genes whose deletion caused a growth rate less than
%CutOffGrowthRate

for a=(1:length(grRateKO))
    if grRateKO(a)<CutOffGrowthRate
        Essentials(b) = Model.genes(a);
        b=b+1;
    end
end

NumberOfEssentialGenes = length(Essentials)
Essentials = Essentials';
Retrieved from "http://2009.igem.org/Team:Alberta/Project/Modeling/FindEssentials"