(define-module
:class "library"
:super-classes ("module")
:documentation "Superclass for all library elements.")
(define-module
:class "sbml-rules"
:super-classes ("module")
:documentation "Rules for SBML models.")
(define-variable
:class "var-rate"
:super-classes ("variable")
:value "100.0d0"
:minimum "0.0d0"
:maximum "100000.0d0"
:unit "1/h")
(define-variable
:class "concentration"
:super-classes ("variable")
:unit "mol/l")
(define-variable
:class "molar-flux"
:super-classes ("variable")
:minimum "0.0d0"
:maximum "1000.0d0"
:unit "mol/s")
(define-variable
:class "mass"
:super-classes ("variable")
:minimum "0.0d0"
:maximum "10000.0d0"
:unit "mol")
(define-variable
:class "scaling-parameter"
:super-classes ("variable")
:documentation "parameter to perform manual numeric scaling of rate-laws
etc. for optimization. It is marked as a help-variable to eventually
remove it by the equation-optimizer in Promot."
:system-theoretic "help"
:value "1.0d0")
(define-variable
:class "kinetic-parameter"
:super-classes ("variable")
:system-theoretic "real-parameter")
(define-variable
:class "stoichiometry"
:super-classes ("variable")
:documentation "Don't know what is the best solution, but these parameters
should not be changed by the user during simulation, therefore I consider them
as helps, which means they can be optimized away."
:system-theoretic "help"
:value "1.0d0"
:unit "1")
(define-variable
:class "sbml-parameter"
:super-classes ("variable")
:documentation "`Parameters' (can be real-parameters, helps, states read from sbml-files")
(define-variable
:class "volume"
:super-classes ("variable")
:minimum "0.0d0"
:maximum "10.0d0"
:unit "l")
(define-variable
:class "var-output"
:super-classes ("variable")
:system-theoretic "output")
(define-module
:class "sbml-model"
:super-classes ("module"))
(define-variable
:class "var-input"
:super-classes ("variable")
:system-theoretic "input")
(define-terminal
:class "tsig"
:super-classes ("terminal")
:icon "icons/term_concentration_in.gif"
:link-color "green"
:variables
(("c"
:is-a "concentration"
:system-theoretic "state"
:value "0.0d0")))
(define-module
:class "teusink2000_glycolysis-rules"
:super-classes ("sbml-rules")
:geometry-width "585"
:geometry-height "200"
:terminals
(("t_adp"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.5")
("t_amp"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.8333333")
("t_atp"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.6666667")
("t_p"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.33333334")
("t_sum_p"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.16666667"))
:equations
(("adp"
:relation "0.0d0 == t_adp.c-(t_sum_p.c-(t_p.c ^ 2*(1-4*root.keqak)+2*t_sum_p.c*t_p.c*(4*root.keqak+-1)+t_sum_p.c ^ 2) ^ 0.5d0)/(1-4*root.keqak)")
("amp"
:relation "0.0d0 == t_amp.c-t_sum_p.c+t_atp.c+t_adp.c")
("atp"
:relation "0.0d0 == t_atp.c-(t_p.c-t_adp.c)/2")))
(define-module
:class "compartment"
:super-classes ("module")
:variables
(("v"
:is-a "volume"
:system-theoretic "real-parameter"
:value "1.0d0")))
(define-module
:class "reaction"
:super-classes ("library")
:documentation "Superclass containing all reaction definitions.")
(define-module
:class "structure"
:super-classes ("reaction")
:documentation "Superclass containing all possible reaction structures."
:variables
(("r")))
(define-terminal
:class "tflux"
:super-classes ("tsig")
:icon "icons/term_react_blue.gif"
:link-color "blue"
:variables
(("r"
:is-a "molar-flux"
:system-theoretic "state"
:positive-vector-direction "out"
:value "0.0d0")))
(define-module
:class "s-1-1"
:super-classes ("structure")
:terminals
(("t_p1"
:is-a "tflux"
:geometry-side "right"
:geometry-position "0.5")
("t_r1"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.5"))
:variables
(("gamma_p1"
:is-a "stoichiometry"
:value "1.0d0")
("gamma_r1"
:is-a "stoichiometry"
:value "-1.0d0"))
:equations
(("e_p1"
:relation "t_p1.r == gamma_p1*r")
("e_r1"
:relation "t_r1.r == gamma_r1*r")))
(define-variable
:class "reaction-rate"
:super-classes ("molar-flux")
:unit "mol/s")
(define-module
:class "kinetics"
:super-classes ("reaction")
:documentation "Superclass containing all kinetics."
:variables
(("r"
:is-a "reaction-rate"
:documentation "reaction-rate"
:system-theoretic "help"
:report "lumped"
:unit "umol/g_dw/h")
("skal"
:is-a "scaling-parameter")))
(define-module
:class "kin-1-1_0_tglco_keqglt_kmgltglci34f"
:super-classes ("s-1-1" "kinetics")
:icon "icons/rev_l"
:variables
(("keqglt"
:is-a "kinetic-parameter")
("kmgltglci"
:is-a "kinetic-parameter")
("kmgltglco"
:is-a "kinetic-parameter")
("vmglt"
:is-a "kinetic-parameter")))
(define-module
:class "storage"
:super-classes ("library")
:documentation "Superclass for species in the network.")
(define-module
:class "species"
:super-classes ("storage")
:documentation "Constant species"
:icon "icons/species"
:parameters
(("cf.r"
:positive-vector-direction "in"
:value "0.0d0"))
:terminals
(("cf"
:is-a "tflux"
:variables
(("c"
:is-eq-to "parent.c"))))
:variables
(("c"
:is-a "concentration"
:system-theoretic "help"
:report "lumped"
:value "n/v")
("n"
:is-a "mass"
:system-theoretic "real-parameter")
("v"
:is-eq-to "parent.v")))
(define-module
:class "dyn-adapter"
:super-classes ("library")
:documentation "Superclass containing adapters. Useful for instance for propagating terminals already connected to other elements.")
(define-variable
:class "var-rate-out"
:super-classes ("var-rate" "var-output")
:absolute-error "1.0d-8")
(define-variable
:class "var-concentration"
:super-classes ("concentration")
:value "1.0d-4"
:minimum "0.0d0"
:maximum "10.0d0")
(define-variable
:class "var-concentration-out"
:super-classes ("var-concentration" "var-output")
:absolute-error "1.0d-8")
(define-variable
:class "var-state-concentration-out"
:super-classes ("var-concentration-out")
:system-theoretic "output")
(define-variable
:class "var-rate-in"
:super-classes ("var-rate" "var-input")
:absolute-error "1.0d-8")
(define-terminal
:class "term-storage-flux"
:super-classes ("terminal")
:documentation "Terminal for a flux in a storage."
:icon "icons/term_stor_blue.gif"
:link-color "blue"
:variables
(("c"
:is-a "var-concentration-out")
("r"
:is-a "var-rate-in"
:positive-vector-direction "in"
:unit "umol/g_dw/h")))
(define-variable
:class "var-concentration-in"
:super-classes ("var-concentration" "var-input")
:absolute-error "1.0d-8")
(define-terminal
:class "term-reaction-flux"
:super-classes ("terminal")
:documentation "Terminal for a flux in a reaction."
:icon "icons/term_react_blue.gif"
:link-color "blue"
:variables
(("c"
:is-a "var-concentration-in"
:unit "umol/g_dw")
("r"
:is-a "var-rate-out"
:positive-vector-direction "out"
:unit "umol/g_dw/h")))
(define-module
:class "adapter-flux"
:super-classes ("dyn-adapter")
:documentation "Flux -> Flux"
:properties
(("abstract"
:value "no"))
:icon "icons/adapter.gif"
:terminals
(("in"
:is-a "term-reaction-flux"
:geometry-side "LEFT"
:geometry-position "0.47333333"
:variables
(("r"
:is-eq-to "parent.r")))
("out"
:is-a "term-storage-flux"
:geometry-side "RIGHT"
:geometry-position "0.47333333"
:variables
(("c"
:is-eq-to "parent.c"))))
:variables
(("c"
:is-a "var-state-concentration-out")
("r"
:is-a "var-rate-out"
:positive-vector-direction "out"
:value "0.0d0"))
:equations
(("hurg"
:relation "c == in.c")
("rate"
:relation "r == out.r")))
(define-module
:class "s-2-1"
:super-classes ("structure")
:terminals
(("t_p1"
:is-a "tflux"
:geometry-side "right"
:geometry-position "0.5")
("t_r1"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.33333334")
("t_r2"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.6666667"))
:variables
(("gamma_p1"
:is-a "stoichiometry"
:value "1.0d0")
("gamma_r1"
:is-a "stoichiometry"
:value "-1.0d0")
("gamma_r2"
:is-a "stoichiometry"
:value "-1.0d0"))
:equations
(("e_p1"
:relation "t_p1.r == gamma_p1*r")
("e_r1"
:relation "t_r1.r == gamma_r1*r")
("e_r2"
:relation "t_r2.r == gamma_r2*r")))
(define-terminal
:class "tvol"
:super-classes ("terminal")
:link-color "yellow"
:variables
(("v"
:is-a "volume"
:system-theoretic "state")))
(define-module
:class "kin-2-1_0_1_ktrehalose"
:super-classes ("s-2-1" "kinetics")
:icon "icons/irrev_r"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("ktrehalose"
:is-a "kinetic-parameter")))
(define-module
:class "s-3-2"
:super-classes ("structure")
:terminals
(("t_p1"
:is-a "tflux"
:geometry-side "right"
:geometry-position "0.33333334")
("t_p2"
:is-a "tflux"
:geometry-side "right"
:geometry-position "0.6666667")
("t_r1"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.25")
("t_r2"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.5")
("t_r3"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.75"))
:variables
(("gamma_p1"
:is-a "stoichiometry"
:value "1.0d0")
("gamma_p2"
:is-a "stoichiometry"
:value "1.0d0")
("gamma_r1"
:is-a "stoichiometry"
:value "-1.0d0")
("gamma_r2"
:is-a "stoichiometry"
:value "-1.0d0")
("gamma_r3"
:is-a "stoichiometry"
:value "-1.0d0"))
:equations
(("e_p1"
:relation "t_p1.r == gamma_p1*r")
("e_p2"
:relation "t_p2.r == gamma_p2*r")
("e_r1"
:relation "t_r1.r == gamma_r1*r")
("e_r2"
:relation "t_r2.r == gamma_r2*r")
("e_r3"
:relation "t_r3.r == gamma_r3*r")))
(define-module
:class "kin-3-2_0_1_ksucc"
:super-classes ("s-3-2" "kinetics")
:icon "icons/irrev_r"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("ksucc"
:is-a "kinetic-parameter")))
(define-module
:class "s-1-2"
:super-classes ("structure")
:terminals
(("t_p1"
:is-a "tflux"
:geometry-side "right"
:geometry-position "0.33333334")
("t_p2"
:is-a "tflux"
:geometry-side "right"
:geometry-position "0.6666667")
("t_r1"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.5"))
:variables
(("gamma_p1"
:is-a "stoichiometry"
:value "1.0d0")
("gamma_p2"
:is-a "stoichiometry"
:value "1.0d0")
("gamma_r1"
:is-a "stoichiometry"
:value "-1.0d0"))
:equations
(("e_p1"
:relation "t_p1.r == gamma_p1*r")
("e_p2"
:relation "t_p2.r == gamma_p2*r")
("e_r1"
:relation "t_r1.r == gamma_r1*r")))
(define-module
:class "kin-1-2_2_1_kpep_kmpykadp_keqpyk_kmpykpyr_kmpykatpos3"
:super-classes ("s-1-2" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5")
("t_m1"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.33333334")
("t_m2"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.6666667"))
:variables
(("keqpyk"
:is-a "kinetic-parameter")
("kmpykadp"
:is-a "kinetic-parameter")
("kmpykatp"
:is-a "kinetic-parameter")
("kmpykpep"
:is-a "kinetic-parameter")
("kmpykpyr"
:is-a "kinetic-parameter")
("vmpyk"
:is-a "kinetic-parameter")))
(define-module
:class "kin-1-1_0_1_mp3g_keqpgm_kmpgmp2g6me"
:super-classes ("s-1-1" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("keqpgm"
:is-a "kinetic-parameter")
("kmpgmp2g"
:is-a "kinetic-parameter")
("kmpgmp3g"
:is-a "kinetic-parameter")
("vmpgm"
:is-a "kinetic-parameter")))
(define-module
:class "kin-1-2_2_1_kp3g_kmpgkatp_keqpgk_kmpgkbpg_kmpgkadp7d8"
:super-classes ("s-1-2" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5")
("t_m1"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.33333334")
("t_m2"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.6666667"))
:variables
(("keqpgk"
:is-a "kinetic-parameter")
("kmpgkadp"
:is-a "kinetic-parameter")
("kmpgkatp"
:is-a "kinetic-parameter")
("kmpgkbpg"
:is-a "kinetic-parameter")
("kmpgkp3g"
:is-a "kinetic-parameter")
("vmpgk"
:is-a "kinetic-parameter")))
(define-module
:class "kin-1-1_0_1_pgig6p_2_keqpgi_2_kmpgif6p_213yw"
:super-classes ("s-1-1" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("keqpgi_2"
:is-a "kinetic-parameter")
("kmpgif6p_2"
:is-a "kinetic-parameter")
("kmpgig6p_2"
:is-a "kinetic-parameter")
("vmpgi_2"
:is-a "kinetic-parameter")))
(define-module
:class "kin-2-1_3_1_vmpfk"
:super-classes ("s-2-1" "kinetics")
:icon "icons/irrev_r"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5")
("t_m1"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.25")
("t_m2"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.5")
("t_m3"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.75"))
:variables
(("vmpfk"
:is-a "kinetic-parameter")))
(define-module
:class "kin-1-2_0_1_vmpdc_npdc_kmpdcpyr"
:super-classes ("s-1-2" "kinetics")
:icon "icons/irrev_r"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("kmpdcpyr"
:is-a "kinetic-parameter")
("npdc"
:is-a "kinetic-parameter")
("vmpdc"
:is-a "kinetic-parameter")))
(define-module
:class "kin-2-1_0_1_kglycogen_3"
:super-classes ("s-2-1" "kinetics")
:icon "icons/irrev_r"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("kglycogen_3"
:is-a "kinetic-parameter")))
(define-module
:class "kin-2-1_2_1_kglci_kmglkatp_keqglk_kmglkg6p_kmglkadpuxj"
:super-classes ("s-2-1" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5")
("t_m1"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.33333334")
("t_m2"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.6666667"))
:variables
(("keqglk"
:is-a "kinetic-parameter")
("kmglkadp"
:is-a "kinetic-parameter")
("kmglkatp"
:is-a "kinetic-parameter")
("kmglkg6p"
:is-a "kinetic-parameter")
("kmglkglci"
:is-a "kinetic-parameter")
("vmglk"
:is-a "kinetic-parameter")))
(define-module
:class "s-2-2"
:super-classes ("structure")
:terminals
(("t_p1"
:is-a "tflux"
:geometry-side "right"
:geometry-position "0.33333334")
("t_p2"
:is-a "tflux"
:geometry-side "right"
:geometry-position "0.6666667")
("t_r1"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.33333334")
("t_r2"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.6666667"))
:variables
(("gamma_p1"
:is-a "stoichiometry"
:value "1.0d0")
("gamma_p2"
:is-a "stoichiometry"
:value "1.0d0")
("gamma_r1"
:is-a "stoichiometry"
:value "-1.0d0")
("gamma_r2"
:is-a "stoichiometry"
:value "-1.0d0"))
:equations
(("e_p1"
:relation "t_p1.r == gamma_p1*r")
("e_p2"
:relation "t_p2.r == gamma_p2*r")
("e_r1"
:relation "t_r1.r == gamma_r1*r")
("e_r2"
:relation "t_r2.r == gamma_r2*r")))
(define-module
:class "kin-2-2_0_1_mgapdhgap_kmgapdhnad_vmgapdhr_kmgapdhbpg_kmgapdhnadh1e35"
:super-classes ("s-2-2" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("kmgapdhbpg"
:is-a "kinetic-parameter")
("kmgapdhgap"
:is-a "kinetic-parameter")
("kmgapdhnad"
:is-a "kinetic-parameter")
("kmgapdhnadh"
:is-a "kinetic-parameter")
("vmgapdhf"
:is-a "kinetic-parameter")
("vmgapdhr"
:is-a "kinetic-parameter")))
(define-module
:class "kin-2-2_0_1_g3pdhdhap_kmg3pdhnadh_keqg3pdh_kmg3pdhgly_kmg3pdhnad15uq"
:super-classes ("s-2-2" "kinetics")
:icon "icons/irrev_r"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("keqg3pdh"
:is-a "kinetic-parameter")
("kmg3pdhdhap"
:is-a "kinetic-parameter")
("kmg3pdhgly"
:is-a "kinetic-parameter")
("kmg3pdhnad"
:is-a "kinetic-parameter")
("kmg3pdhnadh"
:is-a "kinetic-parameter")
("vmg3pdh"
:is-a "kinetic-parameter")))
(define-module
:class "kin-1-1_0_1_op2g_keqeno_kmenopep1e60"
:super-classes ("s-1-1" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("keqeno"
:is-a "kinetic-parameter")
("kmenop2g"
:is-a "kinetic-parameter")
("kmenopep"
:is-a "kinetic-parameter")
("vmeno"
:is-a "kinetic-parameter")))
(define-module
:class "s-1-_"
:super-classes ("structure")
:terminals
(("t_r1"
:is-a "tflux"
:geometry-side "left"
:geometry-position "0.5"))
:variables
(("gamma_r1"
:is-a "stoichiometry"
:value "-1.0d0"))
:equations
(("e_r1"
:relation "t_r1.r == gamma_r1*r")))
(define-module
:class "kin-1-__1_1_katpase"
:super-classes ("s-1-_" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5")
("t_m1"
:is-a "tsig"
:geometry-side "top"
:geometry-position "0.5"))
:variables
(("katpase"
:is-a "kinetic-parameter")))
(define-module
:class "kin-1-1_0_1_df16p_keqald_kmaldgap_kmalddhap_kmaldgapideg"
:super-classes ("s-1-1" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("keqald"
:is-a "kinetic-parameter")
("kmalddhap"
:is-a "kinetic-parameter")
("kmaldf16p"
:is-a "kinetic-parameter")
("kmaldgap"
:is-a "kinetic-parameter")
("kmaldgapi"
:is-a "kinetic-parameter")
("vmald"
:is-a "kinetic-parameter")))
(define-module
:class "kin-2-2_0_1_hnad_kmadhetoh_keqadh_kmadhnad_kmadhnadh_kiadhnadh_kmadhace_kiadhace_kiadhetohcon"
:super-classes ("s-2-2" "kinetics")
:icon "icons/rev_l"
:terminals
(("t_c1"
:is-a "tvol"
:geometry-side "bottom"
:geometry-position "0.5"))
:variables
(("keqadh"
:is-a "kinetic-parameter")
("kiadhace"
:is-a "kinetic-parameter")
("kiadhetoh"
:is-a "kinetic-parameter")
("kiadhnad"
:is-a "kinetic-parameter")
("kiadhnadh"
:is-a "kinetic-parameter")
("kmadhace"
:is-a "kinetic-parameter")
("kmadhetoh"
:is-a "kinetic-parameter")
("kmadhnad"
:is-a "kinetic-parameter")
("kmadhnadh"
:is-a "kinetic-parameter")
("vmadh"
:is-a "kinetic-parameter")))
(define-module
:class "adapter-vol"
:super-classes ("dyn-adapter"))
(define-module
:class "adapter-volume-pop"
:super-classes ("adapter-vol")
:parent-type "compartment"
:documentation "exports one tvol as volume out of parent"
:properties
(("abstract"
:value "no"))
:icon "icons/adapter.gif"
:terminals
(("out"
:is-a "tvol"
:geometry-side "RIGHT"
:variables
(("v"
:is-eq-to "parent.v"))))
:variables
(("v"
:is-eq-to "parent.v")))
(define-module
:class "help"
:super-classes ("library")
:documentation "Superclass containing accessories for making basic operations on signals or converting terminals.")
(define-module
:class "term-conv"
:super-classes ("help")
:documentation "Superclass containing elements for converting different terminal types.")
(define-terminal
:class "term-concentration-out"
:super-classes ("terminal")
:documentation "Output terminal for a concentration (signal)"
:icon "icons/term_concentration_out.gif"
:link-color "green"
:variables
(("c"
:is-a "var-concentration-out"
:unit "umol/g_dw")))
(define-module
:class "cf2co"
:super-classes ("term-conv")
:documentation "Flux -> Signal. Converts a cf-term to a co-term ."
:properties
(("abstract"
:value "no"))
:icon "icons/adapter"
:terminals
(("cf"
:is-a "term-reaction-flux"
:geometry-side "LEFT"
:geometry-position "0.52666664"
:variables
(("c"
:is-eq-to "parent.c")))
("co"
:is-a "term-concentration-out"
:geometry-side "RIGHT"
:geometry-position "0.52666664"))
:variables
(("c"
:is-a "var-concentration-in"))
:equations
(("free"
:relation "co.c == c")
("rate"
:relation "cf.r == 0.0d0")))
(define-module
:class "speciesb"
:super-classes ("species")
:documentation "Non-constant boundary species, value of c! (not n)
is fixed by some external rule.
The rule is necessary, n is unused."
:icon "icons/speciesb"
:parameters
(("c"
:system-theoretic "state"
:value "0")
("n"
:system-theoretic "help"
:value "c*v")))
(define-module
:class "speciesv"
:super-classes ("species")
:documentation "Species with molar balance, defined by connected reactions.
Concentration calculated using surrounding compartment"
:icon "icons/speciesv"
:parameters
(("n"
:system-theoretic "state"))
:equations
(("balc"
:relation ":diff
Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry.
Teusink,B et al.: Eur J Biochem 2000 Sep;267(17):5313-29.
The model reproduces the steady-state fluxes and metabolite concentrations of the branched model as given in Table 4 of the paper. It is derived from the model on JWS online, but has the ATP consumption in the succinate branch with the same stoichiometrie as in the publication. The model was successfully tested on copasi v.4.4(build 26).
For Vmax values, please note that there is a conversion factor of approx. 270 to convert from U/mg-protein as shown in Table 1 of the paper to mmol/(min*L_cytosol). The equilibrium constant for the ADH reaction in the paper is given for the reverse reaction (Keq = 1.45*10
4). The value used in this model is for the forward reaction: 1/Keq = 6.9*10
-5 .
Vmax parameters values used (in [mM/min] except VmGLT):
The result of the G6P steady state concentration (marked in red) differs slightly from the one given in table 4. of the publication
VmGLT
97.264
mmol/min
VmGLK
226.45
VmPGI
339.667
VmPFK
182.903
VmALD
322.258
VmGAPDH_f
1184.52
VmGAPDH_r
6549.68
VmPGK
1306.45
VmPGM
2525.81
VmENO
365.806
VmPYK
1088.71
VmPDC
174.194
VmG3PDH
70.15
Results for steady state:
Authors of the publication also mentioned a few misprints in the original article:
orig. article
this model
Fluxes[mM/min]
Glucose
88
88
Ethanol
129
129
Glycogen
6
6
Trehalose
4.8
4.8
(G6P flux through trehalose branch)
Glycerol
18.2
18.2
Succinate
3.6
3.6
Conc.[mM]
G6P
1.07
1.03
F6P
0.11
0.11
F1,6P
0.6
0.6
DHAP
0.74
0.74
3PGA
0.36
0.36
2PGA
0.04
0.04
PEP
0.07
0.07
PYR
8.52
8.52
AcAld
0.17
0.17
ATP
2.51
2.51
ADP
1.29
1.29
AMP
0.3
0.3
NAD
1.55
1.55
NADH
0.04
0.04
in the kinetic law for
ADH:
in the kinetic law for
PFK:
To make the model easier to curate, the species
ATP,
ADP and
AMP were added. These are calculated via assignment rules from the active phosphate species,
P, and the sum of all
AXP,
SUM_P.
This model originates from BioModels Database: A Database of Annotated Published Models. It is copyright (c) 2005-2009 The BioModels Team.
For more information see the terms of use.
To cite BioModels Database, please use Le Novère N., Bornstein B., Broicher A., Courtot M., Donizelli M., Dharuri H., Li L., Sauro H., Schilstra M., Shapiro B., Snoep J.L., Hucka M. (2006) BioModels Database: A Free, Centralized Database of Curated, Published, Quantitative Kinetic Models of Biochemical and Cellular Systems Nucleic Acids Res., 34: D689-D691.