:-dynamic(model/1).
:-dynamic(modelname/1).
:-dynamic(m1/1).
:-dynamic(r/2).
:-dynamic(int/1).

model('kholodenko_M_MAPk_cascade').



m(X):-m1(X=_).
m1('M_MKKK'='90').
m1('M_MKKK_P'='10').
m1('M_MKK'='280').
m1('M_MKK_P'='10').
m1('M_MKK_PP'='10').
m1('M_MAPK'='280').
m1('M_MAPK_P'='10').
m1('M_MAPK_PP'='10').

r('R_a_1',['M_MKKK_P'+1 , 'M_MKKK'-1]).
r('R_a_2',['M_MKKK'+1 , 'M_MKKK_P'-1]).
r('R_b_1',['M_MKKK_P'+1 , 'M_MKK_P'+1 , 'M_MKKK_P'-1 , 'M_MKK'-1]).
r('R_b_2',['M_MKKK_P'+1 , 'M_MKK_PP'+1 , 'M_MKKK_P'-1 , 'M_MKK_P'-1]).
r('R_b_4',['M_MKK'+1 , 'M_MKK_P'-1]).
r('R_b_3',['M_MKK_P'+1 , 'M_MKK_PP'-1]).
r('R_c_1',['M_MKK_PP'+1 , 'M_MAPK_P'+1 , 'M_MKK_PP'-1 , 'M_MAPK'-1]).
r('R_c_2',['M_MKK_PP'+1 , 'M_MAPK_PP'+1 , 'M_MKK_PP'-1 , 'M_MAPK_P'-1]).
r('R_c_4',['M_MAPK'+1 , 'M_MAPK_P'-1]).
r('R_c_3',['M_MAPK_P'+1 , 'M_MAPK_PP'-1]).
