 %  DOQCS : http://doqcs.ncbs.res.in/ 
 %  Accession Name = NOS_Iwanaga 
 %  Accession Number = 11 
 %  Transcriber = Sudhir Sivakumaran, NCBS 
 %  Developer = Iwanaga, Yamazaki and Kominami 
 %  Species = Rat; Mammalian 
 %  Tissue = Neuronal; Expression in E.coli 
 %  Cell Compartment = Cytosol 
 %  Notes = This model is based on <a href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=pubmed&dopt=Abstract&list_uids=10600125">Iwanaga T et al. Biochemistry (1999) 38(50):16629-35</a>.The model features the catalytic activity of neuronal Nitric Oxide Synthase expressed in Escherichia coli in a four step reaction mechanism. The rates found in the model are used directly from their paper and have been determined from experiments performed at 25 C. 
 
   % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Tue Mar  1 13:48:24 2005

% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
tspan=[0, 400];
% Molecule names
molnames = {
'1 /kinetics/neuroNOS/nNOS'
'2 /kinetics/neuroNOS/Ca-CaM_nNOS'
'3 /kinetics/neuroNOS/L-Arg'
'4 /kinetics/neuroNOS/OH.Arg'
'5 /kinetics/neuroNOS/nNOS.Cit+NO'
'6 /kinetics/neuroNOS/Cit'
'7 /kinetics/neuroNOS/NO'
'8 /kinetics/neuroNOS/nNOS.OH.Arg'
'9 /kinetics/neuroNOS/nNOS.arg'
'10 /kinetics/neuroNOS/CaM-Ca4'
};

bufnames = {
};

sumtotnames = {
};

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

hold on;
 plot(t,y(:,5));
 plot(t,y(:,9));
 plot(t,y(:,10));
 plot(t,y(:,3));
 plot(t,y(:,3));
 plot(t,y(:,8));
 plot(t,y(:,7));
% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 - 3.25 * y(1) * y(10) + 0.05 * y(2)
 + 3.25 * y(1) * y(10) - 0.05 * y(2) - 0.5 * y(3) * y(2) + 81 * y(9)
 - 0.5 * y(3) * y(2) + 81 * y(9)
 + 4.4 * y(8)
 + 19 * y(8) - 2.7 * y(5)
 + 2.7 * y(5)
 + 2.7 * y(5)
 - 4.4 * y(8) - 19 * y(8) + 25 * y(9)
 + 0.5 * y(3) * y(2) - 81 * y(9) - 25 * y(9)
 - 3.25 * y(1) * y(10) + 0.05 * y(2)
];

% end of model
