- Author:
- Soroush Safaei <ssaf006@aucklanduni.ac.nz>
- Date:
- 2019-03-06 21:20:11+13:00
- Desc:
- fixing the param
- Permanent Source URI:
- https://models.fieldml.org/workspace/4bd/rawfile/866c504b253f511a51609d674ff66944e702a8aa/Mnm1_p7.cellml
<?xml version='1.0'?>
<model name="Neurotransmitter" 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_amount" units_ref="per_amount"/>
<units name="per_sec" units_ref="per_sec"/>
<units name="amount_per_sec" units_ref="amount_per_sec"/>
<units name="per_amount_per_sec" units_ref="per_amount_per_sec"/>
<units name="per_amount2_per_sec" units_ref="per_amount2_per_sec"/>
</import>
<import xlink:href="Neurotransmitter_module_Mnm1.cellml">
<component component_ref="Neuron" name="Neuron"/>
<component component_ref="Astrocyte" name="Astrocyte"/>
</import>
<import xlink:href="Metabolic_module_Mnm1.cellml">
<component component_ref="Blood" name="Blood"/>
</import>
<component name="environment">
<variable name="t" public_interface="out" units="second"/>
</component>
<component name="Mnm1">
<variable name="t" private_interface="out" public_interface="in" units="second"/>
<variable name="stimulus" units="amount"/>
<variable initial_value="0" name="inputDelay" private_interface="out" units="amount"/>
<variable name="input_1" units="amount_per_sec"/>
<variable name="v_ID" units="amount_per_sec"/>
<variable initial_value="0" name="glu" private_interface="out" units="amount"/>
<variable name="input_2" units="amount_per_sec"/>
<variable name="glu_sink" units="amount_per_sec"/>
<variable name="y" units="dimensionless"/>
<variable name="oHb" private_interface="in" units="amount"/>
<variable name="dHb" private_interface="in" units="amount"/>
<variable name="O_2" private_interface="in" units="amount"/>
<variable name="glucose" private_interface="in" units="amount"/>
<variable name="k_flow" private_interface="in" units="per_s"/>
<variable name="kflow_glucose" private_interface="in" units="per_s"/>
<variable name="Stimulating" private_interface="in" units="dimensionless"/>
<variable name="Inhibiting" private_interface="in" units="dimensionless"/>
<!-- parameters-->
<variable initial_value="490.35175" name="k_y" units="amount_per_sec"/>
<variable initial_value="1223.49869" name="k_metabolic" units="per_sec"/>
<variable initial_value="114.703694" name="k_neurotrans" units="per_sec"/>
<variable initial_value="4.58755326" name="sink_Input" units="per_sec"/>
<variable initial_value="1.54547324" name="sink_Glu" units="per_sec"/>
<!-- stimulus input to the metabolic module-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>input_1</ci>
<apply>
<times/>
<ci>k_metabolic</ci>
<ci>stimulus</ci>
</apply>
</apply>
<!-- delay state reaction-->
<apply>
<eq/>
<ci>v_ID</ci>
<apply>
<times/>
<ci>sink_Input</ci>
<ci>inputDelay</ci>
</apply>
</apply>
<!-- stimulus input to the neurotransmitter module-->
<apply>
<eq/>
<ci>input_2</ci>
<apply>
<times/>
<ci>k_neurotrans</ci>
<ci>stimulus</ci>
</apply>
</apply>
<!-- glucose breakdown and reuptake-->
<apply>
<eq/>
<ci>glu_sink</ci>
<apply>
<times/>
<ci>glu</ci>
<ci>sink_Glu</ci>
</apply>
</apply>
<!-- stimulus input signal-->
<apply>
<eq/>
<ci>stimulus</ci>
<piecewise>
<piece>
<cn cellml:units="amount">1</cn>
<apply>
<and/>
<apply>
<gt/>
<ci>t</ci>
<cn cellml:units="second">0.5</cn>
</apply>
<apply>
<lt/>
<ci>t</ci>
<cn cellml:units="second">1.0</cn>
</apply>
</apply>
</piece>
<otherwise>
<cn cellml:units="amount">0</cn>
</otherwise>
</piecewise>
</apply>
<!-- delay state-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>inputDelay</ci>
</apply>
<apply>
<minus/>
<ci>input_1</ci>
<ci>v_ID</ci>
</apply>
</apply>
<!-- glutamate release in the synaptic cleft-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>glu</ci>
</apply>
<apply>
<minus/>
<ci>input_2</ci>
<ci>glu_sink</ci>
</apply>
</apply>
<!-- output signal-->
<apply>
<eq/>
<ci>y</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>k_y</ci>
<ci>oHb</ci>
</apply>
<ci>dHb</ci>
</apply>
</apply>
</math>
</component>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="Mnm1">
<component_ref component="Astrocyte"/>
<component_ref component="Neuron"/>
<component_ref component="Blood"/>
</component_ref>
</group>
<connection>
<map_components component_1="environment" component_2="Mnm1"/>
<map_variables variable_1="t" variable_2="t"/>
</connection>
<connection>
<map_components component_1="Mnm1" component_2="Blood"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="inputDelay" variable_2="inputDelay"/>
<map_variables variable_1="oHb" variable_2="oHb"/>
<map_variables variable_1="dHb" variable_2="dHb"/>
<map_variables variable_1="glucose" variable_2="glucose"/>
<map_variables variable_1="O_2" variable_2="O_2"/>
<map_variables variable_1="k_flow" variable_2="k_flow"/>
<map_variables variable_1="kflow_glucose" variable_2="kflow_glucose"/>
<map_variables variable_1="Stimulating" variable_2="Stimulating"/>
<map_variables variable_1="Inhibiting" variable_2="Inhibiting"/>
</connection>
<connection>
<map_components component_1="Mnm1" component_2="Astrocyte"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="glu" variable_2="glu"/>
</connection>
<connection>
<map_components component_1="Mnm1" component_2="Neuron"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="glu" variable_2="glu"/>
</connection>
<connection>
<map_components component_1="Astrocyte" component_2="Neuron"/>
<map_variables variable_1="NO" variable_2="NO"/>
</connection>
<connection>
<map_components component_1="Astrocyte" component_2="Blood"/>
<map_variables variable_1="HETE" variable_2="HETE"/>
<map_variables variable_1="EET" variable_2="EET"/>
<map_variables variable_1="PG" variable_2="PG"/>
</connection>
<connection>
<map_components component_1="Neuron" component_2="Blood"/>
<map_variables variable_1="cGMP" variable_2="cGMP"/>
</connection>
</model>
