 %  DOQCS :http://doqcs.ncbs.res.in/ 
 %  Accession Name = cAMP_pathway 
 %  Accession Number = 25 
 %  Pathway Name = PKA 
 %  Pathway Number = 135 
 %  Transcriber = Upinder S. Bhalla, NCBS 
 %  Developer = Upinder S. Bhalla, NCBS 
 %  Species = Generic mammalian 
 %  Tissue =  Brain - Neuronal 
 %  Cell Compartment = Cell membrane + Cytosol 
 %  Notes =  General ref: <a href=http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=pubmed&dopt=Abstract&list_uids=2165385&query_hl=7>Taylor et al Ann Rev Biochem 1990 59:971-1005</a> 
 
  % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Fri Feb 11 14:53:04 2005

% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
tspan=[0, 1000];
% Molecule names
molnames = {
'1 /kinetics/PKA/R2C2'
'2 /kinetics/PKA/PKA-inhibitor'
'3 /kinetics/PKA/inhibited-PKA'
'4 /kinetics/PKA/cAMP.R2C2'
'5 /kinetics/PKA/cAMP2.R2C2'
'6 /kinetics/PKA/cAMP3.R2C2'
'7 /kinetics/PKA/cAMP4.R2C2'
'8 /kinetics/PKA/cAMP4.R2C'
'9 /kinetics/PKA/cAMP4.R2'
'10 /kinetics/cAMP'
'11 /kinetics/PKA-active'
};

bufnames = {
};

sumtotnames = {
};

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

hold on;
plot(t,y(:,11));
plot(t,y(:,3));
plot(t,y(:,2));
plot(t,y(:,10));
% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 - 54 * y(1) * y(10) + 33 * y(4)
 - 60 * y(11) * y(2) + 1 * y(3)
 + 60 * y(11) * y(2) - 1 * y(3)
 + 54 * y(1) * y(10) - 33 * y(4) - 54 * y(4) * y(10) + 33 * y(5)
 + 54 * y(4) * y(10) - 33 * y(5) - 75 * y(5) * y(10) + 110 * y(6)
 + 75 * y(5) * y(10) - 110 * y(6) - 75 * y(10) * y(6) + 32.5 * y(7)
 + 75 * y(10) * y(6) - 32.5 * y(7) - 60 * y(7) + 18 * y(11) * y(8)
 + 60 * y(7) - 18 * y(11) * y(8) - 60 * y(8) + 18 * y(11) * y(9)
 + 60 * y(8) - 18 * y(11) * y(9)
 - 54 * y(1) * y(10) + 33 * y(4) - 54 * y(4) * y(10) + 33 * y(5) - 75 * y(5) * y(10) + 110 * y(6) - 75 * y(10) * y(6) + 32.5 * y(7)
 + 60 * y(7) - 18 * y(11) * y(8) + 60 * y(8) - 18 * y(11) * y(9) - 60 * y(11) * y(2) + 1 * y(3)
];

% end of model
