- Author:
- WeiweiAi <wai484@aucklanduni.ac.nz>
- Date:
- 2022-08-15 10:06:08+12:00
- Desc:
- Add Ce_BG for chemical species;
Change C_BG to model the electrical capacitor;
Add the first version of build BG script based on stoichiometry matrices.
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/6bc/rawfile/a1dd1872c24d859834e27b7193a2317e7a4af835/BG_fit/GHK_cal.m
function [P,kf,kr,X]=GHK_cal(T,R,F,z,cXi,cXe,nX,W_i,W_e,G)
%cXi/cXe intracelluar/extracelluar concentration, mM
%W_i/W_e intracelluar/extracelluar volume, pL
%nX, channel numbers
%G, conductance of the open channel, mS
% Define constants
N_A = 6.022e23;
%Nernst potential
E_X = R*T/F*log(cXe/cXi);
E_X_norm = z*F*E_X/(R*T);
G_GHK = 2*G*(1-exp(E_X_norm))/(cXi - cXe*exp(E_X_norm))*R*T/F/(z^2); % Unit mA/mM
P = G_GHK/F * 1e12; % Unit pL/s
X=nX/N_A*1e15; % unit fmol
kf=P/(W_i*X); % s^-1 fmol^1
kr=P/(W_e*X); % s^-1 fmol^1
