 %  DOQCS : http://doqcs.ncbs.res.in/ 
 %  Accession Name = CDK2-CyclinA_activation_improved 
 %  Accession Number = 43 
 %  Transcriber =  Sharat J. Vayttaden, NCBS 
 %  Developer = May C. Morris, Claire Gondeau, John A. Tainer and Gilles Divita 
 %  Species = Human, Xenopus 
 %  Tissue = Expressed in E.coli 
 %  Cell Compartment = Cytosol 
 %  Notes = This version of the model for CDK2/Cyclin A complex activation (expressed in E. coli) is a good match to the data and the simulation matches the graphs Fig.1B and Fig.1B inset. The rate constants though differ from the rate constants in <a href = http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=pubmed&dopt=Abstract&list_uids=11959850>Morris MC. et al. J Biol Chem. (2002) 277(26):23847-53</a>. 
 
   % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Wed Feb  2 16:29:40 2005

% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
tspan=[0, 12];
% Molecule names
molnames = {
'1 /kinetics/CDK2cycA'
'2 /kinetics/CyclinA'
'3 /kinetics/Cdk2'
'4 /kinetics/CDK2cycA*'
'5 /kinetics/basal_fluorescence'
};

bufnames = {
};

sumtotnames = {
'1 /kinetics/total_fluorescence'
};

% Initial concs of 5 molecules
y0(1) = 0;
y0(2) = 0.400007;
y0(3) = 0.1;
y0(4) = 0;
y0(5) = 1.21211;
% Sumtotal eqn for /kinetics/total_fluorescence: sumtot(1) =  (y(1) + y(4) + y(5) )
% Call the ODE solver and display results
[t, y] = ode23s(@f, tspan, y0);

hold on;
plot(t, (y(:,1) + y(:,4) + y(:,5) ));
% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 - 1.1 * y(1) + 0.6 * y(4) + 33 * y(3) * y(2) - 50 * y(1)
 - 33 * y(3) * y(2) + 50 * y(1)
 - 33 * y(3) * y(2) + 50 * y(1)
 + 1.1 * y(1) - 0.6 * y(4)
0
];

% end of model
