Location: BG_GPCR_M2_reduced @ e831afee2206 / parameter_finder / TEMP.cellml.txt

Shelley Fong <s.fong@auckland.ac.nz>
2022-02-18 13:28:26+13:00
Updating way cellml is written
Permanent Source URI:

def model individual_GPCR_M2_reduced as
 def import using "units_and_constants/units_BG.cellml" for
        unit mM using unit mM;
unit fmol using unit fmol;
unit per_fmol using unit per_fmol;
        unit J_per_mol using unit J_per_mol;
unit fmol_per_sec using unit fmol_per_sec;
        unit C_per_mol using unit C_per_mol;
  unit J_per_C using unit J_per_C;
        unit microm3 using unit microm3;
  unit fF using unit fF;
        unit fC using unit fC;
  unit fA using unit fA;
        unit per_second using unit per_second;
  unit millivolt using unit millivolt;
        unit per_sec using unit per_sec;
  unit J_per_K_per_mol using unit J_per_K_per_mol;
        unit fmol_per_L using unit fmol_per_L;
  unit fmol_per_L_per_sec using unit fmol_per_L_per_sec;
        unit per_sec_per_fmol_per_L using unit per_sec_per_fmol_per_L;
  unit uM using unit uM;
        unit mM_per_sec using unit mM_per_sec;
  unit uM_per_sec using unit uM_per_sec;
        unit pL using unit pL;
  unit m_to_u using unit m_to_u;
def import using "units_and_constants/constants_BG.cellml" for
            comp constants using comp constants;

    def comp environment as
    var time: second {pub: out};
    // initial values
var q_RB1_inactive_init: fmol {init: 0.00072162};
var q_L_RB1_inactive_init: fmol {init: 1e-18};
        var q_LM2_init: fmol;
        var q_RM2_init: fmol {init: 1e-18};
        var q_Gi_init: fmol {init: 0.00836};
        var q_RM2_Gi_init: fmol {init: 1e-18};
        var q_L_RM2_init: fmol {init: 1e-18};
        var q_L_RM2_Gi_init: fmol {init: 1e-18};
        var q_Gia_GTP_init: fmol {init: 1e-18};
        var q_Gibetagamma_init: fmol {init: 1e-18};
        var q_Gia_GDP_init: fmol {init: 1e-18};
        var q_GTP_init: fmol {init: 2.2};
        var q_GDP_init: fmol {init: 1.1};
        var q_Pi_init: fmol {init: 570};
        var q_RM2_tag_init: fmol {init: 1e-18};
        var q_L_RM2_tag_init: fmol {init: 1e-18};
        var q_RM2_GRKArr_init: fmol {init: 1e-18};
        var q_L_RM2_GRKArr_init: fmol {init: 1e-18};
        var q_GRKArr_init: fmol {init: 1e-18};
        // stimulus

        var stimSt: second {init: 2e2};
        var stimDur: second {init: 5e1};
        var tR: second {init: 5e1};
        var stimMag: fmol {init: 5e2};
        var stimHolding: fmol {init: 1e-8};
        var m: fmol_per_sec;

        m = stimMag/tR;

        q_LM2_init = sel
            case (time < stimSt) and (time > stimSt-tR):
            case (time >= stimSt) and (time < stimSt+stimDur):
            case (time < stimSt+tR+stimDur) and (time >= stimSt+stimDur):
