http://dx.doi.org/10.1126/science.1184913

SBML model exported from PottersWheel.

% PottersWheel model definition file

function m = BeckerSchilling2010_EpoR_AuxiliaryMode()

m             = pwGetEmptyModel();

%% Meta information

m.ID          = 'BeckerSchilling2010_EpoR_AuxiliaryMode';
m.name        = 'BeckerSchilling2010_EpoR_AuxiliaryModel';
m.description = 'BeckerSchilling2010_EpoR_AuxiliaryModel';
m.authors     = {'Verena Becker',' Marcel Schilling'};
m.dates       = {'2010'};
m.type        = 'PW-2-0-42';

%% X: Dynamic variables
% m = pwAddX(m, ID, startValue, type, minValue, maxValue, unit, compartment, name, description, typeOfStartValue)

m = pwAddX(m, 'EpoR'     ,      76, 'fix'   ,   0, 10000,   [], 'cell', []  , []  , []             , []  , 'protein.generic');
m = pwAddX(m, 'SAv'      , 999.293, 'global', 900,  1100,   [], 'cell', []  , []  , []             , []  , 'protein.generic');
m = pwAddX(m, 'SAv_EpoR' ,       0, 'fix'   ,   0, 10000,   [], 'cell', []  , []  , []             , []  , 'protein.generic');
m = pwAddX(m, 'SAv_EpoRi',       0, 'fix'   ,   0, 10000,   [], 'cell', []  , []  , []             , []  , 'protein.generic');
m = pwAddX(m, 'dSAvi'    ,       0, 'fix'   ,   0, 10000,   [], 'cell', []  , []  , []             , []  , 'protein.generic');
m = pwAddX(m, 'dSAve'    ,       0, 'fix'   ,   0, 10000,   [], 'cell', []  , []  , []             , []  , 'protein.generic');


%% R: Reactions
% m = pwAddR(m, reactants, products, modifiers, type, options, rateSignature, parameters, description, ID, name, fast, compartments, parameterTrunks, designerPropsR, stoichiometry, reversible)

m = pwAddR(m, {            }, {'EpoR'      }, {  }, 'C' , [] , 'k1*k2', {'kt','Bmax_SAv'}, [], 'reaction0001');
m = pwAddR(m, {'EpoR'      }, {            }, {  }, 'MA', [] , []     , {'kt'           }, [], 'reaction0002');
m = pwAddR(m, {'SAv','EpoR'}, {'SAv_EpoR'  }, {  }, 'MA', [] , []     , {'kon_SAv'      }, [], 'reaction0003');
m = pwAddR(m, {'SAv_EpoR'  }, {'SAv','EpoR'}, {  }, 'MA', [] , []     , {'koff_SAv'     }, [], 'reaction0004');
m = pwAddR(m, {'SAv_EpoR'  }, {'SAv_EpoRi' }, {  }, 'MA', [] , []     , {'kt'           }, [], 'reaction0005');
m = pwAddR(m, {'SAv_EpoRi' }, {'SAv'       }, {  }, 'MA', [] , []     , {'kex_SAv'      }, [], 'reaction0006');
m = pwAddR(m, {'SAv_EpoRi' }, {'dSAvi'     }, {  }, 'MA', [] , []     , {'kdi'          }, [], 'reaction0007');
m = pwAddR(m, {'SAv_EpoRi' }, {'dSAve'     }, {  }, 'MA', [] , []     , {'kde'          }, [], 'reaction0008');



%% C: Compartments
% m = pwAddC(m, ID, size,  outside, spatialDimensions, name, unit, constant)

m = pwAddC(m, 'cell', 1);


%% K: Dynamical parameters
% m = pwAddK(m, ID, value, type, minValue, maxValue, unit, name, description)

m = pwAddK(m, 'kt'      , 0.0329366   , 'global', 1e-007, 1000);
m = pwAddK(m, 'Bmax_SAv', 76          , 'fix'   , 61    , 91  );
m = pwAddK(m, 'kon_SAv' , 2.29402e-006, 'global', 1e-007, 1000);
m = pwAddK(m, 'koff_SAv', 0.00679946  , 'global', 1e-007, 1000);
m = pwAddK(m, 'kex_SAv' , 0.01101     , 'global', 1e-007, 1000);
m = pwAddK(m, 'kdi'     , 0.00317871  , 'global', 1e-007, 1000);
m = pwAddK(m, 'kde'     , 0.0164042   , 'global', 1e-007, 1000);


%% Default sampling time points
m.t = 0:3:99;


%% Y: Observables
% m = pwAddY(m, rhs, ID, scalingParameter, errorModel, noiseType, unit, name, description, alternativeIDs, designerProps)

m = pwAddY(m, 'SAv + dSAve'      , 'SAv_extracellular_obs');
m = pwAddY(m, 'SAv_EpoR'         , 'SAv_cellsurface_obs'  );
m = pwAddY(m, 'SAv_EpoRi + dSAvi', 'SAv_intracellular_obs');


%% S: Scaling parameters
% m = pwAddS(m, ID, value, type, minValue, maxValue, unit, name, description)

m = pwAddS(m, 'scale_SAv_extracellular_obs', 1, 'fix', 0, 100);
m = pwAddS(m, 'scale_SAv_cellsurface_obs'  , 1, 'fix', 0, 100);
m = pwAddS(m, 'scale_SAv_intracellular_obs', 1, 'fix', 0, 100);


%% Designer properties (do not modify)
m.designerPropsM = [1 1 1 0 0 0 400 250 600 400 1 1 1 0 0 0 0];