- Author:
- Shelley Fong <sfon036@UoA.auckland.ac.nz>
- Date:
- 2024-11-20 11:06:49+13:00
- Desc:
- Removing inactive R as it is all merged into active R species. fast- and slow- kineticconstant are altered to try make system react to L but to no avail.
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/7dd/rawfile/7d0610d3812fadebe4ee5ed14f15212da960620d/BG_GPCR_M2_reduced.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="individual_GPCR_M2_reduced" 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>
<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"/>
<!-- initial values-->
<!-- var q_RM2_inactive: fmol {init: 0.00072162, pub: out};
var q_L_RM2_inactive: fmol {init: 1e-18, pub: out};-->
<variable name="q_LM2" public_interface="out" units="fmol"/>
<variable name="q_LM2_stim" units="fmol"/>
<variable initial_value="0" name="q_LM2_ode" units="fmol"/>
<variable initial_value="1e-18" name="q_RM2" public_interface="out" units="fmol"/>
<variable initial_value="0.00836" name="q_Gi" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_RM2_Gi" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_L_RM2" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_L_RM2_Gi" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_Gia_GTP" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_Gibetagamma" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_Gia_GDP" public_interface="out" units="fmol"/>
<variable initial_value="2.2" name="q_GTP" public_interface="out" units="fmol"/>
<variable initial_value="1.1" name="q_GDP" public_interface="out" units="fmol"/>
<variable initial_value="570" name="q_Pi" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_RM2_tag" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_L_RM2_tag" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_RM2_GRKArr" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_L_RM2_GRKArr" public_interface="out" units="fmol"/>
<variable initial_value="1e-18" name="q_GRKArr" public_interface="out" units="fmol"/>
<!-- stimulus-->
<variable initial_value="2e1" name="stimSt" units="second"/>
<variable initial_value="2e1" name="stimDur" units="second"/>
<variable initial_value="1e1" name="tR" units="second"/>
<variable initial_value="1e5" name="stimMag" units="fmol"/>
<variable initial_value="1e-8" name="stimHolding" units="fmol"/>
<variable name="m" units="fmol_per_sec"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>m</ci>
<apply>
<divide/>
<ci>stimMag</ci>
<ci>tR</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_LM2_stim</ci>
<piecewise>
<piece>
<apply>
<plus/>
<ci>stimHolding</ci>
<apply>
<times/>
<ci>m</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>time</ci>
<ci>stimSt</ci>
</apply>
<ci>tR</ci>
</apply>
</apply>
</apply>
<apply>
<and/>
<apply>
<lt/>
<ci>time</ci>
<ci>stimSt</ci>
</apply>
<apply>
<gt/>
<ci>time</ci>
<apply>
<minus/>
<ci>stimSt</ci>
<ci>tR</ci>
</apply>
</apply>
</apply>
</piece>
<piece>
<apply>
<plus/>
<ci>stimMag</ci>
<ci>stimHolding</ci>
</apply>
<apply>
<and/>
<apply>
<geq/>
<ci>time</ci>
<ci>stimSt</ci>
</apply>
<apply>
<lt/>
<ci>time</ci>
<apply>
<plus/>
<ci>stimSt</ci>
<ci>stimDur</ci>
</apply>
</apply>
</apply>
</piece>
<piece>
<apply>
<plus/>
<ci>stimHolding</ci>
<apply>
<times/>
<apply>
<minus/>
<ci>m</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>time</ci>
<ci>stimSt</ci>
</apply>
<ci>tR</ci>
</apply>
<ci>stimDur</ci>
</apply>
</apply>
</apply>
<apply>
<and/>
<apply>
<lt/>
<ci>time</ci>
<apply>
<plus/>
<ci>stimSt</ci>
<ci>tR</ci>
<ci>stimDur</ci>
</apply>
</apply>
<apply>
<geq/>
<ci>time</ci>
<apply>
<plus/>
<ci>stimSt</ci>
<ci>stimDur</ci>
</apply>
</apply>
</apply>
</piece>
<otherwise>
<ci>stimHolding</ci>
</otherwise>
</piecewise>
</apply>
</math>
<!-- mass conservation checks-->
<variable name="conserv_R_T" units="fmol"/>
<variable name="conserv_L_T" units="fmol"/>
<variable name="conserv_G_T" units="fmol"/>
<variable name="total_R" units="fmol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>conserv_R_T</ci>
<apply>
<plus/>
<ci>q_RM2</ci>
<ci>q_RM2_Gi</ci>
<ci>q_L_RM2</ci>
<ci>q_L_RM2_Gi</ci>
<ci>q_RM2_tag</ci>
<ci>q_L_RM2_tag</ci>
<ci>q_RM2_GRKArr</ci>
<ci>q_L_RM2_GRKArr</ci>
</apply>
</apply>
<!--+q_RM2_inactive+q_L_RM2_inactive-->
<apply>
<eq/>
<ci>total_R</ci>
<apply>
<plus/>
<ci>q_RM2</ci>
<ci>q_RM2_Gi</ci>
<ci>q_L_RM2</ci>
<ci>q_L_RM2_Gi</ci>
</apply>
</apply>
<!--+q_RM2_inactive+q_L_RM2_inactive-->
<apply>
<eq/>
<ci>conserv_L_T</ci>
<apply>
<plus/>
<ci>q_LM2</ci>
<ci>q_L_RM2</ci>
<ci>q_L_RM2_Gi</ci>
<ci>q_L_RM2_tag</ci>
<ci>q_L_RM2_GRKArr</ci>
</apply>
</apply>
<!--+q_L_RM2_inactive-->
<apply>
<eq/>
<ci>conserv_G_T</ci>
<apply>
<plus/>
<ci>q_Gi</ci>
<ci>q_RM2_Gi</ci>
<ci>q_L_RM2_Gi</ci>
<ci>q_Gia_GTP</ci>
<ci>q_Gia_GDP</ci>
</apply>
</apply>
</math>
<!-- From submodule-->
<!--var v_Rswitch_M2: fmol_per_sec {pub: in};
var v_LRswitch_M2: fmol_per_sec {pub: in};-->
<variable name="v_C_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="v_R_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="v_L_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="v_Act1_Gi" public_interface="in" units="fmol_per_sec"/>
<variable name="v_Act2_Gi" public_interface="in" units="fmol_per_sec"/>
<variable name="v_Hyd_Gi" public_interface="in" units="fmol_per_sec"/>
<variable name="v_Reassoc_Gi" public_interface="in" units="fmol_per_sec"/>
<variable name="v_InternR_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="v_InternLR_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="v_Rsynthesis" units="fmol_per_sec"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_Rsynthesis</ci>
<apply>
<plus/>
<apply>
<plus/>
<ci>v_C_M2</ci>
</apply>
<ci>v_L_M2</ci>
</apply>
</apply>
<!--ode(q_RM2_inactive, time) = -v_Rswitch_M2-v_L_M2;
ode(q_L_RM2_inactive, time) = -v_LRswitch_M2;-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_LM2_ode</ci>
</apply>
<apply>
<minus/>
<ci>v_L_M2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_RM2</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_C_M2</ci>
</apply>
<ci>v_L_M2</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1.0</cn>
<ci>v_Rsynthesis</ci>
</apply>
</apply>
</apply>
<!--+v_Rswitch_M2-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_Gi</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_C_M2</ci>
</apply>
<ci>v_R_M2</ci>
</apply>
<ci>v_Reassoc_Gi</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_RM2_Gi</ci>
</apply>
<apply>
<minus/>
<ci>v_C_M2</ci>
<ci>v_Act1_Gi</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_L_RM2</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_R_M2</ci>
</apply>
<ci>v_L_M2</ci>
</apply>
</apply>
<!--+v_LRswitch_M2-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_L_RM2_Gi</ci>
</apply>
<apply>
<minus/>
<ci>v_R_M2</ci>
<ci>v_Act2_Gi</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_Gia_GTP</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_Act1_Gi</ci>
<ci>v_Act2_Gi</ci>
</apply>
<ci>v_Hyd_Gi</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_Gibetagamma</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_Act1_Gi</ci>
<ci>v_Act2_Gi</ci>
</apply>
<ci>v_Reassoc_Gi</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_Gia_GDP</ci>
</apply>
<apply>
<minus/>
<ci>v_Hyd_Gi</ci>
<ci>v_Reassoc_Gi</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_GTP</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_Act1_Gi</ci>
</apply>
<ci>v_Act2_Gi</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_GDP</ci>
</apply>
<apply>
<plus/>
<ci>v_Act1_Gi</ci>
<ci>v_Act2_Gi</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_Pi</ci>
</apply>
<ci>v_Hyd_Gi</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_RM2_tag</ci>
</apply>
<apply>
<minus/>
<ci>v_Act1_Gi</ci>
<ci>v_InternR_M2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_L_RM2_tag</ci>
</apply>
<apply>
<minus/>
<ci>v_Act2_Gi</ci>
<ci>v_InternLR_M2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_RM2_GRKArr</ci>
</apply>
<ci>v_InternR_M2</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_L_RM2_GRKArr</ci>
</apply>
<ci>v_InternLR_M2</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>q_GRKArr</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_InternR_M2</ci>
</apply>
<ci>v_InternLR_M2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_LM2</ci>
<apply>
<plus/>
<ci>q_LM2_stim</ci>
<ci>q_LM2_ode</ci>
</apply>
</apply>
</math>
</component>
<component name="GPCR_M2_reduced_parameters">
<variable initial_value="61.47" name="kappa_C_M2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="20216.3" name="kappa_R_M2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="34.7851" name="kappa_L_M2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.128845" name="kappa_Act1_Gi" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.0277361" name="kappa_Act2_Gi" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.603302" name="kappa_Hyd_Gi" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.0837786" name="kappa_Reassoc_Gi" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.783492" name="kappa_InternR_M2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="0.16866" name="kappa_InternLR_M2" public_interface="out" units="fmol_per_sec"/>
<variable initial_value="2.64271" name="K_LM2" public_interface="out" units="per_fmol"/>
<variable initial_value="12.5289" name="K_RM2" public_interface="out" units="per_fmol"/>
<variable initial_value="0.00346983" name="K_Gi" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0448642" name="K_RM2_Gi" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0380954" name="K_L_RM2" public_interface="out" units="per_fmol"/>
<variable initial_value="10.4206" name="K_L_RM2_Gi" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0385475" name="K_Gia_GTP" public_interface="out" units="per_fmol"/>
<variable initial_value="0.163786" name="K_Gibetagamma" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0745174" name="K_Gia_GDP" public_interface="out" units="per_fmol"/>
<variable initial_value="0.00730947" name="K_GTP" public_interface="out" units="per_fmol"/>
<variable initial_value="0.000413153" name="K_GDP" public_interface="out" units="per_fmol"/>
<variable initial_value="0.000187971" name="K_Pi" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0212477" name="K_RM2_tag" public_interface="out" units="per_fmol"/>
<variable initial_value="0.0987039" name="K_L_RM2_tag" public_interface="out" units="per_fmol"/>
<variable initial_value="0.000371028" name="K_RM2_GRKArr" public_interface="out" units="per_fmol"/>
<variable initial_value="0.00172357" name="K_L_RM2_GRKArr" public_interface="out" units="per_fmol"/>
<variable initial_value="5.07617" name="K_GRKArr" public_interface="out" units="per_fmol"/>
</component>
<component name="GPCR_M2_reduced">
<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-->
<!--var kappa_Rswitch_M2: fmol_per_sec {pub: in};
var kappa_LRswitch_M2: fmol_per_sec {pub: in};-->
<variable name="kappa_C_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="kappa_R_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="kappa_L_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="kappa_Act1_Gi" public_interface="in" units="fmol_per_sec"/>
<variable name="kappa_Act2_Gi" public_interface="in" units="fmol_per_sec"/>
<variable name="kappa_Hyd_Gi" public_interface="in" units="fmol_per_sec"/>
<variable name="kappa_Reassoc_Gi" public_interface="in" units="fmol_per_sec"/>
<variable name="kappa_InternR_M2" public_interface="in" units="fmol_per_sec"/>
<variable name="kappa_InternLR_M2" public_interface="in" units="fmol_per_sec"/>
<!--var K_RM2_inactive: per_fmol {pub: in};
var K_L_RM2_inactive: per_fmol {pub: in};-->
<variable name="K_LM2" public_interface="in" units="per_fmol"/>
<variable name="K_RM2" public_interface="in" units="per_fmol"/>
<variable name="K_Gi" public_interface="in" units="per_fmol"/>
<variable name="K_RM2_Gi" public_interface="in" units="per_fmol"/>
<variable name="K_L_RM2" public_interface="in" units="per_fmol"/>
<variable name="K_L_RM2_Gi" public_interface="in" units="per_fmol"/>
<variable name="K_Gia_GTP" public_interface="in" units="per_fmol"/>
<variable name="K_Gibetagamma" public_interface="in" units="per_fmol"/>
<variable name="K_Gia_GDP" public_interface="in" units="per_fmol"/>
<variable name="K_GTP" public_interface="in" units="per_fmol"/>
<variable name="K_GDP" public_interface="in" units="per_fmol"/>
<variable name="K_Pi" public_interface="in" units="per_fmol"/>
<variable name="K_RM2_tag" public_interface="in" units="per_fmol"/>
<variable name="K_L_RM2_tag" public_interface="in" units="per_fmol"/>
<variable name="K_RM2_GRKArr" public_interface="in" units="per_fmol"/>
<variable name="K_L_RM2_GRKArr" public_interface="in" units="per_fmol"/>
<variable name="K_GRKArr" public_interface="in" units="per_fmol"/>
<!-- Input from global environment-->
<!--var q_RM2_inactive: fmol {pub: in};
var q_L_RM2_inactive: fmol {pub: in};-->
<variable name="q_LM2" public_interface="in" units="fmol"/>
<variable name="q_RM2" public_interface="in" units="fmol"/>
<variable name="q_Gi" public_interface="in" units="fmol"/>
<variable name="q_RM2_Gi" public_interface="in" units="fmol"/>
<variable name="q_L_RM2" public_interface="in" units="fmol"/>
<variable name="q_L_RM2_Gi" public_interface="in" units="fmol"/>
<variable name="q_Gia_GTP" public_interface="in" units="fmol"/>
<variable name="q_Gibetagamma" public_interface="in" units="fmol"/>
<variable name="q_Gia_GDP" public_interface="in" units="fmol"/>
<variable name="q_GTP" public_interface="in" units="fmol"/>
<variable name="q_GDP" public_interface="in" units="fmol"/>
<variable name="q_Pi" public_interface="in" units="fmol"/>
<variable name="q_RM2_tag" public_interface="in" units="fmol"/>
<variable name="q_L_RM2_tag" public_interface="in" units="fmol"/>
<variable name="q_RM2_GRKArr" public_interface="in" units="fmol"/>
<variable name="q_L_RM2_GRKArr" public_interface="in" units="fmol"/>
<variable name="q_GRKArr" public_interface="in" units="fmol"/>
<!-- Constitutive parameters-->
<!--var mu_RM2_inactive: J_per_mol;
var mu_L_RM2_inactive: J_per_mol;-->
<variable name="mu_LM2" units="J_per_mol"/>
<variable name="mu_RM2" units="J_per_mol"/>
<variable name="mu_Gi" units="J_per_mol"/>
<variable name="mu_RM2_Gi" units="J_per_mol"/>
<variable name="mu_L_RM2" units="J_per_mol"/>
<variable name="mu_L_RM2_Gi" units="J_per_mol"/>
<variable name="mu_Gia_GTP" units="J_per_mol"/>
<variable name="mu_Gibetagamma" units="J_per_mol"/>
<variable name="mu_Gia_GDP" units="J_per_mol"/>
<variable name="mu_GTP" units="J_per_mol"/>
<variable name="mu_GDP" units="J_per_mol"/>
<variable name="mu_Pi" units="J_per_mol"/>
<variable name="mu_RM2_tag" units="J_per_mol"/>
<variable name="mu_L_RM2_tag" units="J_per_mol"/>
<variable name="mu_RM2_GRKArr" units="J_per_mol"/>
<variable name="mu_L_RM2_GRKArr" units="J_per_mol"/>
<variable name="mu_GRKArr" units="J_per_mol"/>
<!--var v_Rswitch_M2: fmol_per_sec {pub: out};
var v_LRswitch_M2: fmol_per_sec {pub: out};-->
<variable name="v_C_M2" public_interface="out" units="fmol_per_sec"/>
<variable name="v_R_M2" public_interface="out" units="fmol_per_sec"/>
<variable name="v_L_M2" public_interface="out" units="fmol_per_sec"/>
<variable name="v_Act1_Gi" public_interface="out" units="fmol_per_sec"/>
<variable name="v_Act2_Gi" public_interface="out" units="fmol_per_sec"/>
<variable name="v_Hyd_Gi" public_interface="out" units="fmol_per_sec"/>
<variable name="v_Reassoc_Gi" public_interface="out" units="fmol_per_sec"/>
<variable name="v_InternR_M2" public_interface="out" units="fmol_per_sec"/>
<variable name="v_InternLR_M2" public_interface="out" units="fmol_per_sec"/>
<!--mu_RM2_inactive = R*T*ln(K_RM2_inactive*q_RM2_inactive);
mu_L_RM2_inactive = R*T*ln(K_L_RM2_inactive*q_L_RM2_inactive);-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu_LM2</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_LM2</ci>
<ci>q_LM2</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_RM2</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_RM2</ci>
<ci>q_RM2</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Gi</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Gi</ci>
<ci>q_Gi</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_RM2_Gi</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_RM2_Gi</ci>
<ci>q_RM2_Gi</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_L_RM2</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_L_RM2</ci>
<ci>q_L_RM2</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_L_RM2_Gi</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_L_RM2_Gi</ci>
<ci>q_L_RM2_Gi</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Gia_GTP</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Gia_GTP</ci>
<ci>q_Gia_GTP</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Gibetagamma</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Gibetagamma</ci>
<ci>q_Gibetagamma</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Gia_GDP</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Gia_GDP</ci>
<ci>q_Gia_GDP</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_GTP</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_GTP</ci>
<ci>q_GTP</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_GDP</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_GDP</ci>
<ci>q_GDP</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Pi</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Pi</ci>
<ci>q_Pi</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_RM2_tag</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_RM2_tag</ci>
<ci>q_RM2_tag</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_L_RM2_tag</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_L_RM2_tag</ci>
<ci>q_L_RM2_tag</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_RM2_GRKArr</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_RM2_GRKArr</ci>
<ci>q_RM2_GRKArr</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_L_RM2_GRKArr</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_L_RM2_GRKArr</ci>
<ci>q_L_RM2_GRKArr</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_GRKArr</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_GRKArr</ci>
<ci>q_GRKArr</ci>
</apply>
</apply>
</apply>
</apply>
<!--v_Rswitch_M2 = kappa_Rswitch_M2*(exp(mu_RM2_inactive/(R*T))-exp(mu_RM2/(R*T)));
v_LRswitch_M2 = kappa_LRswitch_M2*(exp(mu_L_RM2_inactive/(R*T))-exp(mu_L_RM2/(R*T)));-->
<apply>
<eq/>
<ci>v_C_M2</ci>
<apply>
<times/>
<ci>kappa_C_M2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_RM2</ci>
<ci>mu_Gi</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_RM2_Gi</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_R_M2</ci>
<apply>
<times/>
<ci>kappa_R_M2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_L_RM2</ci>
<ci>mu_Gi</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_L_RM2_Gi</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<!--v_L_M2 = kappa_L_M2*(exp((mu_RM2_inactive+mu_LM2)/(R*T))-exp(mu_L_RM2_inactive/(R*T)));-->
<apply>
<eq/>
<ci>v_L_M2</ci>
<apply>
<times/>
<ci>kappa_L_M2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_RM2</ci>
<ci>mu_LM2</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_L_RM2</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Act1_Gi</ci>
<apply>
<times/>
<ci>kappa_Act1_Gi</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_RM2_Gi</ci>
<ci>mu_GTP</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_Gia_GTP</ci>
<ci>mu_Gibetagamma</ci>
<ci>mu_RM2_tag</ci>
<ci>mu_GDP</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Act2_Gi</ci>
<apply>
<times/>
<ci>kappa_Act2_Gi</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_L_RM2_Gi</ci>
<ci>mu_GTP</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_Gia_GTP</ci>
<ci>mu_Gibetagamma</ci>
<ci>mu_L_RM2_tag</ci>
<ci>mu_GDP</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Hyd_Gi</ci>
<apply>
<times/>
<ci>kappa_Hyd_Gi</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_Gia_GTP</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_Gia_GDP</ci>
<ci>mu_Pi</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_Reassoc_Gi</ci>
<apply>
<times/>
<ci>kappa_Reassoc_Gi</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_Gia_GDP</ci>
<ci>mu_Gibetagamma</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_Gi</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_InternR_M2</ci>
<apply>
<times/>
<ci>kappa_InternR_M2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_RM2_tag</ci>
<ci>mu_GRKArr</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_RM2_GRKArr</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_InternLR_M2</ci>
<apply>
<times/>
<ci>kappa_InternLR_M2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<plus/>
<ci>mu_L_RM2_tag</ci>
<ci>mu_GRKArr</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>mu_L_RM2_GRKArr</ci>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<connection>
<map_components component_1="environment" component_2="GPCR_M2_reduced"/>
<map_variables variable_1="time" variable_2="time"/>
<!--vars q_RM2_inactive and q_RM2_inactive;
vars q_L_RM2_inactive and q_L_RM2_inactive;-->
<map_variables variable_1="q_LM2" variable_2="q_LM2"/>
<map_variables variable_1="q_RM2" variable_2="q_RM2"/>
<map_variables variable_1="q_Gi" variable_2="q_Gi"/>
<map_variables variable_1="q_RM2_Gi" variable_2="q_RM2_Gi"/>
<map_variables variable_1="q_L_RM2" variable_2="q_L_RM2"/>
<map_variables variable_1="q_L_RM2_Gi" variable_2="q_L_RM2_Gi"/>
<map_variables variable_1="q_Gia_GTP" variable_2="q_Gia_GTP"/>
<map_variables variable_1="q_Gibetagamma" variable_2="q_Gibetagamma"/>
<map_variables variable_1="q_Gia_GDP" variable_2="q_Gia_GDP"/>
<map_variables variable_1="q_GTP" variable_2="q_GTP"/>
<map_variables variable_1="q_GDP" variable_2="q_GDP"/>
<map_variables variable_1="q_Pi" variable_2="q_Pi"/>
<map_variables variable_1="q_RM2_tag" variable_2="q_RM2_tag"/>
<map_variables variable_1="q_L_RM2_tag" variable_2="q_L_RM2_tag"/>
<map_variables variable_1="q_RM2_GRKArr" variable_2="q_RM2_GRKArr"/>
<map_variables variable_1="q_L_RM2_GRKArr" variable_2="q_L_RM2_GRKArr"/>
<map_variables variable_1="q_GRKArr" variable_2="q_GRKArr"/>
<!--vars v_Rswitch_M2 and v_Rswitch_M2;
vars v_LRswitch_M2 and v_LRswitch_M2;-->
<map_variables variable_1="v_C_M2" variable_2="v_C_M2"/>
<map_variables variable_1="v_R_M2" variable_2="v_R_M2"/>
<map_variables variable_1="v_L_M2" variable_2="v_L_M2"/>
<map_variables variable_1="v_Act1_Gi" variable_2="v_Act1_Gi"/>
<map_variables variable_1="v_Act2_Gi" variable_2="v_Act2_Gi"/>
<map_variables variable_1="v_Hyd_Gi" variable_2="v_Hyd_Gi"/>
<map_variables variable_1="v_Reassoc_Gi" variable_2="v_Reassoc_Gi"/>
<map_variables variable_1="v_InternR_M2" variable_2="v_InternR_M2"/>
<map_variables variable_1="v_InternLR_M2" variable_2="v_InternLR_M2"/>
</connection>
<connection>
<map_components component_1="GPCR_M2_reduced" component_2="GPCR_M2_reduced_parameters"/>
<!--vars kappa_Rswitch_M2 and kappa_Rswitch_M2;
vars kappa_LRswitch_M2 and kappa_LRswitch_M2;-->
<map_variables variable_1="kappa_C_M2" variable_2="kappa_C_M2"/>
<map_variables variable_1="kappa_R_M2" variable_2="kappa_R_M2"/>
<map_variables variable_1="kappa_L_M2" variable_2="kappa_L_M2"/>
<map_variables variable_1="kappa_Act1_Gi" variable_2="kappa_Act1_Gi"/>
<map_variables variable_1="kappa_Act2_Gi" variable_2="kappa_Act2_Gi"/>
<map_variables variable_1="kappa_Hyd_Gi" variable_2="kappa_Hyd_Gi"/>
<map_variables variable_1="kappa_Reassoc_Gi" variable_2="kappa_Reassoc_Gi"/>
<map_variables variable_1="kappa_InternR_M2" variable_2="kappa_InternR_M2"/>
<map_variables variable_1="kappa_InternLR_M2" variable_2="kappa_InternLR_M2"/>
<!--vars K_RM2_inactive and K_RM2_inactive;
vars K_L_RM2_inactive and K_L_RM2_inactive;-->
<map_variables variable_1="K_LM2" variable_2="K_LM2"/>
<map_variables variable_1="K_RM2" variable_2="K_RM2"/>
<map_variables variable_1="K_Gi" variable_2="K_Gi"/>
<map_variables variable_1="K_RM2_Gi" variable_2="K_RM2_Gi"/>
<map_variables variable_1="K_L_RM2" variable_2="K_L_RM2"/>
<map_variables variable_1="K_L_RM2_Gi" variable_2="K_L_RM2_Gi"/>
<map_variables variable_1="K_Gia_GTP" variable_2="K_Gia_GTP"/>
<map_variables variable_1="K_Gibetagamma" variable_2="K_Gibetagamma"/>
<map_variables variable_1="K_Gia_GDP" variable_2="K_Gia_GDP"/>
<map_variables variable_1="K_GTP" variable_2="K_GTP"/>
<map_variables variable_1="K_GDP" variable_2="K_GDP"/>
<map_variables variable_1="K_Pi" variable_2="K_Pi"/>
<map_variables variable_1="K_RM2_tag" variable_2="K_RM2_tag"/>
<map_variables variable_1="K_L_RM2_tag" variable_2="K_L_RM2_tag"/>
<map_variables variable_1="K_RM2_GRKArr" variable_2="K_RM2_GRKArr"/>
<map_variables variable_1="K_L_RM2_GRKArr" variable_2="K_L_RM2_GRKArr"/>
<map_variables variable_1="K_GRKArr" variable_2="K_GRKArr"/>
</connection>
<connection>
<map_components component_1="constants" component_2="GPCR_M2_reduced"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
</connection>
</model>
