begin name
leloup5
end name

begin reactions
reaction0 {1.0}$EXTERNAL = {1.0}Mb
reaction1 {1.0}$EXTERNAL = {1.0}Bc
reaction10 {1.0}$EXTERNAL = {1.0}Pc
reaction11 {1.0}Pcp = {1.0}$EXTERNAL
reaction12 {1.0}Pc = {1.0}Pcp
reaction13 {1.0}Cc + {1.0}Pc = {1.0}PCc
reaction14 {1.0}PCc = {1.0}PCcp
reaction15 {1.0}PCcp = {1.0}$EXTERNAL
reaction16 {1.0}PCc = {1.0}PCn
reaction17 {1.0}PCnp = {1.0}$EXTERNAL
reaction18 {1.0}Bcp = {1.0}$EXTERNAL
reaction19 {1.0}Bnp = {1.0}$EXTERNAL
reaction2 {1.0}Mb = {1.0}$EXTERNAL
reaction20 {1.0}$EXTERNAL = {1.0}Mc
reaction21 {1.0}PCn = {1.0}PCnp
reaction22 {1.0}Mp = {1.0}$EXTERNAL
reaction23 {1.0}PCn + {1.0}Bn = {1.0}In
reaction24 {1.0}Mb = {1.0}$EXTERNAL
reaction25 {1.0}Mc = {1.0}$EXTERNAL
reaction26 {1.0}Mp = {1.0}$EXTERNAL
reaction27 {1.0}Pc = {1.0}$EXTERNAL
reaction28 {1.0}Cc = {1.0}$EXTERNAL
reaction29 {1.0}Pcp = {1.0}$EXTERNAL
reaction3 {1.0}Bc = {1.0}Bcp
reaction30 {1.0}Ccp = {1.0}$EXTERNAL
reaction31 {1.0}PCcp = {1.0}$EXTERNAL
reaction32 {1.0}PCc = {1.0}$EXTERNAL
reaction33 {1.0}PCnp = {1.0}$EXTERNAL
reaction34 {1.0}Bc = {1.0}$EXTERNAL
reaction35 {1.0}Bcp = {1.0}$EXTERNAL
reaction36 {1.0}Bn = {1.0}Bnp
reaction37 {1.0}Bnp = {1.0}$EXTERNAL
reaction38 {1.0}In = {1.0}$EXTERNAL
reaction39 {1.0}In = {1.0}$EXTERNAL
reaction4 {1.0}Bc = {1.0}Bn
reaction40 {1.0}Bn = {1.0}$EXTERNAL
reaction41 {1.0}Bcp = {1.0}Bc
reaction42 {1.0}Bnp = {1.0}Bn
reaction43 {1.0}Ccp = {1.0}Cc
reaction44 {1.0}Pcp = {1.0}Pc
reaction45 {1.0}PCnp = {1.0}PCn
reaction46 {1.0}PCn = {1.0}$EXTERNAL
reaction47 {1.0}PCcp = {1.0}PCc
reaction48 {1.0}$EXTERNAL = {1.0}Mr
reaction49 {1.0}Mr = {1.0}$EXTERNAL
reaction5 {1.0}$EXTERNAL = {1.0}Cc
reaction50 {1.0}Mr = {1.0}$EXTERNAL
reaction51 {1.0}$EXTERNAL = {1.0}Rc
reaction52 {1.0}Rc = {1.0}Rn
reaction53 {1.0}Rc = {1.0}$EXTERNAL
reaction54 {1.0}Rc = {1.0}$EXTERNAL
reaction55 {1.0}Rn = {1.0}$EXTERNAL
reaction56 {1.0}Rn = {1.0}$EXTERNAL
reaction6 {1.0}Mc = {1.0}$EXTERNAL
reaction7 {1.0}Cc = {1.0}Ccp
reaction8 {1.0}Ccp = {1.0}$EXTERNAL
reaction9 {1.0}$EXTERNAL = {1.0}Mp
end reactions

