def model individual_GPCR_B1AR_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_R_inactive_init: fmol {init: 0.0004579000e0}; var q_LR_inactive_init: fmol {init: 1e-18}; var q_LB1_init: fmol; var q_RB1_init: fmol {init: 1e-18}; var q_Gs_init: fmol {init: 0.1455400000}; var q_RB1_Gs_init: fmol {init: 1e-18}; var q_L_RB1_init: fmol {init: 1e-18}; var q_L_RB1_Gs_init: fmol {init: 1e-18}; var q_Gsa_GTP_init: fmol {init: 1e-18}; var q_Gsbetagamma_init: fmol {init: 1e-18}; var q_Gsa_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_RB1_tag_init: fmol {init: 1e-18}; var q_L_RB1_tag_init: fmol {init: 1e-18}; var q_RB1_GRKArr_init: fmol {init: 1e-18}; var q_L_RB1_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_LB1_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_R_inactive: fmol {pub: out}; var q_LR_inactive: fmol {pub: out}; var q_LB1: fmol {pub: out}; var q_RB1: fmol {pub: out}; var q_Gs: fmol {pub: out}; var q_RB1_Gs: fmol {pub: out}; var q_L_RB1: fmol {pub: out}; var q_L_RB1_Gs: fmol {pub: out}; var q_Gsa_GTP: fmol {pub: out}; var q_Gsbetagamma: fmol {pub: out}; var q_Gsa_GDP: fmol {pub: out}; var q_GTP: fmol {pub: out}; var q_GDP: fmol {pub: out}; var q_Pi: fmol {pub: out}; var q_RB1_tag: fmol {pub: out}; var q_L_RB1_tag: fmol {pub: out}; var q_RB1_GRKArr: fmol {pub: out}; var q_L_RB1_GRKArr: fmol {pub: out}; var q_GRKArr: fmol {pub: out}; // From submodule var q_R_inactive_mGPCR_B1AR_reduced: fmol {pub: in}; var q_LR_inactive_mGPCR_B1AR_reduced: fmol {pub: in}; var q_LB1_mGPCR_B1AR_reduced: fmol {pub: in}; var q_RB1_mGPCR_B1AR_reduced: fmol {pub: in}; var q_Gs_mGPCR_B1AR_reduced: fmol {pub: in}; var q_RB1_Gs_mGPCR_B1AR_reduced: fmol {pub: in}; var q_L_RB1_mGPCR_B1AR_reduced: fmol {pub: in}; var q_L_RB1_Gs_mGPCR_B1AR_reduced: fmol {pub: in}; var q_Gsa_GTP_mGPCR_B1AR_reduced: fmol {pub: in}; var q_Gsbetagamma_mGPCR_B1AR_reduced: fmol {pub: in}; var q_Gsa_GDP_mGPCR_B1AR_reduced: fmol {pub: in}; var q_GTP_mGPCR_B1AR_reduced: fmol {pub: in}; var q_GDP_mGPCR_B1AR_reduced: fmol {pub: in}; var q_Pi_mGPCR_B1AR_reduced: fmol {pub: in}; var q_RB1_tag_mGPCR_B1AR_reduced: fmol {pub: in}; var q_L_RB1_tag_mGPCR_B1AR_reduced: fmol {pub: in}; var q_RB1_GRKArr_mGPCR_B1AR_reduced: fmol {pub: in}; var q_L_RB1_GRKArr_mGPCR_B1AR_reduced: fmol {pub: in}; var q_GRKArr_mGPCR_B1AR_reduced: fmol {pub: in}; q_R_inactive = q_R_inactive_mGPCR_B1AR_reduced + q_R_inactive_init; q_LR_inactive = q_LR_inactive_mGPCR_B1AR_reduced + q_LR_inactive_init; q_LB1 = q_LB1_mGPCR_B1AR_reduced + q_LB1_init; q_RB1 = q_RB1_mGPCR_B1AR_reduced + q_RB1_init; q_Gs = q_Gs_mGPCR_B1AR_reduced + q_Gs_init; q_RB1_Gs = q_RB1_Gs_mGPCR_B1AR_reduced + q_RB1_Gs_init; q_L_RB1 = q_L_RB1_mGPCR_B1AR_reduced + q_L_RB1_init; q_L_RB1_Gs = q_L_RB1_Gs_mGPCR_B1AR_reduced + q_L_RB1_Gs_init; q_Gsa_GTP = q_Gsa_GTP_mGPCR_B1AR_reduced + q_Gsa_GTP_init; q_Gsbetagamma = q_Gsbetagamma_mGPCR_B1AR_reduced + q_Gsbetagamma_init; q_Gsa_GDP = q_Gsa_GDP_mGPCR_B1AR_reduced + q_Gsa_GDP_init; q_GTP = q_GTP_mGPCR_B1AR_reduced + q_GTP_init; q_GDP = q_GDP_mGPCR_B1AR_reduced + q_GDP_init; q_Pi = q_Pi_mGPCR_B1AR_reduced + q_Pi_init; q_RB1_tag = q_RB1_tag_mGPCR_B1AR_reduced + q_RB1_tag_init; q_L_RB1_tag = q_L_RB1_tag_mGPCR_B1AR_reduced + q_L_RB1_tag_init; q_RB1_GRKArr = q_RB1_GRKArr_mGPCR_B1AR_reduced + q_RB1_GRKArr_init; q_L_RB1_GRKArr = q_L_RB1_GRKArr_mGPCR_B1AR_reduced + q_L_RB1_GRKArr_init; q_GRKArr = q_GRKArr_mGPCR_B1AR_reduced + q_GRKArr_init; enddef; def comp GPCR_B1AR_reduced_parameters as var kappa_Rswitch: fmol_per_sec {init: 128.165, pub: out}; var kappa_LRswitch: fmol_per_sec {init: 0.00270724, pub: out}; var kappa_C_B1: fmol_per_sec {init: 3051.24, pub: out}; var kappa_R_B1: fmol_per_sec {init: 64451.7, pub: out}; var kappa_L_B1: fmol_per_sec {init: 0.000771563, pub: out}; var kappa_Act1_Gs: fmol_per_sec {init: 0.000901275, pub: out}; var kappa_Act2_Gs: fmol_per_sec {init: 10.133, pub: out}; var kappa_Hyd_Gs: fmol_per_sec {init: 0.651714, pub: out}; var kappa_Reassoc_Gs: fmol_per_sec {init: 0.00238072, pub: out}; var kappa_InternR_B1: fmol_per_sec {init: 0.00136649, pub: out}; var kappa_InternLR_B1: fmol_per_sec {init: 15.3633, pub: out}; var K_R_inactive: per_fmol {init: 2.26815, pub: out}; var K_LR_inactive: per_fmol {init: 107378, pub: out}; var K_LB1: per_fmol {init: 4828.8, pub: out}; var K_RB1: per_fmol {init: 0.0226815, pub: out}; var K_Gs: per_fmol {init: 0.122105, pub: out}; var K_RB1_Gs: per_fmol {init: 3.14397, pub: out}; var K_L_RB1: per_fmol {init: 0.00107378, pub: out}; var K_L_RB1_Gs: per_fmol {init: 0.00027964, pub: out}; var K_Gsa_GTP: per_fmol {init: 0.0356841, pub: out}; var K_Gsbetagamma: per_fmol {init: 28.4425, pub: out}; var K_Gsa_GDP: per_fmol {init: 15.1005, pub: out}; var K_GTP: per_fmol {init: 4.77164, pub: out}; var K_GDP: per_fmol {init: 1.61737e-06, pub: out}; var K_Pi: per_fmol {init: 8.58685e-07, pub: out}; var K_RB1_tag: per_fmol {init: 4.82676, pub: out}; var K_L_RB1_tag: per_fmol {init: 0.000429315, pub: out}; var K_RB1_GRKArr: per_fmol {init: 0.00212734, pub: out}; var K_L_RB1_GRKArr: per_fmol {init: 1.89216e-07, pub: out}; var K_GRKArr: per_fmol {init: 1281.21, pub: out}; enddef; def comp GPCR_B1AR_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_B1: fmol_per_sec {pub: in}; var kappa_R_B1: fmol_per_sec {pub: in}; var kappa_L_B1: fmol_per_sec {pub: in}; var kappa_Act1_Gs: fmol_per_sec {pub: in}; var kappa_Act2_Gs: fmol_per_sec {pub: in}; var kappa_Hyd_Gs: fmol_per_sec {pub: in}; var kappa_Reassoc_Gs: fmol_per_sec {pub: in}; var kappa_InternR_B1: fmol_per_sec {pub: in}; var kappa_InternLR_B1: fmol_per_sec {pub: in}; var K_R_inactive: per_fmol {pub: in}; var K_LR_inactive: per_fmol {pub: in}; var K_LB1: per_fmol {pub: in}; var K_RB1: per_fmol {pub: in}; var K_Gs: per_fmol {pub: in}; var K_RB1_Gs: per_fmol {pub: in}; var K_L_RB1: per_fmol {pub: in}; var K_L_RB1_Gs: per_fmol {pub: in}; var K_Gsa_GTP: per_fmol {pub: in}; var K_Gsbetagamma: per_fmol {pub: in}; var K_Gsa_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_RB1_tag: per_fmol {pub: in}; var K_L_RB1_tag: per_fmol {pub: in}; var K_RB1_GRKArr: per_fmol {pub: in}; var K_L_RB1_GRKArr: per_fmol {pub: in}; var K_GRKArr: per_fmol {pub: in}; // Input from global environment var q_R_inactive_global: fmol {pub: in}; var q_LR_inactive_global: fmol {pub: in}; var q_LB1_global: fmol {pub: in}; var q_RB1_global: fmol {pub: in}; var q_Gs_global: fmol {pub: in}; var q_RB1_Gs_global: fmol {pub: in}; var q_L_RB1_global: fmol {pub: in}; var q_L_RB1_Gs_global: fmol {pub: in}; var q_Gsa_GTP_global: fmol {pub: in}; var q_Gsbetagamma_global: fmol {pub: in}; var q_Gsa_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_RB1_tag_global: fmol {pub: in}; var q_L_RB1_tag_global: fmol {pub: in}; var q_RB1_GRKArr_global: fmol {pub: in}; var q_L_RB1_GRKArr_global: fmol {pub: in}; var q_GRKArr_global: fmol {pub: in}; // Output to global environment var q_R_inactive: fmol {init: 1e-16, pub: out}; var q_LR_inactive: fmol {init: 1e-16, pub: out}; var q_LB1: fmol {init: 1e-16, pub: out}; var q_RB1: fmol {init: 1e-16, pub: out}; var q_Gs: fmol {init: 1e-16, pub: out}; var q_RB1_Gs: fmol {init: 1e-16, pub: out}; var q_L_RB1: fmol {init: 1e-16, pub: out}; var q_L_RB1_Gs: fmol {init: 1e-16, pub: out}; var q_Gsa_GTP: fmol {init: 1e-16, pub: out}; var q_Gsbetagamma: fmol {init: 1e-16, pub: out}; var q_Gsa_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_RB1_tag: fmol {init: 1e-16, pub: out}; var q_L_RB1_tag: fmol {init: 1e-16, pub: out}; var q_RB1_GRKArr: fmol {init: 1e-16, pub: out}; var q_L_RB1_GRKArr: fmol {init: 1e-16, pub: out}; var q_GRKArr: fmol {init: 1e-16, pub: out}; // Constitutive parameters var mu_R_inactive: J_per_mol; var mu_LR_inactive: J_per_mol; var mu_LB1: J_per_mol; var mu_RB1: J_per_mol; var mu_Gs: J_per_mol; var mu_RB1_Gs: J_per_mol; var mu_L_RB1: J_per_mol; var mu_L_RB1_Gs: J_per_mol; var mu_Gsa_GTP: J_per_mol; var mu_Gsbetagamma: J_per_mol; var mu_Gsa_GDP: J_per_mol; var mu_GTP: J_per_mol; var mu_GDP: J_per_mol; var mu_Pi: J_per_mol; var mu_RB1_tag: J_per_mol; var mu_L_RB1_tag: J_per_mol; var mu_RB1_GRKArr: J_per_mol; var mu_L_RB1_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_B1: fmol_per_sec; var v_R_B1: fmol_per_sec; var v_L_B1: fmol_per_sec; var v_Act1_Gs: fmol_per_sec; var v_Act2_Gs: fmol_per_sec; var v_Hyd_Gs: fmol_per_sec; var v_Reassoc_Gs: fmol_per_sec; var v_InternR_B1: fmol_per_sec; var v_InternLR_B1: fmol_per_sec; mu_R_inactive = R*T*ln(K_R_inactive*q_R_inactive_global); mu_LR_inactive = R*T*ln(K_LR_inactive*q_LR_inactive_global); mu_LB1 = R*T*ln(K_LB1*q_LB1_global); mu_RB1 = R*T*ln(K_RB1*q_RB1_global); mu_Gs = R*T*ln(K_Gs*q_Gs_global); mu_RB1_Gs = R*T*ln(K_RB1_Gs*q_RB1_Gs_global); mu_L_RB1 = R*T*ln(K_L_RB1*q_L_RB1_global); mu_L_RB1_Gs = R*T*ln(K_L_RB1_Gs*q_L_RB1_Gs_global); mu_Gsa_GTP = R*T*ln(K_Gsa_GTP*q_Gsa_GTP_global); mu_Gsbetagamma = R*T*ln(K_Gsbetagamma*q_Gsbetagamma_global); mu_Gsa_GDP = R*T*ln(K_Gsa_GDP*q_Gsa_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_RB1_tag = R*T*ln(K_RB1_tag*q_RB1_tag_global); mu_L_RB1_tag = R*T*ln(K_L_RB1_tag*q_L_RB1_tag_global); mu_RB1_GRKArr = R*T*ln(K_RB1_GRKArr*q_RB1_GRKArr_global); mu_L_RB1_GRKArr = R*T*ln(K_L_RB1_GRKArr*q_L_RB1_GRKArr_global); mu_GRKArr = R*T*ln(K_GRKArr*q_GRKArr_global); v_Rswitch = kappa_Rswitch*(exp(mu_RB1_inactive/(R*T))-exp(mu_RB1/(R*T))); v_LRswitch = kappa_LRswitch*(exp(mu_L_RB1_inactive/(R*T))-exp(mu_L_RB1/(R*T))); v_C_B1 = kappa_C_B1*(exp((mu_RB1+mu_Gs)/(R*T))-exp(mu_RB1_Gs/(R*T))); v_R_B1 = kappa_R_B1*(exp((mu_L_RB1+mu_Gs)/(R*T))-exp(mu_L_RB1_Gs/(R*T))); v_L_B1 = kappa_L_B1*(exp((mu_RB1_inactive+mu_LB1)/(R*T))-exp(mu_L_RB1_inactive/(R*T))); v_Act1_Gs = kappa_Act1_Gs*(exp((mu_RB1_Gs+mu_GTP)/(R*T))-exp((mu_Gsa_GTP+mu_Gsbetagamma+mu_RB1_tag+mu_GDP)/(R*T))); v_Act2_Gs = kappa_Act2_Gs*(exp((mu_L_RB1_Gs+mu_GTP)/(R*T))-exp((mu_Gsa_GTP+mu_Gsbetagamma+mu_L_RB1_tag+mu_GDP)/(R*T))); v_Hyd_Gs = kappa_Hyd_Gs*(exp(mu_Gsa_GTP/(R*T))-exp((mu_Gsa_GDP+mu_Pi)/(R*T))); v_Reassoc_Gs = kappa_Reassoc_Gs*(exp((mu_Gsa_GDP+mu_Gsbetagamma)/(R*T))-exp(mu_Gs/(R*T))); v_InternR_B1 = kappa_InternR_B1*(exp((mu_RB1_tag+mu_GRKArr)/(R*T))-exp(mu_RB1_GRKArr/(R*T))); v_InternLR_B1 = kappa_InternLR_B1*(exp((mu_L_RB1_tag+mu_GRKArr)/(R*T))-exp(mu_L_RB1_GRKArr/(R*T))); var v_Rsynthesis: fmol_per_sec; v_Rsynthesis = +v_C_B1+v_L_B1; ode(q_R_inactive, time) = -v_Rswitch-v_L_B1; ode(q_LR_inactive, time) = -v_LRswitch; ode(q_LB1, time) = -v_L_B1; ode(q_RB1, time) = v_Rswitch-v_C_B1+1.0{dimensionless}*v_Rsynthesis; ode(q_Gs, time) = -v_C_B1-v_R_B1+v_Reassoc_Gs; ode(q_RB1_Gs, time) = v_C_B1-v_Act1_Gs; ode(q_L_RB1, time) = v_LRswitch-v_R_B1+v_L_B1; ode(q_L_RB1_Gs, time) = v_R_B1-v_Act2_Gs; ode(q_Gsa_GTP, time) = v_Act1_Gs+v_Act2_Gs-v_Hyd_Gs; ode(q_Gsbetagamma, time) = v_Act1_Gs+v_Act2_Gs-v_Reassoc_Gs; ode(q_Gsa_GDP, time) = v_Hyd_Gs-v_Reassoc_Gs; ode(q_GTP, time) = -v_Act1_Gs-v_Act2_Gs; ode(q_GDP, time) = v_Act1_Gs+v_Act2_Gs; ode(q_Pi, time) = v_Hyd_Gs; ode(q_RB1_tag, time) = v_Act1_Gs-v_InternR_B1; ode(q_L_RB1_tag, time) = v_Act2_Gs-v_InternLR_B1; ode(q_RB1_GRKArr, time) = v_InternR_B1; ode(q_L_RB1_GRKArr, time) = v_InternLR_B1; ode(q_GRKArr, time) = -v_InternR_B1-v_InternLR_B1; enddef; def map between environment and GPCR_B1AR_reduced for vars time and time; vars q_R_inactive_mGPCR_B1AR_reduced and q_R_inactive; vars q_R_inactive and q_R_inactive_global; vars q_LR_inactive_mGPCR_B1AR_reduced and q_LR_inactive; vars q_LR_inactive and q_LR_inactive_global; vars q_LB1_mGPCR_B1AR_reduced and q_LB1; vars q_LB1 and q_LB1_global; vars q_RB1_mGPCR_B1AR_reduced and q_RB1; vars q_RB1 and q_RB1_global; vars q_Gs_mGPCR_B1AR_reduced and q_Gs; vars q_Gs and q_Gs_global; vars q_RB1_Gs_mGPCR_B1AR_reduced and q_RB1_Gs; vars q_RB1_Gs and q_RB1_Gs_global; vars q_L_RB1_mGPCR_B1AR_reduced and q_L_RB1; vars q_L_RB1 and q_L_RB1_global; vars q_L_RB1_Gs_mGPCR_B1AR_reduced and q_L_RB1_Gs; vars q_L_RB1_Gs and q_L_RB1_Gs_global; vars q_Gsa_GTP_mGPCR_B1AR_reduced and q_Gsa_GTP; vars q_Gsa_GTP and q_Gsa_GTP_global; vars q_Gsbetagamma_mGPCR_B1AR_reduced and q_Gsbetagamma; vars q_Gsbetagamma and q_Gsbetagamma_global; vars q_Gsa_GDP_mGPCR_B1AR_reduced and q_Gsa_GDP; vars q_Gsa_GDP and q_Gsa_GDP_global; vars q_GTP_mGPCR_B1AR_reduced and q_GTP; vars q_GTP and q_GTP_global; vars q_GDP_mGPCR_B1AR_reduced and q_GDP; vars q_GDP and q_GDP_global; vars q_Pi_mGPCR_B1AR_reduced and q_Pi; vars q_Pi and q_Pi_global; vars q_RB1_tag_mGPCR_B1AR_reduced and q_RB1_tag; vars q_RB1_tag and q_RB1_tag_global; vars q_L_RB1_tag_mGPCR_B1AR_reduced and q_L_RB1_tag; vars q_L_RB1_tag and q_L_RB1_tag_global; vars q_RB1_GRKArr_mGPCR_B1AR_reduced and q_RB1_GRKArr; vars q_RB1_GRKArr and q_RB1_GRKArr_global; vars q_L_RB1_GRKArr_mGPCR_B1AR_reduced and q_L_RB1_GRKArr; vars q_L_RB1_GRKArr and q_L_RB1_GRKArr_global; vars q_GRKArr_mGPCR_B1AR_reduced and q_GRKArr; vars q_GRKArr and q_GRKArr_global; enddef; def map between GPCR_B1AR_reduced and GPCR_B1AR_reduced_parameters for vars kappa_Rswitch and kappa_Rswitch; vars kappa_LRswitch and kappa_LRswitch; vars kappa_C_B1 and kappa_C_B1; vars kappa_R_B1 and kappa_R_B1; vars kappa_L_B1 and kappa_L_B1; vars kappa_Act1_Gs and kappa_Act1_Gs; vars kappa_Act2_Gs and kappa_Act2_Gs; vars kappa_Hyd_Gs and kappa_Hyd_Gs; vars kappa_Reassoc_Gs and kappa_Reassoc_Gs; vars kappa_InternR_B1 and kappa_InternR_B1; vars kappa_InternLR_B1 and kappa_InternLR_B1; vars K_R_inactive and K_R_inactive; vars K_LR_inactive and K_LR_inactive; vars K_LB1 and K_LB1; vars K_RB1 and K_RB1; vars K_Gs and K_Gs; vars K_RB1_Gs and K_RB1_Gs; vars K_L_RB1 and K_L_RB1; vars K_L_RB1_Gs and K_L_RB1_Gs; vars K_Gsa_GTP and K_Gsa_GTP; vars K_Gsbetagamma and K_Gsbetagamma; vars K_Gsa_GDP and K_Gsa_GDP; vars K_GTP and K_GTP; vars K_GDP and K_GDP; vars K_Pi and K_Pi; vars K_RB1_tag and K_RB1_tag; vars K_L_RB1_tag and K_L_RB1_tag; vars K_RB1_GRKArr and K_RB1_GRKArr; vars K_L_RB1_GRKArr and K_L_RB1_GRKArr; vars K_GRKArr and K_GRKArr; enddef; def map between constants and GPCR_B1AR_reduced for vars R and R; vars T and T; enddef; enddef;