// Global value
var q_RM2_inactive: fmol {pub: out};
var q_L_RM2_inactive: fmol {pub: out};
var q_LM2: fmol {pub: out};
var q_RM2: fmol {pub: out};
var q_Gi: fmol {pub: out};
var q_RM2_Gi: fmol {pub: out};
var q_L_RM2: fmol {pub: out};
var q_L_RM2_Gi: fmol {pub: out};
var q_Gia_GTP: fmol {pub: out};
var q_Gibetagamma: fmol {pub: out};
var q_Gia_GDP: fmol {pub: out};
var q_GTP: fmol {pub: out};
var q_GDP: fmol {pub: out};
var q_Pi: fmol {pub: out};
var q_RM2_tag: fmol {pub: out};
var q_L_RM2_tag: fmol {pub: out};
var q_RM2_GRKArr: fmol {pub: out};
var q_L_RM2_GRKArr: fmol {pub: out};
var q_GRKArr: fmol {pub: out};
// From submodule
var q_RM2_inactive_mGPCR_M2_reduced: fmol {pub: in};
var q_L_RM2_inactive_mGPCR_M2_reduced: fmol {pub: in};
var q_LM2_mGPCR_M2_reduced: fmol {pub: in};
var q_RM2_mGPCR_M2_reduced: fmol {pub: in};
var q_Gi_mGPCR_M2_reduced: fmol {pub: in};
var q_RM2_Gi_mGPCR_M2_reduced: fmol {pub: in};
var q_L_RM2_mGPCR_M2_reduced: fmol {pub: in};
var q_L_RM2_Gi_mGPCR_M2_reduced: fmol {pub: in};
var q_Gia_GTP_mGPCR_M2_reduced: fmol {pub: in};
var q_Gibetagamma_mGPCR_M2_reduced: fmol {pub: in};
var q_Gia_GDP_mGPCR_M2_reduced: fmol {pub: in};
var q_GTP_mGPCR_M2_reduced: fmol {pub: in};
var q_GDP_mGPCR_M2_reduced: fmol {pub: in};
var q_Pi_mGPCR_M2_reduced: fmol {pub: in};
var q_RM2_tag_mGPCR_M2_reduced: fmol {pub: in};
var q_L_RM2_tag_mGPCR_M2_reduced: fmol {pub: in};
var q_RM2_GRKArr_mGPCR_M2_reduced: fmol {pub: in};
var q_L_RM2_GRKArr_mGPCR_M2_reduced: fmol {pub: in};
var q_GRKArr_mGPCR_M2_reduced: fmol {pub: in};
q_RM2_inactive = q_RM2_inactive_mGPCR_M2_reduced + q_RM2_inactive_init;
q_L_RM2_inactive = q_L_RM2_inactive_mGPCR_M2_reduced + q_L_RM2_inactive_init;
q_LM2 = q_LM2_mGPCR_M2_reduced + q_LM2_init;
q_RM2 = q_RM2_mGPCR_M2_reduced + q_RM2_init;
q_Gi = q_Gi_mGPCR_M2_reduced + q_Gi_init;
q_RM2_Gi = q_RM2_Gi_mGPCR_M2_reduced + q_RM2_Gi_init;
q_L_RM2 = q_L_RM2_mGPCR_M2_reduced + q_L_RM2_init;
q_L_RM2_Gi = q_L_RM2_Gi_mGPCR_M2_reduced + q_L_RM2_Gi_init;
q_Gia_GTP = q_Gia_GTP_mGPCR_M2_reduced + q_Gia_GTP_init;
q_Gibetagamma = q_Gibetagamma_mGPCR_M2_reduced + q_Gibetagamma_init;
q_Gia_GDP = q_Gia_GDP_mGPCR_M2_reduced + q_Gia_GDP_init;
q_GTP = q_GTP_mGPCR_M2_reduced + q_GTP_init;
q_GDP = q_GDP_mGPCR_M2_reduced + q_GDP_init;
q_Pi = q_Pi_mGPCR_M2_reduced + q_Pi_init;
q_RM2_tag = q_RM2_tag_mGPCR_M2_reduced + q_RM2_tag_init;
q_L_RM2_tag = q_L_RM2_tag_mGPCR_M2_reduced + q_L_RM2_tag_init;
q_RM2_GRKArr = q_RM2_GRKArr_mGPCR_M2_reduced + q_RM2_GRKArr_init;
q_L_RM2_GRKArr = q_L_RM2_GRKArr_mGPCR_M2_reduced + q_L_RM2_GRKArr_init;
q_GRKArr = q_GRKArr_mGPCR_M2_reduced + q_GRKArr_init;

