Location: Hodgkin & Huxley (1952) model @ bd02dedb491f / sed-ml / scripts_BG / BG_paras_Na.m

Author:

WeiweiAi <wai484@aucklanduni.ac.nz>

Date:

2021-11-24 17:49:12+13:00

Desc:

Add BG params related scripts and matrix

Permanent Source URI:

http://models.cellml.org/workspace/64f/rawfile/bd02dedb491faf9960ccd23496c19ccd3e339db1/sed-ml/scripts_BG/BG_paras_Na.m

clear
%% Define volumes (unit pL)
W_i = 38;
W_e = 5.182;

opts = detectImportOptions('Na_forward_matrix.xlsx');
preview('Na_forward_matrix.xlsx',opts)
opts.SelectedVariableNames = [2:12]; 
opts.DataRange = '2:11';
N_f = readmatrix('Na_forward_matrix.xlsx',opts);
opts.SelectedVariableNames = 'species';
opts.DataRange = '2:11';
Kname=readmatrix('Na_forward_matrix.xlsx',opts);
opts.SelectedVariableNames = [2:12]; 
kappa_name=opts.SelectedVariableNames';

opts = detectImportOptions('Na_reverse_matrix.xlsx');
opts.SelectedVariableNames = [2:12]; 
opts.DataRange = '2:11';
N_r = readmatrix('Na_reverse_matrix.xlsx',opts);

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

N = N_r - N_f;
N_T = N_rT - N_fT;

num_cols = size(N,2);
I = eye(num_cols);
zerofill = zeros(1,num_cols);

N_c=[1; -1; 0;0;0;0;0;0;0;0;];
K_c=[W_e/W_i;];
N_c_T=transpose(N_c);

M = [I N_fT; I N_rT;zerofill N_c_T;];

%% Import gate transition parameters
% Na channel
load(['Na_GHK_paras.mat']);
load(['Na_m_parameters.mat']);
alpha_m0_bg = params_vec(1)*1e3; % unit s^-1
beta_m0_bg = params_vec(3)*1e3; % unit s^-1
load(['Na_h_parameters.mat']);
alpha_h0_bg = params_vec(1)*1e3; % unit s^-1
beta_h0_bg = params_vec(3)*1e3; % unit s^-1

kf = [kf_GHK; ... % R_GHK
    3*alpha_m0_bg; ... % Rm10
    2*alpha_m0_bg; ... % Rm20
    alpha_m0_bg; ... % Rm30
    3*alpha_m0_bg; ... % Rm11
    2*alpha_m0_bg; ... % Rm21
    alpha_m0_bg; ... % Rm31    
    alpha_h0_bg; ... % Rh0
    alpha_h0_bg; ... % Rh1
    alpha_h0_bg; ... % Rh2
    alpha_h0_bg; ... % Rh3
    ]; 
    
kr = [kr_GHK; ... % R_GHK
    1*beta_m0_bg; ... % Rm10
    2*beta_m0_bg; ... % Rm20
    3*beta_m0_bg; ... % Rm30
    1*beta_m0_bg; ... % Rm11
    2*beta_m0_bg; ... % Rm21
    3*beta_m0_bg; ... % Rm31    
    beta_h0_bg; ... % Rh0
    beta_h0_bg; ... % Rh1
    beta_h0_bg; ... % Rh2
    beta_h0_bg; ... % Rh3
    ];

k = [kf;kr;K_c];
W = [ones(size(N,2),1);W_i;W_e;ones(8,1);];

lambdaW = exp(pinv(M)*log(k));
lambda = lambdaW./W;
kappa = lambda(1:size(N,2));
K = lambda(size(N,2)+1:end);
save(['BG_paras_Na.mat'],'kappa','K','Kname','kappa_name'); %

fileID = fopen('BG_paras_Na_CellML.txt','w');
var_kappa=num2cell(kappa);
data_kappa=[kappa_name,var_kappa]';
var_K=num2cell(K);
data_K=[Kname,var_K]';
fprintf(fileID,'var kappa_%s: fmol_per_sec {init: %g};\n',data_kappa{:});
fprintf(fileID,'var K_%s: per_fmol {init: %g};\n',data_K{:});
fclose(fileID);
%% Checks
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));