 %  DOQCS :http://doqcs.ncbs.res.in/ 
 %  Accession Name = PKA_2003 
 %  Accession Number = 47 
 %  Pathway Name = Gs 
 %  Pathway Number = 198 
 %  Transcriber = Upinder S. Bhalla, NCBS 
 %  Developer = Upinder S. Bhalla, NCBS 
 %  Species = Generic mammalian 
 %  Tissue =  Generic 
 %  Cell Compartment = Cytosol + Cell membrane 
 %  Notes =  We assume GTP is present in fixed amounts, so we leave it out of the explicit equations in this model. Normally we would expect it to associate along with the G-Receptor-ligand complex. Most info is from Berstein et al JBC 267:12 8081-8088 1992 Structure of rec activation of Gq from Fay et al Biochem 30 5066-5075 1991
 
  % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Fri Feb 11 15:53:33 2005

% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
tspan=[0, 200];
% Molecule names
molnames = {
'1 /kinetics/Gs/R'
'2 /kinetics/Gs/L.R'
'3 /kinetics/Gs/GDP.Gabc'
'4 /kinetics/Gs/L.R.GDP.Gabc'
'5 /kinetics/Gs/GDP.Ga'
'6 /kinetics/Gs/Gbg'
'7 /kinetics/Gs/R.GDP.Gabc'
'8 /kinetics/GTP.Ga'
};

bufnames = {
'1 /kinetics/Gs/L'
};

sumtotnames = {
};

% Initial concs of 8 molecules
y0(1) = 0.0833333;
y0(2) = 0;
y0(3) = 1;
y0(4) = 0;
y0(5) = 0;
y0(6) = 0;
y0(7) = 0;
y0(8) = 0;
buff(1) = 1;
% Call the ODE solver and display results
[t, y] = ode23s(@f, tspan, y0);

hold on;
plot(t,y(:,4));
plot(t,y(:,2));
plot(t,y(:,6));
plot(t,y(:,8));
plot(t,y(:,5));
plot(t,y(:,7));
% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 - 0.100002 * y(1) * buff(1) + 0.1 * y(2) - 0.199998 * y(3) * y(1) + 0.1 * y(7)
 + 0.100002 * y(1) * buff(1) - 0.1 * y(2) - 10.0002 * y(3) * y(2) + 0.1 * y(4) + 0.025 * y(4)
 + 6 * y(5) * y(6) - 10.0002 * y(3) * y(2) + 0.1 * y(4) - 0.199998 * y(3) * y(1) + 0.1 * y(7)
 + 10.0002 * y(3) * y(2) - 0.1 * y(4) + 4.99998 * buff(1) * y(7) - 0.1 * y(4) - 0.025 * y(4)
 + 0.066667 * y(8) - 6 * y(5) * y(6)
 - 6 * y(5) * y(6) + 0.025 * y(4)
 - 4.99998 * buff(1) * y(7) + 0.1 * y(4) + 0.199998 * y(3) * y(1) - 0.1 * y(7)
 - 0.066667 * y(8) + 0.025 * y(4)
];

% end of model
