Location: 12 L Platform 1 model codes @ 4c1ab73f48d7 / VSM_ftu / yang_2003.cellml

Author:
aram148 <a.rampadarath@auckland.ac.nz>
Date:
2021-11-04 16:02:42+13:00
Desc:
Updated USMC-Bursztyn model
Permanent Source URI:
http://models.cellml.org/workspace/6b0/rawfile/4c1ab73f48d7150c2094cf8017425482a1594ccf/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>