Location: BG_cAMP @ b22d4b553a41 / BG_cAMP.cellml

Author:
Shelley Fong <s.fong@auckland.ac.nz>
Date:
2021-05-31 18:20:34+12:00
Desc:
increased guess of forward kinetic rate constant in all reactions
Permanent Source URI:
https://staging.physiomeproject.org/workspace/674/rawfile/b22d4b553a41f8ed56c3e4a790af928d6bbe8213/BG_cAMP.cellml

<?xml version='1.0' encoding='UTF-8'?>
<!-- BG model of cAMP, based on saucerman2003, v1
 20may21 initialising file 
 24may21: basal reaction with AC. 
 given Km and kcat for kinetic versions of enzyme reaction, recreate BG parameters by setting k1+ = 1, k2+ = 0, K_E = 1 to be able to solve for K_S, K_C, kappa1 and kappa2 
 See Pan thesis section 4.2 -->
<model name="BG_cAMP" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
    <units name="per_sec">
        <unit exponent="-1" units="second"/>
    </units>
    <units name="J_per_K_per_mol">
        <unit units="joule"/>
        <unit exponent="-1" units="kelvin"/>
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="C_per_mol">
        <unit units="coulomb"/>
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="fmol">
        <unit prefix="femto" units="mole"/>
    </units>
    <units name="fmol_per_L">
        <unit prefix="femto" units="mole"/>
        <unit exponent="-1" units="litre"/>
    </units>
    <units name="fmol_per_L_per_sec">
        <unit units="fmol_per_L"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_sec_per_fmol_per_L">
        <unit exponent="-1" units="second"/>
        <unit exponent="-1" units="fmol_per_L"/>
    </units>
    <units name="per_fmol">
        <unit exponent="-1" units="fmol"/>
    </units>
    <units name="fmol_per_sec">
        <unit units="fmol"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="J_per_mol">
        <unit units="joule"/>
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="mM">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="litre"/>
    </units>
    <units name="pL">
        <unit prefix="pico" units="litre"/>
    </units>
    <component name="environment">
        <variable initial_value="0" name="t" units="second"/>
        <!-- Physical parameters-->
        <variable initial_value="8.31" name="R" units="J_per_K_per_mol"/>
        <variable initial_value="310" name="T" units="kelvin"/>
        <variable initial_value="96485" name="F" units="C_per_mol"/>
        <!-- kinetic parameters - for comparing. (basal reaction)-->
        <variable initial_value="1.03e0" name="Km_basal_mM" units="mM"/>
        <variable initial_value="1.03e12" name="Km_basal" units="fmol_per_L"/>
        <variable initial_value="1e-12" name="k1p_basal" units="per_sec_per_fmol_per_L"/>
        <variable initial_value="0.83" name="k1m_basal" units="per_sec"/>
        <variable initial_value="0.2" name="k2p_basal" units="per_sec"/>
        <variable initial_value="2.41e-13" name="k2m_basal" units="per_sec_per_fmol_per_L"/>
        <variable initial_value="49.7e-6" name="AC_tot" units="mM"/>
        <variable initial_value="0.00263157894736842" name="ATP_kinetic" units="mM"/>
        <!-- kinetic states-->
        <variable initial_value="0.002631579" name="cAMP_kinetic" units="mM"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>cAMP_kinetic</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>k2p_basal</ci>
                        <ci>AC_tot</ci>
                        <ci>ATP_kinetic</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>Km_basal_mM</ci>
                        <ci>ATP_kinetic</ci>
                    </apply>
                </apply>
            </apply>
        </math>
        <!-- Bond graph parameters, found using pinv of kinetic parameters-->
        <variable initial_value="5.40E-07" name="K_ATP" units="per_fmol"/>
        <variable initial_value="5.40E-07" name="K_cAMP" units="per_fmol"/>
        <variable initial_value="0.021284929" name="K_AC" units="per_fmol"/>
        <variable initial_value="11840.92476" name="K_AC_ATP" units="per_fmol"/>
        <variable initial_value="51435.17655" name="K_Gs_AC" units="per_fmol"/>
        <variable initial_value="8.75E+09" name="K_Gs_AC_ATP" units="per_fmol"/>
        <variable initial_value="265.9246991" name="K_FSK_AC" units="per_fmol"/>
        <variable initial_value="1.24E+08" name="K_FSK_AC_ATP" units="per_fmol"/>
        <variable initial_value="64.60624363" name="K_PDE" units="per_fmol"/>
        <variable initial_value="45362.1928" name="K_PDE_cAMP" units="per_fmol"/>
        <variable initial_value="5.40E-07" name="K_5AMP" units="per_fmol"/>
        <variable initial_value="0.009277763" name="K_IBMX" units="per_fmol"/>
        <variable initial_value="1.80E+10" name="K_PDEinh" units="per_fmol"/>
        <variable initial_value="0.006041267" name="K_Gs" units="per_fmol"/>
        <variable initial_value="2.84E-07" name="K_FSK" units="per_fmol"/>
        <variable initial_value="8.6986E+15" name="kappa_1a" units="fmol_per_sec"/>
        <variable initial_value="1.6891E-05" name="kappa_1b" units="fmol_per_sec"/>
        <variable initial_value="3.5997E+09" name="kappa_2a" units="fmol_per_sec"/>
        <variable initial_value="9.7134E-10" name="kappa_2b" units="fmol_per_sec"/>
        <variable initial_value="6.9625E+11" name="kappa_3a" units="fmol_per_sec"/>
        <variable initial_value="8.0959E-25" name="kappa_3b" units="fmol_per_sec"/>
        <variable initial_value="2.8658E+12" name="kappa_4a" units="fmol_per_sec"/>
        <variable initial_value="1.1022E-04" name="kappa_4b" units="fmol_per_sec"/>
        <variable initial_value="1.6683E+08" name="kappa_5" units="fmol_per_sec"/>
        <variable initial_value="7.7768E+11" name="kappa_6" units="fmol_per_sec"/>
        <variable initial_value="1.6546E+16" name="kappa_7" units="fmol_per_sec"/>
        <!-- TEMPLATES
        var ATP: per_fmol {init: };
        var cAMP: per_fmol {init: };
        var AC: per_fmol {init: };
        var ATP_AC: per_fmol {init: };
        var GsaGTPAC: per_fmol {init: };
        var GsaGTPAC_ATP: per_fmol {init: };
        var FSKAC: per_fmol {init: };
        var FSKAC_ATP: per_fmol {init: };
        var PDE: per_fmol {init: };
        var PDE_cAMP: per_fmol {init: };
        var 5AMP: per_fmol {init: };
        var IBMX: per_fmol {init: };
        var PDEinh: per_fmol {init: };
        var GsaGTP: per_fmol {init: };
        var FSK: per_fmol {init: };
        var kappa_1a: fmol_per_sec {init: };
        var kappa_1b: fmol_per_sec {init: };
        var kappa_2a: fmol_per_sec {init: };
        var kappa_2b: fmol_per_sec {init: };
        var kappa_3a: fmol_per_sec {init: };
        var kappa_3b: fmol_per_sec {init: };
        var kappa_4a: fmol_per_sec {init: };
        var kappa_4b: fmol_per_sec {init: };
        var kappa_5: fmol_per_sec {init: };
        var kappa_6: fmol_per_sec {init: };
        var kappa_7: fmol_per_sec {init: };-->
        <!-- States (and initial conditions)-->
        <!-- ORIGINAL-->
        <variable initial_value="1.90E+02" name="q_ATP" units="fmol"/>
        <variable initial_value="1.89E-03" name="q_AC" units="fmol"/>
        <variable initial_value="3.21E-02" name="q_cAMP" units="fmol"/>
        <variable initial_value="1e0" name="q_AC_ATP" units="fmol"/>
        <variable initial_value="1e-4" name="q_FSK" units="fmol"/>
        <variable initial_value="1e0" name="q_FSK_AC" units="fmol"/>
        <variable initial_value="1e0" name="q_FSK_AC_ATP" units="fmol"/>
        <variable initial_value="9.52E-04" name="q_Gs" units="fmol"/>
        <variable initial_value="1e0" name="q_Gs_AC" units="fmol"/>
        <variable initial_value="1e0" name="q_Gs_AC_ATP" units="fmol"/>
        <variable initial_value="1.48E-03" name="q_PDE" units="fmol"/>
        <variable initial_value="1e0" name="q_PDEinh" units="fmol"/>
        <variable initial_value="1e0" name="q_PDE_cAMP" units="fmol"/>
        <variable initial_value="3.80E-02" name="q_IBMX" units="fmol"/>
        <variable initial_value="1e0" name="q_5AMP" units="fmol"/>
        <!-- Concentrations - for plotting-->
        <variable name="c_ATP" units="mM"/>
        <variable name="c_cAMP" units="mM"/>
        <variable name="c_FSK" units="mM"/>
        <variable name="c_Gs" units="mM"/>
        <variable name="c_PDE" units="mM"/>
        <variable name="c_IBMX" units="mM"/>
        <variable name="c_5AMP" units="mM"/>
        <variable name="c_AC" units="mM"/>
        <variable name="c_AC_ATP" units="mM"/>
        <variable initial_value="38.0" name="vol" units="pL"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>c_ATP</ci>
                <apply>
                    <divide/>
                    <ci>q_ATP</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>c_cAMP</ci>
                <apply>
                    <divide/>
                    <ci>q_cAMP</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>c_FSK</ci>
                <apply>
                    <divide/>
                    <ci>q_FSK</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>c_IBMX</ci>
                <apply>
                    <divide/>
                    <ci>q_IBMX</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>c_Gs</ci>
                <apply>
                    <divide/>
                    <ci>q_Gs</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>c_PDE</ci>
                <apply>
                    <divide/>
                    <ci>q_PDE</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>c_5AMP</ci>
                <apply>
                    <divide/>
                    <ci>q_5AMP</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>c_AC</ci>
                <apply>
                    <divide/>
                    <ci>q_AC</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>c_AC_ATP</ci>
                <apply>
                    <divide/>
                    <ci>q_AC_ATP</ci>
                    <ci>vol</ci>
                </apply>
            </apply>
        </math>
        <!-- Bond variables-->
        <variable name="v1a" units="fmol_per_sec"/>
        <variable name="v1b" units="fmol_per_sec"/>
        <variable name="v2a" units="fmol_per_sec"/>
        <variable name="v2b" units="fmol_per_sec"/>
        <variable name="v3a" units="fmol_per_sec"/>
        <variable name="v3b" units="fmol_per_sec"/>
        <variable name="v4a" units="fmol_per_sec"/>
        <variable name="v4b" units="fmol_per_sec"/>
        <variable name="v5" units="fmol_per_sec"/>
        <variable name="v6" units="fmol_per_sec"/>
        <variable name="v7" units="fmol_per_sec"/>
        <!-- Constitutive equations-->
        <variable name="mu_ATP" units="J_per_mol"/>
        <variable name="mu_AC" units="J_per_mol"/>
        <variable name="mu_cAMP" units="J_per_mol"/>
        <variable name="mu_AC_ATP" units="J_per_mol"/>
        <variable name="mu_FSK" units="J_per_mol"/>
        <variable name="mu_FSK_AC" units="J_per_mol"/>
        <variable name="mu_FSK_AC_ATP" units="J_per_mol"/>
        <variable name="mu_Gs" units="J_per_mol"/>
        <variable name="mu_Gs_AC" units="J_per_mol"/>
        <variable name="mu_Gs_AC_ATP" units="J_per_mol"/>
        <variable name="mu_PDE" units="J_per_mol"/>
        <variable name="mu_PDEinh" units="J_per_mol"/>
        <variable name="mu_PDE_cAMP" units="J_per_mol"/>
        <variable name="mu_IBMX" units="J_per_mol"/>
        <variable name="mu_5AMP" units="J_per_mol"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>mu_ATP</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_ATP</ci>
                            <ci>q_ATP</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_AC</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_AC</ci>
                            <ci>q_AC</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_cAMP</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_cAMP</ci>
                            <ci>q_cAMP</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_AC_ATP</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_AC_ATP</ci>
                            <ci>q_AC_ATP</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_FSK_AC</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_FSK_AC</ci>
                            <ci>q_FSK_AC</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_FSK</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_FSK</ci>
                            <ci>q_FSK</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_FSK_AC_ATP</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_FSK_AC_ATP</ci>
                            <ci>q_FSK_AC_ATP</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_Gs_AC_ATP</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_Gs_AC_ATP</ci>
                            <ci>q_Gs_AC_ATP</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_Gs_AC</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_Gs_AC</ci>
                            <ci>q_Gs_AC</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_Gs</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_Gs</ci>
                            <ci>q_Gs</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_PDE_cAMP</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_PDE_cAMP</ci>
                            <ci>q_PDE_cAMP</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_PDE</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_PDE</ci>
                            <ci>q_PDE</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_IBMX</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_IBMX</ci>
                            <ci>q_IBMX</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_PDEinh</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_PDEinh</ci>
                            <ci>q_PDEinh</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>mu_5AMP</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <ci>T</ci>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <ci>K_5AMP</ci>
                            <ci>q_5AMP</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v1a</ci>
                <apply>
                    <times/>
                    <ci>kappa_1a</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_AC</ci>
                                    <ci>mu_ATP</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_AC_ATP</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v1b</ci>
                <apply>
                    <times/>
                    <ci>kappa_1b</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_AC_ATP</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_AC</ci>
                                    <ci>mu_cAMP</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v2a</ci>
                <apply>
                    <times/>
                    <ci>kappa_2a</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_Gs_AC</ci>
                                    <ci>mu_ATP</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_Gs_AC_ATP</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v2b</ci>
                <apply>
                    <times/>
                    <ci>kappa_2b</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_Gs_AC_ATP</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_Gs_AC</ci>
                                    <ci>mu_cAMP</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v3a</ci>
                <apply>
                    <times/>
                    <ci>kappa_3a</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_FSK_AC</ci>
                                    <ci>mu_ATP</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_FSK_AC_ATP</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v3b</ci>
                <apply>
                    <times/>
                    <ci>kappa_3b</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_FSK_AC_ATP</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_FSK_AC</ci>
                                    <ci>mu_cAMP</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v4a</ci>
                <apply>
                    <times/>
                    <ci>kappa_4a</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_PDE</ci>
                                    <ci>mu_cAMP</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_PDE_cAMP</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v4b</ci>
                <apply>
                    <times/>
                    <ci>kappa_4b</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_PDE_cAMP</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_PDE</ci>
                                    <ci>mu_5AMP</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v5</ci>
                <apply>
                    <times/>
                    <ci>kappa_5</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_PDE</ci>
                                    <ci>mu_IBMX</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_PDEinh</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v6</ci>
                <apply>
                    <times/>
                    <ci>kappa_6</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_AC</ci>
                                    <ci>mu_Gs</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_Gs_AC</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v7</ci>
                <apply>
                    <times/>
                    <ci>kappa_7</ci>
                    <apply>
                        <minus/>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>mu_FSK</ci>
                                    <ci>mu_AC</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <ci>mu_FSK_AC</ci>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_ATP</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <ci>v1a</ci>
                        </apply>
                        <ci>v3a</ci>
                    </apply>
                    <ci>v2a</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_AC</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <ci>v1b</ci>
                            <ci>v1a</ci>
                        </apply>
                        <ci>v6</ci>
                    </apply>
                    <ci>v7</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_AC_ATP</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v1a</ci>
                    <ci>v1b</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_cAMP</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <ci>v1b</ci>
                        <ci>v3b</ci>
                        <ci>v2b</ci>
                    </apply>
                    <ci>v4a</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_FSK</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v7</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_FSK_AC</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <ci>v7</ci>
                        <ci>v3b</ci>
                    </apply>
                    <ci>v3a</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_FSK_AC_ATP</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v3a</ci>
                    <ci>v3b</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_Gs</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v6</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_Gs_AC</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v6</ci>
                        <ci>v2a</ci>
                    </apply>
                    <ci>v2b</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_Gs_AC_ATP</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v2a</ci>
                    <ci>v2b</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_PDE_cAMP</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v4a</ci>
                    <ci>v4b</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_PDE</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <minus/>
                        <ci>v4b</ci>
                        <ci>v4a</ci>
                    </apply>
                    <ci>v5</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_IBMX</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v5</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_PDEinh</ci>
                </apply>
                <ci>v5</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_5AMP</ci>
                </apply>
                <ci>v4b</ci>
            </apply>
        </math>
    </component>
</model>