- Author:
- ssaf006 <ssaf006@aucklanduni.ac.nz>
- Date:
- 2019-06-06 20:47:37+12:00
- Desc:
- fixed the nvc model
- Permanent Source URI:
- http://models.cellml.org/workspace/4bd/rawfile/4e32e1754f465d7374b6a52d1371fba0a346c5c4/bondgraph/Metabolic_module_Mnm1.cellml
<?xml version='1.0'?>
<!--parameter set p9-->
<model name="Metabolic_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"/>
<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>
<component name="Blood">
<variable name="t" public_interface="in" units="second"/>
<variable name="inputDelay" public_interface="in" units="amount"/>
<variable initial_value="0" name="oHb" public_interface="out" units="amount"/>
<variable initial_value="0" name="dHb" public_interface="out" units="amount"/>
<variable initial_value="0" name="O_2" public_interface="out" units="amount"/>
<variable initial_value="0" name="glucose" public_interface="out" units="amount"/>
<variable name="k_flow" public_interface="out" units="per_sec"/>
<variable initial_value="10" name="G_body" units="amount"/>
<variable initial_value="10" name="O_2_body" units="amount"/>
<variable initial_value="10" name="oHb_body" units="amount"/>
<variable initial_value="10" name="dHb_body" units="amount"/>
<variable name="v1_f" units="amount_per_sec"/>
<variable name="v1_b" units="amount_per_sec"/>
<variable name="vin_oHb" units="amount_per_sec"/>
<variable name="vout_oHb" units="amount_per_sec"/>
<variable name="vin_dHb" units="amount_per_sec"/>
<variable name="vout_dHb" units="amount_per_sec"/>
<variable name="vin_G" units="amount_per_sec"/>
<variable name="vout_G" units="amount_per_sec"/>
<variable name="vin_O2" units="amount_per_sec"/>
<variable name="vout_O2" units="amount_per_sec"/>
<variable name="v_basal" units="amount_per_sec"/>
<variable name="v_stim" units="amount_per_sec"/>
<variable name="HETE" public_interface="in" units="amount"/>
<variable name="EET" public_interface="in" units="amount"/>
<variable name="PG" public_interface="in" units="amount"/>
<variable name="cGMP" public_interface="in" units="amount"/>
<variable name="Stimulating" public_interface="out" units="per_sec"/>
<variable name="Inhibiting" public_interface="out" units="per_sec"/>
<variable name="Act" units="per_sec"/>
<!-- parameters-->
<variable initial_value="1104.22097" name="k1_f" public_interface="out" units="per_sec"/>
<variable initial_value="177.331822" name="k1_b" public_interface="out" units="per_amount_per_sec"/>
<variable initial_value="1122.55652" name="k_basal" public_interface="out" units="per_amount_per_sec"/>
<variable initial_value="91.3332081" name="kflow_glucose" public_interface="out" units="per_sec"/>
<variable initial_value="4.33455353" name="proportion_1" public_interface="out" units="dimensionless"/>
<variable initial_value="2.41884183" name="proportion_2" public_interface="out" units="dimensionless"/>
<variable initial_value="40.8202198" name="b_1" public_interface="out" units="per_amount_per_sec"/>
<variable initial_value="16.5399601" name="b_2" public_interface="out" units="per_amount_per_sec"/>
<variable initial_value="18.6398884" name="b_3" public_interface="out" units="per_amount_per_sec"/>
<variable initial_value="4.83994729" name="b_4" public_interface="out" units="per_amount_per_sec"/>
<variable initial_value="1e-18" name="k_i" public_interface="out" units="per_amount2_per_sec"/>
<!-- vasodilation-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Stimulating</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>b_1</ci>
<ci>cGMP</ci>
</apply>
<apply>
<times/>
<ci>b_2</ci>
<ci>PG</ci>
</apply>
<apply>
<times/>
<ci>b_3</ci>
<ci>EET</ci>
</apply>
</apply>
</apply>
<!-- vasoconstriction-->
<apply>
<eq/>
<ci>Inhibiting</ci>
<apply>
<times/>
<ci>b_4</ci>
<ci>HETE</ci>
</apply>
</apply>
<!-- signal substance effect on blood flow-->
<apply>
<eq/>
<ci>Act</ci>
<apply>
<minus/>
<ci>Stimulating</ci>
<ci>Inhibiting</ci>
</apply>
</apply>
<!-- blood flow-->
<apply>
<eq/>
<ci>k_flow</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>kflow_glucose</ci>
<ci>Stimulating</ci>
</apply>
<ci>Inhibiting</ci>
</apply>
</apply>
<!-- rate of releasing oxyhemoglobin into oxygen and deoxyhemoglobin-->
<apply>
<eq/>
<ci>v1_f</ci>
<apply>
<times/>
<ci>k1_f</ci>
<ci>oHb</ci>
</apply>
</apply>
<!-- rate of binding oxygen and deoxyhemoglobin into oxyhemoglobin-->
<apply>
<eq/>
<ci>v1_b</ci>
<apply>
<times/>
<ci>k1_b</ci>
<ci>dHb</ci>
<ci>O_2</ci>
</apply>
</apply>
<!-- oxyhemoglobin influx-->
<apply>
<eq/>
<ci>vin_oHb</ci>
<apply>
<times/>
<ci>oHb_body</ci>
<ci>k_flow</ci>
</apply>
</apply>
<!-- oxyhemoglobin outflux-->
<apply>
<eq/>
<ci>vout_oHb</ci>
<apply>
<times/>
<ci>oHb</ci>
<ci>k_flow</ci>
</apply>
</apply>
<!-- deoxyhemoglobin influx-->
<apply>
<eq/>
<ci>vin_dHb</ci>
<apply>
<times/>
<ci>dHb_body</ci>
<ci>k_flow</ci>
</apply>
</apply>
<!-- deoxyhemoglobin outflux-->
<apply>
<eq/>
<ci>vout_dHb</ci>
<apply>
<times/>
<ci>dHb</ci>
<ci>k_flow</ci>
</apply>
</apply>
<!-- glucose influx-->
<apply>
<eq/>
<ci>vin_G</ci>
<apply>
<times/>
<ci>G_body</ci>
<ci>k_flow</ci>
</apply>
</apply>
<!-- glucose outflux-->
<apply>
<eq/>
<ci>vout_G</ci>
<apply>
<times/>
<ci>glucose</ci>
<ci>k_flow</ci>
</apply>
</apply>
<!-- oxygen influx-->
<apply>
<eq/>
<ci>vin_O2</ci>
<apply>
<times/>
<ci>O_2_body</ci>
<ci>k_flow</ci>
</apply>
</apply>
<!-- oxygen outflux-->
<apply>
<eq/>
<ci>vout_O2</ci>
<apply>
<times/>
<ci>O_2</ci>
<ci>k_flow</ci>
</apply>
</apply>
<!-- basal metabolism-->
<apply>
<eq/>
<ci>v_basal</ci>
<apply>
<times/>
<ci>k_basal</ci>
<apply>
<power/>
<ci>O_2</ci>
<ci>proportion_1</ci>
</apply>
<ci>glucose</ci>
</apply>
</apply>
<!-- metabolism during stimulation-->
<apply>
<eq/>
<ci>v_stim</ci>
<apply>
<times/>
<ci>inputDelay</ci>
<apply>
<power/>
<ci>O_2</ci>
<ci>proportion_2</ci>
</apply>
<ci>glucose</ci>
</apply>
</apply>
<!-- change in oxyhemoglobin level-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>oHb</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v1_b</ci>
<ci>v1_f</ci>
</apply>
<ci>vin_oHb</ci>
</apply>
<ci>vout_oHb</ci>
</apply>
</apply>
<!-- change in deoxyhemoglobin level-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>dHb</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v1_f</ci>
<ci>v1_b</ci>
</apply>
<ci>vin_dHb</ci>
</apply>
<ci>vout_dHb</ci>
</apply>
</apply>
<!-- change in oxygen level-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>O_2</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v1_f</ci>
<ci>v1_b</ci>
</apply>
<apply>
<times/>
<ci>v_basal</ci>
<ci>proportion_1</ci>
</apply>
</apply>
<apply>
<times/>
<ci>v_stim</ci>
<ci>proportion_2</ci>
</apply>
</apply>
<ci>vin_O2</ci>
</apply>
<ci>vout_O2</ci>
</apply>
</apply>
<!-- change in glucose level-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>glucose</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>vin_G</ci>
<ci>vout_G</ci>
</apply>
<ci>v_basal</ci>
</apply>
<ci>v_stim</ci>
</apply>
</apply>
</math>
</component>
</model>