def comp GPCR_M2_reduced_parameters as
var kappa_Rswitch: fmol_per_sec {init: 77.8337, pub: out};
var kappa_LRswitch: fmol_per_sec {init: 0.00321972, pub: out};
var kappa_C_M2: fmol_per_sec {init: 4441.42, pub: out};
var kappa_R_M2: fmol_per_sec {init: 183727, pub: out};
var kappa_L_M2: fmol_per_sec {init: 0.0354169, pub: out};
var kappa_Act1_Gi: fmol_per_sec {init: 0.00199813, pub: out};
var kappa_Act2_Gi: fmol_per_sec {init: 1.08204, pub: out};
var kappa_Hyd_Gi: fmol_per_sec {init: 0.604042, pub: out};
var kappa_Reassoc_Gi: fmol_per_sec {init: 0.0057063, pub: out};
var kappa_InternR_M2: fmol_per_sec {init: 0.00590542, pub: out};
var kappa_InternLR_M2: fmol_per_sec {init: 3.19794, pub: out};
var K_RM2_inactive: per_fmol {init: 3.73486, pub: out};
var K_L_RM2_inactive: per_fmol {init: 90286.7, pub: out};
var K_LM2: per_fmol {init: 63.8849, pub: out};
var K_RM2: per_fmol {init: 0.0373486, pub: out};
var K_Gi: per_fmol {init: 0.0509433, pub: out};
var K_RM2_Gi: per_fmol {init: 1.96354, pub: out};
var K_L_RM2: per_fmol {init: 0.000902867, pub: out};
var K_L_RM2_Gi: per_fmol {init: 0.00362594, pub: out};
var K_Gia_GTP: per_fmol {init: 0.0385003, pub: out};
var K_Gibetagamma: per_fmol {init: 19.8214, pub: out};
var K_Gia_GDP: per_fmol {init: 9.04022, pub: out};
var K_GTP: per_fmol {init: 0.538466, pub: out};
var K_GDP: per_fmol {init: 3.39308e-06, pub: out};
var K_Pi: per_fmol {init: 1.54752e-06, pub: out};
var K_RM2_tag: per_fmol {init: 1.38022, pub: out};
var K_L_RM2_tag: per_fmol {init: 0.00254876, pub: out};
var K_RM2_GRKArr: per_fmol {init: 0.000492256, pub: out};
var K_L_RM2_GRKArr: per_fmol {init: 9.09015e-07, pub: out};
var K_GRKArr: per_fmol {init: 1036.77, pub: out};
def comp GPCR_M2_reduced as
        var time: second {pub: in};
        var R: J_per_K_per_mol {pub: in};
        var T: kelvin {pub: in};
        // parameters
