Location: BG_TCC @ ef16f59ea55f / tcc_1.cellml
- Author:
- Shelley Fong <sfon036@UoA.auckland.ac.nz>
- Date:
- 2024-12-13 14:16:35+13:00
- Desc:
- Updating params to match LRd cell vol
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/831/rawfile/ef16f59ea55fa2a5373cdc4c4cb63686afd04653/tcc_1.cellml
<?xml version='1.0'?>
<!-- transient outward 1-->
<model name="t_type_ca_1" 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">
<import xlink:href="../Units_BG.cellml">
<units name="mM" units_ref="mM"/>
<units name="fmol" units_ref="fmol"/>
<units name="per_fmol" units_ref="per_fmol"/>
<units name="J_per_mol" units_ref="J_per_mol"/>
<units name="fmol_per_s" units_ref="fmol_per_sec"/>
<units name="C_per_mol" units_ref="C_per_mol"/>
<units name="J_per_C" units_ref="J_per_C"/>
<units name="microm3" units_ref="microm3"/>
<units name="fF" units_ref="fF"/>
<units name="fC" units_ref="fC"/>
<units name="fA" units_ref="fA"/>
</import>
<import xlink:href="../Units_kinetic.cellml">
<units name="microA_per_microF" units_ref="microA_per_microF"/>
<units name="millivolt" units_ref="millivolt"/>
<units name="per_second" units_ref="per_second"/>
<units name="coulomb_per_s" units_ref="coulomb_per_s"/>
</import>
<component name="membrane">
<variable name="time" public_interface="out" units="second"/>
<variable initial_value="1.0e6" name="C_m" public_interface="out" units="fF"/>
<variable initial_value="0.04653" name="qCao" public_interface="out" units="fmol"/>
<variable initial_value="0.000001551" name="qCai" public_interface="out" units="fmol"/>
<variable initial_value="1e4" name="vol_o" public_interface="out" units="microm3"/>
<variable name="q_membrane" public_interface="out" units="fC"/>
<variable name="Vmem" units="volt"/>
<!---9E-11-->
<variable initial_value="2494.34" name="RT" public_interface="out" units="J_per_mol"/>
<variable initial_value="2" name="z" public_interface="out" units="dimensionless"/>
<variable initial_value="96500" name="F" public_interface="out" units="C_per_mol"/>
<variable name="v_tcc_GHK" public_interface="in" units="fmol_per_s"/>
<variable name="v_E" public_interface="in" units="fA"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>qCai</ci>
</apply>
<apply>
<minus/>
<ci>v_tcc_GHK</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>qCao</ci>
</apply>
<ci>v_tcc_GHK</ci>
</apply>
<!--ode(q_membrane, time) = -v_E;-->
<apply>
<eq/>
<ci>q_membrane</ci>
<apply>
<plus/>
<apply>
<minus/>
<cn cellml:units="fC" type="e-notation">10<sep/>4</cn>
</apply>
<apply>
<times/>
<apply>
<floor/>
<apply>
<divide/>
<ci>time</ci>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<cn cellml:units="dimensionless" type="e-notation">2<sep/>3</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Vmem</ci>
<apply>
<divide/>
<ci>q_membrane</ci>
<ci>C_m</ci>
</apply>
</apply>
</math>
</component>
<component name="tcc_channel">
<variable name="t" private_interface="out" public_interface="in" units="second"/>
<variable name="RT" private_interface="out" public_interface="in" units="J_per_mol"/>
<variable name="z" public_interface="in" units="dimensionless"/>
<variable name="F" private_interface="out" public_interface="in" units="C_per_mol"/>
<variable name="C_m" public_interface="in" units="fF"/>
<!-- Parameters-->
<!-- overall ion permeability (GHK) - but you can choose the channel density (embedded in kappa_Ca)-->
<!--var kappa_Ca: fmol_per_s {init: 26249.98853};-->
<variable initial_value="26249.98853e-3" name="kappa_Ca" units="fmol_per_s"/>
<variable initial_value="162.0184821" name="K_Cao" units="per_fmol"/>
<variable initial_value="162.0184821" name="K_Cai" units="per_fmol"/>
<variable name="Cao" units="mM"/>
<variable name="vol_o" public_interface="in" units="microm3"/>
<variable name="q_s21" private_interface="in" units="fmol"/>
<!-- state variables-->
<variable name="qCai" public_interface="in" units="fmol"/>
<variable name="qCao" public_interface="in" units="fmol"/>
<variable name="q_E" public_interface="in" units="fC"/>
<variable name="i_BG_in_tcc_out" units="microA_per_microF"/>
<variable name="v_tcc_GHK" public_interface="out" units="fmol_per_s"/>
<variable name="v_E" public_interface="out" units="fA"/>
<variable name="u_Cai" units="J_per_mol"/>
<variable name="u_Cao" units="J_per_mol"/>
<variable name="u_E" private_interface="out" units="J_per_C"/>
<variable name="U_mem" units="J_per_mol"/>
<variable name="V_BG_mV" units="millivolt"/>
<variable name="Af_Ca" units="J_per_mol"/>
<variable name="Ar_Ca" units="J_per_mol"/>
<variable name="mu_s21" private_interface="in" units="J_per_mol"/>
<!-- Conservation laws-->
<!-- Constitutive relations - Convention: positive flow from inside to outside-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Cao</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>-3</cn>
<ci>qCao</ci>
</apply>
<ci>vol_o</ci>
</apply>
</apply>
<!-- Am 01 the same as U_mem-->
<apply>
<eq/>
<ci>U_mem</ci>
<apply>
<times/>
<ci>F</ci>
<ci>u_E</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_Cai</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Cai</ci>
<ci>qCai</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>u_Cao</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Cao</ci>
<ci>qCao</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_Ca</ci>
<apply>
<plus/>
<ci>mu_s21</ci>
<ci>u_Cai</ci>
<ci>U_mem</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_Ca</ci>
<apply>
<plus/>
<ci>mu_s21</ci>
<ci>u_Cao</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_tcc_GHK</ci>
<piecewise>
<piece>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>kappa_Ca</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_Ca</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_Ca</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>U_mem</ci>
<cn cellml:units="J_per_mol">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>kappa_Ca</ci>
<ci>U_mem</ci>
</apply>
<ci>RT</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>U_mem</ci>
<ci>RT</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_Ca</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_Ca</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>u_E</ci>
<apply>
<divide/>
<ci>q_E</ci>
<ci>C_m</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_E</ci>
<apply>
<times/>
<ci>z</ci>
<ci>F</ci>
<ci>v_tcc_GHK</ci>
</apply>
</apply>
<!-- for comparison, calculate the change in concentration AND current as well-->
<apply>
<eq/>
<ci>i_BG_in_tcc_out</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>v_tcc_GHK</ci>
<ci>F</ci>
</apply>
<ci>C_m</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_BG_mV</ci>
<apply>
<times/>
<ci>u_E</ci>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>3</cn>
</apply>
</apply>
</math>
</component>
<component name="tcc_gate">
<variable name="t" public_interface="in" units="second"/>
<variable name="V" public_interface="in" units="volt"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="RT" public_interface="in" units="J_per_mol"/>
<!-- channel states transitions-->
<variable name="v_s00" units="fmol_per_s"/>
<variable name="v_s10" units="fmol_per_s"/>
<variable name="v_s20" units="fmol_per_s"/>
<variable name="v_s01" units="fmol_per_s"/>
<variable name="v_s11" units="fmol_per_s"/>
<variable name="v_s21" units="fmol_per_s"/>
<!-- gate variable transitions-->
<variable name="v_b00" units="fmol_per_s"/>
<variable name="v_b10" units="fmol_per_s"/>
<variable name="v_b01" units="fmol_per_s"/>
<variable name="v_b11" units="fmol_per_s"/>
<variable name="v_g00" units="fmol_per_s"/>
<variable name="v_g10" units="fmol_per_s"/>
<variable name="v_g20" units="fmol_per_s"/>
<!-- potential of states-->
<variable name="mu_s00" units="J_per_mol"/>
<variable name="mu_s10" units="J_per_mol"/>
<variable name="mu_s20" units="J_per_mol"/>
<variable name="mu_s01" units="J_per_mol"/>
<variable name="mu_s11" units="J_per_mol"/>
<variable name="mu_s21" public_interface="out" units="J_per_mol"/>
<!-- forward/reverse affinities at each reaction-->
<variable name="Af_b00" units="J_per_mol"/>
<variable name="Ar_b00" units="J_per_mol"/>
<variable name="Af_b10" units="J_per_mol"/>
<variable name="Ar_b10" units="J_per_mol"/>
<variable name="Af_b01" units="J_per_mol"/>
<variable name="Ar_b01" units="J_per_mol"/>
<variable name="Af_b11" units="J_per_mol"/>
<variable name="Ar_b11" units="J_per_mol"/>
<variable name="Af_g00" units="J_per_mol"/>
<variable name="Ar_g00" units="J_per_mol"/>
<variable name="Af_g10" units="J_per_mol"/>
<variable name="Ar_g10" units="J_per_mol"/>
<variable name="Af_g20" units="J_per_mol"/>
<variable name="Ar_g20" units="J_per_mol"/>
<!-- gating variables-->
<variable initial_value="59.49036954" name="kappa_b00" units="fmol_per_s"/>
<variable initial_value="148.0184892" name="kappa_b10" units="fmol_per_s"/>
<variable initial_value="2.147331948" name="kappa_b01" units="fmol_per_s"/>
<variable initial_value="5.342794695" name="kappa_b11" units="fmol_per_s"/>
<variable initial_value="0.387904368" name="kappa_g00" units="fmol_per_s"/>
<variable initial_value="1.930296249" name="kappa_g10" units="fmol_per_s"/>
<variable initial_value="2.401393173" name="kappa_g20" units="fmol_per_s"/>
<!-- -->
<variable initial_value="4.48207646" name="K_00" units="per_fmol"/>
<variable initial_value="0.900699589" name="K_10" units="per_fmol"/>
<variable initial_value="0.72400349" name="K_20" units="per_fmol"/>
<variable initial_value="124.1728766" name="K_01" units="per_fmol"/>
<variable initial_value="24.95326885" name="K_11" units="per_fmol"/>
<variable initial_value="20.05802374" name="K_21" units="per_fmol"/>
<!-- -->
<variable initial_value="0.455486889" name="z_bf" units="dimensionless"/>
<variable initial_value="-0.203216832" name="z_br" units="dimensionless"/>
<variable initial_value="-0.461802026" name="z_gf" units="dimensionless"/>
<variable initial_value="0.381995996" name="z_gr" units="dimensionless"/>
<!-- BG constants : -->
<variable initial_value="0.999995" name="q_s00" units="fmol"/>
<variable initial_value="1e-6" name="q_s10" units="fmol"/>
<variable initial_value="1e-6" name="q_s20" units="fmol"/>
<variable initial_value="1e-6" name="q_s01" units="fmol"/>
<variable initial_value="1e-6" name="q_s11" units="fmol"/>
<variable initial_value="1e-6" name="q_s21" public_interface="out" units="fmol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Af_b00</ci>
<apply>
<plus/>
<ci>mu_s00</ci>
<apply>
<times/>
<ci>z_bf</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_b00</ci>
<apply>
<plus/>
<ci>mu_s10</ci>
<apply>
<times/>
<ci>z_br</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_b10</ci>
<apply>
<plus/>
<ci>mu_s10</ci>
<apply>
<times/>
<ci>z_bf</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_b10</ci>
<apply>
<plus/>
<ci>mu_s20</ci>
<apply>
<times/>
<ci>z_br</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_b01</ci>
<apply>
<plus/>
<ci>mu_s01</ci>
<apply>
<times/>
<ci>z_bf</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_b01</ci>
<apply>
<plus/>
<ci>mu_s11</ci>
<apply>
<times/>
<ci>z_br</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_b11</ci>
<apply>
<plus/>
<ci>mu_s11</ci>
<apply>
<times/>
<ci>z_bf</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_b11</ci>
<apply>
<plus/>
<ci>mu_s21</ci>
<apply>
<times/>
<ci>z_br</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_g00</ci>
<apply>
<plus/>
<ci>mu_s00</ci>
<apply>
<times/>
<ci>z_gf</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_g00</ci>
<apply>
<plus/>
<ci>mu_s01</ci>
<apply>
<times/>
<ci>z_gr</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_g10</ci>
<apply>
<plus/>
<ci>mu_s10</ci>
<apply>
<times/>
<ci>z_gf</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_g10</ci>
<apply>
<plus/>
<ci>mu_s11</ci>
<apply>
<times/>
<ci>z_gr</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_g20</ci>
<apply>
<plus/>
<ci>mu_s20</ci>
<apply>
<times/>
<ci>z_gf</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_g20</ci>
<apply>
<plus/>
<ci>mu_s21</ci>
<apply>
<times/>
<ci>z_gr</ci>
<ci>F</ci>
<ci>V</ci>
</apply>
</apply>
</apply>
<!-- -->
<apply>
<eq/>
<ci>v_b00</ci>
<apply>
<times/>
<ci>kappa_b00</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_b00</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_b00</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_b10</ci>
<apply>
<times/>
<ci>kappa_b10</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_b10</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_b10</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_b01</ci>
<apply>
<times/>
<ci>kappa_b01</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_b01</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_b01</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_b11</ci>
<apply>
<times/>
<ci>kappa_b11</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_b11</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_b11</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_g00</ci>
<apply>
<times/>
<ci>kappa_g00</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_g00</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_g00</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_g10</ci>
<apply>
<times/>
<ci>kappa_g10</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_g10</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_g10</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_g20</ci>
<apply>
<times/>
<ci>kappa_g20</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_g20</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_g20</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- gate variable velocities-->
<apply>
<eq/>
<ci>v_s00</ci>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_b00</ci>
</apply>
<ci>v_g00</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_s10</ci>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_b00</ci>
<ci>v_b10</ci>
</apply>
<ci>v_g10</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_s20</ci>
<apply>
<minus/>
<ci>v_b10</ci>
<ci>v_g20</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_s01</ci>
<apply>
<minus/>
<ci>v_g00</ci>
<ci>v_b01</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_s11</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_b01</ci>
<ci>v_b11</ci>
</apply>
<ci>v_g10</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_s21</ci>
<apply>
<plus/>
<ci>v_b11</ci>
<ci>v_g20</ci>
</apply>
</apply>
<!-- gate variable potentials-->
<apply>
<eq/>
<ci>mu_s00</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_00</ci>
<ci>q_s00</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_s10</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_10</ci>
<ci>q_s10</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_s20</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_20</ci>
<ci>q_s20</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_s01</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_01</ci>
<ci>q_s01</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_s11</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_11</ci>
<ci>q_s11</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_s21</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_21</ci>
<ci>q_s21</ci>
</apply>
</apply>
</apply>
</apply>
<!-- rate of change of state variables-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_s00</ci>
</apply>
<ci>v_s00</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_s10</ci>
</apply>
<ci>v_s10</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_s20</ci>
</apply>
<ci>v_s20</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_s01</ci>
</apply>
<ci>v_s01</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_s11</ci>
</apply>
<ci>v_s11</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_s21</ci>
</apply>
<ci>v_s21</ci>
</apply>
</math>
</component>
<connection>
<map_components component_1="membrane" component_2="tcc_channel"/>
<map_variables variable_1="time" variable_2="t"/>
<map_variables variable_1="RT" variable_2="RT"/>
<map_variables variable_1="z" variable_2="z"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="qCai" variable_2="qCai"/>
<map_variables variable_1="qCao" variable_2="qCao"/>
<map_variables variable_1="q_membrane" variable_2="q_E"/>
<map_variables variable_1="C_m" variable_2="C_m"/>
<map_variables variable_1="vol_o" variable_2="vol_o"/>
<map_variables variable_1="v_tcc_GHK" variable_2="v_tcc_GHK"/>
<map_variables variable_1="v_E" variable_2="v_E"/>
</connection>
<connection>
<map_components component_1="tcc_channel" component_2="tcc_gate"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="RT" variable_2="RT"/>
<map_variables variable_1="u_E" variable_2="V"/>
<map_variables variable_1="mu_s21" variable_2="mu_s21"/>
<map_variables variable_1="q_s21" variable_2="q_s21"/>
</connection>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="tcc_channel">
<component_ref component="tcc_gate"/>
</component_ref>
</group>
</model>