begin rate equations
reaction0 = cell*function0[reaction0vsb, reaction0K, reaction0m, Rn[t]]
reaction1 = cell*function1[reaction1k, Mb[t]]
reaction10 = cell*function1[reaction10k, Mp[t]]
reaction11 = cell*function2[reaction11V, Pcp[t], reaction11Km]
reaction12 = cell*function2[reaction12V, Pc[t], reaction12Km]
reaction13 = cell*(reaction13k1*Cc[t]*Pc[t] - reaction13k2*PCc[t])
reaction14 = cell*function2[reaction14V, PCc[t], reaction14Km]
reaction15 = cell*function2[reaction15V, PCcp[t], reaction15Km]
reaction16 = cell*(reaction16k1*PCc[t] - reaction16k2*PCn[t])
reaction17 = cell*reaction17k1*PCnp[t]
reaction18 = cell*reaction18k1*Bcp[t]
reaction19 = cell*reaction19k1*Bnp[t]
reaction2 = cell*reaction2k1*Mb[t]
reaction20 = cell*function3[reaction20Vs, Bn[t], reaction20n, reaction20K]
reaction21 = cell*function2[reaction21V, PCn[t], reaction21Km]
reaction22 = cell*reaction22k1*Mp[t]
reaction23 = cell*(-(reaction23k2*In[t]) + reaction23k1*Bn[t]*PCn[t])
reaction24 = cell*function2[reaction24V, Mb[t], reaction24Km]
reaction25 = cell*function2[reaction25V, Mc[t], reaction25Km]
reaction26 = cell*function2[reaction26V, Mp[t], reaction26Km]
reaction27 = cell*reaction27k1*Pc[t]
reaction28 = cell*reaction28k1*Cc[t]
reaction29 = cell*reaction29k1*Pcp[t]
reaction3 = cell*function2[reaction3V, Bc[t], reaction3Km]
reaction30 = cell*reaction30k1*Ccp[t]
reaction31 = cell*reaction31k1*PCcp[t]
reaction32 = cell*reaction32k1*PCc[t]
reaction33 = cell*function2[reaction33V, PCnp[t], reaction33Km]
reaction34 = cell*reaction34k1*Bc[t]
reaction35 = cell*function2[reaction35V, Bcp[t], reaction35Km]
reaction36 = cell*function2[reaction36V, Bn[t], reaction36Km]
reaction37 = cell*function2[reaction37V, Bnp[t], reaction37Km]
reaction38 = cell*reaction38k1*In[t]
reaction39 = cell*function2[reaction39V, In[t], reaction39Km]
reaction4 = cell*(reaction4k1*Bc[t] - reaction4k2*Bn[t])
reaction40 = cell*reaction40k1*Bn[t]
reaction41 = cell*function2[reaction41V, Bcp[t], reaction41Km]
reaction42 = cell*function2[reaction42V, Bnp[t], reaction42Km]
reaction43 = cell*function2[reaction43V, Ccp[t], reaction43Km]
reaction44 = cell*function2[reaction44V, Pcp[t], reaction44Km]
reaction45 = cell*function2[reaction45V, PCnp[t], reaction45Km]
reaction46 = cell*reaction46k1*PCn[t]
reaction47 = cell*function2[reaction47V, PCcp[t], reaction47Km]
reaction48 = cell*function3[reaction48Vs, Bn[t], reaction48n, reaction48K]
reaction49 = cell*reaction49k1*Mr[t]
reaction5 = cell*function1[reaction5k, Mc[t]]
reaction50 = cell*function2[reaction50V, Mr[t], reaction50Km]
reaction51 = cell*function1[reaction51k, Mr[t]]
reaction52 = cell*(reaction52k1*Rc[t] - reaction52k2*Rn[t])
reaction53 = cell*function2[reaction53V, Rc[t], reaction53Km]
reaction54 = cell*reaction54k1*Rc[t]
reaction55 = cell*function2[reaction55V, Rn[t], reaction55Km]
reaction56 = cell*reaction56k1*Rn[t]
reaction6 = cell*reaction6k1*Mc[t]
reaction7 = cell*function2[reaction7V, Cc[t], reaction7Km]
reaction8 = cell*function2[reaction8V, Ccp[t], reaction8Km]
reaction9 = cell*function3[parameter0000097, Bn[t], reaction9n, reaction9K]
end rate equations

