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; enddef; def import using "units_and_constants/constants_BG.cellml" for comp constants using comp constants; enddef; 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): stimHolding+m*(time-stimSt+tR); case (time >= stimSt) and (time < stimSt+stimDur): stimMag+stimHolding; case (time < stimSt+tR+stimDur) and (time >= stimSt+stimDur): stimHolding+-m*(time-stimSt-tR-stimDur); otherwise: stimHolding; endsel; // 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; enddef; 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}; enddef; 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; enddef; 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; enddef; 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; enddef; def map between constants and GPCR_M2_reduced for vars R and R; vars T and T; enddef; enddef;