Location: NVC @ 778f87972db2 / Neurotransmitter_module_Mnm1.cellml

Author:
Soroush Safaei <ssaf006@aucklanduni.ac.nz>
Date:
2018-03-12 20:05:34+13:00
Desc:
convert to bg
Permanent Source URI:
https://models.fieldml.org/workspace/4bd/rawfile/778f87972db2ce4cd94b41a154e877ade7bdb3dc/Neurotransmitter_module_Mnm1.cellml

<?xml version='1.0'?>
<model name="Neurotransmitter_module" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
    <import xlink:href="Units.cellml">
        <units name="amount" units_ref="amount"/>
        <units name="per_sec" units_ref="per_sec"/>
        <units name="amount_per_sec" units_ref="amount_per_sec"/>
    </import>
    <component name="Neuron">
        <variable name="t" public_interface="in" units="second"/>
        <variable name="glu" public_interface="in" units="amount"/>
        <variable initial_value="0" name="Ca_Neuro" units="amount"/>
        <variable initial_value="0" name="NO" public_interface="out" units="amount"/>
        <variable initial_value="0" name="cGMP" public_interface="out" units="amount"/>
        <variable initial_value="0" name="S1_cGMP" units="amount"/>
        <variable initial_value="0" name="S2_cGMP" units="amount"/>
        <variable initial_value="0" name="S3_cGMP" units="amount"/>
        <variable name="Glutamate_N" units="amount_per_sec"/>
        <variable name="calcium_Neuro1" units="amount_per_sec"/>
        <variable name="calcium_Neuro2" units="amount_per_sec"/>
        <variable name="NO_cGMP" units="amount_per_sec"/>
        <variable name="sink_NO" units="amount_per_sec"/>
        <variable name="v1_cGMP" units="amount_per_sec"/>
        <variable name="v2_cGMP" units="amount_per_sec"/>
        <variable name="v3_cGMP" units="amount_per_sec"/>
        <variable name="cGMP_sink" units="amount_per_sec"/>
        <!-- parameters-->
        <variable initial_value="7.39467132" name="k_3" public_interface="out" units="per_sec"/>
        <variable initial_value="0.17179358" name="k_10" public_interface="out" units="per_sec"/>
        <variable initial_value="5.88395274" name="k_11" public_interface="out" units="per_sec"/>
        <variable initial_value="350.783338" name="k_21" public_interface="out" units="per_sec"/>
        <variable initial_value="1609.96028" name="k_22" public_interface="out" units="per_sec"/>
        <variable initial_value="538.467776" name="k_23" public_interface="out" units="per_sec"/>
        <variable initial_value="9279.94197" name="Ca_Nbas" public_interface="out" units="amount_per_sec"/>
        <variable initial_value="3.88188977" name="NO_sink" public_interface="out" units="per_sec"/>
        <variable initial_value="3.21617846" name="sink_N" public_interface="out" units="per_sec"/>
        <variable initial_value="49.1725342" name="sink_C" public_interface="out" units="per_sec"/>
        <!-- calcium influx in the neuron-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>Glutamate_N</ci>
                <apply>
                    <times/>
                    <ci>k_3</ci>
                    <ci>glu</ci>
                </apply>
            </apply>
            <!-- calcium outflux in the neuron-->
            <apply>
                <eq/>
                <ci>calcium_Neuro1</ci>
                <apply>
                    <times/>
                    <ci>Ca_Neuro</ci>
                    <ci>sink_N</ci>
                </apply>
            </apply>
            <!-- calcium induced NO-->
            <apply>
                <eq/>
                <ci>calcium_Neuro2</ci>
                <apply>
                    <times/>
                    <ci>k_10</ci>
                    <ci>Ca_Neuro</ci>
                </apply>
            </apply>
            <!-- NO induced cGMP-->
            <apply>
                <eq/>
                <ci>NO_cGMP</ci>
                <apply>
                    <times/>
                    <ci>k_11</ci>
                    <ci>NO</ci>
                </apply>
            </apply>
            <!-- NO breakdown-->
            <apply>
                <eq/>
                <ci>sink_NO</ci>
                <apply>
                    <times/>
                    <ci>NO_sink</ci>
                    <ci>NO</ci>
                </apply>
            </apply>
            <!-- delay state reactions-->
            <apply>
                <eq/>
                <ci>v1_cGMP</ci>
                <apply>
                    <times/>
                    <ci>k_21</ci>
                    <ci>S1_cGMP</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v2_cGMP</ci>
                <apply>
                    <times/>
                    <ci>k_22</ci>
                    <ci>S2_cGMP</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v3_cGMP</ci>
                <apply>
                    <times/>
                    <ci>k_23</ci>
                    <ci>S3_cGMP</ci>
                </apply>
            </apply>
            <!-- cGMP breakdown-->
            <apply>
                <eq/>
                <ci>cGMP_sink</ci>
                <apply>
                    <times/>
                    <ci>cGMP</ci>
                    <ci>sink_C</ci>
                </apply>
            </apply>
            <!-- calcium influx in the neuron-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>Ca_Neuro</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>Glutamate_N</ci>
                        <ci>calcium_Neuro1</ci>
                    </apply>
                    <ci>Ca_Nbas</ci>
                </apply>
            </apply>
            <!-- change in NO level-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>NO</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>calcium_Neuro2</ci>
                    <ci>sink_NO</ci>
                </apply>
            </apply>
            <!-- cGMP effecting the blood vesels-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>cGMP</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v3_cGMP</ci>
                    <ci>cGMP_sink</ci>
                </apply>
            </apply>
            <!-- delay states-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S1_cGMP</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>NO_cGMP</ci>
                    <ci>v1_cGMP</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S2_cGMP</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v1_cGMP</ci>
                    <ci>v2_cGMP</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S3_cGMP</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v2_cGMP</ci>
                    <ci>v3_cGMP</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Astrocyte">
        <variable name="t" public_interface="in" units="second"/>
        <variable name="glu" public_interface="in" units="amount"/>
        <variable initial_value="0" name="Ca_Astro" units="amount"/>
        <variable initial_value="0" name="AA" units="amount"/>
        <variable name="NO" public_interface="in" units="amount"/>
        <variable initial_value="0" name="HETE" public_interface="out" units="amount"/>
        <variable initial_value="0" name="S1_HETE" units="amount"/>
        <variable initial_value="0" name="S2_HETE" units="amount"/>
        <variable initial_value="0" name="S3_HETE" units="amount"/>
        <variable initial_value="0" name="PG" public_interface="out" units="amount"/>
        <variable initial_value="0" name="S1_PG" units="amount"/>
        <variable initial_value="0" name="S2_PG" units="amount"/>
        <variable initial_value="0" name="S3_PG" units="amount"/>
        <variable initial_value="0" name="EET" public_interface="out" units="amount"/>
        <variable initial_value="0" name="S1_EET" units="amount"/>
        <variable initial_value="0" name="S2_EET" units="amount"/>
        <variable initial_value="0" name="S3_EET" units="amount"/>
        <variable name="Glutamate_A" units="amount_per_sec"/>
        <variable name="calcium_Astro1" units="amount_per_sec"/>
        <variable name="calcium_Astro2" units="amount_per_sec"/>
        <variable name="AA_HETE" units="amount_per_sec"/>
        <variable name="AA_PG" units="amount_per_sec"/>
        <variable name="AA_EET" units="amount_per_sec"/>
        <variable name="v1_HETE" units="amount_per_sec"/>
        <variable name="v2_HETE" units="amount_per_sec"/>
        <variable name="v3_HETE" units="amount_per_sec"/>
        <variable name="HETE_sink" units="amount_per_sec"/>
        <variable name="v1_PG" units="amount_per_sec"/>
        <variable name="v2_PG" units="amount_per_sec"/>
        <variable name="v3_PG" units="amount_per_sec"/>
        <variable name="PG_sink" units="amount_per_sec"/>
        <variable name="v1_EET" units="amount_per_sec"/>
        <variable name="v2_EET" units="amount_per_sec"/>
        <variable name="v3_EET" units="amount_per_sec"/>
        <variable name="EET_sink" units="amount_per_sec"/>
        <!-- parameters-->
        <variable initial_value="1.18018801" name="k_2" public_interface="out" units="per_sec"/>
        <variable initial_value="4.14315769" name="k_4" public_interface="out" units="per_amount_per_sec"/>
        <variable initial_value="2.00837352" name="k_5" public_interface="out" units="amount_per_sec"/>
        <variable initial_value="0.888485394" name="k_6" public_interface="out" units="per_sec"/>
        <variable initial_value="3.58168609" name="k_7" public_interface="out" units="amount_per_sec"/>
        <variable initial_value="0.349135552" name="k_8" public_interface="out" units="dimensionless"/>
        <variable initial_value="1.0227789" name="k_9" public_interface="out" units="dimensionless"/>
        <variable initial_value="49.704919" name="k_12" public_interface="out" units="per_sec"/>
        <variable initial_value="0.564723084" name="k_13" public_interface="out" units="per_sec"/>
        <variable initial_value="0.920793674" name="k_14" public_interface="out" units="per_sec"/>
        <variable initial_value="485.388306" name="k_15" public_interface="out" units="per_sec"/>
        <variable initial_value="3481.42693" name="k_16" public_interface="out" units="per_sec"/>
        <variable initial_value="15170.2739" name="k_17" public_interface="out" units="per_sec"/>
        <variable initial_value="0.90356415" name="k_18" public_interface="out" units="per_sec"/>
        <variable initial_value="1.05200861" name="k_19" public_interface="out" units="per_sec"/>
        <variable initial_value="1.12487644" name="k_20" public_interface="out" units="per_sec"/>
        <variable initial_value="8.85798932" name="km_1" public_interface="out" units="amount"/>
        <variable initial_value="12.7144158" name="km_2" public_interface="out" units="amount"/>
        <variable initial_value="7.88859454" name="Ca_Abas" public_interface="out" units="amount_per_sec"/>
        <variable initial_value="0.910278754" name="sink_A" public_interface="out" units="per_sec"/>
        <variable initial_value="0.656835507" name="sink_H" public_interface="out" units="per_sec"/>
        <variable initial_value="19985.5644" name="sink_P" public_interface="out" units="per_sec"/>
        <variable initial_value="1.11950323" name="sink_E" public_interface="out" units="per_sec"/>
        <variable initial_value="37.96142" name="PL" public_interface="out" units="amount"/>
        <!-- calcium influx in the astrocyte-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>Glutamate_A</ci>
                <apply>
                    <times/>
                    <ci>k_2</ci>
                    <ci>glu</ci>
                </apply>
            </apply>
            <!-- calcium outflux in the astrocyte-->
            <apply>
                <eq/>
                <ci>calcium_Astro1</ci>
                <apply>
                    <times/>
                    <ci>Ca_Astro</ci>
                    <ci>sink_A</ci>
                </apply>
            </apply>
            <!-- calcium induced AA-->
            <apply>
                <eq/>
                <ci>calcium_Astro2</ci>
                <apply>
                    <times/>
                    <ci>PL</ci>
                    <ci>Ca_Astro</ci>
                    <ci>k_4</ci>
                </apply>
            </apply>
            <!-- AA turning into HETE-->
            <apply>
                <eq/>
                <ci>AA_HETE</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>k_5</ci>
                        <ci>AA</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>km_1</ci>
                        <apply>
                            <times/>
                            <ci>k_8</ci>
                            <ci>NO</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!-- AA turning into PG-->
            <apply>
                <eq/>
                <ci>AA_PG</ci>
                <apply>
                    <times/>
                    <ci>k_6</ci>
                    <ci>AA</ci>
                </apply>
            </apply>
            <!-- AA turning into EET-->
            <apply>
                <eq/>
                <ci>AA_EET</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>k_7</ci>
                        <ci>AA</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>km_2</ci>
                        <apply>
                            <times/>
                            <ci>k_9</ci>
                            <ci>NO</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!-- delay state reactions-->
            <apply>
                <eq/>
                <ci>v1_HETE</ci>
                <apply>
                    <times/>
                    <ci>k_12</ci>
                    <ci>S1_HETE</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v2_HETE</ci>
                <apply>
                    <times/>
                    <ci>k_13</ci>
                    <ci>S2_HETE</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v3_HETE</ci>
                <apply>
                    <times/>
                    <ci>k_14</ci>
                    <ci>S3_HETE</ci>
                </apply>
            </apply>
            <!-- HETE breakdown-->
            <apply>
                <eq/>
                <ci>HETE_sink</ci>
                <apply>
                    <times/>
                    <ci>HETE</ci>
                    <ci>sink_H</ci>
                </apply>
            </apply>
            <!-- delay state reactions-->
            <apply>
                <eq/>
                <ci>v1_PG</ci>
                <apply>
                    <times/>
                    <ci>k_15</ci>
                    <ci>S1_PG</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v2_PG</ci>
                <apply>
                    <times/>
                    <ci>k_16</ci>
                    <ci>S2_PG</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v3_PG</ci>
                <apply>
                    <times/>
                    <ci>k_17</ci>
                    <ci>S3_PG</ci>
                </apply>
            </apply>
            <!-- PG breakdown-->
            <apply>
                <eq/>
                <ci>PG_sink</ci>
                <apply>
                    <times/>
                    <ci>PG</ci>
                    <ci>sink_P</ci>
                </apply>
            </apply>
            <!-- delay state reactions-->
            <apply>
                <eq/>
                <ci>v1_EET</ci>
                <apply>
                    <times/>
                    <ci>k_18</ci>
                    <ci>S1_EET</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v2_EET</ci>
                <apply>
                    <times/>
                    <ci>k_19</ci>
                    <ci>S2_EET</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>v3_EET</ci>
                <apply>
                    <times/>
                    <ci>k_20</ci>
                    <ci>S3_EET</ci>
                </apply>
            </apply>
            <!-- EET breakdown-->
            <apply>
                <eq/>
                <ci>EET_sink</ci>
                <apply>
                    <times/>
                    <ci>EET</ci>
                    <ci>sink_E</ci>
                </apply>
            </apply>
            <!-- calcium influx in the astrocyte-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>Ca_Astro</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>Glutamate_A</ci>
                        <ci>calcium_Astro1</ci>
                    </apply>
                    <ci>Ca_Abas</ci>
                </apply>
            </apply>
            <!-- change in AA level-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>AA</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>calcium_Astro2</ci>
                    <apply>
                        <plus/>
                        <ci>AA_HETE</ci>
                        <ci>AA_PG</ci>
                        <ci>AA_EET</ci>
                    </apply>
                </apply>
            </apply>
            <!-- HETE effecting the blood vessels-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>HETE</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v3_HETE</ci>
                    <ci>HETE_sink</ci>
                </apply>
            </apply>
            <!-- delay states-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S1_HETE</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>AA_HETE</ci>
                    <ci>v1_HETE</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S2_HETE</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v1_HETE</ci>
                    <ci>v2_HETE</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S3_HETE</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v2_HETE</ci>
                    <ci>v3_HETE</ci>
                </apply>
            </apply>
            <!-- PG effecting blood vessels-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>PG</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v3_PG</ci>
                    <ci>PG_sink</ci>
                </apply>
            </apply>
            <!-- delay states-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S1_PG</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>AA_PG</ci>
                    <ci>v1_PG</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S2_PG</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v1_PG</ci>
                    <ci>v2_PG</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S3_PG</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v2_PG</ci>
                    <ci>v3_PG</ci>
                </apply>
            </apply>
            <!-- EET effeting the blood vessels-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>EET</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v3_EET</ci>
                    <ci>EET_sink</ci>
                </apply>
            </apply>
            <!-- delay states-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S1_EET</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>AA_EET</ci>
                    <ci>v1_EET</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S2_EET</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v1_EET</ci>
                    <ci>v2_EET</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>S3_EET</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v2_EET</ci>
                    <ci>v3_EET</ci>
                </apply>
            </apply>
        </math>
    </component>
</model>