Location: BG_Ca_leak @ 3027d789e7b5 / matlab_parameter_fitting / K_kappa_ABmassAction.m

Author:

Shelley Fong <s.fong@auckland.ac.nz>

Date:

2022-04-28 11:07:01+12:00

Desc:

Updating with Kernik stuff

Permanent Source URI:

https://models.physiomeproject.org/workspace/835/rawfile/3027d789e7b520983200af71015fd49f6c3851eb/matlab_parameter_fitting/K_kappa_ABmassAction.m

% finding K, kappa for simple mass-action reaction of Ca leaking from SR
% into cytosol.
% Not gated, and not voltage dependent.


W_i = 3.71; %25.8;
W_SR = W_i * 0.06; %W_i * 0.088235;

%% stoichiometric matrices
% include kf and kr in matrices
N_f = [1; 0]; % [Ca_SR, Cai]
N_r = [0; 1];

N_fT = transpose(N_f);
N_rT = transpose(N_r);
N = N_r - N_f;

num_cols = size(N,2); % number of reactions
I = eye(num_cols); 

M = [I N_fT; I N_rT];

% k+ = k- = K_leak (on ClancyRudy original model)
K_leak = 1.6e-6; %0.5833333;
kf = [K_leak]; 
kr = [K_leak];

k = [kf;kr];
W = [1;W_SR;W_i];

% solve lambda
lambdaW = exp(pinv(M)*log(k));
lambda = lambdaW./W;
kappa = lambda(1:size(N,2));
K = lambda(size(N,2)+1:end);

save(['K_kappa.mat'],'kappa','K');

% Checks for kappa and K
N_rref = rref(N);
R_mat = null(N,'r');

K_eq = kf./kr;
zero_est = R_mat'*K_eq;

k_est = exp(M*log(lambdaW));
diff = sum(abs((k-k_est)./k));

%% display text for cellml
rname = {'leak'};
Kname = {'CaSR','Cai'};
for ik = 1:length(kappa)
    fprintf('var kappa_%s: fmol_per_sec {init: %g, pub: out};\n', rname{ik}, kappa(ik));
end
for ik = 1:length(K)
    fprintf('var K_%s: per_fmol {init: %g, pub: out};\n', Kname{ik}, K(ik));
end
if true
    % component import parameters
    for ik = 1:length(kappa)
        fprintf('var kappa_%s: fmol_per_sec {pub: in};\n', rname{ik});
    end
    for ik = 1:length(K)
        fprintf('var K_%s: per_fmol {pub: in};\n', Kname{ik});
    end
    
    % mapping
    for ik = 1:length(kappa)
        fprintf('vars kappa_%s and kappa_%s;\n', rname{ik}, rname{ik});
    end
    for ik = 1:length(K)
        fprintf('vars K_%s and K_%s;\n', Kname{ik}, Kname{ik});
    end    
end