- Author:
- Shelley Fong <s.fong@auckland.ac.nz>
- Date:
- 2021-11-23 15:24:25+13:00
- Desc:
- Fitting exponential curve to k_off and adding it as product of reaction
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/7a9/rawfile/de8455407814dd6856b754511afd2eb172ef2001/BG_TRPN_tension_koff.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="individual_TRPN" 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_and_constants/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_sec" 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"/>
<units name="per_second" units_ref="per_second"/>
<units name="millivolt" units_ref="millivolt"/>
<units name="per_sec" units_ref="per_sec"/>
<units name="J_per_K_per_mol" units_ref="J_per_K_per_mol"/>
<units name="fmol_per_L" units_ref="fmol_per_L"/>
<units name="fmol_per_L_per_sec" units_ref="fmol_per_L_per_sec"/>
<units name="per_sec_per_fmol_per_L" units_ref="per_sec_per_fmol_per_L"/>
<units name="uM" units_ref="uM"/>
<units name="mM_per_sec" units_ref="mM_per_sec"/>
<units name="uM_per_sec" units_ref="uM_per_sec"/>
<units name="pL" units_ref="pL"/>
<units name="m_to_u" units_ref="m_to_u"/>
</import>
<units name="mm">
<unit prefix="milli" units="metre"/>
</units>
<units name="mm2">
<unit exponent="2" units="mm"/>
</units>
<units name="N_per_mm2">
<unit units="newton"/>
<unit exponent="-1" units="mm2"/>
</units>
<units name="mm3_per_mol">
<unit exponent="3" units="mm"/>
<unit exponent="-1" units="mole"/>
</units>
<units name="J_per_mm3">
<unit units="joule"/>
<unit exponent="-3" units="mm"/>
</units>
<units name="mil_to_1">
<unit units="second"/>
<unit exponent="-1" prefix="milli" units="second"/>
</units>
<import xlink:href="units_and_constants/constants_BG.cellml">
<component component_ref="constants" name="constants"/>
</import>
<component name="environment">
<variable name="time" public_interface="out" units="second"/>
<!-- Tension-->
<variable name="tension" public_interface="out" units="N_per_mm2"/>
<variable initial_value="7.85e-5" name="Area" public_interface="out" units="mm2"/>
<variable initial_value="1" name="strainlength" public_interface="out" units="mm"/>
<variable initial_value="1e-4" name="q_XB" public_interface="out" units="mole"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>tension</ci>
<piecewise>
<piece>
<cn cellml:units="N_per_mm2">1</cn>
<apply>
<gt/>
<ci>time</ci>
<cn cellml:units="second">1</cn>
</apply>
</piece>
<otherwise>
<cn cellml:units="N_per_mm2">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
<!-- initial values-->
<variable initial_value="6.82e-1" name="q_Cai_init" units="fmol"/>
<variable initial_value="2.57" name="q_TRPN_init" units="fmol"/>
<variable initial_value="1e-16" name="q_Ca_TRPN_init" units="fmol"/>
<!-- Global value-->
<variable name="q_Cai" public_interface="out" units="fmol"/>
<variable name="q_TRPN" public_interface="out" units="fmol"/>
<variable name="q_Ca_TRPN" public_interface="out" units="fmol"/>
<!-- From submodule-->
<variable name="q_Cai_mTRPN" public_interface="in" units="fmol"/>
<variable name="q_TRPN_mTRPN" public_interface="in" units="fmol"/>
<variable name="q_Ca_TRPN_mTRPN" public_interface="in" units="fmol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_Cai</ci>
<apply>
<plus/>
<ci>q_Cai_mTRPN</ci>
<ci>q_Cai_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_TRPN</ci>
<apply>
<plus/>
<ci>q_TRPN_mTRPN</ci>
<ci>q_TRPN_init</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Ca_TRPN</ci>
<apply>
<plus/>
<ci>q_Ca_TRPN_mTRPN</ci>
<ci>q_Ca_TRPN_init</ci>
</apply>
</apply>
</math>
</component>
<component name="TRPN_parameters">
<variable initial_value="83.2553" name="kappa_R_TRPNCa" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="1.00748" name="K_Cai" public_interface="out" units="per_fmol"/>
<variable initial_value="1.00748" name="K_TRPN" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0698328" name="K_Ca_TRPN" public_interface="out" units="per_fmol"/>
</component>
<component name="TRPN">
<variable name="time" public_interface="in" units="second"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<!-- parameters-->
<variable name="kappa_R_TRPNCa" public_interface="in" units="fmol_per_sec"/>
<variable name="K_Cai" public_interface="in" units="per_fmol"/>
<variable name="K_TRPN" public_interface="in" units="per_fmol"/>
<variable name="K_Ca_TRPN" public_interface="in" units="per_fmol"/>
<!-- Input from global environment-->
<variable name="q_Cai_global" public_interface="in" units="fmol"/>
<variable name="q_TRPN_global" public_interface="in" units="fmol"/>
<variable name="q_Ca_TRPN_global" public_interface="in" units="fmol"/>
<!-- Output to global environment-->
<variable initial_value="1e-16" name="q_Cai" public_interface="out" units="fmol"/>
<variable initial_value="1e-16" name="q_TRPN" public_interface="out" units="fmol"/>
<variable initial_value="1e-16" name="q_Ca_TRPN" public_interface="out" units="fmol"/>
<!-- Tension and transformer-->
<variable name="tension" public_interface="in" units="N_per_mm2"/>
<variable name="Area" public_interface="in" units="mm2"/>
<variable name="strainlength" public_interface="in" units="mm"/>
<variable name="q_XB" public_interface="in" units="mole"/>
<variable initial_value="1e3" name="m" units="mil_to_1"/>
<variable name="conv" units="mm3_per_mol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>conv</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>Area</ci>
<ci>strainlength</ci>
</apply>
<ci>q_XB</ci>
</apply>
</apply>
</math>
<!-- exponential curve fit to k_off-->
<variable name="mu_TF" units="J_per_mol"/>
<variable initial_value="1.0247382029957026" name="A" units="J_per_mm3"/>
<variable initial_value="0.011770022274421879" name="k" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu_TF</ci>
<apply>
<times/>
<ci>conv</ci>
<ci>A</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>k</ci>
</apply>
<ci>tension</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Constitutive parameters-->
<variable name="mu_Cai" units="J_per_mol"/>
<variable name="mu_TRPN" units="J_per_mol"/>
<variable name="mu_Ca_TRPN" units="J_per_mol"/>
<variable name="v_R_TRPNCa" units="fmol_per_sec"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu_Cai</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Cai</ci>
<ci>q_Cai_global</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_TRPN</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_TRPN</ci>
<ci>q_TRPN_global</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Ca_TRPN</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Ca_TRPN</ci>
<ci>q_Ca_TRPN_global</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_R_TRPNCa</ci>
<apply>
<times/>
<ci>kappa_R_TRPNCa</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_Cai</ci>
<ci>mu_TRPN</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_Ca_TRPN</ci>
<ci>mu_TF</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_Cai</ci>
</apply>
<apply>
<minus/>
<ci>v_R_TRPNCa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_TRPN</ci>
</apply>
<apply>
<minus/>
<ci>v_R_TRPNCa</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_Ca_TRPN</ci>
</apply>
<ci>v_R_TRPNCa</ci>
</apply>
</math>
</component>
<connection>
<map_components component_1="environment" component_2="TRPN"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="q_Cai_mTRPN" variable_2="q_Cai"/>
<map_variables variable_1="q_Cai" variable_2="q_Cai_global"/>
<map_variables variable_1="q_TRPN_mTRPN" variable_2="q_TRPN"/>
<map_variables variable_1="q_TRPN" variable_2="q_TRPN_global"/>
<map_variables variable_1="q_Ca_TRPN_mTRPN" variable_2="q_Ca_TRPN"/>
<map_variables variable_1="q_Ca_TRPN" variable_2="q_Ca_TRPN_global"/>
<map_variables variable_1="tension" variable_2="tension"/>
<map_variables variable_1="Area" variable_2="Area"/>
<map_variables variable_1="strainlength" variable_2="strainlength"/>
<map_variables variable_1="q_XB" variable_2="q_XB"/>
</connection>
<connection>
<map_components component_1="TRPN" component_2="TRPN_parameters"/>
<map_variables variable_1="kappa_R_TRPNCa" variable_2="kappa_R_TRPNCa"/>
<map_variables variable_1="K_Cai" variable_2="K_Cai"/>
<map_variables variable_1="K_TRPN" variable_2="K_TRPN"/>
<map_variables variable_1="K_Ca_TRPN" variable_2="K_Ca_TRPN"/>
</connection>
<connection>
<map_components component_1="constants" component_2="TRPN"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
</connection>
</model>