 %  DOQCS : http://doqcs.ncbs.res.in/ 
 %  Accession Name = RasGTPase 
 %  Accession Number = 71 
 %  Transcriber = Sharat J. Vayttaden and Pragati Jain, NCBS 
 %  Developer = Robert A. Phillips, Jackie L. Hunter, John F. Eccleston, and Martin R. Webb 
 %  Species = General Mammalian 
 %  Tissue = E.coli Expression system 
 %  Cell Compartment = Cytosol 
 %  Notes =  Ras is an important regulator of cell growth in all eukaryotic cells. The model represent hydrolysis of active Ras-bound GTP to give inactive Ras-bound GDP catalyzed by GTPase activating proteins i.e NF1. The inactive Ras-bound GDP turns signalling protein off. 
 
   % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Fri Apr 14 16:38:38 2006
% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
global t;
global y;

% Define Tables here
tspan=[0, 0.6];
% Molecule names
molnames = {
'1 /kinetics/Pi'
'2 /kinetics/RasGTP'
'3 /kinetics/RasGTP-NF1'
'4 /kinetics/RasGTP-NF1*'
'5 /kinetics/RasGDP-NF1_Pi'
'6 /kinetics/RasGDP_NF1'
'7 /kinetics/RasGDP'
'8 /kinetics/NF1'
};

bufnames = {
};

sumtotnames = {
};

% Initial concs of 8 molecules
y0(1) = 0;
y0(2) = 1;
y0(3) = 0;
y0(4) = 0;
y0(5) = 0;
y0(6) = 0;
y0(7) = 0;
y0(8) = 9.99936;

% Call the ODE solver and display results
[t, y] = ode23s(@f, tspan, y0);

hold on;
do_plots;

function do_plots()
global t;
global y;
plot(t,y(:,4));
plot(t,y(:,2));
plot(t,y(:,1));

% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 + 40 * y(5) - 0.000107999 * y(1) * y(6)
 - 1.2 * y(2) * y(8) + 6.36 * y(3)
 + 1.2 * y(2) * y(8) - 6.36 * y(3) - 418 * y(3) + 5.5 * y(4)
 + 418 * y(3) - 5.5 * y(4) - 19.5 * y(4) + 0.22 * y(5)
 + 19.5 * y(4) - 0.22 * y(5) - 40 * y(5) + 0.000107999 * y(1) * y(6)
 + 40 * y(5) - 0.000107999 * y(1) * y(6) - 46.5 * y(6) + 1.2 * y(7) * y(8)
 + 46.5 * y(6) - 1.2 * y(7) * y(8)
 + 46.5 * y(6) - 1.2 * y(7) * y(8) - 1.2 * y(2) * y(8) + 6.36 * y(3)
];

% end of model