var kappa_Rswitch: fmol_per_sec {pub: in};
var kappa_LRswitch: fmol_per_sec {pub: in};
var kappa_C_M2: fmol_per_sec {pub: in};
var kappa_R_M2: fmol_per_sec {pub: in};
var kappa_L_M2: fmol_per_sec {pub: in};
var kappa_Act1_Gi: fmol_per_sec {pub: in};
var kappa_Act2_Gi: fmol_per_sec {pub: in};
var kappa_Hyd_Gi: fmol_per_sec {pub: in};
var kappa_Reassoc_Gi: fmol_per_sec {pub: in};
var kappa_InternR_M2: fmol_per_sec {pub: in};
var kappa_InternLR_M2: fmol_per_sec {pub: in};
var K_RM2_inactive: per_fmol {pub: in};
var K_L_RM2_inactive: per_fmol {pub: in};
var K_LM2: per_fmol {pub: in};
var K_RM2: per_fmol {pub: in};
var K_Gi: per_fmol {pub: in};
var K_RM2_Gi: per_fmol {pub: in};
var K_L_RM2: per_fmol {pub: in};
var K_L_RM2_Gi: per_fmol {pub: in};
var K_Gia_GTP: per_fmol {pub: in};
var K_Gibetagamma: per_fmol {pub: in};
var K_Gia_GDP: per_fmol {pub: in};
var K_GTP: per_fmol {pub: in};
var K_GDP: per_fmol {pub: in};
var K_Pi: per_fmol {pub: in};
var K_RM2_tag: per_fmol {pub: in};
var K_L_RM2_tag: per_fmol {pub: in};
var K_RM2_GRKArr: per_fmol {pub: in};
var K_L_RM2_GRKArr: per_fmol {pub: in};
var K_GRKArr: per_fmol {pub: in};
// Input from global environment
var q_RM2_inactive_global: fmol {pub: in};
var q_L_RM2_inactive_global: fmol {pub: in};
var q_LM2_global: fmol {pub: in};
var q_RM2_global: fmol {pub: in};
var q_Gi_global: fmol {pub: in};
var q_RM2_Gi_global: fmol {pub: in};
var q_L_RM2_global: fmol {pub: in};
var q_L_RM2_Gi_global: fmol {pub: in};
var q_Gia_GTP_global: fmol {pub: in};
var q_Gibetagamma_global: fmol {pub: in};
var q_Gia_GDP_global: fmol {pub: in};
var q_GTP_global: fmol {pub: in};
var q_GDP_global: fmol {pub: in};
var q_Pi_global: fmol {pub: in};
var q_RM2_tag_global: fmol {pub: in};
var q_L_RM2_tag_global: fmol {pub: in};
var q_RM2_GRKArr_global: fmol {pub: in};
var q_L_RM2_GRKArr_global: fmol {pub: in};
var q_GRKArr_global: fmol {pub: in};
// Output to global environment
var q_RM2_inactive: fmol {init: 1e-16, pub: out};
var q_L_RM2_inactive: fmol {init: 1e-16, pub: out};
var q_LM2: fmol {init: 1e-16, pub: out};
var q_RM2: fmol {init: 1e-16, pub: out};
var q_Gi: fmol {init: 1e-16, pub: out};
var q_RM2_Gi: fmol {init: 1e-16, pub: out};
var q_L_RM2: fmol {init: 1e-16, pub: out};
var q_L_RM2_Gi: fmol {init: 1e-16, pub: out};
var q_Gia_GTP: fmol {init: 1e-16, pub: out};
var q_Gibetagamma: fmol {init: 1e-16, pub: out};
var q_Gia_GDP: fmol {init: 1e-16, pub: out};
var q_GTP: fmol {init: 1e-16, pub: out};
var q_GDP: fmol {init: 1e-16, pub: out};
var q_Pi: fmol {init: 1e-16, pub: out};
var q_RM2_tag: fmol {init: 1e-16, pub: out};
var q_L_RM2_tag: fmol {init: 1e-16, pub: out};
var q_RM2_GRKArr: fmol {init: 1e-16, pub: out};
var q_L_RM2_GRKArr: fmol {init: 1e-16, pub: out};
var q_GRKArr: fmol {init: 1e-16, pub: out};
// Constitutive parameters
var mu_RM2_inactive: J_per_mol;
var mu_L_RM2_inactive: J_per_mol;
var mu_LM2: J_per_mol;
var mu_RM2: J_per_mol;
var mu_Gi: J_per_mol;
var mu_RM2_Gi: J_per_mol;
var mu_L_RM2: J_per_mol;
var mu_L_RM2_Gi: J_per_mol;
var mu_Gia_GTP: J_per_mol;
var mu_Gibetagamma: J_per_mol;
var mu_Gia_GDP: J_per_mol;
var mu_GTP: J_per_mol;
var mu_GDP: J_per_mol;
var mu_Pi: J_per_mol;
var mu_RM2_tag: J_per_mol;
var mu_L_RM2_tag: J_per_mol;
var mu_RM2_GRKArr: J_per_mol;
var mu_L_RM2_GRKArr: J_per_mol;
var mu_GRKArr: J_per_mol;
var v_Rswitch: fmol_per_sec;
var v_LRswitch: fmol_per_sec;
var v_C_M2: fmol_per_sec;
var v_R_M2: fmol_per_sec;
var v_L_M2: fmol_per_sec;
var v_Act1_Gi: fmol_per_sec;
var v_Act2_Gi: fmol_per_sec;
var v_Hyd_Gi: fmol_per_sec;
var v_Reassoc_Gi: fmol_per_sec;
var v_InternR_M2: fmol_per_sec;
var v_InternLR_M2: fmol_per_sec;
mu_RM2_inactive = R*T*ln(K_RM2_inactive*q_RM2_inactive_global);
mu_L_RM2_inactive = R*T*ln(K_L_RM2_inactive*q_L_RM2_inactive_global);
mu_LM2 = R*T*ln(K_LM2*q_LM2_global);
mu_RM2 = R*T*ln(K_RM2*q_RM2_global);
mu_Gi = R*T*ln(K_Gi*q_Gi_global);
mu_RM2_Gi = R*T*ln(K_RM2_Gi*q_RM2_Gi_global);
mu_L_RM2 = R*T*ln(K_L_RM2*q_L_RM2_global);
mu_L_RM2_Gi = R*T*ln(K_L_RM2_Gi*q_L_RM2_Gi_global);
mu_Gia_GTP = R*T*ln(K_Gia_GTP*q_Gia_GTP_global);
mu_Gibetagamma = R*T*ln(K_Gibetagamma*q_Gibetagamma_global);
mu_Gia_GDP = R*T*ln(K_Gia_GDP*q_Gia_GDP_global);
mu_GTP = R*T*ln(K_GTP*q_GTP_global);
mu_GDP = R*T*ln(K_GDP*q_GDP_global);
mu_Pi = R*T*ln(K_Pi*q_Pi_global);
mu_RM2_tag = R*T*ln(K_RM2_tag*q_RM2_tag_global);
mu_L_RM2_tag = R*T*ln(K_L_RM2_tag*q_L_RM2_tag_global);
mu_RM2_GRKArr = R*T*ln(K_RM2_GRKArr*q_RM2_GRKArr_global);
mu_L_RM2_GRKArr = R*T*ln(K_L_RM2_GRKArr*q_L_RM2_GRKArr_global);
mu_GRKArr = R*T*ln(K_GRKArr*q_GRKArr_global);
v_Rswitch = kappa_Rswitch*(exp(mu_RM2_inactive/(R*T))-exp(mu_RM2/(R*T)));
v_LRswitch = kappa_LRswitch*(exp(mu_L_RM2_inactive/(R*T))-exp(mu_L_RM2/(R*T)));
        v_C_M2 = kappa_C_M2*(exp((mu_RM2+mu_Gi)/(R*T))-exp(mu_RM2_Gi/(R*T)));
        v_R_M2 = kappa_R_M2*(exp((mu_L_RM2+mu_Gi)/(R*T))-exp(mu_L_RM2_Gi/(R*T)));
        v_L_M2 = kappa_L_M2*(exp((mu_RM2_inactive+mu_LM2)/(R*T))-exp(mu_L_RM2_inactive/(R*T)));
        v_Act1_Gi = kappa_Act1_Gi*(exp((mu_RM2_Gi+mu_GTP)/(R*T))-exp((mu_Gia_GTP+mu_Gibetagamma+mu_RM2_tag+mu_GDP)/(R*T)));
        v_Act2_Gi = kappa_Act2_Gi*(exp((mu_L_RM2_Gi+mu_GTP)/(R*T))-exp((mu_Gia_GTP+mu_Gibetagamma+mu_L_RM2_tag+mu_GDP)/(R*T)));
        v_Hyd_Gi = kappa_Hyd_Gi*(exp(mu_Gia_GTP/(R*T))-exp((mu_Gia_GDP+mu_Pi)/(R*T)));
        v_Reassoc_Gi = kappa_Reassoc_Gi*(exp((mu_Gia_GDP+mu_Gibetagamma)/(R*T))-exp(mu_Gi/(R*T)));
        v_InternR_M2 = kappa_InternR_M2*(exp((mu_RM2_tag+mu_GRKArr)/(R*T))-exp(mu_RM2_GRKArr/(R*T)));
        v_InternLR_M2 = kappa_InternLR_M2*(exp((mu_L_RM2_tag+mu_GRKArr)/(R*T))-exp(mu_L_RM2_GRKArr/(R*T)));
        var v_Rsynthesis: fmol_per_sec;

        v_Rsynthesis = +v_C_M2+v_L_M2;

        ode(q_RB1_inactive, time) = -v_Rswitch-v_L_M2;
        ode(q_L_RB1_inactive, time) = -v_LRswitch;
        ode(q_LM2, time) = -v_L_M2;
        ode(q_RM2, time) = v_Rswitch-v_C_M2+1.0{dimensionless}*v_Rsynthesis;
        ode(q_Gi, time) = -v_C_M2-v_R_M2+v_Reassoc_Gi;
        ode(q_RM2_Gi, time) = v_C_M2-v_Act1_Gi;
        ode(q_L_RM2, time) = v_LRswitch-v_R_M2+v_L_M2;
        ode(q_L_RM2_Gi, time) = v_R_M2-v_Act2_Gi;
        ode(q_Gia_GTP, time) = v_Act1_Gi+v_Act2_Gi-v_Hyd_Gi;
        ode(q_Gibetagamma, time) = v_Act1_Gi+v_Act2_Gi-v_Reassoc_Gi;
        ode(q_Gia_GDP, time) = v_Hyd_Gi-v_Reassoc_Gi;
        ode(q_GTP, time) = -v_Act1_Gi-v_Act2_Gi;
        ode(q_GDP, time) = v_Act1_Gi+v_Act2_Gi;
        ode(q_Pi, time) = v_Hyd_Gi;
        ode(q_RM2_tag, time) = v_Act1_Gi-v_InternR_M2;
        ode(q_L_RM2_tag, time) = v_Act2_Gi-v_InternLR_M2;
        ode(q_RM2_GRKArr, time) = v_InternR_M2;
        ode(q_L_RM2_GRKArr, time) = v_InternLR_M2;
        ode(q_GRKArr, time) = -v_InternR_M2-v_InternLR_M2;