begin parameters
reaction0vsb = 1.8
reaction0K = 1.0
reaction0m = 2.0
reaction1k = 0.32
reaction2k1 = 0.02
reaction3V = 1.4
reaction3Km = 1.006
reaction4k1 = 0.8
reaction4k2 = 0.4
reaction5k = 3.2
reaction6k1 = 0.02
reaction7V = 1.2
reaction7Km = 1.006
reaction8V = 1.4
reaction8Km = 0.3
reaction9n = 2.0
reaction9K = 0.6
reaction10k = 1.2
reaction11V = 3.4
reaction11Km = 0.3
reaction12V = 9.6
reaction12Km = 1.006
reaction13k1 = 0.8
reaction13k2 = 0.4
reaction14V = 2.4
reaction14Km = 1.006
reaction15V = 1.4
reaction15Km = 0.3
reaction16k1 = 0.8
reaction16k2 = 0.4
reaction17k1 = 0.02
reaction18k1 = 0.02
reaction19k1 = 0.02
reaction20Vs = 2.2
reaction20n = 2.0
reaction20K = 0.6
reaction21V = 2.4
reaction21Km = 1.006
reaction22k1 = 0.02
reaction23k1 = 1.0
reaction23k2 = 0.2
reaction24V = 1.3
reaction24Km = 0.4
reaction25V = 2.0
reaction25Km = 0.4
reaction26V = 2.2
reaction26Km = 0.3
reaction27k1 = 0.02
reaction28k1 = 0.02
reaction29k1 = 0.02
reaction40k1 = 0.02
reaction41V = 0.2
reaction41Km = 0.1
reaction42V = 0.4
reaction42Km = 0.1
reaction43V = 0.2
reaction43Km = 0.1
reaction44V = 0.6
reaction44Km = 0.1
reaction45V = 0.2
reaction45Km = 0.1
reaction46k1 = 0.02
reaction47V = 0.2
reaction47Km = 0.1
reaction48Vs = 1.6
reaction48n = 2.0
reaction48K = 0.6
reaction49k1 = 0.02
reaction50V = 1.6
reaction50Km = 0.4
reaction51k = 1.7
reaction30k1 = 0.02
reaction31k1 = 0.02
reaction32k1 = 0.02
reaction33V = 1.4
reaction33Km = 0.3
reaction34k1 = 0.02
reaction35V = 3.0
reaction35Km = 0.3
reaction36V = 1.4
reaction36Km = 1.006
reaction37V = 3.0
reaction37Km = 0.3
reaction38k1 = 0.02
reaction39V = 1.6
reaction39Km = 0.3
reaction52k1 = 0.8
reaction52k2 = 0.4
reaction53V = 4.4
reaction53Km = 0.3
reaction54k1 = 0.02
reaction55V = 0.8
reaction55Km = 0.3
reaction56k1 = 0.02
cell = 1.0
end parameters

begin initial conditions
Bc[0] = Bci
Bcp[0] = Bcpi
Bn[0] = Bni
Bnp[0] = Bnpi
Cc[0] = Cci
Ccp[0] = Ccpi
In[0] = Ini
Mb[0] = Mbi
Mc[0] = Mci
Mp[0] = Mpi
Mr[0] = Mri
PCc[0] = PCci
PCcp[0] = PCcpi
PCn[0] = PCni
PCnp[0] = PCnpi
Pc[0] = Pci
Pcp[0] = Pcpi
Rc[0] = Rci
Rn[0] = Rni
end initial conditions

begin initial values
Bci = 0.0
Bcpi = 0.0
Bni = 0.0
Bnpi = 0.0
Cci = 0.0
Ccpi = 0.0
Ini = 0.0
Mbi = 0.0
Mci = 0.0
Mpi = 0.0
Mri = 0.0
PCci = 0.0
PCcpi = 0.0
PCni = 0.0
PCnpi = 0.0
Pci = 0.0
Pcpi = 0.0
Rci = 0.0
Rni = 0.0
end initial values

begin assignment rules
parameter0000097 := 2.4 + 0.6000000000000001*parameter0000096
parameter0000096 := function5[12, t]
end assignment rules

begin function definitions
function2[V_,substrate_,Km_] := (substrate*V)/(Km + substrate)
function3[Vs_,B_,n_,K_] := (B^n*Vs)/(B^n + K^n)
function5[length_,tt_] := Ceiling[Sin[0.001 + (Pi*tt)/length]/2]
function0[vsb_,K_,m_,Bn_] := (K^m*vsb)/(Bn^m + K^m)
function1[k_,mRNA_] := k*mRNA
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units