- Author:
- aram148 <a.rampadarath@auckland.ac.nz>
- Date:
- 2021-12-17 11:46:34+13:00
- Desc:
- Added the Matlab 7 airway code
- Permanent Source URI:
- https://models.cellml.org/workspace/6b0/rawfile/9ac5ec0fe3f8ae60e6020b99146573defeb5a5e9/VSM_ftu/yang_2003.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="yang_2003" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<units name="per_nM_per_sec">
<unit exponent="-1" prefix="nano" units="mole"/>
<unit exponent="-1" units="second"/>
</units>
<units name="uM_per_sec">
<unit prefix="micro" units="mole"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_sec">
<unit exponent="-1" units="second"/>
</units>
<units name="uM">
<unit prefix="micro" units="mole"/>
</units>
<units name="nM">
<unit prefix="nano" units="mole"/>
</units>
<!-- def comp environment as
var time: second {pub: out};
enddef;-->
<component name="sGC">
<variable initial_value="1" name="Eb" public_interface="out" units="dimensionless"/>
<variable initial_value="0" name="E6c" public_interface="out" units="dimensionless"/>
<variable name="E5c" public_interface="out" units="dimensionless"/>
<variable initial_value="220" name="NO" public_interface="out" units="nM"/>
<variable initial_value="2.0" name="k_one" units="per_nM_per_sec"/>
<variable initial_value="100" name="km1" units="per_sec"/>
<variable initial_value="0.1" name="k_two" units="per_sec"/>
<variable initial_value="0.003" name="k_three" units="per_nM_per_sec"/>
<variable name="k_four" units="per_sec"/>
<variable initial_value="0" name="cGMP" public_interface="out" units="uM"/>
<variable name="vd" public_interface="out" units="dimensionless"/>
<variable name="vp" public_interface="out" units="dimensionless"/>
<variable initial_value="0.8520" name="Vm_sgc" units="dimensionless"/>
<variable name="Vm_pde" units="dimensionless"/>
<variable initial_value="2.0" name="Km_pde" units="uM"/>
<variable initial_value="0.0195" name="K_pde" units="per_sec"/>
<variable initial_value="0.01" name="k_dno" units="per_sec"/>
<variable initial_value="220" name="J_no" units="nM"/>
<variable initial_value="2" name="m" units="dimensionless"/>
<variable initial_value="0.011" name="K_four" units="dimensionless"/>
<variable name="time" public_interface="in" units="second"/>
<variable initial_value="0.0086" name="k_mclpb" units="per_sec"/>
<variable initial_value="0.0327" name="k_mlcpc" units="per_sec"/>
<variable initial_value="5.5" name="K_mlcpm" units="uM"/>
<variable initial_value="2" name="n_mlcoh" units="dimensionless"/>
<variable name="k2y" public_interface="out" units="per_sec"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Eb</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<ci>k_one</ci>
</apply>
<ci>Eb</ci>
<ci>NO</ci>
</apply>
<apply>
<times/>
<ci>km1</ci>
<ci>E6c</ci>
</apply>
<apply>
<times/>
<ci>k_four</ci>
<ci>E5c</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>E6c</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>k_one</ci>
<ci>Eb</ci>
<ci>NO</ci>
</apply>
<apply>
<times/>
<ci>km1</ci>
<ci>E6c</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_two</ci>
<ci>E6c</ci>
</apply>
</apply>
<apply>
<times/>
<ci>k_three</ci>
<ci>E6c</ci>
<ci>NO</ci>
</apply>
</apply>
</apply>
<!-- ode(E5c, time) = k_3*E6c*NO+k_2*E6c-k4*E5c;-->
<apply>
<eq/>
<ci>E5c</ci>
<apply>
<minus/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>Eb</ci>
</apply>
<ci>E6c</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>cGMP</ci>
</apply>
<apply>
<minus/>
<ci>vd</ci>
<ci>vp</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>NO</ci>
</apply>
<apply>
<minus/>
<ci>J_no</ci>
<apply>
<times/>
<ci>k_dno</ci>
<ci>NO</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>vd</ci>
<apply>
<times/>
<ci>Vm_sgc</ci>
<ci>E5c</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>vp</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>cGMP</ci>
<ci>Vm_pde</ci>
</apply>
<apply>
<plus/>
<ci>Km_pde</ci>
<ci>cGMP</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Vm_pde</ci>
<apply>
<times/>
<ci>K_pde</ci>
<ci>cGMP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_four</ci>
<apply>
<times/>
<ci>K_four</ci>
<apply>
<power/>
<ci>cGMP</ci>
<ci>m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k2y</ci>
<apply>
<plus/>
<ci>k_mclpb</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>k_mlcpc</ci>
<apply>
<power/>
<ci>cGMP</ci>
<ci>n_mlcoh</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>cGMP</ci>
<ci>n_mlcoh</ci>
</apply>
<apply>
<power/>
<ci>K_mlcpm</ci>
<ci>n_mlcoh</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<!-- def map between environment and sGC for
vars time and time;
enddef;-->
</model>
