- Author:
- nima <nafs080@aucklanduni.ac.nz>
- Date:
- 2022-04-05 19:17:27+12:00
- Desc:
- CellML and SEDML Files
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/841/rawfile/2fed4cbe51bfb5b3c2a6b10546351fd3c2daf026/Figure02.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="Nima_2018" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#">
<!-- electroneutrality maintained also apical and basolateral concentration considered as constant without any equations
Basolateral concentration are constant and they are same as high glucose DMEM cell culture
No conditon for glucose concentration in the lumen
Output of this model then fitted to experimental data
water transport added to the curent model-->
<units name="per_second">
<unit exponent="-1" units="second"/>
</units>
<units name="per_second2">
<unit exponent="-1" units="second"/>
</units>
<units name="Mm2">
<unit units="M"/>
<unit units="m2"/>
</units>
<!-- def unit per_second5 as
unit second {expo: -5};
enddef;-->
<units name="volt_per_second">
<unit units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_meter">
<unit exponent="-1" units="meter"/>
</units>
<units name="per_volt_per_second">
<unit exponent="-1" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt4_per_second">
<unit exponent="-4" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt3_per_second">
<unit exponent="-3" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt2_per_second">
<unit exponent="-2" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_mol">
<unit exponent="-1" units="mole"/>
</units>
<units name="umol">
<unit prefix="micro" units="mole"/>
</units>
<units name="per_umol">
<unit exponent="-1" units="umol"/>
</units>
<units name="M">
<unit units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="per_M">
<unit exponent="-1" units="M"/>
</units>
<units name="m2">
<unit exponent="2" units="metre"/>
</units>
<units name="cm">
<unit prefix="centi" units="metre"/>
</units>
<units name="um2">
<unit exponent="2" prefix="micro" units="metre"/>
</units>
<units name="um">
<unit prefix="micro" units="meter"/>
</units>
<units name="um3">
<unit exponent="3" prefix="micro" units="metre"/>
</units>
<units name="m3">
<unit exponent="3" units="metre"/>
</units>
<units name="m3_per_mole">
<unit units="m3"/>
<unit exponent="-1" units="mole"/>
</units>
<units name="m3_per_second">
<unit units="m3"/>
<unit units="second"/>
</units>
<!-- def unit litre as
unit m3 {pref: milli};
enddef;-->
<units name="mol_per_um2">
<unit units="mole"/>
<unit exponent="-1" units="um2"/>
</units>
<units name="mM">
<unit prefix="milli" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="mV">
<unit prefix="milli" units="volt"/>
</units>
<units name="J_per_K_per_mol">
<unit exponent="-1" units="mole"/>
<unit units="joule"/>
<unit exponent="-1" units="kelvin"/>
</units>
<units name="pA_per_cm2">
<unit prefix="pico" units="ampere"/>
<unit exponent="-2" prefix="centi" units="metre"/>
</units>
<units name="uA">
<unit prefix="micro" units="ampere"/>
</units>
<units name="uA_per_um2">
<unit units="uA"/>
<unit exponent="-1" units="um2"/>
</units>
<units name="uA_per_A">
<unit units="uA"/>
<unit exponent="-1" units="ampere"/>
</units>
<units name="umol_per_mol">
<unit units="umol"/>
<unit exponent="-1" units="mole"/>
</units>
<units name="C_per_mol">
<unit exponent="-1" units="mole"/>
<unit units="coulomb"/>
</units>
<units name="per_M_per_second">
<unit exponent="-1" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M2_per_second">
<unit exponent="-2" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M4_per_second2">
<unit exponent="-4" units="M"/>
<unit exponent="-2" units="second"/>
</units>
<units name="per_M4_per_second">
<unit exponent="-4" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="uM">
<unit prefix="micro" units="M"/>
</units>
<units name="M2">
<unit exponent="2" units="M"/>
</units>
<units name="M3">
<unit exponent="3" units="M"/>
</units>
<units name="M3_per_second">
<unit exponent="3" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M3_per_second5">
<unit exponent="-3" units="M"/>
<unit exponent="-5" units="second"/>
</units>
<units name="M_per_second">
<unit units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="cm2">
<unit exponent="2" prefix="centi" units="metre"/>
</units>
<units name="cm2_per_s">
<unit units="cm2"/>
<unit exponent="-1" units="second"/>
</units>
<units name="mole_per_m2">
<unit units="mole"/>
<unit exponent="-1" units="m2"/>
</units>
<!-- def unit uM_per_cm2_per_s as
unit uM;
unit cm2 {expo: -1};
unit second {expo: -1};
enddef;-->
<units name="Farad">
<unit exponent="-1" units="volt"/>
<unit units="coulomb"/>
</units>
<units name="uF">
<unit prefix="micro" units="Farad"/>
</units>
<units name="uF_per_cm2">
<unit prefix="micro" units="Farad"/>
<unit exponent="-1" units="cm2"/>
</units>
<units name="umol_per_cm2_per_s">
<unit prefix="micro" units="mole"/>
<unit exponent="-1" units="cm2"/>
<unit exponent="-1" units="second"/>
</units>
<units name="umol_per_s">
<unit prefix="micro" units="mole"/>
<unit exponent="-1" units="second"/>
</units>
<units name="uA_per_cm2">
<unit units="uA"/>
<unit exponent="-2" prefix="centi" units="metre"/>
</units>
<units name="m_per_s">
<unit units="meter"/>
<unit exponent="-1" units="second"/>
</units>
<units name="m3_mole_per_litre_per_umol">
<unit units="m3"/>
<unit units="mole"/>
<unit exponent="-1" units="litre"/>
<unit exponent="-1" units="umol"/>
</units>
<units name="litre_umol_per_mol_per_m3">
<unit units="litre"/>
<unit units="umol"/>
<unit exponent="-1" units="mole"/>
<unit exponent="-1" units="m3"/>
</units>
<units name="S">
<unit units="siemens"/>
</units>
<units name="S_per_m2">
<unit units="siemens"/>
<unit exponent="-1" units="m2"/>
</units>
<units name="per_m2">
<unit exponent="-1" units="m2"/>
</units>
<units name="cm_per_s">
<unit units="cm"/>
<unit exponent="-1" units="second"/>
</units>
<component name="parameters">
<variable initial_value="-1" name="z_Cl" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="z_Na" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="z_k" public_interface="out" units="dimensionless"/>
<variable initial_value="96485" name="F" public_interface="out" units="C_per_mol"/>
<variable initial_value="8.314" name="R" public_interface="out" units="J_per_K_per_mol"/>
<variable initial_value="310" name="T" public_interface="out" units="kelvin"/>
<variable initial_value="1e-05" name="capacitance" public_interface="out" units="uF"/>
<variable initial_value="7e-10" name="A_membrane" public_interface="out" units="m2"/>
<!--A_membrane = v_cell/2.5e-6{metre};-->
<variable name="v_cell" public_interface="in" units="m3"/>
<variable initial_value="0" name="P_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_k" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Cl" public_interface="out" units="m_per_s"/>
<variable initial_value="1e-08" name="P_P_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="9e-07" name="P_P_k" public_interface="out" units="m_per_s"/>
<variable initial_value="1e-08" name="P_P_Cl" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_k" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_Cl" public_interface="out" units="m_per_s"/>
<variable name="time" public_interface="out" units="second"/>
<variable initial_value="14000" name="k0_12" public_interface="out" units="per_M2_per_second"/>
<variable initial_value="300" name="k0_21" public_interface="out" units="per_second"/>
<variable initial_value="15000" name="k0_23" public_interface="out" units="per_M_per_second"/>
<variable initial_value="20" name="k0_32" public_interface="out" units="per_second"/>
<variable initial_value="0.2" name="k0_25" public_interface="out" units="per_second"/>
<variable initial_value="0.01" name="k0_52" public_interface="out" units="per_second"/>
<variable initial_value="50" name="k0_34" public_interface="out" units="per_second"/>
<variable initial_value="50" name="k0_43" public_interface="out" units="per_second"/>
<variable initial_value="6" name="k0_45" public_interface="out" units="per_second"/>
<variable initial_value="450" name="k0_54" public_interface="out" units="per_M_per_second"/>
<variable initial_value="10" name="k0_56" public_interface="out" units="per_second"/>
<variable initial_value="200" name="k0_65" public_interface="out" units="per_M2_per_second"/>
<variable initial_value="25" name="k0_61" public_interface="out" units="per_second"/>
<variable initial_value="600" name="k0_16" public_interface="out" units="per_second"/>
<variable initial_value="0.7" name="delta" public_interface="out" units="dimensionless"/>
<variable initial_value="0.3" name="alpha_p" public_interface="out" units="dimensionless"/>
<variable initial_value="0" name="alpha_pp" public_interface="out" units="dimensionless"/>
<!--ar N_C: dimensionless {init: 6e10};
var N_Avo: per_mol {init: 6.022e23};-->
<!--var C_T: umol {pub: out};-->
<variable initial_value="2" name="n" public_interface="out" units="dimensionless"/>
<variable initial_value="-2" name="z_c" public_interface="out" units="dimensionless"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="mu" public_interface="out" units="dimensionless"/>
<!--alpha_p = 1{dimensionless}-delta-alpha_pp;-->
<!--C_T = 1e6{umol_per_mol}*N_C/N_Avo;-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Apical_concentrations">
<variable initial_value="0.15" name="Na_m" public_interface="out" units="M"/>
<variable initial_value="0.0054" name="K_m" public_interface="out" units="M"/>
<variable name="glucose_m" public_interface="out" units="M"/>
<variable initial_value="0.15" name="Cl_m" public_interface="out" units="M"/>
<variable name="time" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>glucose_m</ci>
<piecewise>
<piece>
<cn cellml:units="M">0</cn>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<cn cellml:units="second">0</cn>
</apply>
<apply>
<lt/>
<ci>time</ci>
<cn cellml:units="second">5000</cn>
</apply>
</apply>
</piece>
<piece>
<cn cellml:units="M">0.05</cn>
<apply>
<geq/>
<ci>time</ci>
<cn cellml:units="second">5000</cn>
</apply>
</piece>
</piecewise>
</apply>
</math>
</component>
<component name="Apical_voltage">
<variable initial_value="-0.046" name="v_mc" public_interface="out" units="volt"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="capacitance" public_interface="in" units="uF"/>
<!-- var I_D_Na: uA {pub: in};
var I_D_K: uA {pub: in};
var I_D_Cl: uA {pub: in};-->
<variable name="I_D_P_Na" public_interface="in" units="uA"/>
<variable name="I_D_P_K" public_interface="in" units="uA"/>
<variable name="I_D_P_Cl" public_interface="in" units="uA"/>
<variable name="I_NHE3_H" public_interface="in" units="uA"/>
<variable name="I_NHE3_Na" public_interface="in" units="uA"/>
<variable name="I_BK" public_interface="in" units="uA"/>
<!-- var I_AE1_HCO3: uA {pub: in};
var I_AE1_Cl: uA {pub: in};-->
<variable name="I_ENaC" public_interface="in" units="uA"/>
<!--var I_NaCl: uA {pub: in};-->
<variable name="I_Na_SGLT" public_interface="in" units="uA"/>
<variable name="I_CFTR" public_interface="in" units="uA"/>
<!--var I_CFTR_HCO3: uA {pub: in};-->
<!--var I_Gl_SGLT: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_mc</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_Na_SGLT</ci>
</apply>
<ci>I_BK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_CFTR</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_ENaC</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<plus/>
<ci>I_D_P_Na</ci>
</apply>
<ci>I_D_P_K</ci>
<ci>I_D_P_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<ci>capacitance</ci>
</apply>
</apply>
</math>
</component>
<component name="Cell_concentration">
<variable initial_value="0.045" name="Na_i" public_interface="out" units="M"/>
<!--var Na_m: M {pub: in};
var Na_s: M {pub: in};-->
<variable initial_value="0.009" name="glucose_i" public_interface="out" units="M"/>
<variable initial_value="0.105" name="K_i" public_interface="out" units="M"/>
<variable initial_value="6e-05" name="Ca_i" public_interface="out" units="M"/>
<variable initial_value="0.057" name="Cl_i" public_interface="out" units="M"/>
<variable name="HCO3_i" public_interface="out" units="M"/>
<variable name="H_int" public_interface="out" units="M"/>
<variable name="HCO3_m" public_interface="out" units="M"/>
<variable name="HCO3_s" public_interface="out" units="M"/>
<variable initial_value="6.3e-08" name="H_Ext" public_interface="out" units="M"/>
<variable initial_value="6.3e-08" name="H_s" public_interface="out" units="M"/>
<variable initial_value="2e-15" name="v_cell" public_interface="out" units="m3"/>
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="J_Na_SGLT" public_interface="in" units="umol_per_s"/>
<variable name="J_Gl_SGLT" public_interface="in" units="umol_per_s"/>
<variable name="J_Na_NaK" public_interface="in" units="umol_per_s"/>
<variable name="J_K_NaK" public_interface="in" units="umol_per_s"/>
<variable name="J_NBC_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_NBC_HCO3" public_interface="in" units="umol_per_s"/>
<!--var J_NaCl: umol_per_s {pub: in};-->
<variable name="J_NHE3_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_NHE3_H" public_interface="in" units="umol_per_s"/>
<!-- var J_D_Na: umol_per_s {pub: in};
var J_D_Bl_Na: umol_per_s {pub: in};-->
<variable name="J_GLUT" public_interface="in" units="umol_per_s"/>
<variable name="J_A_GLUT" public_interface="in" units="umol_per_s"/>
<!-- var J_D_K: umol_per_s {pub: in};
var J_D_Bl_K: umol_per_s {pub: in};
var J_D_Cl: umol_per_s {pub: in};
var J_D_Bl_Cl: umol_per_s {pub: in};-->
<variable name="J_AE1_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_AE1_HCO3" public_interface="in" units="umol_per_s"/>
<variable name="J_CFTR" public_interface="in" units="umol_per_s"/>
<variable name="J_CLC2" public_interface="in" units="umol_per_s"/>
<!--var J_CFTR_HCO3: umol_per_s {pub: in};-->
<variable name="J_ENaC" public_interface="in" units="umol_per_s"/>
<variable name="J_NaKCC" public_interface="in" units="umol_per_s"/>
<variable name="J_IK" public_interface="in" units="umol_per_s"/>
<variable name="J_BK" public_interface="in" units="umol_per_s"/>
<variable name="J_D_P_Na" public_interface="in" units="umol_per_s"/>
<!-- var J_buf_C: M_per_second {pub: in};
var J_buf_A: M_per_second {pub: in};-->
<variable name="time" public_interface="in" units="second"/>
<variable initial_value="7.2" name="pH_int" public_interface="out" units="per_second"/>
<variable initial_value="7.2" name="pH_Ext" public_interface="out" units="per_second"/>
<variable initial_value="7.2" name="pH_s" public_interface="out" units="per_second"/>
<variable initial_value="1680000" name="rho" public_interface="out" units="per_meter"/>
<variable name="beta" public_interface="out" units="M"/>
<variable initial_value="6.4" name="p_K" public_interface="out" units="dimensionless"/>
<variable initial_value="40" name="P_CO2" public_interface="out" units="dimensionless"/>
<variable initial_value="1e-05" name="s" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="theta_1" units="dimensionless"/>
<variable initial_value="1" name="theta_2" units="dimensionless"/>
<variable initial_value="1" name="theta_3" units="dimensionless"/>
<variable initial_value="1" name="theta_4" units="dimensionless"/>
<variable initial_value="1" name="theta_5" units="dimensionless"/>
<variable initial_value="1" name="theta_6" units="dimensionless"/>
<variable initial_value="1" name="theta_7" units="dimensionless"/>
<variable initial_value="1" name="theta_8" units="dimensionless"/>
<variable initial_value="1" name="theta_9" units="dimensionless"/>
<variable initial_value="1" name="theta_10" units="dimensionless"/>
<variable initial_value="1" name="theta_11" units="dimensionless"/>
<variable initial_value="1" name="theta_12" units="dimensionless"/>
<variable initial_value="1" name="theta_13" units="dimensionless"/>
<variable initial_value="1" name="theta_14" units="dimensionless"/>
<variable initial_value="1" name="theta_15" units="dimensionless"/>
<variable initial_value="1" name="theta_16" units="dimensionless"/>
<variable initial_value="1" name="theta_17" units="dimensionless"/>
<variable initial_value="1" name="theta_18" units="dimensionless"/>
<variable initial_value="1" name="theta_19" units="dimensionless"/>
<variable initial_value="1" name="theta_20" units="dimensionless"/>
<variable initial_value="1" name="theta_21" units="dimensionless"/>
<variable initial_value="1" name="theta_22" units="dimensionless"/>
<variable initial_value="1" name="theta_23" units="dimensionless"/>
<variable initial_value="1" name="theta_24" units="dimensionless"/>
<variable initial_value="1" name="theta_25" units="dimensionless"/>
<variable initial_value="1" name="theta_26" units="dimensionless"/>
<variable initial_value="1" name="theta_27" units="dimensionless"/>
<variable initial_value="1.8e-05" name="v_w" public_interface="out" units="m3_per_mole"/>
<variable name="J_w_A" public_interface="out" units="umol_per_s"/>
<variable name="J_w_B" public_interface="out" units="umol_per_s"/>
<variable name="J_w_P" public_interface="out" units="umol_per_s"/>
<variable name="J_D_P_Gl" public_interface="out" units="umol_per_s"/>
<variable initial_value="6e-06" name="L_A" public_interface="out" units="m_per_s"/>
<variable initial_value="6e-06" name="L_B" public_interface="out" units="m_per_s"/>
<variable initial_value="4e-07" name="L_P" public_interface="out" units="m_per_s"/>
<variable initial_value="3e-10" name="L_P_Gl" public_interface="out" units="m_per_s"/>
<variable name="NII" public_interface="out" units="M"/>
<variable name="chi_i" public_interface="out" units="M"/>
<variable initial_value="0.001" name="psi" public_interface="out" units="M"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="glucose_s" public_interface="in" units="M"/>
<variable name="Cl_s" public_interface="in" units="M"/>
<!--var dif: dimensionless {pub: out};-->
<!-- var alpha: dimensionless {init: 0.1, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>NII</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>-3</cn>
<ci>chi_i</ci>
</apply>
<ci>v_cell</ci>
</apply>
</apply>
</math>
<!--chi_i = 1e-3{dimensionless}*(NII * v_cell);-->
<variable initial_value="1" name="alpha_1" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_2" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_3" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_4" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_5" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_6" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_7" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_8" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_9" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_10" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_11" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_12" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_13" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_14" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="alpha_15" public_interface="out" units="dimensionless"/>
<variable name="OS_A" public_interface="out" units="M"/>
<variable name="OS_C" public_interface="out" units="M"/>
<variable name="OS_B" public_interface="out" units="M"/>
<variable name="A" public_interface="out" units="M"/>
<variable name="B" public_interface="out" units="M"/>
<variable initial_value="0" name="Gl_net_P" public_interface="out" units="M"/>
<variable initial_value="0" name="Gl_net_A" public_interface="out" units="M"/>
<variable initial_value="0" name="Gl_net_B" public_interface="out" units="M"/>
<variable name="Net_gl" public_interface="out" units="M"/>
<variable initial_value="0" name="Cell_electroneutrality" public_interface="out" units="M"/>
<variable initial_value="0" name="Blood_electroneutrality" public_interface="out" units="M"/>
<variable initial_value="0" name="Lumen_electroneutrality" public_interface="out" units="M"/>
<variable name="J_A_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_B_Na" public_interface="out" units="umol_per_s"/>
<!--var Concentration_Na: M {pub: out};-->
<!-- Concentration_Na = Na_m+Na_i+Na_s;-->
<!--var K_g: per_second {init: 0.0015, pub: out};-->
<!--dif = v_w*(J_w_B-J_w_A);-->
<!--ode(dif,time) = 150{dimensionless};-->
<!--var K_g: per_second {init: 0.0015, pub: out};-->
<!-- Lumen_electroneutrality = -(Cl_m+HCO3_m)+Na_m+K_m;
Cell_electroneutrality = -(Cl_i+HCO3_i+NII)+Na_i+K_i;-->
<!-- Blood_electroneutrality = -(Cl_s+HCO3_s)+Na_s+K_s;-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>HCO3_s</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>Na_s</ci>
<ci>K_s</ci>
</apply>
<ci>Blood_electroneutrality</ci>
</apply>
<ci>Cl_s</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>HCO3_m</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>Na_m</ci>
<ci>K_m</ci>
</apply>
<ci>Lumen_electroneutrality</ci>
</apply>
<ci>Cl_m</ci>
</apply>
</apply>
<!--NII = Na_i+K_i-Cell_electroneutrality-Cl_i-HCO3_i;-->
<apply>
<eq/>
<ci>J_A_Na</ci>
<apply>
<plus/>
<ci>J_Na_SGLT</ci>
<ci>J_ENaC</ci>
<ci>J_NHE3_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_B_Na</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>J_Na_NaK</ci>
</apply>
<ci>J_NaKCC</ci>
<ci>J_NBC_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>chi_i</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>Na_i</ci>
<ci>K_i</ci>
</apply>
<ci>Cell_electroneutrality</ci>
</apply>
<ci>Cl_i</ci>
</apply>
<ci>HCO3_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>A</ci>
<ci>OS_C</ci>
</apply>
<apply>
<eq/>
<ci>B</ci>
<apply>
<divide/>
<apply>
<plus/>
<ci>OS_A</ci>
<ci>OS_B</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>OS_A</ci>
<apply>
<plus/>
<ci>Cl_m</ci>
<ci>HCO3_m</ci>
<ci>Na_m</ci>
<ci>K_m</ci>
<ci>glucose_m</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>OS_C</ci>
<apply>
<plus/>
<ci>Cl_i</ci>
<ci>HCO3_i</ci>
<ci>Na_i</ci>
<ci>K_i</ci>
<ci>chi_i</ci>
<ci>glucose_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>OS_B</ci>
<apply>
<plus/>
<ci>Cl_s</ci>
<ci>HCO3_s</ci>
<ci>Na_s</ci>
<ci>K_s</ci>
<ci>glucose_s</ci>
</apply>
</apply>
<!--psi = 2{dimensionless}*(Cl_i+HCO3_i+Na_i+K_i+chi_c+glucose_i)-(Cl_s+HCO3_s+Na_s+K_s+glucose_s)-(Cl_m+HCO3_m+Na_m+K_m+glucose_m);-->
<apply>
<eq/>
<ci>beta</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2.3</cn>
<ci>HCO3_i</ci>
<ci>A_membrane</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>theta_1</ci>
<ci>J_Na_SGLT</ci>
</apply>
<apply>
<times/>
<ci>theta_2</ci>
<ci>J_Na_NaK</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_15</ci>
<ci>J_ENaC</ci>
</apply>
<apply>
<times/>
<ci>theta_16</ci>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_19</ci>
<ci>J_NBC_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_3</ci>
<ci>J_NHE3_Na</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_cell</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>glucose_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_26</ci>
<cn cellml:units="dimensionless">1</cn>
<ci>J_A_GLUT</ci>
</apply>
<apply>
<times/>
<ci>theta_6</ci>
<ci>J_Gl_SGLT</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_7</ci>
<cn cellml:units="dimensionless">1</cn>
<ci>J_GLUT</ci>
</apply>
</apply>
<apply>
<times/>
<ci>glucose_i</ci>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_cell</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_8</ci>
<ci>J_K_NaK</ci>
</apply>
<apply>
<times/>
<ci>theta_17</ci>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_27</ci>
<ci>J_IK</ci>
</apply>
<ci>J_BK</ci>
</apply>
<apply>
<times/>
<ci>K_i</ci>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_cell</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Cl_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<apply>
<minus/>
<ci>theta_11</ci>
</apply>
<ci>J_CFTR</ci>
</apply>
<ci>J_CLC2</ci>
</apply>
<apply>
<times/>
<ci>theta_18</ci>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_14</ci>
<ci>J_AE1_Cl</ci>
</apply>
</apply>
<apply>
<times/>
<ci>Cl_i</ci>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_cell</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<!--ode(HCO3_i, time) = 1e-9{m3_mole_per_litre_per_umol}/v_cell*theta_20*J_NBC_HCO3+theta_25*J_AE1_HCO3+theta_21*J_buf_C;
psi = (J_w_B - (ode(v_cell,time)/(1e-6{umole_per_mole}*v_w)) - 1e9{dimensionless}*L_A*A_membrane*(alpha_1*Na_m+alpha_2*Cl_m+alpha_3*glucose_m+alpha_4*K_m+alpha_5*HCO3_m-(alpha_6*Na_i+alpha_7*Cl_i+alpha_8*glucose_i+alpha_9*K_i+alpha_10*HCO3_i+chi_c)))/(1e9{dimensionless}*L_A*A_membrane);-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>pH_int</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>rho</ci>
</apply>
<ci>beta</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci>theta_20</ci>
<ci>J_NBC_HCO3</ci>
</apply>
<apply>
<times/>
<ci>theta_22</ci>
<ci>J_NHE3_H</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_25</ci>
<ci>J_AE1_HCO3</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_cell</ci>
</apply>
<apply>
<times/>
<cn cellml:units="umole_per_mole" type="e-notation">1<sep/>-6</cn>
<ci>v_w</ci>
<apply>
<minus/>
<ci>J_w_B</ci>
<ci>J_w_A</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_w_A</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<ci>L_A</ci>
<ci>A_membrane</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>Na_m</ci>
</apply>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>Cl_m</ci>
</apply>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>glucose_m</ci>
</apply>
<apply>
<times/>
<ci>alpha_4</ci>
<ci>K_m</ci>
</apply>
<apply>
<times/>
<ci>alpha_5</ci>
<ci>HCO3_m</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_6</ci>
<ci>Na_i</ci>
</apply>
<apply>
<times/>
<ci>alpha_7</ci>
<ci>Cl_i</ci>
</apply>
<apply>
<times/>
<ci>alpha_8</ci>
<ci>glucose_i</ci>
</apply>
<apply>
<times/>
<ci>alpha_9</ci>
<ci>K_i</ci>
</apply>
<apply>
<times/>
<ci>alpha_10</ci>
<ci>HCO3_i</ci>
</apply>
<ci>chi_i</ci>
</apply>
</apply>
<ci>psi</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_w_B</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<ci>L_B</ci>
<ci>A_membrane</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_6</ci>
<ci>Na_i</ci>
</apply>
<apply>
<times/>
<ci>alpha_7</ci>
<ci>Cl_i</ci>
</apply>
<apply>
<times/>
<ci>alpha_8</ci>
<ci>glucose_i</ci>
</apply>
<apply>
<times/>
<ci>alpha_9</ci>
<ci>K_i</ci>
</apply>
<apply>
<times/>
<ci>alpha_10</ci>
<ci>HCO3_i</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_11</ci>
<ci>Na_s</ci>
</apply>
<apply>
<times/>
<ci>alpha_12</ci>
<ci>Cl_s</ci>
</apply>
<apply>
<times/>
<ci>alpha_13</ci>
<ci>glucose_s</ci>
</apply>
<apply>
<times/>
<ci>alpha_14</ci>
<ci>K_s</ci>
</apply>
<apply>
<times/>
<ci>alpha_15</ci>
<ci>HCO3_s</ci>
</apply>
</apply>
</apply>
<ci>chi_i</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_w_P</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<ci>L_P</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_1</ci>
<ci>Na_m</ci>
</apply>
<apply>
<times/>
<ci>alpha_2</ci>
<ci>Cl_m</ci>
</apply>
<apply>
<times/>
<ci>alpha_3</ci>
<ci>glucose_m</ci>
</apply>
<apply>
<times/>
<ci>alpha_4</ci>
<ci>K_m</ci>
</apply>
<apply>
<times/>
<ci>alpha_5</ci>
<ci>HCO3_m</ci>
</apply>
<ci>psi</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>alpha_11</ci>
<ci>Na_s</ci>
</apply>
<apply>
<times/>
<ci>alpha_12</ci>
<ci>Cl_s</ci>
</apply>
<apply>
<times/>
<ci>alpha_13</ci>
<ci>glucose_s</ci>
</apply>
<apply>
<times/>
<ci>alpha_14</ci>
<ci>K_s</ci>
</apply>
<apply>
<times/>
<ci>alpha_15</ci>
<ci>HCO3_s</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_Gl</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<ci>L_P_Gl</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>glucose_m</ci>
<ci>glucose_s</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Gl_net_A</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_26</ci>
<ci>J_A_GLUT</ci>
</apply>
<apply>
<times/>
<ci>theta_6</ci>
<ci>J_Gl_SGLT</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>glucose_i</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>3</cn>
<ci>v_cell</ci>
</apply>
</apply>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>-6</cn>
<ci>v_w</ci>
<apply>
<minus/>
<ci>J_w_B</ci>
<ci>J_w_A</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Gl_net_B</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<minus/>
<ci>theta_7</ci>
</apply>
<ci>J_GLUT</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<ci>glucose_i</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>3</cn>
<ci>v_cell</ci>
</apply>
</apply>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>-6</cn>
<ci>v_w</ci>
<apply>
<minus/>
<ci>J_w_B</ci>
<ci>J_w_A</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Net_gl</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>J_A_GLUT</ci>
<ci>J_Gl_SGLT</ci>
</apply>
<ci>J_GLUT</ci>
</apply>
</apply>
<!--J_D_P_Gl = 1e9{dimensionless}*alpha*L_P*A_membrane*(OS_A-OS_B);
J_D_P_Gl = alpha*J_w_P;-->
<!--ode(Gl_net_P, time) = J_D_P_Gl;-->
<apply>
<eq/>
<ci>H_int</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_int</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>HCO3_i</ci>
<apply>
<times/>
<ci>s</ci>
<ci>P_CO2</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_int</ci>
<ci>p_K</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Blood_concentrations">
<variable initial_value="0.15" name="Cl_s" public_interface="out" units="M"/>
<variable initial_value="0.0054" name="K_s" public_interface="out" units="M"/>
<variable initial_value="0.15" name="Na_s" public_interface="out" units="M"/>
<variable initial_value="0.004" name="glucose_s" public_interface="out" units="M"/>
<variable initial_value="8e-17" name="Q_in" public_interface="out" units="m3_per_second"/>
<variable name="Q_out" public_interface="out" units="m3_per_second"/>
<variable initial_value="1e-16" name="v_B" public_interface="out" units="m3"/>
<variable initial_value="0.005" name="glucose_in" public_interface="out" units="M"/>
<variable initial_value="0.15" name="Cl_in" public_interface="out" units="M"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="J_GLUT" public_interface="in" units="umol_per_s"/>
<!-- var Number: dimensionless {init: 4000000000};-->
<variable initial_value="0.15" name="Na_in" public_interface="out" units="M"/>
<variable name="J_Na_NaK" public_interface="in" units="umol_per_s"/>
<variable name="J_NBC_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_NaKCC" public_interface="in" units="umol_per_s"/>
<variable name="J_D_P_Na" public_interface="in" units="umol_per_s"/>
<variable initial_value="0.0054" name="K_in" public_interface="out" units="M"/>
<!--var K_out: M {init: 0.0055, pub: out};-->
<variable name="J_K_NaK" public_interface="in" units="umol_per_s"/>
<variable name="J_D_P_K" public_interface="in" units="umol_per_s"/>
<variable name="J_IK" public_interface="in" units="umol_per_s"/>
<variable name="J_CLC2" public_interface="in" units="umol_per_s"/>
<variable name="J_D_P_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_D_P_Gl" public_interface="in" units="umol_per_s"/>
<variable initial_value="1.8e-05" name="v_w1" public_interface="out" units="m3_per_mole"/>
<variable name="J_w_B" public_interface="in" units="umol_per_s"/>
<variable name="J_w_P" public_interface="in" units="umol_per_s"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Q_out</ci>
<apply>
<minus/>
<ci>Q_in</ci>
<apply>
<times/>
<cn cellml:units="umole_per_mole" type="e-notation">1<sep/>-6</cn>
<ci>v_w1</ci>
<apply>
<plus/>
<apply>
<plus/>
<ci>J_w_P</ci>
</apply>
<ci>J_w_B</ci>
</apply>
</apply>
</apply>
</apply>
<!--ode(v_B, time) = 1e-6{umole_per_mole}*v_w1*(-J_w_P-J_w_B);-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>glucose_s</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>v_B</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>Q_in</ci>
<ci>glucose_in</ci>
</apply>
<apply>
<times/>
<ci>Q_out</ci>
<ci>glucose_s</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="mole_per_umol" type="e-notation">1<sep/>-9</cn>
<apply>
<plus/>
<ci>J_GLUT</ci>
<ci>J_D_P_Gl</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_s</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>v_B</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>Q_in</ci>
<ci>Na_in</ci>
</apply>
<apply>
<times/>
<ci>Q_out</ci>
<ci>Na_s</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="mole_per_umol" type="e-notation">1<sep/>-9</cn>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_Na_NaK</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_NBC_Na</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_NaKCC</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_D_P_Na</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_s</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>v_B</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>Q_in</ci>
<ci>K_in</ci>
</apply>
<apply>
<times/>
<ci>Q_out</ci>
<ci>K_s</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="mole_per_umol" type="e-notation">1<sep/>-9</cn>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>J_K_NaK</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_NaKCC</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_D_P_K</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_IK</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Cl_s</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>v_B</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>Q_in</ci>
<ci>Cl_in</ci>
</apply>
<apply>
<times/>
<ci>Q_out</ci>
<ci>Cl_s</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="mole_per_umol" type="e-notation">1<sep/>-9</cn>
<apply>
<minus/>
<apply>
<minus/>
<ci>J_CLC2</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NaKCC</ci>
</apply>
</apply>
<ci>J_D_P_Cl</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Basol_voltage">
<variable initial_value="-0.05" name="v_sc" public_interface="out" units="volt"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="capacitance" public_interface="in" units="uF"/>
<variable name="I_Na_NaK" public_interface="in" units="uA"/>
<variable name="I_K_NaK" public_interface="in" units="uA"/>
<variable name="I_D_P_Na" public_interface="in" units="uA"/>
<variable name="I_D_P_K" public_interface="in" units="uA"/>
<variable name="I_D_P_Cl" public_interface="in" units="uA"/>
<variable name="I_NBC_Na" public_interface="in" units="uA"/>
<variable name="I_NBC_HCO3" public_interface="in" units="uA"/>
<variable name="I_IK" public_interface="in" units="uA"/>
<variable name="I_CLC2" public_interface="in" units="uA"/>
<!--var I_GLUT: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_sc</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_CLC2</ci>
<apply>
<minus/>
<ci>I_Na_NaK</ci>
</apply>
<ci>I_K_NaK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_NBC_Na</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_NBC_HCO3</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_Na</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_K</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_Cl</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_IK</ci>
</apply>
</apply>
</apply>
<ci>capacitance</ci>
</apply>
</apply>
</math>
</component>
<component name="rate_constants">
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<!--var k_12: per_second {pub: out};-->
<variable name="ks_12" public_interface="out" units="per_M2_per_second"/>
<variable name="k_21" public_interface="out" units="per_second"/>
<!--var k_23: per_second {pub: out};-->
<variable name="k_32" public_interface="out" units="per_second"/>
<variable name="k_25" public_interface="out" units="per_second"/>
<variable name="k_52" public_interface="out" units="per_second"/>
<variable name="k_34" public_interface="out" units="per_second"/>
<variable name="k_43" public_interface="out" units="per_second"/>
<variable name="k_45" public_interface="out" units="per_second"/>
<variable name="k_54" public_interface="out" units="per_second"/>
<variable name="k_56" public_interface="out" units="per_second"/>
<variable name="k_65" public_interface="out" units="per_second"/>
<variable name="k_61" public_interface="out" units="per_second"/>
<variable name="k_16" public_interface="out" units="per_second"/>
<variable name="k0_12" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_21" public_interface="in" units="per_second"/>
<!--var k0_23: per_M_per_second {pub: in};-->
<variable name="k0_32" public_interface="in" units="per_second"/>
<variable name="k0_25" public_interface="in" units="per_second"/>
<variable name="k0_52" public_interface="in" units="per_second"/>
<variable name="k0_34" public_interface="in" units="per_second"/>
<variable name="k0_43" public_interface="in" units="per_second"/>
<variable name="k0_45" public_interface="in" units="per_second"/>
<variable name="k0_54" public_interface="in" units="per_M_per_second"/>
<!--var k0_25: per_second {pub: in};
var k0_52: per_second {pub: in};-->
<variable name="k0_56" public_interface="in" units="per_second"/>
<variable name="k0_65" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_61" public_interface="in" units="per_second"/>
<variable name="k0_16" public_interface="in" units="per_second"/>
<!--var k0_54_temp: per_M_per_second {pub: out};-->
<!--var k_52_temp: per_second {pub: out};-->
<variable name="delta" public_interface="in" units="dimensionless"/>
<variable name="n" public_interface="in" units="dimensionless"/>
<variable name="z_c" public_interface="in" units="dimensionless"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="alpha_p" public_interface="in" units="dimensionless"/>
<!--var alpha_pp: dimensionless {pub: in};-->
<!--var C_T: umol {pub: in};-->
<variable name="mu" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>ks_12</ci>
<apply>
<times/>
<ci>k0_12</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>n</ci>
</apply>
<ci>alpha_p</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<!-- k_12 = ks_12*pow(Na_m, n);-->
<!--k_12 = k0_12*exp(-alpha_p*mu)*sqr(Na_m);-->
<apply>
<eq/>
<ci>k_21</ci>
<apply>
<times/>
<ci>k0_21</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>n</ci>
<ci>z_Na</ci>
<ci>alpha_p</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<!--k_23 = k0_23*glucose_m;-->
<apply>
<eq/>
<ci>k_32</ci>
<ci>k0_32</ci>
</apply>
<apply>
<eq/>
<ci>k_34</ci>
<ci>k0_34</ci>
</apply>
<apply>
<eq/>
<ci>k_43</ci>
<ci>k0_43</ci>
</apply>
<apply>
<eq/>
<ci>k_45</ci>
<ci>k0_45</ci>
</apply>
<apply>
<eq/>
<ci>k_54</ci>
<apply>
<times/>
<ci>k0_54</ci>
<ci>glucose_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_25</ci>
<ci>k0_25</ci>
</apply>
<!--k_52_temp = k0_52*exp((z_c+n)*delta*mu/2{dimensionless});-->
<apply>
<eq/>
<ci>k_56</ci>
<ci>k0_56</ci>
</apply>
<apply>
<eq/>
<ci>k_65</ci>
<apply>
<times/>
<ci>k0_65</ci>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_61</ci>
<apply>
<times/>
<ci>k0_61</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_c</ci>
<ci>delta</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_16</ci>
<apply>
<times/>
<ci>k0_16</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>z_c</ci>
</apply>
<ci>delta</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_52</ci>
<ci>k0_52</ci>
</apply>
</math>
<!--k0_54_temp = k0_23*k_34*k_45/(k_43*k_32);-->
</component>
<component name="phenomonological_constants">
<variable name="epsilon" public_interface="out" units="per_second"/>
<variable name="lambda" public_interface="out" units="per_M3_per_second5"/>
<variable name="chi" public_interface="out" units="M"/>
<variable name="alpha" public_interface="out" units="M3"/>
<variable name="beta" public_interface="out" units="M2"/>
<variable name="gamma" public_interface="out" units="M3_per_second"/>
<variable name="phi" public_interface="out" units="M_per_second"/>
<!--var k_12: per_second {pub: in};-->
<variable name="k_21" public_interface="in" units="per_second"/>
<!--var k_23: per_second {pub: in};-->
<variable name="k_32" public_interface="in" units="per_second"/>
<variable name="k_34" public_interface="in" units="per_second"/>
<variable name="k_43" public_interface="in" units="per_second"/>
<variable name="k_45" public_interface="in" units="per_second"/>
<variable name="k_54" public_interface="in" units="per_second"/>
<variable name="k_25" public_interface="in" units="per_second"/>
<variable name="k_52" public_interface="in" units="per_second"/>
<variable name="k_56" public_interface="in" units="per_second"/>
<variable name="k_65" public_interface="in" units="per_second"/>
<variable name="k_61" public_interface="in" units="per_second"/>
<variable name="k_16" public_interface="in" units="per_second"/>
<variable name="ks_12" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_23" public_interface="in" units="per_M_per_second"/>
<!--var C_T: umol {pub: in};-->
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<variable initial_value="60000000" name="n_SGLT" public_interface="out" units="dimensionless"/>
<variable name="J_Na_SGLT" public_interface="out" units="umol_per_s"/>
<variable name="I_Na_SGLT" public_interface="out" units="uA"/>
<variable name="J_Gl_SGLT" public_interface="out" units="umol_per_s"/>
<variable name="J_GL_SGLT1" public_interface="out" units="umol_per_s"/>
<!--var I_Gl_SGLT: uA {pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>lambda</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>chi</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>k0_23</ci>
<ci>k_16</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k0_23</ci>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k0_23</ci>
<ci>k_16</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k0_23</ci>
<ci>k_16</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k0_23</ci>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k0_23</ci>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k0_23</ci>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_65</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>phi</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>lambda</ci>
</apply>
<ci>ks_12</ci>
<ci>k_56</ci>
<ci>k_61</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>epsilon</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>lambda</ci>
</apply>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Na_SGLT</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>n_SGLT</ci>
</apply>
<cn cellml:units="per_umol" type="e-notation">6.023<sep/>17</cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>epsilon</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>glucose_m</ci>
</apply>
<apply>
<times/>
<ci>phi</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<ci>gamma</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>alpha</ci>
<apply>
<times/>
<ci>beta</ci>
<ci>glucose_m</ci>
</apply>
<apply>
<times/>
<ci>chi</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>glucose_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Na_SGLT</ci>
<apply>
<times/>
<ci>J_Na_SGLT</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Gl_SGLT</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>J_Na_SGLT</ci>
</apply>
</apply>
<!--to change the unit to picomol/hour-->
<apply>
<eq/>
<ci>J_GL_SGLT1</ci>
<apply>
<times/>
<ci>J_Gl_SGLT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
</math>
<!--I_Gl_SGLT = -J_Gl_SGLT*F;-->
</component>
<component name="CFTR">
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="V_A_Cl" public_interface="out" units="volt"/>
<!--var V_A_HCO3: volt {pub: out};-->
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_CFTR" public_interface="out" units="umol_per_s"/>
<variable name="J_CFTR1" public_interface="out" units="umol_per_s"/>
<!--var J_CFTR_HCO3: umol_per_s {pub: out};-->
<variable name="I_CFTR" public_interface="out" units="uA"/>
<!--var I_CFTR_HCO3: uA {pub: out};-->
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="1.2" name="G_CFTR" public_interface="out" units="S_per_m2"/>
<!-- var HCO3_i: M {pub: in};
var HCO3_m: M {pub: in};
var z_HCO3: dimensionless {init: -1, pub: out};-->
<!-- alpha is multipled by area in order to have a dimensionless unit-->
<!--var alpha_CFTR: per_m2 {init: 0.1, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_CFTR</ci>
<apply>
<divide/>
<ci>I_CFTR</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_CFTR1</ci>
<apply>
<times/>
<ci>J_CFTR</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Cl_m</ci>
<ci>Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_CFTR</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_CFTR</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_Cl</ci>
</apply>
</apply>
</apply>
</math>
<!-- V_A_HCO3 = R*T/(z_HCO3*F)*ln(HCO3_m/HCO3_i);
I_CFTR_HCO3 = 0.25e6{dimensionless}*G_CFTR*A_membrane*(v_mc-V_A_HCO3);
J_CFTR_HCO3 = I_CFTR_HCO3/F;-->
</component>
<component name="CLC2">
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="Cl_s" public_interface="in" units="M"/>
<variable name="V_A_Cl" public_interface="out" units="volt"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="J_CLC2" public_interface="out" units="umol_per_s"/>
<variable name="I_CLC2" public_interface="out" units="uA"/>
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="1" name="G_CLC2" public_interface="out" units="S_per_m2"/>
<!-- alpha is multipled by area in order to have a dimensionless unit-->
<!--var alpha_CFTR: per_m2 {init: 0.1, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_CLC2</ci>
<apply>
<divide/>
<ci>I_CLC2</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Cl_s</ci>
<ci>Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_CLC2</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_CLC2</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_sc</ci>
<ci>V_A_Cl</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="ENaC">
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="V_A_Na" public_interface="out" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_ENaC" public_interface="out" units="umol_per_s"/>
<variable name="J_ENaC1" public_interface="out" units="umol_per_s"/>
<variable name="I_ENaC" public_interface="out" units="uA"/>
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="1.6" name="G_ENaC" public_interface="out" units="S_per_m2"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<!--var alpha_ENaC: per_m2 {init: 0.1, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_ENaC</ci>
<apply>
<divide/>
<ci>I_ENaC</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_ENaC1</ci>
<apply>
<times/>
<ci>J_ENaC</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Na_m</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_ENaC</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
</apply>
<ci>G_ENaC</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="NaKCC">
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="S_C" public_interface="out" units="M"/>
<variable name="J_NaKCC" public_interface="out" units="umol_per_s"/>
<variable name="I_NaKCC" public_interface="out" units="uA"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="u_NaKCC" public_interface="out" units="per_second"/>
<variable initial_value="157.55" name="gamma_1" public_interface="out" units="per_second2"/>
<variable initial_value="20000000" name="gamma_2" public_interface="out" units="per_M4_per_second2"/>
<variable initial_value="1.03" name="gamma_3" public_interface="out" units="per_second"/>
<variable initial_value="1380000" name="gamma_4" public_interface="out" units="per_M4_per_second"/>
<variable initial_value="1e-09" name="alpha_NaKCC" public_interface="out" units="mole_per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>S_C</ci>
<apply>
<times/>
<ci>Na_i</ci>
<ci>Cl_i</ci>
<ci>Cl_i</ci>
<ci>K_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_NaKCC</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>gamma_1</ci>
<apply>
<times/>
<ci>gamma_2</ci>
<ci>S_C</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>gamma_3</ci>
<apply>
<times/>
<ci>gamma_4</ci>
<ci>S_C</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NaKCC</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<ci>alpha_NaKCC</ci>
<ci>u_NaKCC</ci>
<ci>A_membrane</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NaKCC</ci>
<apply>
<times/>
<ci>J_NaKCC</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="A_GLUT2">
<variable initial_value="1000" name="K_12" public_interface="out" units="per_M_per_second"/>
<variable initial_value="12.6" name="K_21" public_interface="out" units="per_second"/>
<variable initial_value="1350" name="K_23" public_interface="out" units="per_second"/>
<variable initial_value="90" name="K_32" public_interface="out" units="per_second"/>
<variable initial_value="12.6" name="K_34" public_interface="out" units="per_second"/>
<variable initial_value="1000" name="K_43" public_interface="out" units="per_M_per_second"/>
<variable initial_value="0.8" name="K_41" public_interface="out" units="per_second"/>
<variable initial_value="12" name="K_14" public_interface="out" units="per_second"/>
<!-- var C_1: M {init: 1e-3, pub: out};
var C_2: M {init: 0, pub: out};-->
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<!-- var CS_1: M {pub: out};
var CS_2: M {pub: out};-->
<variable name="K_S" public_interface="out" units="M"/>
<variable name="J_A_GLUT" public_interface="out" units="umol_per_s"/>
<variable name="C_tot" public_interface="out" units="umol"/>
<variable initial_value="2e8" name="n_GLUT" public_interface="out" units="dimensionless"/>
<variable name="D" public_interface="out" units="per_second"/>
<variable name="J_A_GLUT1" public_interface="out" units="umol_per_s"/>
<!--var time: second {init: 0, pub: out};-->
<variable name="v_cell" public_interface="in" units="m3"/>
<!--var I_GLUT: uA {pub: out};-->
<!--var n_GLUT: dimensionless {init: 13.2e7};
var F: C_per_mol {pub: in};-->
<!-- D = sel
case (glucose_m > 0.02{M}):
(1{dimensionless}+glucose_i/K_S)*(K_32*glucose_m/K_S+K_41)+(1{dimensionless}+glucose_m/K_S)*(K_23*glucose_i/K_S+K_14);
case glucose_m < 0.02{M}:
1e100{per_second};
endsel;-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>C_tot</ci>
<apply>
<divide/>
<ci>n_GLUT</ci>
<cn cellml:units="per_umole" type="e-notation">6.022<sep/>17</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>K_S</ci>
<apply>
<divide/>
<ci>K_21</ci>
<ci>K_12</ci>
</apply>
</apply>
<!-- CS_1 = C_tot* ((K_32*(S_2)/K_S) + K_41)*(S_1/K_S);
CS_2 = C_tot* ((K_23*(S_1)/K_S) + K_14)*(S_2/K_S);-->
<apply>
<eq/>
<ci>J_A_GLUT</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>C_tot</ci>
<ci>D</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>K_14</ci>
<ci>glucose_m</ci>
</apply>
<ci>K_S</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>K_41</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_A_GLUT1</ci>
<apply>
<times/>
<ci>J_A_GLUT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>D</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_m</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_14</ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_i</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>glucose_m</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_41</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="GLUT2">
<variable initial_value="1000" name="K_12" public_interface="out" units="per_M_per_second"/>
<variable initial_value="12.6" name="K_21" public_interface="out" units="per_second"/>
<variable initial_value="1350" name="K_23" public_interface="out" units="per_second"/>
<variable initial_value="90" name="K_32" public_interface="out" units="per_second"/>
<variable initial_value="12.6" name="K_34" public_interface="out" units="per_second"/>
<variable initial_value="1000" name="K_43" public_interface="out" units="per_M_per_second"/>
<variable initial_value="0.8" name="K_41" public_interface="out" units="per_second"/>
<variable initial_value="12" name="K_14" public_interface="out" units="per_second"/>
<!-- var C_1: M {init: 1e-3, pub: out};
var C_2: M {init: 0, pub: out};-->
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_s" public_interface="in" units="M"/>
<!-- var CS_1: M {pub: out};
var CS_2: M {pub: out};-->
<variable name="K_S" public_interface="out" units="M"/>
<variable name="J_GLUT" public_interface="out" units="umol_per_s"/>
<variable name="C_tot" public_interface="out" units="umol"/>
<variable initial_value="2e8" name="n_GLUT" public_interface="out" units="dimensionless"/>
<variable name="D" public_interface="out" units="per_second"/>
<variable name="J_GLUT1" public_interface="out" units="umol_per_s"/>
<!--var time: second {init: 0, pub: out};-->
<variable name="v_cell" public_interface="in" units="m3"/>
<!--var I_GLUT: uA {pub: out};-->
<!--var n_GLUT: dimensionless {init: 13.2e7};
var F: C_per_mol {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>C_tot</ci>
<apply>
<divide/>
<ci>n_GLUT</ci>
<cn cellml:units="per_umole" type="e-notation">6.022<sep/>17</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>K_S</ci>
<apply>
<divide/>
<ci>K_21</ci>
<ci>K_12</ci>
</apply>
</apply>
<!-- CS_1 = C_tot* ((K_32*(S_2)/K_S) + K_41)*(S_1/K_S);
CS_2 = C_tot* ((K_23*(S_1)/K_S) + K_14)*(S_2/K_S);-->
<apply>
<eq/>
<ci>J_GLUT</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>C_tot</ci>
<ci>D</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>K_41</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>K_14</ci>
<ci>glucose_s</ci>
</apply>
<ci>K_S</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_GLUT1</ci>
<apply>
<times/>
<ci>J_GLUT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>D</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_i</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>glucose_s</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_41</ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_s</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_14</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="NaK_channel">
<!-- var time: second {pub: in};-->
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="TO_NaK" public_interface="out" units="per_second"/>
<variable initial_value="5000000" name="n_NaK" public_interface="out" units="dimensionless"/>
<variable initial_value="0.01" name="K_Na_i_half" public_interface="out" units="M"/>
<variable initial_value="0.0015" name="K_K_s_half" public_interface="out" units="M"/>
<variable name="I_Na_NaK" public_interface="out" units="uA"/>
<variable name="I_K_NaK" public_interface="out" units="uA"/>
<variable name="J_NaK" public_interface="out" units="umol_per_s"/>
<variable name="J_Na_NaK" public_interface="out" units="umol_per_s"/>
<variable name="J_K_NaK" public_interface="out" units="umol_per_s"/>
<!--var J_K_NaK: umol_per_s {pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_K_NaK</ci>
<apply>
<times/>
<ci>J_K_NaK</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Na_NaK</ci>
<apply>
<times/>
<ci>J_Na_NaK</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NaK</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>TO_NaK</ci>
<ci>n_NaK</ci>
</apply>
<cn cellml:units="per_umol" type="e-notation">6.023<sep/>17</cn>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
<cn cellml:units="dimensionless">1.13</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_Na_i_half</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1.3</cn>
<ci>K_s</ci>
</apply>
<apply>
<plus/>
<ci>K_K_s_half</ci>
<ci>K_s</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>TO_NaK</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_volt4_per_second" type="e-notation">5.46<sep/>+4</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt3_per_second" type="e-notation">2.43<sep/>+2</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt2_per_second" type="e-notation">2<sep/>+3</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt_per_second">160</cn>
<ci>v_sc</ci>
</apply>
<cn cellml:units="per_second">51</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Na_NaK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>J_NaK</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_K_NaK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NaK</ci>
</apply>
</apply>
</math>
<!--J_Na_NaK = -108e8{dimensionless}*J_NaK;-->
<!--J_K_NaK = -72e8{dimensionless}*J_NaK;-->
</component>
<component name="NBC">
<variable name="time" public_interface="in" units="second"/>
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="HCO3_i" public_interface="in" units="M"/>
<variable name="HCO3_s" public_interface="in" units="M"/>
<variable initial_value="0.08" name="g_NBC" public_interface="out" units="S_per_m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<!--kelvin {init: 300};-->
<variable initial_value="2" name="n" units="dimensionless"/>
<!--var eps: M {init: 1e-12};-->
<variable name="E_NBC" public_interface="out" units="volt"/>
<variable name="J_NBC_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_NBC_HCO3" public_interface="out" units="umol_per_s"/>
<variable name="I_NBC_Na" public_interface="out" units="uA"/>
<variable name="I_NBC_HCO3" public_interface="out" units="uA"/>
<variable initial_value="1e-18" name="eps" units="M"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_NBC</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>n</ci>
</apply>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>Na_s</ci>
<ci>Na_i</ci>
</apply>
<apply>
<power/>
<ci>HCO3_s</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<power/>
<ci>HCO3_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NBC_Na</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<ci>g_NBC</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_sc</ci>
<ci>E_NBC</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NBC_HCO3</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3.0</cn>
<ci>I_NBC_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NBC_Na</ci>
<apply>
<divide/>
<ci>I_NBC_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NBC_HCO3</ci>
<apply>
<divide/>
<ci>I_NBC_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="NHE3">
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="H_Ext" public_interface="in" units="M"/>
<variable name="H_int" public_interface="in" units="M"/>
<variable initial_value="0" name="NH4_ext" public_interface="out" units="M"/>
<variable initial_value="0" name="NH4_int" public_interface="out" units="M"/>
<variable name="J_NHE3_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_H1" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_H" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_NH4" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_Na_Max" public_interface="out" units="umol_per_s"/>
<variable name="I_NHE3_H" public_interface="out" units="uA"/>
<variable name="I_NHE3_Na" public_interface="out" units="uA"/>
<variable initial_value="0.001" name="x_T" units="M"/>
<variable name="sigma" units="per_second"/>
<variable name="P_Na" units="per_second"/>
<variable name="P_H" units="per_second"/>
<variable name="P_NH4" units="per_second"/>
<variable initial_value="1.6" name="P0_Na" units="per_second"/>
<variable initial_value="0.48" name="P0_H" units="per_second"/>
<variable initial_value="1.6" name="P0_NH4" units="per_second"/>
<variable initial_value="0.03" name="K_Na" units="M"/>
<variable initial_value="7.2e-08" name="K_H" units="M"/>
<variable initial_value="0.027" name="K_NH4" units="M"/>
<variable initial_value="1e-09" name="K_I" units="M"/>
<variable initial_value="0" name="f_m" units="dimensionless"/>
<variable initial_value="2" name="f_M" units="dimensionless"/>
<variable name="alpha_ext_Na" units="dimensionless"/>
<variable name="alpha_int_Na" units="dimensionless"/>
<variable name="beta_ext_H" units="dimensionless"/>
<variable name="beta_int_H" units="dimensionless"/>
<variable name="gamma_ext_NH4" units="dimensionless"/>
<variable name="gamma_int_NH4" units="dimensionless"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<!-- var pH_Ext: dimensionless;
var pH_int: dimensionless;
var p_K: dimensionless {init: 6.1, pub: out};
var P_CO2: dimensionless {init: 40, pub: out};
var s: dimensionless {init: 0.03e-3, pub: out};
var HCO3_i: M {pub: in};
var HCO3_m: M {pub: in};-->
<!-- pH_Ext = p_K+log(HCO3_m/(s*P_CO2));
pH_int = p_K+log(HCO3_i/(s*P_CO2));-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_ext_Na</ci>
<apply>
<divide/>
<ci>Na_m</ci>
<ci>K_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_int_Na</ci>
<apply>
<divide/>
<ci>Na_i</ci>
<ci>K_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_ext_H</ci>
<apply>
<divide/>
<ci>H_Ext</ci>
<ci>K_H</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_int_H</ci>
<apply>
<divide/>
<ci>H_int</ci>
<ci>K_H</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_ext_NH4</ci>
<apply>
<divide/>
<ci>NH4_ext</ci>
<ci>K_NH4</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_int_NH4</ci>
<apply>
<divide/>
<ci>NH4_int</ci>
<ci>K_NH4</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_Na</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_Na</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_H</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_H</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_NH4</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_NH4</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>alpha_ext_Na</ci>
<ci>beta_ext_H</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>alpha_int_Na</ci>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>beta_int_H</ci>
</apply>
<apply>
<times/>
<ci>P_NH4</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>alpha_int_Na</ci>
<ci>beta_int_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>alpha_ext_Na</ci>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>beta_ext_H</ci>
</apply>
<apply>
<times/>
<ci>P_NH4</ci>
<ci>gamma_ext_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
</apply>
<ci>v_cell</ci>
<ci>x_T</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_H</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>beta_ext_H</ci>
</apply>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>beta_int_H</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_H</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
</apply>
<ci>v_cell</ci>
<ci>x_T</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_H</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>beta_int_H</ci>
</apply>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>beta_ext_H</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>beta_int_H</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<times/>
<ci>beta_ext_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_H1</ci>
<apply>
<times/>
<ci>J_NHE3_H</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_NH4</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>gamma_ext_NH4</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>beta_ext_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<times/>
<ci>gamma_ext_NH4</ci>
<ci>beta_int_H</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_Na_Max</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<ci>P_Na</ci>
<ci>P_H</ci>
</apply>
<apply>
<plus/>
<ci>P_Na</ci>
<ci>P_H</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NHE3_H</ci>
<apply>
<times/>
<ci>J_NHE3_H</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NHE3_Na</ci>
<apply>
<times/>
<ci>J_NHE3_Na</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="AE1">
<!--var A_membrane: m2 {pub: in};-->
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="HCO3_i" public_interface="in" units="M"/>
<variable name="HCO3_m" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="J_AE1_HCO3" public_interface="out" units="umol_per_s"/>
<!--var J_AE1_HCO31: umol_per_s {pub: out};-->
<variable name="J_AE1_Cl" public_interface="out" units="umol_per_s"/>
<!--var J_AE1_Cl1: umol_per_s {pub: out};-->
<variable name="I_AE1_HCO3" public_interface="out" units="uA"/>
<variable name="I_AE1_Cl" public_interface="out" units="uA"/>
<variable name="J_HCO3_influx" public_interface="out" units="umol_per_s"/>
<variable name="J_Cl_influx" public_interface="out" units="umol_per_s"/>
<variable initial_value="0.198" name="K_HCO3_m" units="M"/>
<variable initial_value="0.198" name="K_HCO3_i" units="M"/>
<variable initial_value="0.05" name="K_Cl_s" units="M"/>
<variable initial_value="0.05" name="K_Cl_i" units="M"/>
<variable initial_value="1247" name="P_HCO3_m" units="per_second"/>
<variable initial_value="135" name="P_HCO3_i" units="per_second"/>
<variable initial_value="562" name="P_Cl_s" units="per_second"/>
<variable initial_value="61" name="P_Cl_i" units="per_second"/>
<variable name="beta_ext" units="dimensionless"/>
<variable name="beta_int" units="dimensionless"/>
<variable name="gamma_ext" units="dimensionless"/>
<variable name="gamma_int" units="dimensionless"/>
<variable name="sigma" units="per_second"/>
<variable initial_value="0.001" name="x_Tmax" units="M"/>
<variable initial_value="0.172" name="K_I" units="M"/>
<variable name="x_T" units="M"/>
<variable name="x_ext" units="M"/>
<variable name="x_int" units="M"/>
<variable name="Jo_bm" units="umol_per_s"/>
<variable name="Ji_bm" units="umol_per_s"/>
<variable name="Js_bm" units="umol_per_s"/>
<variable name="Jo_cm" units="umol_per_s"/>
<variable name="Ji_cm" units="umol_per_s"/>
<variable name="Js_cm" units="umol_per_s"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>x_T</ci>
<apply>
<divide/>
<ci>x_Tmax</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>HCO3_i</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_ext</ci>
<apply>
<divide/>
<ci>HCO3_m</ci>
<ci>K_HCO3_m</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_int</ci>
<apply>
<divide/>
<ci>HCO3_i</ci>
<ci>K_HCO3_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_ext</ci>
<apply>
<divide/>
<ci>Cl_m</ci>
<ci>K_Cl_s</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_int</ci>
<apply>
<divide/>
<ci>Cl_i</ci>
<ci>K_Cl_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>beta_ext</ci>
<ci>gamma_ext</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>beta_int</ci>
<ci>gamma_int</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_AE1_HCO3</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<ci>v_cell</ci>
<ci>x_Tmax</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<!--J_AE1_HCO31 = J_AE1_HCO3*36e8{dimensionless};
J_AE1_Cl1 = -J_AE1_HCO31;-->
<apply>
<eq/>
<ci>J_AE1_Cl</ci>
<ci>J_AE1_HCO3</ci>
</apply>
<apply>
<eq/>
<ci>I_AE1_HCO3</ci>
<apply>
<times/>
<ci>J_AE1_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_AE1_Cl</ci>
<apply>
<times/>
<ci>J_AE1_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>x_ext</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
<ci>sigma</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>x_int</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
</apply>
</apply>
<ci>sigma</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_HCO3_influx</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Cl_influx</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Jo_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
<apply>
<divide/>
<ci>K_HCO3_i</ci>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>HCO3_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ji_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
<apply>
<divide/>
<ci>K_HCO3_m</ci>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>HCO3_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Js_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Jo_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
<apply>
<divide/>
<ci>K_Cl_i</ci>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ji_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
<apply>
<divide/>
<ci>K_Cl_s</ci>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>Cl_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Js_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<!-- def comp Hydration as
var time: second {pub: in};
var v_cell: m3 {pub: in};
var CO2_i: M {init: 1e-3, pub: out};
var CO2_s: M {init: 1.2e-3, pub: out};
var CO2_m: M {init: 1.2e-3, pub: out};
var H_int: M {pub: in};
var H_Ext: M {pub: in};
var HCO3_i: M {pub: in};
var HCO3_m: M {pub: in};
var p_CO2: per_second {init: 5e-6, pub: out};
var k_f: per_second {init: 0.03, pub: out};
var k_r: per_M_per_second {init: 20e-3, pub: out};
var J_CDF_A: umol_per_s {pub: out};
var J_CDF_B: umol_per_s {pub: out};
var J_buf_C: umol_per_s {pub: out};
var J_buf_C1: umol_per_s {pub: out};
var J_buf_A: umol_per_s {pub: out};
var J_D_P_CO2: umol_per_s {pub: in};-->
<!-- J_CDF_A = -1e9{dimensionless}*v_cell*p_CO2*(CO2_i-CO2_m);
J_CDF_B = 1e9{dimensionless}*v_cell*p_CO2*(CO2_i-CO2_s);
J_buf_C = 1e9{dimensionless}*v_cell*(k_f*CO2_i-k_r*HCO3_i*H_int);
J_buf_C1 = J_buf_C*36e8{dimensionless};
J_buf_A = 1e9{dimensionless}*v_cell*(k_f*CO2_m-k_r*HCO3_m*H_Ext);
ode(CO2_i, time) = -J_buf_C+J_CDF_A+J_CDF_B;
ode(CO2_m, time) = -J_buf_A-J_CDF_A-0{dimensionless}*J_D_P_CO2;
enddef;-->
<component name="IK">
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Ca_i" public_interface="in" units="M"/>
<variable name="V_A_IK" public_interface="out" units="volt"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="J_IK" public_interface="out" units="umol_per_s"/>
<variable name="I_IK" public_interface="out" units="uA"/>
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="2.1" name="G_IK" public_interface="out" units="S_per_m2"/>
<!--var P_IK: dimensionless {pub: out};-->
<variable initial_value="0.0003" name="K_s_IK" public_interface="out" units="M"/>
<variable initial_value="1.7" name="eta" public_interface="out" units="dimensionless"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<!--var alpha_IK: per_m2 {init: 1, pub: out};-->
<!--P_IK = 1{dimensionless}/(1{dimensionless}+pow(K_s_IK/Ca_i, eta));-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_A_IK</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_s</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_IK</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
</apply>
<ci>G_IK</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_sc</ci>
<ci>V_A_IK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_IK</ci>
<apply>
<divide/>
<ci>I_IK</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="BK">
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="Ca_i" public_interface="in" units="M"/>
<variable name="V_A_BK" public_interface="out" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_BK" public_interface="out" units="umol_per_s"/>
<variable name="I_BK" public_interface="out" units="uA"/>
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="0.8" name="G_BK" public_interface="out" units="S_per_m2"/>
<variable name="P_BK" public_interface="out" units="dimensionless"/>
<variable initial_value="0.0003" name="K_s_BK" public_interface="out" units="M"/>
<variable initial_value="1.7" name="eta" public_interface="out" units="dimensionless"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<!--var alpha_BK: per_m2 {init: 0.1, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>P_BK</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>K_s_BK</ci>
<ci>Ca_i</ci>
</apply>
<ci>eta</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_BK</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_m</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_BK</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
</apply>
<ci>G_BK</ci>
<ci>P_BK</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_BK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_BK</ci>
<apply>
<divide/>
<ci>I_BK</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<!-- def comp Apical_Diffusion as
var v_mc: volt {pub: in};
var Na_m: M {pub: in};
var K_m: M {pub: in};
var Cl_m: M {pub: in};
var Na_i: M {pub: in};
var K_i: M {pub: in};
var Cl_i: M {pub: in};
var F: C_per_mol {pub: in};
var z_Na: dimensionless {pub: in};
var z_k: dimensionless {pub: in};
var z_Cl: dimensionless {pub: in};
var R: J_per_K_per_mol {pub: in};
var A_membrane: m2 {pub: in};
var T: kelvin {pub: in};
var P_Na: m_per_s {pub: in};
var P_k: m_per_s {pub: in};
var P_Cl: m_per_s {pub: in};
var J_D_Na: umol_per_s {pub: out};
var J_D_NA: umol_per_s {pub: out};
var I_D_Na: uA {pub: out};
var J_D_K: umol_per_s {pub: out};
var I_D_K: uA {pub: out};
var J_D_Cl: umol_per_s {pub: out};
var I_D_Cl: uA {pub: out};
var J_D_Apical: umol_per_s {pub: out};-->
<!-- J_D_Na = -1e9{litre_umol_per_mol_per_m3}*P_Na*F*z_Na*v_mc/(R*T)*(Na_m-Na_i*exp(z_Na*F*v_mc/(R*T)))/(1{dimensionless}-exp(z_Na*F*v_mc/(R*T)))*A_membrane;
J_D_NA = 36e8{dimensionless}*J_D_Na;
I_D_Na = J_D_Na*F;
J_D_K = -1e9{litre_umol_per_mol_per_m3}*P_k*F*z_k*v_mc/(R*T)*(K_m-K_i*exp(z_k*F*v_mc/(R*T)))/(1{dimensionless}-exp(z_k*F*v_mc/(R*T)))*A_membrane;
I_D_K = J_D_K*F;
J_D_Cl = -1e9{litre_umol_per_mol_per_m3}*P_Cl*F*z_Cl*v_mc/(R*T)*(Cl_m-Cl_i*exp(z_Cl*F*v_mc/(R*T)))/(1{dimensionless}-exp(z_Cl*F*v_mc/(R*T)))*A_membrane;
I_D_Cl = J_D_Cl*F;
J_D_Apical = J_D_Na+J_D_K+J_D_Cl;
enddef;-->
<!-- def comp Basolateral_Diffusion as
var v_sc: volt {pub: in};
var Na_s: M {pub: in};
var K_s: M {pub: in};
var Cl_s: M {pub: in};
var Na_i: M {pub: in};
var K_i: M {pub: in};
var Cl_i: M {pub: in};
var F: C_per_mol {pub: in};
var z_Na: dimensionless {pub: in};
var z_k: dimensionless {pub: in};
var z_Cl: dimensionless {pub: in};
var R: J_per_K_per_mol {pub: in};
var A_Basol: m2 {pub: in};
var T: kelvin {pub: in};
var P_Bl_Na: m_per_s {pub: in};
var P_Bl_k: m_per_s {pub: in};
var P_Bl_Cl: m_per_s {pub: in};
var J_D_Bl_Na: umol_per_s {pub: out};
var I_D_Bl_Na: uA {pub: out};
var J_D_Bl_K: umol_per_s {pub: out};
var I_D_Bl_K: uA {pub: out};
var J_D_Bl_Cl: umol_per_s {pub: out};
var I_D_Bl_Cl: uA {pub: out};
var J_D_Basolateral: umol_per_s {pub: out};-->
<!-- J_D_Bl_Na = -1e9{litre_umol_per_mol_per_m3}*P_Bl_Na*F*z_Na*v_sc/(R*T)*(Na_s-Na_i*exp(z_Na*F*v_sc/(R*T)))/(1{dimensionless}-exp(z_Na*F*v_sc/(R*T)))*A_Basol;
I_D_Bl_Na = J_D_Bl_Na*F;
J_D_Bl_K = -1e9{litre_umol_per_mol_per_m3}*P_Bl_k*F*z_k*v_sc/(R*T)*(K_s-K_i*exp(z_k*F*v_sc/(R*T)))/(1{dimensionless}-exp(z_k*F*v_sc/(R*T)))*A_Basol;
I_D_Bl_K = J_D_Bl_K*F;
J_D_Bl_Cl = -1e9{litre_umol_per_mol_per_m3}*P_Bl_Cl*F*z_Cl*v_sc/(R*T)*(Cl_s-Cl_i*exp(z_Cl*F*v_sc/(R*T)))/(1{dimensionless}-exp(z_Cl*F*v_sc/(R*T)))*A_Basol;
I_D_Bl_Cl = J_D_Bl_Cl*F;
J_D_Basolateral = J_D_Bl_Na+J_D_Bl_K+J_D_Bl_Cl;
enddef;-->
<component name="Paracellular_voltage">
<variable name="v_ms" public_interface="out" units="volt"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<!-- var time: second {pub: in};
var capacitance: uF {pub: in};
var I_D_P_Na: uA {pub: in};
var I_D_P_K: uA {pub: in};
var I_D_P_Cl: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ms</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>v_sc</ci>
</apply>
</apply>
</math>
<!--ode(v_ms, time) = (I_D_P_Na+I_D_P_K+I_D_P_Cl)/capacitance;-->
</component>
<component name="Paracellular_Diffusion">
<variable name="v_ms" public_interface="in" units="volt"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Cl_s" public_interface="in" units="M"/>
<!-- var CO2_m: M {pub: in};
var CO2_i: M {pub: in};-->
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="A_membrane" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_P_Na" public_interface="in" units="m_per_s"/>
<variable name="P_P_k" public_interface="in" units="m_per_s"/>
<variable name="P_P_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_P_Na" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_Na" public_interface="out" units="uA"/>
<variable name="J_D_P_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_K" public_interface="out" units="uA"/>
<variable name="J_D_P_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Paracellular" public_interface="out" units="umol_per_s"/>
<variable name="V_P_Na" public_interface="out" units="volt"/>
<variable name="V_P_K" public_interface="out" units="volt"/>
<variable name="V_P_Cl" public_interface="out" units="volt"/>
<variable initial_value="4" name="G_P_Na" public_interface="out" units="S_per_m2"/>
<variable initial_value="1" name="G_P_K" public_interface="out" units="S_per_m2"/>
<variable initial_value="4" name="G_P_Cl" public_interface="out" units="S_per_m2"/>
<!-- var I_D_P_CO2: uA {pub: out};
var J_D_P_CO2: umol_per_s {pub: out};
var G_P_CO2: S_per_m2 {init: 7, pub: out};
var v_P_CO2: volt {pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_P_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Na_m</ci>
<ci>Na_s</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_Na</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_P_Na</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_ms</ci>
<ci>V_P_Na</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_Na</ci>
<apply>
<divide/>
<ci>I_D_P_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_P_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_m</ci>
<ci>K_s</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_K</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_P_K</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_ms</ci>
<ci>V_P_K</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_K</ci>
<apply>
<divide/>
<ci>I_D_P_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_P_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Cl_m</ci>
<ci>Cl_s</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_Cl</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_P_Cl</ci>
<ci>A_membrane</ci>
<apply>
<minus/>
<ci>v_ms</ci>
<ci>V_P_Cl</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_Cl</ci>
<apply>
<divide/>
<ci>I_D_P_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<!-- v_P_CO2 = R*T/F*ln(CO2_m/CO2_i);
I_D_P_CO2 = 1e6{dimensionless}*G_P_CO2*A_membrane*(v_ms-v_P_CO2);
J_D_P_CO2 = I_D_P_CO2/F;-->
<apply>
<eq/>
<ci>J_D_Paracellular</ci>
<apply>
<plus/>
<ci>J_D_P_Na</ci>
<ci>J_D_P_K</ci>
<ci>J_D_P_Cl</ci>
</apply>
</apply>
</math>
</component>
<!-- def comp interstitium_concentration as
var glucose_s: M {init: 0, pub: out};-->
<!-- //var J_GLUT: umol_per_s {pub: in};
//var J_G_Blood: umol_per_s {pub: in};
//var v_cell: m3 {pub: in};
//var time: second {pub: in};-->
<!-- //ode(glucose_s, time) = 1e-9{m3_mole_per_litre_per_umol}/v_cell*(-J_GLUT-J_G_Blood);
enddef;-->
<!-- def comp Blood_concentration as
var glucose_s: M {pub: in};
var glucose_B: M {init: 0.012, pub: out};
var J_G_Blood: umol_per_s {pub: out};
var K: dimensionless {init: 1e-6, pub: out};-->
<!-- //var v_cell: m3 {pub: in};-->
<!-- var time: second {pub: in};-->
<!-- J_G_Blood = K*(glucose_s-glucose_B);
ode(glucose_B, time) = 0{dimensionless};
enddef;-->
<connection>
<map_components component_1="rate_constants" component_2="phenomonological_constants"/>
<map_variables variable_1="ks_12" variable_2="ks_12"/>
<!--vars k_12 and k_12;-->
<map_variables variable_1="k_21" variable_2="k_21"/>
<!--vars k_23 and k_23;-->
<map_variables variable_1="k_32" variable_2="k_32"/>
<map_variables variable_1="k_25" variable_2="k_25"/>
<map_variables variable_1="k_52" variable_2="k_52"/>
<map_variables variable_1="k_34" variable_2="k_34"/>
<map_variables variable_1="k_43" variable_2="k_43"/>
<map_variables variable_1="k_45" variable_2="k_45"/>
<map_variables variable_1="k_54" variable_2="k_54"/>
<map_variables variable_1="k_56" variable_2="k_56"/>
<map_variables variable_1="k_65" variable_2="k_65"/>
<map_variables variable_1="k_61" variable_2="k_61"/>
<map_variables variable_1="k_16" variable_2="k_16"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Cell_concentration"/>
<map_variables variable_1="J_Na_SGLT" variable_2="J_Na_SGLT"/>
<map_variables variable_1="J_Gl_SGLT" variable_2="J_Gl_SGLT"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="parameters"/>
<map_variables variable_1="k0_23" variable_2="k0_23"/>
<!--vars C_T and C_T;-->
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Apical_voltage"/>
<map_variables variable_1="I_Na_SGLT" variable_2="I_Na_SGLT"/>
<!--vars I_Gl_SGLT and I_Gl_SGLT;-->
</connection>
<connection>
<map_components component_1="rate_constants" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
</connection>
<connection>
<map_components component_1="rate_constants" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
</connection>
<connection>
<map_components component_1="rate_constants" component_2="parameters"/>
<map_variables variable_1="k0_12" variable_2="k0_12"/>
<map_variables variable_1="k0_21" variable_2="k0_21"/>
<!--vars k0_23 and k0_23;-->
<map_variables variable_1="k0_32" variable_2="k0_32"/>
<map_variables variable_1="k0_25" variable_2="k0_25"/>
<map_variables variable_1="k0_52" variable_2="k0_52"/>
<map_variables variable_1="k0_34" variable_2="k0_34"/>
<map_variables variable_1="k0_43" variable_2="k0_43"/>
<map_variables variable_1="k0_45" variable_2="k0_45"/>
<map_variables variable_1="k0_54" variable_2="k0_54"/>
<map_variables variable_1="k0_56" variable_2="k0_56"/>
<map_variables variable_1="k0_65" variable_2="k0_65"/>
<map_variables variable_1="k0_61" variable_2="k0_61"/>
<map_variables variable_1="k0_16" variable_2="k0_16"/>
<map_variables variable_1="delta" variable_2="delta"/>
<map_variables variable_1="alpha_p" variable_2="alpha_p"/>
<!--vars alpha_pp and alpha_pp;-->
<!--vars C_T and C_T;-->
<map_variables variable_1="n" variable_2="n"/>
<map_variables variable_1="z_c" variable_2="z_c"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="mu" variable_2="mu"/>
</connection>
<connection>
<map_components component_1="CFTR" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
</connection>
<!-- def map between CFTR_B and parameters for
vars F and F;
vars z_Cl and z_Cl;
vars R and R;
vars T and T;
vars A_membrane and A_membrane;
enddef;-->
<connection>
<map_components component_1="CFTR" component_2="Cell_concentration"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<!-- vars HCO3_i and HCO3_i;
vars HCO3_m and HCO3_m;-->
<map_variables variable_1="J_CFTR" variable_2="J_CFTR"/>
<!-- vars J_CFTR_HCO3 and J_CFTR_HCO3;-->
</connection>
<!-- def map between CFTR_B and Cell_concentration for
vars Cl_i and Cl_i;-->
<!-- // vars HCO3_i and HCO3_i;
// vars HCO3_m and HCO3_m;-->
<!-- vars J_CFTR_B and J_CFTR_B;-->
<!-- // vars J_CFTR_HCO3 and J_CFTR_HCO3;
enddef;-->
<connection>
<map_components component_1="CFTR" component_2="Apical_concentrations"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
</connection>
<!-- def map between CFTR_B and Blood_concentrations for
vars Cl_s and Cl_s;
vars J_CFTR_B and J_CFTR_B;
enddef;-->
<connection>
<map_components component_1="CFTR" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_CFTR" variable_2="I_CFTR"/>
<!--vars I_CFTR_HCO3 and I_CFTR_HCO3;-->
</connection>
<connection>
<map_components component_1="CLC2" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
</connection>
<connection>
<map_components component_1="CLC2" component_2="Cell_concentration"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="J_CLC2" variable_2="J_CLC2"/>
</connection>
<connection>
<map_components component_1="CLC2" component_2="Blood_concentrations"/>
<map_variables variable_1="Cl_s" variable_2="Cl_s"/>
<map_variables variable_1="J_CLC2" variable_2="J_CLC2"/>
</connection>
<connection>
<map_components component_1="CLC2" component_2="Basol_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
<map_variables variable_1="I_CLC2" variable_2="I_CLC2"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="J_ENaC" variable_2="J_ENaC"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_ENaC" variable_2="I_ENaC"/>
</connection>
<connection>
<map_components component_1="IK" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Cell_concentration"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="J_IK" variable_2="J_IK"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Blood_concentrations"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
<map_variables variable_1="J_IK" variable_2="J_IK"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Basol_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
<map_variables variable_1="I_IK" variable_2="I_IK"/>
</connection>
<connection>
<map_components component_1="BK" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Cell_concentration"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="J_BK" variable_2="J_BK"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Apical_concentrations"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_BK" variable_2="I_BK"/>
</connection>
<connection>
<map_components component_1="NaKCC" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="J_NaKCC" variable_2="J_NaKCC"/>
</connection>
<connection>
<map_components component_1="NaKCC" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="GLUT2"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
<map_variables variable_1="J_GLUT" variable_2="J_GLUT"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Blood_concentrations" component_2="GLUT2"/>
<map_variables variable_1="glucose_s" variable_2="glucose_s"/>
<map_variables variable_1="J_GLUT" variable_2="J_GLUT"/>
</connection>
<!-- def map between Blood_concentrations and Cell_concentration for
vars Na_s and Na_s;
vars glucose_s and glucose_s;-->
<!-- //vars v_B and v_B;-->
<!-- vars J_D_P_Gl and J_D_P_Gl;
vars J_w_B and J_w_B;
vars J_w_P and J_w_P;-->
<!-- //vars v_w and v_w;
enddef;-->
<!-- def map between GLUT2 and Basol_voltage for
vars I_GLUT and I_GLUT;
enddef;-->
<!-- def map between Cell_concentration and Na_Cl for
vars Na_i and Na_i;
vars Cl_i and Cl_i;
vars J_NaCl and J_NaCl;
enddef;-->
<connection>
<map_components component_1="Cell_concentration" component_2="NHE3"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="H_int" variable_2="H_int"/>
<map_variables variable_1="H_Ext" variable_2="H_Ext"/>
<!-- vars HCO3_i and HCO3_i;
vars HCO3_m and HCO3_m;-->
<!-- vars Cl_i and Cl_i;-->
<map_variables variable_1="J_NHE3_Na" variable_2="J_NHE3_Na"/>
<map_variables variable_1="J_NHE3_H" variable_2="J_NHE3_H"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NaK_channel"/>
<map_variables variable_1="F" variable_2="F"/>
<!-- vars time and time;-->
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="NaK_channel"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="J_Na_NaK" variable_2="J_Na_NaK"/>
<map_variables variable_1="J_K_NaK" variable_2="J_K_NaK"/>
</connection>
<connection>
<map_components component_1="Blood_concentrations" component_2="NaK_channel"/>
<map_variables variable_1="J_Na_NaK" variable_2="J_Na_NaK"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
<map_variables variable_1="J_K_NaK" variable_2="J_K_NaK"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="NaK_channel"/>
<map_variables variable_1="I_Na_NaK" variable_2="I_Na_NaK"/>
<map_variables variable_1="I_K_NaK" variable_2="I_K_NaK"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="NBC"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<!-- vars Cl_i and Cl_i;-->
<map_variables variable_1="J_NBC_Na" variable_2="J_NBC_Na"/>
<map_variables variable_1="J_NBC_HCO3" variable_2="J_NBC_HCO3"/>
<map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
<map_variables variable_1="HCO3_s" variable_2="HCO3_s"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NBC"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="NBC"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
<map_variables variable_1="I_NBC_HCO3" variable_2="I_NBC_HCO3"/>
<map_variables variable_1="I_NBC_Na" variable_2="I_NBC_Na"/>
</connection>
<connection>
<map_components component_1="Blood_concentrations" component_2="NBC"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
<map_variables variable_1="J_NBC_Na" variable_2="J_NBC_Na"/>
</connection>
<connection>
<map_components component_1="Blood_concentrations" component_2="NaKCC"/>
<map_variables variable_1="J_NaKCC" variable_2="J_NaKCC"/>
</connection>
<!-- def map between Cell_concentration and NBC for
vars Na_i and Na_i;
enddef;-->
<!-- def map between Cell_concentration and Apical_Diffusion for
vars Na_i and Na_i;
vars K_i and K_i;
vars Cl_i and Cl_i;
vars J_D_Na and J_D_Na;
vars J_D_K and J_D_K;
vars J_D_Cl and J_D_Cl;
enddef;-->
<connection>
<map_components component_1="Apical_concentrations" component_2="Cell_concentration"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
</connection>
<!-- def map between Apical_concentrations and Apical_Diffusion for
vars Na_m and Na_m;
vars K_m and K_m;
vars Cl_m and Cl_m;-->
<!-- //vars F and F;
enddef;-->
<!-- def map between parameters and Apical_Diffusion for
vars F and F;
vars z_Na and z_Na;
vars z_k and z_k;
vars z_Cl and z_Cl;
vars R and R;
vars T and T;
vars A_membrane and A_membrane;
vars P_Na and P_Na;
vars P_k and P_k;
vars P_Cl and P_Cl;
enddef;-->
<!-- def map between parameters and Cell_concentration for
vars A_membrane and A_membrane;
enddef;-->
<!-- def map between parameters and Apical_voltage for
vars time and time;
enddef;-->
<!-- def map between Apical_Diffusion and Apical_voltage for
vars I_D_Na and I_D_Na;
vars I_D_K and I_D_K;
vars I_D_Cl and I_D_Cl;
vars v_mc and v_mc;
enddef;-->
<connection>
<map_components component_1="parameters" component_2="Apical_voltage"/>
<map_variables variable_1="capacitance" variable_2="capacitance"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
</connection>
<!-- def map between Cell_concentration and Basolateral_Diffusion for
vars Na_i and Na_i;
vars K_i and K_i;
vars Cl_i and Cl_i;
vars J_D_Bl_Na and J_D_Bl_Na;
vars J_D_Bl_K and J_D_Bl_K;
vars J_D_Bl_Cl and J_D_Bl_Cl;
enddef;-->
<connection>
<map_components component_1="Blood_concentrations" component_2="Cell_concentration"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
<map_variables variable_1="Cl_s" variable_2="Cl_s"/>
<map_variables variable_1="glucose_s" variable_2="glucose_s"/>
<map_variables variable_1="J_D_P_Gl" variable_2="J_D_P_Gl"/>
<map_variables variable_1="J_w_B" variable_2="J_w_B"/>
<map_variables variable_1="J_w_P" variable_2="J_w_P"/>
</connection>
<!-- def map between Blood_concentrations and Basolateral_Diffusion for
vars Na_s and Na_s;
vars K_s and K_s;
vars Cl_s and Cl_s;-->
<!-- // vars J_D_Bl_Na and J_D_Bl_Na;
// vars J_D_Bl_K and J_D_Bl_K;
// vars J_D_Bl_Cl and J_D_Bl_Cl;-->
<!-- //vars F and F;
enddef;-->
<!-- def map between parameters and Basolateral_Diffusion for
vars F and F;
vars z_Na and z_Na;
vars z_k and z_k;
vars z_Cl and z_Cl;
vars R and R;
vars T and T;
vars A_Basol and A_Basol;
vars P_Bl_Na and P_Bl_Na;
vars P_Bl_k and P_Bl_k;
vars P_Bl_Cl and P_Bl_Cl;
enddef;-->
<connection>
<map_components component_1="parameters" component_2="Cell_concentration"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NHE3"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
<!-- vars time and time;
vars R and R;
vars T and T;-->
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="NHE3"/>
<!--vars v_sc and v_sc;-->
<map_variables variable_1="I_NHE3_H" variable_2="I_NHE3_H"/>
<map_variables variable_1="I_NHE3_Na" variable_2="I_NHE3_Na"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="NHE3"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
</connection>
<!-- def map between AE1 and NBC for
vars HCO3_s and HCO3_s;
vars HCO3_i and HCO3_i;
enddef; -->
<connection>
<map_components component_1="parameters" component_2="AE1"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
<!--vars A_membrane and A_membrane;-->
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="AE1"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
<map_variables variable_1="HCO3_m" variable_2="HCO3_m"/>
<map_variables variable_1="J_AE1_Cl" variable_2="J_AE1_Cl"/>
<map_variables variable_1="J_AE1_HCO3" variable_2="J_AE1_HCO3"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="AE1"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
</connection>
<!-- def map between Apical_voltage and AE1 for
//vars v_sc and v_sc;-->
<!-- vars I_AE1_Cl and I_AE1_Cl;
vars I_AE1_HCO3 and I_AE1_HCO3;
enddef;-->
<!-- def map between Basolateral_concentrations and Basol_voltage for
vars time and time;
enddef;-->
<!-- def map between Basolateral_Diffusion and Basol_voltage for
vars I_D_Bl_Na and I_D_Bl_Na;
vars I_D_Bl_K and I_D_Bl_K;
vars I_D_Bl_Cl and I_D_Bl_Cl;
vars v_sc and v_sc;
enddef;-->
<connection>
<map_components component_1="parameters" component_2="Basol_voltage"/>
<map_variables variable_1="capacitance" variable_2="capacitance"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
<!--vars F and F;-->
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="parameters"/>
<map_variables variable_1="time" variable_2="time"/>
<!--vars F and F;-->
</connection>
<connection>
<map_components component_1="parameters" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_membrane" variable_2="A_membrane"/>
<map_variables variable_1="P_P_Na" variable_2="P_P_Na"/>
<map_variables variable_1="P_P_k" variable_2="P_P_k"/>
<map_variables variable_1="P_P_Cl" variable_2="P_P_Cl"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="J_D_P_Na" variable_2="J_D_P_Na"/>
</connection>
<!-- def map between Apical_concentrations and Paracellular_voltage for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Paracellular_Diffusion" component_2="Paracellular_voltage"/>
<!-- vars I_D_P_Na and I_D_P_Na;
vars I_D_P_K and I_D_P_K;
vars I_D_P_Cl and I_D_P_Cl;-->
<map_variables variable_1="v_ms" variable_2="v_ms"/>
</connection>
<!-- def map between parameters and Paracellular_voltage for
vars capacitance and capacitance;
vars time and time;
enddef;-->
<connection>
<map_components component_1="Basol_voltage" component_2="Paracellular_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="Paracellular_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="I_D_P_Na" variable_2="I_D_P_Na"/>
<map_variables variable_1="I_D_P_K" variable_2="I_D_P_K"/>
<map_variables variable_1="I_D_P_Cl" variable_2="I_D_P_Cl"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="I_D_P_Na" variable_2="I_D_P_Na"/>
<map_variables variable_1="I_D_P_K" variable_2="I_D_P_K"/>
<map_variables variable_1="I_D_P_Cl" variable_2="I_D_P_Cl"/>
</connection>
<connection>
<map_components component_1="Blood_concentrations" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
<map_variables variable_1="Cl_s" variable_2="Cl_s"/>
<map_variables variable_1="J_D_P_Na" variable_2="J_D_P_Na"/>
<map_variables variable_1="J_D_P_K" variable_2="J_D_P_K"/>
<map_variables variable_1="J_D_P_Cl" variable_2="J_D_P_Cl"/>
<!--vars F and F;-->
</connection>
<connection>
<map_components component_1="parameters" component_2="Blood_concentrations"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="A_GLUT2"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
<map_variables variable_1="J_A_GLUT" variable_2="J_A_GLUT"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="A_GLUT2"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
<!--vars J_GLUT and J_GLUT;-->
</connection>
<!-- def map between Cell_concentration and Hydration for
vars H_int and H_int;
vars HCO3_i and HCO3_i;
vars J_buf_C and J_buf_C;
vars H_Ext and H_Ext;
vars HCO3_m and HCO3_m;
vars J_buf_A and J_buf_A;
enddef;-->
<!-- def map between parameters and Hydration for
vars time and time;
vars v_cell and v_cell;
enddef;-->
<!-- def map between Paracellular_Diffusion and Hydration for
vars CO2_i and CO2_i;
vars CO2_m and CO2_m;
vars J_D_P_CO2 and J_D_P_CO2;
enddef;-->
</model>