def map between environment and GPCR_M2_reduced for
vars time and time;
vars q_RM2_inactive_mGPCR_M2_reduced and q_RM2_inactive;
vars q_RM2_inactive and q_RM2_inactive_global;
vars q_L_RM2_inactive_mGPCR_M2_reduced and q_L_RM2_inactive;
vars q_L_RM2_inactive and q_L_RM2_inactive_global;
vars q_LM2_mGPCR_M2_reduced and q_LM2;
vars q_LM2 and q_LM2_global;
vars q_RM2_mGPCR_M2_reduced and q_RM2;
vars q_RM2 and q_RM2_global;
vars q_Gi_mGPCR_M2_reduced and q_Gi;
vars q_Gi and q_Gi_global;
vars q_RM2_Gi_mGPCR_M2_reduced and q_RM2_Gi;
vars q_RM2_Gi and q_RM2_Gi_global;
vars q_L_RM2_mGPCR_M2_reduced and q_L_RM2;
vars q_L_RM2 and q_L_RM2_global;
vars q_L_RM2_Gi_mGPCR_M2_reduced and q_L_RM2_Gi;
vars q_L_RM2_Gi and q_L_RM2_Gi_global;
vars q_Gia_GTP_mGPCR_M2_reduced and q_Gia_GTP;
vars q_Gia_GTP and q_Gia_GTP_global;
vars q_Gibetagamma_mGPCR_M2_reduced and q_Gibetagamma;
vars q_Gibetagamma and q_Gibetagamma_global;
vars q_Gia_GDP_mGPCR_M2_reduced and q_Gia_GDP;
vars q_Gia_GDP and q_Gia_GDP_global;
vars q_GTP_mGPCR_M2_reduced and q_GTP;
vars q_GTP and q_GTP_global;
vars q_GDP_mGPCR_M2_reduced and q_GDP;
vars q_GDP and q_GDP_global;
vars q_Pi_mGPCR_M2_reduced and q_Pi;
vars q_Pi and q_Pi_global;
vars q_RM2_tag_mGPCR_M2_reduced and q_RM2_tag;
vars q_RM2_tag and q_RM2_tag_global;
vars q_L_RM2_tag_mGPCR_M2_reduced and q_L_RM2_tag;
vars q_L_RM2_tag and q_L_RM2_tag_global;
vars q_RM2_GRKArr_mGPCR_M2_reduced and q_RM2_GRKArr;
vars q_RM2_GRKArr and q_RM2_GRKArr_global;
vars q_L_RM2_GRKArr_mGPCR_M2_reduced and q_L_RM2_GRKArr;
vars q_L_RM2_GRKArr and q_L_RM2_GRKArr_global;
vars q_GRKArr_mGPCR_M2_reduced and q_GRKArr;
vars q_GRKArr and q_GRKArr_global;
def map between GPCR_M2_reduced and GPCR_M2_reduced_parameters for
vars kappa_Rswitch and kappa_Rswitch;
vars kappa_LRswitch and kappa_LRswitch;
vars kappa_C_M2 and kappa_C_M2;
vars kappa_R_M2 and kappa_R_M2;
vars kappa_L_M2 and kappa_L_M2;
vars kappa_Act1_Gi and kappa_Act1_Gi;
vars kappa_Act2_Gi and kappa_Act2_Gi;
vars kappa_Hyd_Gi and kappa_Hyd_Gi;
vars kappa_Reassoc_Gi and kappa_Reassoc_Gi;
vars kappa_InternR_M2 and kappa_InternR_M2;
vars kappa_InternLR_M2 and kappa_InternLR_M2;
vars K_RM2_inactive and K_RM2_inactive;
vars K_L_RM2_inactive and K_L_RM2_inactive;
vars K_LM2 and K_LM2;
vars K_RM2 and K_RM2;
vars K_Gi and K_Gi;
vars K_RM2_Gi and K_RM2_Gi;
vars K_L_RM2 and K_L_RM2;
vars K_L_RM2_Gi and K_L_RM2_Gi;
vars K_Gia_GTP and K_Gia_GTP;
vars K_Gibetagamma and K_Gibetagamma;
vars K_Gia_GDP and K_Gia_GDP;
vars K_GTP and K_GTP;
vars K_GDP and K_GDP;
vars K_Pi and K_Pi;
vars K_RM2_tag and K_RM2_tag;
vars K_L_RM2_tag and K_L_RM2_tag;
vars K_RM2_GRKArr and K_RM2_GRKArr;
vars K_L_RM2_GRKArr and K_L_RM2_GRKArr;
vars K_GRKArr and K_GRKArr;
def map between constants and GPCR_M2_reduced for
vars R and R;
 vars T and T;