 %  DOQCS : http://doqcs.ncbs.res.in/ 
 %  Accession Name = CondorelliGC 
 %  Accession Number = 19 
 %  Transcriber = Sudhir Sivakumaran, NCBS 
 %  Developer = Condorelli P and George SC. 
 %  Species = Generic mammalian 
 %  Tissue = Pulmonary and vascular smooth muscle 
 %  Cell Compartment = Cytosol 
 %  Notes = This model features the observations of <a href=http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=pubmed&dopt=Abstract&list_uids=11325714>Condorelli P, George SC. Biophys J. (2001) 80(5):2110-9</a>. They propose a in vivo mechanism for the binding of free NO to sGC based on in vitro data within both pulmonary and vascular smooth muscle. All the parameters used have been obtained from reported results based on experiments. One important conclusion they make is that control of sGC in vivo is most likely to be ultra-sensitive and that activation in vivo occurs at lower NO concentrations than previously reported. 
 
   % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Tue Mar  1 14:31:40 2005

% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
tspan=[0, 1000];
% Molecule names
molnames = {
'1 /kinetics/condorelliGCmod/NO_sGCpart_act'
'2 /kinetics/condorelliGCmod/NO_sGCfull_act'
'3 /kinetics/condorelliGCmod/NO_sGCfull_act/full_act_sGC/full_act_sGC_cplx'
'4 /kinetics/condorelliGCmod/GTP'
'5 /kinetics/condorelliGCmod/cGMP'
'6 /kinetics/condorelliGCmod/NO_a'
'7 /kinetics/condorelliGCmod/sGC_basal'
'8 /kinetics/condorelliGCmod/NO'
};

bufnames = {
};

sumtotnames = {
};

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

hold on;
plot(t,y(:,1));
plot(t,y(:,2));
plot(t,y(:,7));
plot(t,y(:,5));
plot(t,y(:,6));
% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 + 980 * y(7) * y(8) - 350 * y(1) - 0.0609 * y(1) + 0.014 * y(2) - 1.68 * y(1) * y(6)
 + 0.0609 * y(1) - 0.014 * y(2) + 1.68 * y(1) * y(6) - 3.4483 * y(4) * y(2) + 160 * y(3) + 40 * y(3) - 0.07 * y(2)
 + 3.4483 * y(4) * y(2) - 160 * y(3) - 40 * y(3)
 - 3.4483 * y(4) * y(2) + 160 * y(3)
 + 40 * y(3)
 - 1.68 * y(1) * y(6)
 - 980 * y(7) * y(8) + 350 * y(1) + 0.07 * y(2)
 - 980 * y(7) * y(8) + 350 * y(1) + 0.07 * y(2)
];

% end of model
