- Author:
- Soroush Safaei <ssaf006@aucklanduni.ac.nz>
- Date:
- 2018-07-29 23:30:09+12:00
- Desc:
- modular
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/436/rawfile/a3e40926098096b3269538f412265b538889ec5e/main.cellml
<?xml version='1.0'?>
<model name="Glucose_Uptake" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
<!-- 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-->
<import xlink:href="Units.cellml">
<units name="per_s" units_ref="per_s"/>
<units name="per_s2" units_ref="per_s2"/>
<units name="V_per_s" units_ref="V_per_s"/>
<units name="per_m" units_ref="per_m"/>
<units name="per_V_s" units_ref="per_V_s"/>
<units name="per_V4_s" units_ref="per_V4_s"/>
<units name="per_V3_s" units_ref="per_V3_s"/>
<units name="per_V2_s" units_ref="per_V2_s"/>
<units name="per_mol" units_ref="per_mol"/>
<units name="umol" units_ref="umol"/>
<units name="per_umol" units_ref="per_umol"/>
<units name="M" units_ref="M"/>
<units name="per_M" units_ref="per_M"/>
<units name="m2" units_ref="m2"/>
<units name="cm" units_ref="cm"/>
<units name="um2" units_ref="um2"/>
<units name="um" units_ref="um"/>
<units name="um3" units_ref="um3"/>
<units name="m3" units_ref="m3"/>
<units name="mol_per_um2" units_ref="mol_per_um2"/>
<units name="mM" units_ref="mM"/>
<units name="mV" units_ref="mV"/>
<units name="J_per_mol_K" units_ref="J_per_mol_K"/>
<units name="pA_per_cm2" units_ref="pA_per_cm2"/>
<units name="uA" units_ref="uA"/>
<units name="uA_per_um2" units_ref="uA_per_um2"/>
<units name="uA_per_A" units_ref="uA_per_A"/>
<units name="umol_per_mol" units_ref="umol_per_mol"/>
<units name="C_per_mol" units_ref="C_per_mol"/>
<units name="per_M_s" units_ref="per_M_s"/>
<units name="per_M2_s" units_ref="per_M2_s"/>
<units name="per_M4_s2" units_ref="per_M4_s2"/>
<units name="per_M4_s" units_ref="per_M4_s"/>
<units name="uM" units_ref="uM"/>
<units name="M2" units_ref="M2"/>
<units name="M3" units_ref="M3"/>
<units name="M3_per_s" units_ref="M3_per_s"/>
<units name="per_M3_s5" units_ref="per_M3_s5"/>
<units name="M_per_s" units_ref="M_per_s"/>
<units name="cm2" units_ref="cm2"/>
<units name="cm2_per_s" units_ref="cm2_per_s"/>
<units name="mol_per_m2" units_ref="mol_per_m2"/>
<units name="Farad" units_ref="Farad"/>
<units name="uF" units_ref="uF"/>
<units name="uF_per_cm2" units_ref="uF_per_cm2"/>
<units name="umol_per_cm2_s" units_ref="umol_per_cm2_s"/>
<units name="umol_per_s" units_ref="umol_per_s"/>
<units name="uA_per_cm2" units_ref="uA_per_cm2"/>
<units name="m_per_s" units_ref="m_per_s"/>
<units name="m3_mol_per_l_umol" units_ref="m3_mol_per_l_umol"/>
<units name="l_umol_per_mol_m3" units_ref="l_umol_per_mol_m3"/>
<units name="S" units_ref="S"/>
<units name="S_per_m2" units_ref="S_per_m2"/>
<units name="per_m2" units_ref="per_m2"/>
<units name="cm_per_s" units_ref="cm_per_s"/>
</import>
<import xlink:href="SGLT1.cellml">
<component component_ref="SGLT1" name="SGLT1"/>
</import>
<import xlink:href="GLUT2.cellml">
<component component_ref="GLUT2" name="GLUT2"/>
<component component_ref="GLUT2_Apical" name="GLUT2_Apical"/>
</import>
<import xlink:href="ENaC.cellml">
<component component_ref="ENaC" name="ENaC"/>
</import>
<import xlink:href="NaK.cellml">
<component component_ref="NaK_Pump" name="NaK_Pump"/>
</import>
<import xlink:href="NBC.cellml">
<component component_ref="NBC" name="NBC"/>
</import>
<import xlink:href="NKCC.cellml">
<component component_ref="NKCC" name="NKCC"/>
</import>
<import xlink:href="IK.cellml">
<component component_ref="IK" name="IK"/>
</import>
<import xlink:href="BK.cellml">
<component component_ref="BK" name="BK"/>
</import>
<import xlink:href="AE1.cellml">
<component component_ref="AE1" name="AE1"/>
</import>
<import xlink:href="NHE3.cellml">
<component component_ref="NHE3" name="NHE3"/>
</import>
<import xlink:href="CFTR.cellml">
<component component_ref="CFTR" name="CFTR"/>
</import>
<import xlink:href="Diffusion.cellml">
<component component_ref="Diffusion" name="Diffusion"/>
</import>
<component name="Apical_Membrane">
<variable initial_value="0.05" name="q_Gl_a" public_interface="out" units="M"/>
<variable initial_value="0.105" name="q_Na_a" public_interface="out" units="M"/>
<variable initial_value="0.106" name="q_Cl_a" public_interface="out" units="M"/>
<variable initial_value="0.004" name="q_K_a" public_interface="out" units="M"/>
<variable initial_value="-0.046" name="V_a" public_interface="out" units="volt"/>
<variable name="I_D_Na" public_interface="in" units="uA"/>
<variable name="I_D_K" public_interface="in" units="uA"/>
<variable name="I_D_Cl" 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_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"/>
<variable name="I_ENaC" public_interface="in" units="uA"/>
<variable name="I_Na_SGLT" public_interface="in" units="uA"/>
<variable name="I_CFTR" public_interface="in" units="uA"/>
<variable name="C" public_interface="in" units="uF"/>
<variable name="t" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>V_a</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_D_Na</ci>
<ci>I_D_K</ci>
</apply>
<ci>I_D_Cl</ci>
</apply>
<ci>I_Na_SGLT</ci>
</apply>
<ci>I_BK</ci>
</apply>
<ci>I_CFTR</ci>
<ci>I_ENaC</ci>
<ci>I_D_P_Na</ci>
<ci>I_D_P_K</ci>
</apply>
<ci>I_D_P_Cl</ci>
</apply>
<ci>C</ci>
</apply>
</apply>
</math>
</component>
<component name="Epithelial_Cell">
<variable initial_value="0.01" name="q_Gl_c" public_interface="out" units="M"/>
<variable initial_value="0.047" name="q_Na_c" public_interface="out" units="M"/>
<variable initial_value="0.057" name="q_Cl_c" public_interface="out" units="M"/>
<variable initial_value="0.129" name="q_K_c" public_interface="out" units="M"/>
<variable initial_value="6e-5" name="q_Ca_c" public_interface="out" units="M"/>
<variable name="q_H_int" public_interface="out" units="M"/>
<variable name="q_H_ext" public_interface="out" units="M"/>
<variable name="q_H_b" public_interface="out" units="M"/>
<variable initial_value="0.003" name="q_HCO3_a" public_interface="out" units="M"/>
<variable name="q_HCO3_c" public_interface="out" units="M"/>
<variable initial_value="0.041" name="q_HCO3_b" public_interface="out" units="M"/>
<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"/>
<variable name="J_NHE3_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_NHE3_H" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_GLUT" public_interface="in" units="umol_per_s"/>
<variable name="J_A_GLUT" public_interface="in" units="umol_per_s"/>
<variable name="J_D_K" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_K" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_Cl" public_interface="in" units="umol_per_s"/>
<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_CFTR_HCO3" public_interface="in" units="umol_per_s"/>
<variable name="J_ENaC" public_interface="in" units="umol_per_s"/>
<variable name="J_NKCC" 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"/>
<variable initial_value="7.2" name="pq_H_int" public_interface="out" units="per_s"/>
<variable initial_value="7.3" name="pq_H_ext" public_interface="out" units="per_s"/>
<variable initial_value="7.3" name="pq_H_b" public_interface="out" units="per_s"/>
<variable initial_value="1680000" name="rho" public_interface="out" units="per_m"/>
<variable name="beta" public_interface="out" units="M"/>
<variable initial_value="6.1" name="p_K" public_interface="out" units="dimensionless"/>
<variable initial_value="40" name="P_CO2" public_interface="out" units="dimensionless"/>
<variable initial_value="1e-5" 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 name="K_Concentration" public_interface="out" units="M"/>
<variable name="Na_Concentration" public_interface="out" units="M"/>
<variable name="Cl_Concentration" public_interface="out" units="M"/>
<variable name="Total_Gl_Flux_Apical" public_interface="out" units="umol_per_s"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="t" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Total_Gl_Flux_Apical</ci>
<apply>
<plus/>
<ci>J_A_GLUT</ci>
<ci>J_Gl_SGLT</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Na_Concentration</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_15</ci>
<ci>J_ENaC</ci>
</apply>
<apply>
<times/>
<ci>theta_3</ci>
<ci>J_NHE3_Na</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>K_Concentration</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>J_K_NaK</ci>
<ci>J_NKCC</ci>
</apply>
<ci>J_D_K</ci>
<ci>J_D_Bl_K</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0</cn>
<ci>J_IK</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Cl_Concentration</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<apply>
<minus/>
<apply>
<times/>
<ci>theta_11</ci>
<ci>J_CFTR</ci>
</apply>
<apply>
<times/>
<ci>theta_14</ci>
<ci>J_AE1_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2.3</cn>
<ci>q_HCO3_c</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Na_c</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mol_per_l_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<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_NKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_19</ci>
<ci>J_NBC_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_4</ci>
<ci>J_D_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_5</ci>
<ci>J_D_Bl_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_3</ci>
<ci>J_NHE3_Na</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Gl_c</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mol_per_l_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<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>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_K_c</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mol_per_l_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_8</ci>
<ci>J_K_NaK</ci>
</apply>
<apply>
<times/>
<ci>theta_17</ci>
<ci>J_NKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_9</ci>
<ci>J_D_K</ci>
</apply>
<apply>
<times/>
<ci>theta_10</ci>
<ci>J_D_Bl_K</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_27</ci>
<ci>J_IK</ci>
</apply>
</apply>
<ci>J_BK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Cl_c</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mol_per_l_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_11</ci>
<apply>
<minus/>
<ci>J_CFTR</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_18</ci>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_12</ci>
<ci>J_D_Cl</ci>
</apply>
<apply>
<times/>
<ci>theta_13</ci>
<ci>J_D_Bl_Cl</ci>
</apply>
<apply>
<times/>
<ci>theta_14</ci>
<ci>J_AE1_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>pq_H_int</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="m3_mol_per_l_umol" type="e-notation">1<sep/>-9</cn>
<ci>rho</ci>
</apply>
<apply>
<times/>
<ci>beta</ci>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<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/>
<ci>q_H_int</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pq_H_int</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_HCO3_c</ci>
<apply>
<times/>
<ci>s</ci>
<ci>P_CO2</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pq_H_int</ci>
<ci>p_K</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H_ext</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pq_H_ext</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H_b</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pq_H_b</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Basolateral_Membrane">
<variable initial_value="0" name="q_Gl_b" public_interface="out" units="M"/>
<variable initial_value="0.155" name="q_Na_b" public_interface="out" units="M"/>
<variable initial_value="0.119" name="q_Cl_b" public_interface="out" units="M"/>
<variable initial_value="0.0054" name="q_K_b" public_interface="out" units="M"/>
<variable name="I_D_Bl_Na" public_interface="in" units="uA"/>
<variable name="I_D_Bl_K" public_interface="in" units="uA"/>
<variable name="I_D_Bl_Cl" public_interface="in" units="uA"/>
<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="C" public_interface="in" units="uF"/>
<variable initial_value="-0.05" name="V_b" public_interface="out" units="volt"/>
<variable name="t" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>V_b</ci>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_D_Bl_Na</ci>
<ci>I_D_Bl_K</ci>
</apply>
<ci>I_D_Bl_Cl</ci>
</apply>
<ci>I_Na_NaK</ci>
</apply>
<ci>I_K_NaK</ci>
<ci>I_NBC_Na</ci>
</apply>
<ci>I_NBC_HCO3</ci>
</apply>
<ci>I_D_P_Na</ci>
</apply>
<ci>I_D_P_K</ci>
</apply>
<ci>I_D_P_Cl</ci>
</apply>
<ci>I_IK</ci>
</apply>
<ci>C</ci>
</apply>
</apply>
</math>
</component>
<component name="Blood">
<variable name="q_Gl_b" public_interface="in" units="M"/>
<variable initial_value="0.012" name="q_Gl_blood" public_interface="out" units="M"/>
<variable name="J_Gl_blood" public_interface="out" units="umol_per_s"/>
<variable initial_value="1e-6" name="K" public_interface="out" units="dimensionless"/>
<variable name="t" public_interface="in" units="second"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_Gl_blood</ci>
<apply>
<times/>
<ci>K</ci>
<apply>
<minus/>
<ci>q_Gl_b</ci>
<ci>q_Gl_blood</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Gl_blood</ci>
</apply>
<cn cellml:units="dimensionless">0</cn>
</apply>
</math>
</component>
<!-- Apical_Membrane-->
<!-- def comp SGLT1 as
var epsilon: per_s {pub: out};
var lambda: per_M3_s5 {pub: out};
var chi: M {pub: out};
var alpha: M3 {pub: out};
var beta: M2 {pub: out};
var gamma: M3_per_s {pub: out};
var phi: M_per_s {pub: out};
var k_12: per_s {pub: in};
var k_21: per_s {pub: in};
var k_32: per_s {pub: in};
var k_34: per_s {pub: in};
var k_43: per_s {pub: in};
var k_45: per_s {pub: in};
var k_54: per_s {pub: in};
var k_25: per_s {pub: in};
var k_52: per_s {pub: in};
var k_56: per_s {pub: in};
var k_65: per_s {pub: in};
var k_61: per_s {pub: in};
var k_16: per_s {pub: in};
var ks_12: per_M2_s {pub: in};
var k0_23: per_M_s {pub: in};
var F: C_per_mol {pub: in};
var q_Na_a: M {pub: in};
var q_Gl_a: M {pub: in};
var n_SGLT: dimensionless {init: 180000000, pub: out};
var J_Na_SGLT: umol_per_s {pub: out};
var I_Na_SGLT: uA {pub: out};
var J_Gl_SGLT: umol_per_s {pub: out};
var J_GL_SGLT1: umol_per_s {pub: out};-->
<!-- lambda = ks_12*k0_23*(k_43*k_54*k_65+k_34*k_54*k_65+k_45*k_52*k_65+k_43*k_52*k_65+k_34*k_52*k_65+k_34*k_45*k_65+k_45*k_56*k_61+k_43*k_56*k_61+k_34*k_56*k_61+k_43*k_54*k_61+k_34*k_54*k_61+k_45*k_52*k_61+k_43*k_52*k_61+k_34*k_52*k_61+k_34*k_45*k_61+k_34*k_45*k_56);
chi = 1{dimensionless}/lambda*ks_12*(k_32*k_43*k_54*k_65+k_25*k_43*k_54*k_65+k_25*k_34*k_54*k_65+k_25*k_32*k_54*k_65+k_34*k_45*k_52*k_65+k_32*k_45*k_52*k_65+k_32*k_43*k_52*k_65+k_25*k_34*k_45*k_65+k_25*k_32*k_45*k_65+k_25*k_32*k_43*k_65+k_34*k_45*k_56*k_61+k_32*k_45*k_56*k_61+k_32*k_43*k_56*k_61+k_32*k_43*k_54*k_61+k_25*k_43*k_54*k_61+k_25*k_34*k_54*k_61+k_25*k_32*k_54*k_61+k_34*k_45*k_52*k_61+k_32*k_45*k_52*k_61+k_32*k_43*k_52*k_61+k_25*k_34*k_45*k_61+k_25*k_32*k_45*k_61+k_25*k_32*k_43*k_61+k_25*k_34*k_45*k_56+k_25*k_32*k_45*k_56+k_25*k_32*k_43*k_56);
beta = 1{dimensionless}/lambda*k0_23*(k_16*k_43*k_54*k_65+k_16*k_34*k_54*k_65+k_16*k_45*k_52*k_65+k_16*k_43*k_52*k_65+k_16*k_34*k_52*k_65+k_16*k_34*k_45*k_65+k_34*k_45*k_56*k_61+k_16*k_34*k_45*k_56);
alpha = 1{dimensionless}/lambda*(k_21*k_32*k_43*k_54*k_65+k_16*k_32*k_43*k_54*k_65+k_16*k_25*k_43*k_54*k_65+k_16*k_21*k_43*k_54*k_65+k_16*k_25*k_34*k_54*k_65+k_16*k_21*k_34*k_54*k_65+k_16*k_25*k_32*k_54*k_65+k_16*k_21*k_32*k_54*k_65+k_21*k_34*k_45*k_52*k_65+k_16*k_34*k_45*k_52*k_65+k_21*k_32*k_45*k_52*k_65+k_16*k_32*k_45*k_52*k_65+k_21*k_32*k_43*k_52*k_65+k_16*k_32*k_43*k_52*k_65+k_16*k_25*k_34*k_45*k_65+k_16*k_21*k_34*k_45*k_65+k_16*k_25*k_32*k_45*k_65+k_16*k_21*k_32*k_45*k_65+k_16*k_25*k_32*k_43*k_65+k_16*k_21*k_32*k_43*k_65+k_25*k_34*k_45*k_56*k_61+k_21*k_34*k_45*k_56*k_61+k_25*k_32*k_45*k_56*k_61+k_21*k_32*k_45*k_56*k_61+k_25*k_32*k_43*k_56*k_61+k_21*k_32*k_43*k_56*k_61+k_21*k_32*k_43*k_54*k_61+k_21*k_34*k_45*k_52*k_61+k_21*k_32*k_45*k_52*k_61+k_21*k_32*k_43*k_52*k_61+k_16*k_25*k_34*k_45*k_56+k_16*k_21*k_34*k_45*k_56+k_16*k_25*k_32*k_45*k_56+k_16*k_21*k_32*k_45*k_56+k_16*k_25*k_32*k_43*k_56+k_16*k_21*k_32*k_43*k_56+k_16*k_21*k_32*k_43*k_54+k_16*k_21*k_34*k_45*k_52+k_16*k_21*k_32*k_45*k_52+k_16*k_21*k_32*k_43*k_52);
gamma = 1{dimensionless}/lambda*k_16*k_21*k_65*(k_32*k_43*k_54+k_34*k_45*k_52+k_32*k_45*k_52+k_32*k_43*k_52);
phi = -1{dimensionless}/lambda*ks_12*k_56*k_61*(k_25*k_34*k_45+k_25*k_32*k_45+k_25*k_32*k_43);
epsilon = -1{dimensionless}/lambda*ks_12*k0_23*k_34*k_45*k_56*k_61;
J_Na_SGLT = -2{dimensionless}*n_SGLT/6.023e17{per_umol}*(epsilon*sqr(q_Na_a)*q_Gl_a+phi*sqr(q_Na_a)+gamma)/(alpha+beta*q_Gl_a+chi*sqr(q_Na_a)+sqr(q_Na_a)*q_Gl_a);
I_Na_SGLT = J_Na_SGLT*F;
J_Gl_SGLT = 1{dimensionless}/2{dimensionless}*J_Na_SGLT;-->
<!-- //to change the unit to picomol/hour-->
<!-- J_GL_SGLT1 = J_Gl_SGLT*36e8{dimensionless};
enddef;-->
<!-- def comp GLUT2_Apical as
var K_12: per_M_s {init: 10, pub: out};
var K_21: per_s {init: 12.6, pub: out};
var K_23: per_s {init: 500, pub: out};
var K_32: per_s {init: 60, pub: out};
var K_34: per_s {init: 12.6, pub: out};
var K_43: per_M_s {init: 1, pub: out};
var K_41: per_s {init: 0.8, pub: out};
var K_14: per_s {init: 12, pub: out};
var q_Gl_c: M {pub: in};
var q_Gl_a: M {pub: in};
var K_S: M {pub: out};
var J_A_GLUT: umol_per_s {pub: out};
var C_tot: mole {pub: out};
var n_GLUT: dimensionless {init: 90000000, pub: out};
var D: per_s {pub: out};
var J_A_GLUT1: umol_per_s {pub: out};
var v_cell: m3 {pub: in};-->
<!-- C_tot = n_GLUT/6.022e23{per_mole};
K_S = K_21/K_12;
J_A_GLUT = 1e6{dimensionless}*C_tot/D*(K_23*K_41*q_Gl_a/K_S-K_32*K_14*q_Gl_c/K_S);
J_A_GLUT1 = J_A_GLUT*36e8{dimensionless};
D = (1{dimensionless}+q_Gl_a/K_S)*(K_32*q_Gl_c/K_S+K_41)+(1{dimensionless}+q_Gl_c/K_S)*(K_23*q_Gl_a/K_S+K_14);
enddef;-->
<!-- def comp NHE3 as
var A_Apical: m2 {pub: in};
var v_cell: m3 {pub: in};
var q_Na_a: M {pub: in};
var q_Na_c: M {pub: in};
var q_H_ext: M {pub: in};
var q_H_int: M {pub: in};
var NH4_ext: M {init: 0, pub: out};
var NH4_int: M {init: 0, pub: out};
var J_NHE3_Na: umol_per_s {pub: out};
var J_NHE3_H1: umol_per_s {pub: out};
var J_NHE3_H: umol_per_s {pub: out};
var J_NHE3_NH4: umol_per_s {pub: out};
var J_NHE3_Na_Max: umol_per_s {pub: out};
var I_NHE3_H: uA {pub: out};
var I_NHE3_Na: uA {pub: out};
var x_T: M {init: 0.002};
var sigma: per_s;
var P_Na: per_s;
var P_H: per_s;
var P_NH4: per_s;
var P0_Na: per_s {init: 16};
var P0_H: per_s {init: 4.8};
var P0_NH4: per_s {init: 16};
var K_Na: M {init: 0.03};
var K_H: M {init: 7.2e-8};
var K_NH4: M {init: 0.027};
var K_I: M {init: 1e-9};
var f_m: dimensionless {init: 0};
var f_M: dimensionless {init: 2};
var alpha_ext_Na: dimensionless;
var alpha_int_Na: dimensionless;
var beta_ext_H: dimensionless;
var beta_int_H: dimensionless;
var gamma_ext_NH4: dimensionless;
var gamma_int_NH4: dimensionless;
var F: C_per_mol {pub: in};-->
<!-- alpha_ext_Na = q_Na_a/K_Na;
alpha_int_Na = q_Na_c/K_Na;
beta_ext_H = q_H_ext/K_H;
beta_int_H = q_H_int/K_H;
gamma_ext_NH4 = NH4_ext/K_NH4;
gamma_int_NH4 = NH4_int/K_NH4;
P_Na = P0_Na*(f_M*q_H_int+f_m*K_I)/(q_H_int+K_I);
P_H = P0_H*(f_M*q_H_int+f_m*K_I)/(q_H_int+K_I);
P_NH4 = P0_NH4*(f_M*q_H_int+f_m*K_I)/(q_H_int+K_I);
sigma = (1{dimensionless}+alpha_ext_Na+beta_ext_H+gamma_ext_NH4)*(P_Na*alpha_int_Na+P_H*beta_int_H+P_NH4*gamma_int_NH4)+(1{dimensionless}+alpha_int_Na+beta_int_H+gamma_int_NH4)*(P_Na*alpha_ext_Na+P_H*beta_ext_H+P_NH4*gamma_ext_NH4);
J_NHE3_Na = -1e9{dimensionless}*v_cell*x_T/sigma*(P_Na*P_H*(alpha_int_Na*beta_ext_H-alpha_ext_Na*beta_int_H)+P_Na*P_NH4*(alpha_int_Na*gamma_ext_NH4-alpha_ext_Na*gamma_int_NH4));
J_NHE3_H = 1e9{dimensionless}*v_cell*x_T/sigma*(P_Na*P_H*(alpha_ext_Na*beta_int_H-alpha_int_Na*beta_ext_H)+P_H*P_NH4*(beta_int_H*gamma_ext_NH4-beta_ext_H*gamma_int_NH4));
J_NHE3_H1 = J_NHE3_H*36e8{dimensionless};
J_NHE3_NH4 = x_T/sigma*(P_Na*P_NH4*(alpha_ext_Na*gamma_int_NH4-alpha_int_Na*gamma_ext_NH4)+P_H*P_NH4*(beta_ext_H*gamma_int_NH4-gamma_ext_NH4*beta_int_H));
J_NHE3_Na_Max = x_T*P_Na*P_H/(P_Na+P_H);
I_NHE3_H = J_NHE3_H*F;
I_NHE3_Na = J_NHE3_Na*F;
enddef;-->
<!-- def comp AE1 as
var A_Apical: m2 {pub: in};
var v_cell: m3 {pub: in};
var q_HCO3_c: M {pub: in};
var q_HCO3_a: M {pub: in};
var q_Cl_c: M {pub: in};
var q_Cl_a: M {pub: in};
var J_AE1_HCO3: umol_per_s {pub: out};
var J_AE1_Cl: umol_per_s {pub: out};
var I_AE1_HCO3: uA {pub: out};
var I_AE1_Cl: uA {pub: out};
var J_q_HCO3_cnflux: umol_per_s {pub: out};
var J_q_Cl_cnflux: umol_per_s {pub: out};
var K_q_HCO3_a: M {init: 0.198};
var K_q_HCO3_c: M {init: 0.198};
var K_q_Cl_b: M {init: 0.05};
var K_q_Cl_c: M {init: 0.05};
var P_q_HCO3_a: per_s {init: 1247};
var P_q_HCO3_c: per_s {init: 135};
var P_q_Cl_b: per_s {init: 562};
var P_q_Cl_c: per_s {init: 61};
var beta_ext: dimensionless;
var beta_int: dimensionless;
var gamma_ext: dimensionless;
var gamma_int: dimensionless;
var sigma: per_s;
var x_Tmax: M {init: 0.001};
var K_I: M {init: 0.172};
var x_T: M;
var x_ext: M;
var x_int: M;
var Jo_bm: umol_per_s;
var Ji_bm: umol_per_s;
var Js_bm: umol_per_s;
var Jo_cm: umol_per_s;
var Ji_cm: umol_per_s;
var Js_cm: umol_per_s;
var F: C_per_mol {pub: in};-->
<!-- x_T = x_Tmax/(1{dimensionless}+q_HCO3_c/K_I);
beta_ext = q_HCO3_a/K_q_HCO3_a;
beta_int = q_HCO3_c/K_q_HCO3_c;
gamma_ext = q_Cl_a/K_q_Cl_b;
gamma_int = q_Cl_c/K_q_Cl_c;
sigma = (1{dimensionless}+beta_ext+gamma_ext)*(P_q_HCO3_c*beta_int+P_q_Cl_c*gamma_int)+(1{dimensionless}+beta_int+gamma_int)*(P_q_HCO3_a*beta_ext+P_q_Cl_b*gamma_ext);
J_AE1_HCO3 = 1e9{dimensionless}*v_cell*x_Tmax/sigma*(P_q_HCO3_c*beta_int*P_q_Cl_b*gamma_ext-P_q_HCO3_a*beta_ext*P_q_Cl_c*gamma_int);
J_AE1_Cl = J_AE1_HCO3;
I_AE1_HCO3 = J_AE1_HCO3*F;
I_AE1_Cl = J_AE1_Cl*F;
x_ext = x_T*(P_q_HCO3_c*beta_int+P_q_Cl_c*gamma_int)/sigma;
x_int = x_T*(P_q_HCO3_a*beta_ext+P_q_Cl_b*gamma_ext)/sigma;
J_q_HCO3_cnflux = x_T/sigma*P_q_HCO3_a*beta_ext*(P_q_HCO3_c*beta_int+P_q_Cl_c*gamma_int);
J_q_Cl_cnflux = x_T/sigma*P_q_Cl_b*gamma_ext*(P_q_HCO3_c*beta_int+P_q_Cl_c*gamma_int);
Jo_bm = pow(1{dimensionless}/x_T*(1{dimensionless}/P_q_HCO3_a+1{dimensionless}/P_q_HCO3_c+K_q_HCO3_c/(P_q_HCO3_c*q_HCO3_c)), -1{dimensionless});
Ji_bm = pow(1{dimensionless}/x_T*(1{dimensionless}/P_q_HCO3_a+1{dimensionless}/P_q_HCO3_c+K_q_HCO3_a/(P_q_HCO3_a*q_HCO3_a)), -1{dimensionless});
Js_bm = pow(1{dimensionless}/x_T*(1{dimensionless}/P_q_HCO3_a+1{dimensionless}/P_q_HCO3_c), -1{dimensionless});
Jo_cm = pow(1{dimensionless}/x_T*(1{dimensionless}/P_q_Cl_b+1{dimensionless}/P_q_Cl_c+K_q_Cl_c/(P_q_Cl_c*q_Cl_c)), -1{dimensionless});
Ji_cm = pow(1{dimensionless}/x_T*(1{dimensionless}/P_q_Cl_b+1{dimensionless}/P_q_Cl_c+K_q_Cl_b/(P_q_Cl_b*q_Cl_a)), -1{dimensionless});
Js_cm = pow(1{dimensionless}/x_T*(1{dimensionless}/P_q_Cl_b+1{dimensionless}/P_q_Cl_c), -1{dimensionless});
enddef;-->
<!-- def comp CFTR as
var z_Cl: dimensionless {pub: in};
var R: J_per_mol_K {pub: in};
var T: kelvin {pub: in};
var q_Cl_c: M {pub: in};
var q_Cl_a: M {pub: in};
var V_A_Cl: volt {pub: out};
var V_A_HCO3: volt {pub: out};
var V_a: volt {pub: in};
var J_CFTR: umol_per_s {pub: out};
var J_CFTR1: umol_per_s {pub: out};
var J_CFTR_HCO3: umol_per_s {pub: out};
var I_CFTR: uA {pub: out};
var I_CFTR_HCO3: uA {pub: out};
var A_Apical: m2 {pub: in};
var F: C_per_mol {pub: in};
var G_CFTR: S {init: 17, pub: out};
var q_HCO3_c: M {pub: in};
var q_HCO3_a: 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};-->
<!-- J_CFTR = I_CFTR/F;
J_CFTR1 = J_CFTR*36e8{dimensionless};
V_A_Cl = R*T/(z_Cl*F)*ln(q_Cl_a/q_Cl_c);
I_CFTR = -1e6{dimensionless}*G_CFTR*alpha_CFTR*A_Apical*(V_a-V_A_Cl);
V_A_HCO3 = R*T/(z_HCO3*F)*ln(q_HCO3_a/q_HCO3_c);
I_CFTR_HCO3 = 0.25e6{dimensionless}*G_CFTR*alpha_CFTR*A_Apical*(V_a-V_A_HCO3);
J_CFTR_HCO3 = I_CFTR_HCO3/F;
enddef;-->
<!-- def comp ENaC as
var z_Na: dimensionless {pub: in};
var R: J_per_mol_K {pub: in};
var T: kelvin {pub: in};
var q_Na_c: M {pub: in};
var q_Na_a: M {pub: in};
var V_A_Na: volt {pub: out};
var V_a: volt {pub: in};
var J_ENaC: umol_per_s {pub: out};
var J_ENaC1: umol_per_s {pub: out};
var I_ENaC: uA {pub: out};
var A_Apical: m2 {pub: in};
var F: C_per_mol {pub: in};
var G_ENaC: S {init: 9, pub: out};-->
<!-- // 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};-->
<!-- J_ENaC = I_ENaC/F;
J_ENaC1 = J_ENaC*36e8{dimensionless};
V_A_Na = R*T/(z_Na*F)*ln(q_Na_a/q_Na_c);
I_ENaC = -1e+6{dimensionless}*G_ENaC*alpha_ENaC*A_Apical*(V_a-V_A_Na);
enddef;-->
<!-- def comp BK as
var z_k: dimensionless {pub: in};
var R: J_per_mol_K {pub: in};
var T: kelvin {pub: in};
var q_K_c: M {pub: in};
var q_K_a: M {pub: in};
var q_Ca_c: M {pub: in};
var V_A_BK: volt {pub: out};
var V_a: volt {pub: in};
var J_BK: umol_per_s {pub: out};
var I_BK: uA {pub: out};
var A_Apical: m2 {pub: in};
var F: C_per_mol {pub: in};
var G_BK: S {init: 30, pub: out};
var P_BK: dimensionless {pub: out};
var q_K_b_BK: M {init: 0.0003, pub: out};
var eta: dimensionless {init: 1.7, pub: out};-->
<!-- // 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};-->
<!-- P_BK = 1{dimensionless}/(1{dimensionless}+pow(q_K_b_BK/q_Ca_c, eta));
V_A_BK = R*T/(z_k*F)*ln(q_K_a/q_K_c);
I_BK = 1e+6{dimensionless}*G_BK*alpha_BK*P_BK*A_Apical*(V_a-V_A_BK);
J_BK = I_BK/F;
enddef;-->
<!-- // Basolateral_Membrane-->
<!-- def comp GLUT2 as
var K_12: per_M_s {init: 10, pub: out};
var K_21: per_s {init: 12.6, pub: out};
var K_23: per_s {init: 500, pub: out};
var K_32: per_s {init: 60, pub: out};
var K_34: per_s {init: 12.6, pub: out};
var K_43: per_M_s {init: 1, pub: out};
var K_41: per_s {init: 0.8, pub: out};
var K_14: per_s {init: 12, pub: out};
var q_Gl_c: M {pub: in};
var q_Gl_b: M {pub: in};
var K_S: M {pub: out};
var J_GLUT: umol_per_s {pub: out};
var C_tot: mole {pub: out};
var n_GLUT: dimensionless {init: 15000000, pub: out};
var D: per_s {pub: out};
var J_GLUT1: umol_per_s {pub: out};
var v_cell: m3 {pub: in};-->
<!-- C_tot = n_GLUT/6.022e23{per_mole};
K_S = K_21/K_12;
J_GLUT = 1e6{dimensionless}*C_tot/D*(K_23*K_41*q_Gl_c/K_S-K_32*K_14*q_Gl_b/K_S);
J_GLUT1 = J_GLUT*36e8{dimensionless};
D = (1{dimensionless}+q_Gl_c/K_S)*(K_32*q_Gl_b/K_S+K_41)+(1{dimensionless}+q_Gl_b/K_S)*(K_23*q_Gl_c/K_S+K_14);
enddef;-->
<!-- def comp NaK_Pump as
var t: second {pub: in};
var V_b: volt {pub: in};
var F: C_per_mol {pub: in};
var q_K_b: M {pub: in};
var q_Na_c: M {pub: in};
var TO_NaK: per_s {pub: out};
var n_NaK: dimensionless {init: 5000000, pub: out};
var K_q_Na_c_half: M {init: 0.01, pub: out};
var K_q_K_b_half: M {init: 0.0015, pub: out};
var I_Na_NaK: uA {pub: out};
var I_K_NaK: uA {pub: out};
var J_NaK: umol_per_s {pub: out};
var J_Na_NaK: umol_per_s {pub: out};
var J_K_NaK: umol_per_s {pub: out};-->
<!-- I_K_NaK = J_K_NaK*F;
I_Na_NaK = J_Na_NaK*F;
J_NaK = TO_NaK*n_NaK/6.023e17{per_umol}*pow(q_Na_c, 1.36{dimensionless})*1.13{dimensionless}/(pow(K_q_Na_c_half, 1.36{dimensionless})+pow(q_Na_c, 1.36{dimensionless}))*1.3{dimensionless}*q_K_b/(K_q_K_b_half+q_K_b);
TO_NaK = 5.46e+4{per_V4_s}*pow(V_b, 4{dimensionless})-2.43e+2{per_V3_s}*pow(V_b, 3{dimensionless})-2e+3{per_V2_s}*sqr(V_b)+160{per_V_s}*V_b+51{per_s};
J_Na_NaK = 3{dimensionless}*J_NaK;
J_K_NaK = 2{dimensionless}*J_NaK;
enddef;-->
<!-- def comp NBC as
var t: second {pub: in};
var A_Apical: m2 {pub: in};
var V_b: volt {pub: in};
var v_cell: m3 {pub: in};
var q_Na_c: M {pub: in};
var q_Na_b: M {pub: in};
var q_HCO3_c: M {pub: in};
var q_HCO3_b: M {init: 0.003, pub: out};
var g_NBC: S_per_m2 {init: 0.01, pub: out};
var F: C_per_mol {pub: in};
var R: J_per_mol_K {pub: in};
var T: kelvin {pub: in};
var n: dimensionless {init: 2};
var E_NBC: volt {pub: out};
var J_NBC_Na: umol_per_s {pub: out};
var J_NBC_HCO3: umol_per_s {pub: out};
var I_NBC_Na: uA {pub: out};
var I_NBC_HCO3: uA {pub: out};
var eps: M {init: 1e-18};-->
<!-- E_NBC = R*T/((1{dimensionless}-n)*F)*ln(q_Na_b/q_Na_c*sqr(q_HCO3_b/q_HCO3_c));
J_NBC_Na = 1e6{dimensionless}*g_NBC*A_Apical*(E_NBC-V_b)/F;
J_NBC_HCO3 = 3.0{dimensionless}*J_NBC_Na;
I_NBC_Na = J_NBC_Na*F;
I_NBC_HCO3 = J_NBC_HCO3*F;
enddef;-->
<!-- def comp NKCC as
var A_Apical: m2 {pub: in};
var S_C: M {pub: out};
var J_NKCC: umol_per_s {pub: out};
var I_NKCC: uA {pub: out};
var q_Na_c: M {pub: in};
var q_K_c: M {pub: in};
var q_Cl_c: M {pub: in};
var F: C_per_mol {pub: in};
var u_NKCC: per_s {pub: out};
var gamma_1: per_s2 {init: 157.55, pub: out};
var gamma_2: per_M4_s2 {init: 20000000, pub: out};
var gamma_3: per_s {init: 1.03, pub: out};
var gamma_4: per_M4_s {init: 1380000, pub: out};
var alpha_NKCC: mol_per_m2 {init: 1e-9, pub: out};-->
<!-- S_C = q_Na_c*q_Cl_c*q_Cl_c*q_K_c;
u_NKCC = (gamma_1-gamma_2*S_C)/(gamma_3+gamma_4*S_C);
J_NKCC = -1e6{dimensionless}*alpha_NKCC*u_NKCC*A_Apical;
I_NKCC = J_NKCC*F;
enddef;-->
<!-- def comp IK as
var z_k: dimensionless {pub: in};
var R: J_per_mol_K {pub: in};
var T: kelvin {pub: in};
var q_K_c: M {pub: in};
var q_K_b: M {pub: in};
var q_Ca_c: M {pub: in};
var V_A_IK: volt {pub: out};
var V_b: volt {pub: in};
var J_IK: umol_per_s {pub: out};
var I_IK: uA {pub: out};
var A_Apical: m2 {pub: in};
var F: C_per_mol {pub: in};
var G_IK: S {init: 24, pub: out};
var P_IK: dimensionless {pub: out};
var q_K_b_IK: M {init: 0.0003, pub: out};
var eta: dimensionless {init: 1.7, pub: out};-->
<!-- // 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(q_K_b_IK/q_Ca_c, eta));
V_A_IK = R*T/(z_k*F)*ln(q_K_b/q_K_c);
I_IK = 1e+6{dimensionless}*G_IK*alpha_IK*P_IK*A_Apical*(V_b-V_A_IK);
J_IK = I_IK/F;
enddef;-->
<component name="Apical_Diffusion">
<variable name="V_a" public_interface="in" units="volt"/>
<variable name="q_Na_a" public_interface="in" units="M"/>
<variable name="q_K_a" public_interface="in" units="M"/>
<variable name="q_Cl_a" public_interface="in" units="M"/>
<variable name="q_Na_c" public_interface="in" units="M"/>
<variable name="q_K_c" public_interface="in" units="M"/>
<variable name="q_Cl_c" public_interface="in" units="M"/>
<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_mol_K"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_Na" public_interface="in" units="m_per_s"/>
<variable name="P_k" public_interface="in" units="m_per_s"/>
<variable name="P_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_D_NA" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Na" public_interface="out" units="uA"/>
<variable name="J_D_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_K" public_interface="out" units="uA"/>
<variable name="J_D_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Apical" public_interface="out" units="umol_per_s"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_D_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Na</ci>
<ci>F</ci>
<ci>z_Na</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_Na_a</ci>
<apply>
<times/>
<ci>q_Na_c</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_NA</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
<ci>J_D_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Na</ci>
<apply>
<times/>
<ci>J_D_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_k</ci>
<ci>F</ci>
<ci>z_k</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_K_a</ci>
<apply>
<times/>
<ci>q_K_c</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_K</ci>
<apply>
<times/>
<ci>J_D_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Cl</ci>
<ci>F</ci>
<ci>z_Cl</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_Cl_a</ci>
<apply>
<times/>
<ci>q_Cl_c</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Cl</ci>
<apply>
<times/>
<ci>J_D_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Apical</ci>
<apply>
<plus/>
<ci>J_D_Na</ci>
<ci>J_D_K</ci>
<ci>J_D_Cl</ci>
</apply>
</apply>
</math>
</component>
<component name="Basolateral_Diffusion">
<variable name="V_b" public_interface="in" units="volt"/>
<variable name="q_Na_b" public_interface="in" units="M"/>
<variable name="q_K_b" public_interface="in" units="M"/>
<variable name="q_Cl_b" public_interface="in" units="M"/>
<variable name="q_Na_c" public_interface="in" units="M"/>
<variable name="q_K_c" public_interface="in" units="M"/>
<variable name="q_Cl_c" public_interface="in" units="M"/>
<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_mol_K"/>
<variable name="A_Basol" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_Bl_Na" public_interface="in" units="m_per_s"/>
<variable name="P_Bl_k" public_interface="in" units="m_per_s"/>
<variable name="P_Bl_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_Bl_Na" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_Na" public_interface="out" units="uA"/>
<variable name="J_D_Bl_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_K" public_interface="out" units="uA"/>
<variable name="J_D_Bl_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Basolateral" public_interface="out" units="umol_per_s"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_D_Bl_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_Na</ci>
<ci>F</ci>
<ci>z_Na</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_Na_b</ci>
<apply>
<times/>
<ci>q_Na_c</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_Na</ci>
<apply>
<times/>
<ci>J_D_Bl_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Bl_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_k</ci>
<ci>F</ci>
<ci>z_k</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_K_b</ci>
<apply>
<times/>
<ci>q_K_c</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_K</ci>
<apply>
<times/>
<ci>J_D_Bl_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Bl_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_Cl</ci>
<ci>F</ci>
<ci>z_Cl</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_Cl_b</ci>
<apply>
<times/>
<ci>q_Cl_c</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>V_b</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_Cl</ci>
<apply>
<times/>
<ci>J_D_Bl_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Basolateral</ci>
<apply>
<plus/>
<ci>J_D_Bl_Na</ci>
<ci>J_D_Bl_K</ci>
<ci>J_D_Bl_Cl</ci>
</apply>
</apply>
</math>
</component>
<component name="Paracellular_Pathway">
<variable name="V_ab" public_interface="out" units="volt"/>
<variable name="q_Na_a" public_interface="in" units="M"/>
<variable name="q_K_a" public_interface="in" units="M"/>
<variable name="q_Cl_a" public_interface="in" units="M"/>
<variable name="q_Na_b" public_interface="in" units="M"/>
<variable name="q_K_b" public_interface="in" units="M"/>
<variable name="q_Cl_b" public_interface="in" units="M"/>
<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_mol_K"/>
<variable name="A_Apical" 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_b" public_interface="in" units="volt"/>
<variable name="V_a" public_interface="in" units="volt"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_ab</ci>
<apply>
<minus/>
<ci>V_a</ci>
<ci>V_b</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_P_Na</ci>
<ci>F</ci>
<ci>z_Na</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_Na_a</ci>
<apply>
<times/>
<ci>q_Na_b</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_Na</ci>
<apply>
<times/>
<ci>J_D_P_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_P_k</ci>
<ci>F</ci>
<ci>z_k</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_K_a</ci>
<apply>
<times/>
<ci>q_K_b</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_K</ci>
<apply>
<times/>
<ci>J_D_P_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="l_umol_per_mol_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_P_Cl</ci>
<ci>F</ci>
<ci>z_Cl</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>q_Cl_a</ci>
<apply>
<times/>
<ci>q_Cl_b</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>V_ab</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_Cl</ci>
<apply>
<times/>
<ci>J_D_P_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<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>
<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_mol_K"/>
<variable initial_value="310" name="T" public_interface="out" units="kelvin"/>
<variable initial_value="1e-5" name="C" public_interface="out" units="uF"/>
<variable initial_value="8.4e-10" name="A_Apical" public_interface="out" units="m2"/>
<variable initial_value="5e-10" name="A_Basol" public_interface="out" units="m2"/>
<variable initial_value="5e-16" name="v_cell" public_interface="out" 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-6" name="P_P_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="1e-6" name="P_P_k" public_interface="out" units="m_per_s"/>
<variable initial_value="3e-8" 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="t" public_interface="out" units="second"/>
<variable initial_value="14000" name="k0_12" public_interface="out" units="per_M2_s"/>
<variable initial_value="300" name="k0_21" public_interface="out" units="per_s"/>
<variable initial_value="15000" name="k0_23" public_interface="out" units="per_M_s"/>
<variable initial_value="20" name="k0_32" public_interface="out" units="per_s"/>
<variable initial_value="0.1" name="k0_25" public_interface="out" units="per_s"/>
<variable initial_value="0.01" name="k0_52" public_interface="out" units="per_s"/>
<variable initial_value="50" name="k0_34" public_interface="out" units="per_s"/>
<variable initial_value="50" name="k0_43" public_interface="out" units="per_s"/>
<variable initial_value="6" name="k0_45" public_interface="out" units="per_s"/>
<variable initial_value="400" name="k0_54" public_interface="out" units="per_M_s"/>
<variable initial_value="10" name="k0_56" public_interface="out" units="per_s"/>
<variable initial_value="200" name="k0_65" public_interface="out" units="per_M2_s"/>
<variable initial_value="25" name="k0_61" public_interface="out" units="per_s"/>
<variable initial_value="600" name="k0_16" public_interface="out" units="per_s"/>
<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"/>
<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_a" public_interface="in" units="volt"/>
<variable name="mu" public_interface="out" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>F</ci>
<ci>V_a</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Rate_Constants">
<variable name="q_Na_c" public_interface="in" units="M"/>
<variable name="q_Na_a" public_interface="in" units="M"/>
<variable name="q_Gl_c" public_interface="in" units="M"/>
<variable name="q_Gl_a" public_interface="in" units="M"/>
<variable name="k_12" public_interface="out" units="per_s"/>
<variable name="ks_12" public_interface="out" units="per_M2_s"/>
<variable name="k_21" public_interface="out" units="per_s"/>
<variable name="k_32" public_interface="out" units="per_s"/>
<variable name="k_25" public_interface="out" units="per_s"/>
<variable name="k_52" public_interface="out" units="per_s"/>
<variable name="k_34" public_interface="out" units="per_s"/>
<variable name="k_43" public_interface="out" units="per_s"/>
<variable name="k_45" public_interface="out" units="per_s"/>
<variable name="k_54" public_interface="out" units="per_s"/>
<variable name="k_56" public_interface="out" units="per_s"/>
<variable name="k_65" public_interface="out" units="per_s"/>
<variable name="k_61" public_interface="out" units="per_s"/>
<variable name="k_16" public_interface="out" units="per_s"/>
<variable name="k0_12" public_interface="in" units="per_M2_s"/>
<variable name="k0_21" public_interface="in" units="per_s"/>
<variable name="k0_32" public_interface="in" units="per_s"/>
<variable name="k0_25" public_interface="in" units="per_s"/>
<variable name="k0_52" public_interface="in" units="per_s"/>
<variable name="k0_34" public_interface="in" units="per_s"/>
<variable name="k0_43" public_interface="in" units="per_s"/>
<variable name="k0_45" public_interface="in" units="per_s"/>
<variable name="k0_54" public_interface="in" units="per_M_s"/>
<variable name="k0_56" public_interface="in" units="per_s"/>
<variable name="k0_65" public_interface="in" units="per_M2_s"/>
<variable name="k0_61" public_interface="in" units="per_s"/>
<variable name="k0_16" public_interface="in" units="per_s"/>
<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"/>
<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>
<apply>
<eq/>
<ci>k_12</ci>
<apply>
<times/>
<ci>k0_12</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>alpha_p</ci>
</apply>
<ci>mu</ci>
</apply>
</apply>
<apply>
<power/>
<ci>q_Na_a</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<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>
<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>q_Gl_c</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_25</ci>
<ci>k0_25</ci>
</apply>
<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>q_Na_c</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>
</component>
<!-- Mapping Apical Membrane-->
<connection>
<map_components component_1="Apical_Membrane" component_2="SGLT1"/>
<map_variables variable_1="q_Na_a" variable_2="q_Na_a"/>
<map_variables variable_1="q_Gl_a" variable_2="q_Gl_a"/>
<map_variables variable_1="I_Na_SGLT" variable_2="I_Na_SGLT"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="GLUT2_Apical"/>
<map_variables variable_1="q_Gl_a" variable_2="q_Gl_a"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="NHE3"/>
<map_variables variable_1="I_NHE3_H" variable_2="I_NHE3_H"/>
<map_variables variable_1="I_NHE3_Na" variable_2="I_NHE3_Na"/>
<map_variables variable_1="q_Na_a" variable_2="q_Na_a"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="AE1"/>
<map_variables variable_1="q_Cl_a" variable_2="q_Cl_a"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="CFTR"/>
<map_variables variable_1="q_Cl_a" variable_2="q_Cl_a"/>
<map_variables variable_1="V_a" variable_2="V_a"/>
<map_variables variable_1="I_CFTR" variable_2="I_CFTR"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="ENaC"/>
<map_variables variable_1="q_Na_a" variable_2="q_Na_a"/>
<map_variables variable_1="V_a" variable_2="V_a"/>
<map_variables variable_1="I_ENaC" variable_2="I_ENaC"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="BK"/>
<map_variables variable_1="q_K_a" variable_2="q_K_a"/>
<map_variables variable_1="V_a" variable_2="V_a"/>
<map_variables variable_1="I_BK" variable_2="I_BK"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="Apical_Diffusion"/>
<map_variables variable_1="I_D_Na" variable_2="I_D_Na"/>
<map_variables variable_1="I_D_K" variable_2="I_D_K"/>
<map_variables variable_1="I_D_Cl" variable_2="I_D_Cl"/>
<map_variables variable_1="V_a" variable_2="V_a"/>
<map_variables variable_1="q_Na_a" variable_2="q_Na_a"/>
<map_variables variable_1="q_K_a" variable_2="q_K_a"/>
<map_variables variable_1="q_Cl_a" variable_2="q_Cl_a"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="Paracellular_Pathway"/>
<map_variables variable_1="q_Na_a" variable_2="q_Na_a"/>
<map_variables variable_1="q_K_a" variable_2="q_K_a"/>
<map_variables variable_1="q_Cl_a" variable_2="q_Cl_a"/>
<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"/>
<map_variables variable_1="V_a" variable_2="V_a"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="Parameters"/>
<map_variables variable_1="C" variable_2="C"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="V_a" variable_2="V_a"/>
</connection>
<connection>
<map_components component_1="Apical_Membrane" component_2="Rate_Constants"/>
<map_variables variable_1="q_Na_a" variable_2="q_Na_a"/>
<map_variables variable_1="q_Gl_a" variable_2="q_Gl_a"/>
</connection>
<!-- Mapping Epithelial Cell-->
<connection>
<map_components component_1="Epithelial_Cell" component_2="SGLT1"/>
<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="Epithelial_Cell" component_2="GLUT2"/>
<map_variables variable_1="q_Gl_c" variable_2="q_Gl_c"/>
<map_variables variable_1="J_GLUT" variable_2="J_GLUT"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="GLUT2_Apical"/>
<map_variables variable_1="q_Gl_c" variable_2="q_Gl_c"/>
<map_variables variable_1="J_A_GLUT" variable_2="J_A_GLUT"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="NHE3"/>
<map_variables variable_1="q_Na_c" variable_2="q_Na_c"/>
<map_variables variable_1="q_H_int" variable_2="q_H_int"/>
<map_variables variable_1="q_H_ext" variable_2="q_H_ext"/>
<map_variables variable_1="J_NHE3_Na" variable_2="J_NHE3_Na"/>
<map_variables variable_1="J_NHE3_H" variable_2="J_NHE3_H"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="AE1"/>
<map_variables variable_1="q_Cl_c" variable_2="q_Cl_c"/>
<map_variables variable_1="q_HCO3_c" variable_2="q_HCO3_c"/>
<map_variables variable_1="q_HCO3_a" variable_2="q_HCO3_a"/>
<map_variables variable_1="J_AE1_Cl" variable_2="J_AE1_Cl"/>
<map_variables variable_1="J_AE1_HCO3" variable_2="J_AE1_HCO3"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="NaK_Pump"/>
<map_variables variable_1="q_Na_c" variable_2="q_Na_c"/>
<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="Epithelial_Cell" component_2="NBC"/>
<map_variables variable_1="q_Na_c" variable_2="q_Na_c"/>
<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="q_HCO3_c" variable_2="q_HCO3_c"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="CFTR"/>
<map_variables variable_1="q_Cl_c" variable_2="q_Cl_c"/>
<map_variables variable_1="q_HCO3_c" variable_2="q_HCO3_c"/>
<map_variables variable_1="q_HCO3_a" variable_2="q_HCO3_a"/>
<map_variables variable_1="J_CFTR" variable_2="J_CFTR"/>
<map_variables variable_1="J_CFTR_HCO3" variable_2="J_CFTR_HCO3"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="ENaC"/>
<map_variables variable_1="q_Na_c" variable_2="q_Na_c"/>
<map_variables variable_1="J_ENaC" variable_2="J_ENaC"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="NKCC"/>
<map_variables variable_1="q_Na_c" variable_2="q_Na_c"/>
<map_variables variable_1="q_Cl_c" variable_2="q_Cl_c"/>
<map_variables variable_1="q_K_c" variable_2="q_K_c"/>
<map_variables variable_1="J_NKCC" variable_2="J_NKCC"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="IK"/>
<map_variables variable_1="q_K_c" variable_2="q_K_c"/>
<map_variables variable_1="q_Ca_c" variable_2="q_Ca_c"/>
<map_variables variable_1="J_IK" variable_2="J_IK"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="BK"/>
<map_variables variable_1="q_K_c" variable_2="q_K_c"/>
<map_variables variable_1="q_Ca_c" variable_2="q_Ca_c"/>
<map_variables variable_1="J_BK" variable_2="J_BK"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="Apical_Diffusion"/>
<map_variables variable_1="q_Na_c" variable_2="q_Na_c"/>
<map_variables variable_1="q_K_c" variable_2="q_K_c"/>
<map_variables variable_1="q_Cl_c" variable_2="q_Cl_c"/>
<map_variables variable_1="J_D_Na" variable_2="J_D_Na"/>
<map_variables variable_1="J_D_K" variable_2="J_D_K"/>
<map_variables variable_1="J_D_Cl" variable_2="J_D_Cl"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="Basolateral_Diffusion"/>
<map_variables variable_1="q_Na_c" variable_2="q_Na_c"/>
<map_variables variable_1="q_K_c" variable_2="q_K_c"/>
<map_variables variable_1="q_Cl_c" variable_2="q_Cl_c"/>
<map_variables variable_1="J_D_Bl_Na" variable_2="J_D_Bl_Na"/>
<map_variables variable_1="J_D_Bl_K" variable_2="J_D_Bl_K"/>
<map_variables variable_1="J_D_Bl_Cl" variable_2="J_D_Bl_Cl"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="Paracellular_Pathway"/>
<map_variables variable_1="J_D_P_Na" variable_2="J_D_P_Na"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="Parameters"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="Epithelial_Cell" component_2="Rate_Constants"/>
<map_variables variable_1="q_Na_c" variable_2="q_Na_c"/>
<map_variables variable_1="q_Gl_c" variable_2="q_Gl_c"/>
</connection>
<!-- Mapping Basolateral Membrane-->
<connection>
<map_components component_1="Basolateral_Membrane" component_2="GLUT2"/>
<map_variables variable_1="q_Gl_b" variable_2="q_Gl_b"/>
</connection>
<connection>
<map_components component_1="Basolateral_Membrane" component_2="NaK_Pump"/>
<map_variables variable_1="q_K_b" variable_2="q_K_b"/>
<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_b" variable_2="V_b"/>
</connection>
<connection>
<map_components component_1="Basolateral_Membrane" component_2="NBC"/>
<map_variables variable_1="V_b" variable_2="V_b"/>
<map_variables variable_1="I_NBC_HCO3" variable_2="I_NBC_HCO3"/>
<map_variables variable_1="I_NBC_Na" variable_2="I_NBC_Na"/>
<map_variables variable_1="q_Na_b" variable_2="q_Na_b"/>
</connection>
<connection>
<map_components component_1="Basolateral_Membrane" component_2="IK"/>
<map_variables variable_1="q_K_b" variable_2="q_K_b"/>
<map_variables variable_1="V_b" variable_2="V_b"/>
<map_variables variable_1="I_IK" variable_2="I_IK"/>
</connection>
<connection>
<map_components component_1="Basolateral_Membrane" component_2="Paracellular_Pathway"/>
<map_variables variable_1="V_b" variable_2="V_b"/>
<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"/>
<map_variables variable_1="q_Na_b" variable_2="q_Na_b"/>
<map_variables variable_1="q_K_b" variable_2="q_K_b"/>
<map_variables variable_1="q_Cl_b" variable_2="q_Cl_b"/>
</connection>
<connection>
<map_components component_1="Basolateral_Membrane" component_2="Basolateral_Diffusion"/>
<map_variables variable_1="q_Na_b" variable_2="q_Na_b"/>
<map_variables variable_1="q_K_b" variable_2="q_K_b"/>
<map_variables variable_1="q_Cl_b" variable_2="q_Cl_b"/>
<map_variables variable_1="I_D_Bl_Na" variable_2="I_D_Bl_Na"/>
<map_variables variable_1="I_D_Bl_K" variable_2="I_D_Bl_K"/>
<map_variables variable_1="I_D_Bl_Cl" variable_2="I_D_Bl_Cl"/>
<map_variables variable_1="V_b" variable_2="V_b"/>
</connection>
<connection>
<map_components component_1="Basolateral_Membrane" component_2="Blood"/>
<map_variables variable_1="q_Gl_b" variable_2="q_Gl_b"/>
</connection>
<connection>
<map_components component_1="Basolateral_Membrane" component_2="Parameters"/>
<map_variables variable_1="C" variable_2="C"/>
<map_variables variable_1="t" variable_2="t"/>
</connection>
<!-- Mapping Parameters-->
<connection>
<map_components component_1="Parameters" component_2="SGLT1"/>
<map_variables variable_1="k0_23" variable_2="k0_23"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="NHE3"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="AE1"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="GLUT2"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="NaK_Pump"/>
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="NBC"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="BK"/>
<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_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="NKCC"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="Apical_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_Apical" variable_2="A_Apical"/>
<map_variables variable_1="P_Na" variable_2="P_Na"/>
<map_variables variable_1="P_k" variable_2="P_k"/>
<map_variables variable_1="P_Cl" variable_2="P_Cl"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="Basolateral_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_Basol" variable_2="A_Basol"/>
<map_variables variable_1="P_Bl_Na" variable_2="P_Bl_Na"/>
<map_variables variable_1="P_Bl_k" variable_2="P_Bl_k"/>
<map_variables variable_1="P_Bl_Cl" variable_2="P_Bl_Cl"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="Paracellular_Pathway"/>
<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_Apical" variable_2="A_Apical"/>
<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="Parameters" component_2="Blood"/>
<map_variables variable_1="t" variable_2="t"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="CFTR"/>
<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_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="ENaC"/>
<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_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="IK"/>
<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_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="GLUT2_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Parameters" component_2="Rate_Constants"/>
<map_variables variable_1="k0_12" variable_2="k0_12"/>
<map_variables variable_1="k0_21" variable_2="k0_21"/>
<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"/>
<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="Rate_Constants" component_2="SGLT1"/>
<map_variables variable_1="ks_12" variable_2="ks_12"/>
<map_variables variable_1="k_12" variable_2="k_12"/>
<map_variables variable_1="k_21" variable_2="k_21"/>
<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>
</model>
