Cellerator Examples

This page provides an examples of how a Cellerator Notebook looks, translated into html format, so that Mathematica is not required. It demonstrates several reaction to differential equation translations, a time course prediction, and an input-output relation

<<cellerator.m

Cellerator™ 1.5.0 (03-Jan-2005) loaded 13-Jan-2005 18:23 using Mathematica Version 5.1 for Mac OS X (October 25, 2004)

complexLeft = "$" ; complexRight = "$" ; dash = "$" ;

Creation, Annihilation, and Conversion Reactions

des = interpret[{{∅A, k1}, {AB, k2}, {B∅, k3}, {P⇄Q}, {R + S⇄Q}}]

{{A^′[t] k1 - k2 A[t], B^′[t] k2 A[t] - k3 B[t], P^′[t] &# ... 004 Q[t] - k003 R[t] S[t], S^′[t] k004 Q[t] - k003 R[t] S[t]}, {A, B, P, Q, R, S}}

Rate constants are automatically generated when they are omitted:

TableForm[des[[1]]]

A^′[t] k1 - k2 A[t]
B^′[t] k2 A[t] - k3 B[t]
P^′[t]  -k001 P[t] + k002 Q[t]
Q^′[t] k001 P[t] - k002 Q[t] - k004 Q[t] + k003 R[t] S[t]
R^′[t] k004 Q[t] - k003 R[t] S[t]
S^′[t] k004 Q[t] - k003 R[t] S[t]

Unidirectional Catalytic Reaction with intermediate complex formation

lowLevelReactionTable[{{Overscript[S⇄P, e]}}]

e+S→$S$e$
$S$e$→e+S
$S$e$→e+P
e+P→$S$e$

The following arrow is slightly different:

lowLevelReactionTable[{{Overscript[S⇌P, e]}}]

e+S→$S$e$
$S$e$→e+S
$S$e$→$P$e$
$P$e$→$S$e$
$P$e$→e+P
e+P→$P$e$

interpret[{{Overscript[S⇄P, e], a, d, k}}]

{{e^′[t]  -a e[t] S[t] + d $S$e$[t] + k $S$e$[t], P^′[t] k $S$e$ ... + d $S$e$[t], $S$e$^′[t] a e[t] S[t] - d $S$e$[t] - k $S$e$[t]}, {e, P, S, $S$e$}}

interpret[{{Overscript[S⇌P, e], a, d, k}}]

{{e^′[t]  -a e[t] S[t] + d $S$e$[t], P^′[t] 0, S^′[t] > ... e$[t], $S$e$^′[t] a e[t] S[t] - d $S$e$[t] - k $S$e$[t]}, {e, P, S, $P$e$, $S$e$}}

Bidirectional Catalytic Reaction with intermediate complex formation

complexLeft = "" ; complexRight = "" ; dash = "" ;

lowLevelReactionTable[{{Underoverscript[S⇄P, Er, arg3]}}]

Ef+S→SEf
SEf→Ef+S
SEf→Ef+P
Ef+P→SEf
Er+P→PEr
PEr→Er+P
PEr→Er+S
Er+S→PEr

interpret[{{Underoverscript[S⇄P, Er, arg3], af, df, kf, ar, dr, kr}}]

{{Ef^′[t]  -af Ef[t] S[t] + df SEf[t] + kf SEf[t], Er^′[t]  -ar  ... SEf[t], SEf^′[t] af Ef[t] S[t] - df SEf[t] - kf SEf[t]}, {Ef, Er, P, PEr, S, SEf}}

Catalytic Reaction without intermediate complex formation

lowLevelReactionTable[{{Overscript[AB, C]}}]

A+C→B+C

interpret[{{Overscript[AB, C], k}}]

{{A^′[t]  -k A[t] C[t], B^′[t] k A[t] C[t], C^′[t] 0}, {A, B, C}}

Catalytic Reaction - Michaelis Menten

interpret[{{Overscript[A⟹B, e], MM[v, Kd]}}]

{{A^′[t]  -(Kd A[t] e[t])/(v + A[t]), B^′[t]  (Kd A[t] e[t])/(v + A[t]), e^′[t] 0}, {A, B, e}}

interpret[{{A⟹B, MM[v, Kd]}}]

{{A^′[t]  -(Kd A[t])/(v + A[t]), B^′[t]  (Kd A[t])/(v + A[t])}, {A, B}}

interpret[{{Overscript[A⟹B, e], MM[a, d, k]}}]

{{A^′[t]  -(k A[t] e[t])/((d + k)/a + A[t]), B^′[t]  (k A[t] e[t])/((d + k)/a + A[t]), e^′[t] 0}, {A, B, e}}

interpret[{{Underoverscript[A⟺B, Er, arg3], MM[{vf, Kf}, {vr, Kr}]}}]

{{A^′[t]  -(Kf A[t] Ef[t])/(vf + A[t]) + (Kr B[t] Er[t])/(vr + B[t]), B^′ ...  (Kr B[t] Er[t])/(vr + B[t]), Ef^′[t] 0, Er^′[t] 0}, {A, B, Ef, Er}}

Catalytic Reaction - Hill function

interpret[{{Overscript[A↦B, C], vmax V, nhill n, khalf K}}]

{{A^′[t]  -(V A[t]^n C[t])/(K^n + A[t]^n), B^′[t]  (V A[t]^n C[t])/(K^n + A[t]^n), C^′[t] 0}, {A, B, C}}

interpret[{{Overscript[A↦∅, C], vmax V, nhill n, khalf K}}]

{{A^′[t]  -(V A[t]^n C[t])/(K^n + A[t]^n), C^′[t] 0}, {A, C}}

Annihilation by Hill Function

interpret[{{Overscript[A↦∅, ∅], vmax V, nhill n, khalf K}}]

{{A^′[t]  -(V A[t]^n)/(K^n + A[t]^n)}, {A}}

S-System

interpret[{ {va↦v1, SSystem[taua1, kp, km, ca1p, ca1m]}, {vb↦v1, SSystem ... b1, kp3, km3, cp3, cm3]},  {vb↦v2, SSystem[tau2, kp2, km2, cb2p, cb2m]} }]

{{v1^′[t]  (-km va[t]^ca1m vb[t]^cb1m vc[t]^cm3 + kp va[t]^ca1p vb[t]^cb1p vc[ ... a^′[t] 0, vb^′[t] 0, vc^′[t] 0}, {v1, v2, va, vb, vc}}

Diffusion

interpret[{{AB + C, k}, {A▽, dA}, {B▽, dB * DiagonalMatrix[{dB1, dB2, dB3}]}, {C▽, Array[dC, {3, 3}]} }]

{{A^(1, 0, 0, 0)[t, x, y, z]  -k A[t, x, y, z] + dA A^(0, 0, 0, 2)[t, x, y, z] + dA  ... , y, z] + dC[2, 1] C^(0, 1, 1, 0)[t, x, y, z] + dC[1, 1] C^(0, 2, 0, 0)[t, x, y, z]}, {A, B, C}}

Cascades

no catalyst

lowLevelReactionTable[{{PQRS}}]

P→Q
Q→R
R→S

single catalyst

lowLevelReactionTable[{{Overscript[P⇄Q⇄R⇄T, S]}}]

P+S→PS
PS→P+S
PS→Q+S
Q+S→PS
Q+S→QS
QS→Q+S
QS→R+S
R+S→QS
R+S→RS
RS→R+S
RS→S+T
S+T→RS

different catalysts

lowLevelReactionTable[{{Overscript[P⇄Q⇄R, {S1, S2}]}}]

P+S1→PS1
PS1→P+S1
PS1→Q+S1
Q+S1→PS1
Q+S2→QS2
QS2→Q+S2
QS2→R+S2
R+S2→QS2

s = interpret[{{Underoverscript[P⇄Q⇄R, revCat, arg3]}}] ;

TableForm[First[s]]

forCat^′[t]  -k005 forCat[t] P[t] + k006 PforCat[t] + k007 PforCat[t] - k005 forCat[t] Q[t] + k006 QforCat[t] + k007 QforCat[t]
P^′[t]  -k005 forCat[t] P[t] + k006 PforCat[t] + k010 QrevCat[t]
PforCat^′[t] k005 forCat[t] P[t] - k006 PforCat[t] - k007 PforCat[t]
Q^′[t] k007 PforCat[t] - k005 forCat[t] Q[t] + k006 QforCat[t] + k009 QrevCat[t] - k008 Q[t] revCat[t] + k010 RrevCat[t]
QforCat^′[t] k005 forCat[t] Q[t] - k006 QforCat[t] - k007 QforCat[t]
QrevCat^′[t]  -k009 QrevCat[t] - k010 QrevCat[t] + k008 Q[t] revCat[t]
R^′[t] k007 QforCat[t] - k008 R[t] revCat[t] + k009 RrevCat[t]
revCat^′[t] k009 QrevCat[t] + k010 QrevCat[t] - k008 Q[t] revCat[t] - k008 R[t] revCat[t] + k009 RrevCat[t] + k010 RrevCat[t]
RrevCat^′[t] k008 R[t] revCat[t] - k009 RrevCat[t] - k010 RrevCat[t]

Transcriptional and Hill Reactions

interpret[{{A1↦B, type hill, vmax V1, khalf K, nhill n}}]

{{A1^′[t] 0, B^′[t]  (V1 A1[t]^n)/(K^n + A1[t]^n)}, {A1, B}}

cooperative activation with Hill model

TableForm[First[interpret[{ {A1↦B, type  hill, nhill n1, khalf ... B↦A1, type hill, vmax V1, khalf K, nhill n4} }] ]]

A1^′[t]  (V1 B[t]^n4)/(K^n4 + B[t]^n4)
A2^′[t] 0
A3^′[t] 0
B^′[t]  (v1 A1[t] + v2 A2[t] + v3 A3[t])/(1 + v1 A1[t] + v2 A2[t] + v3 A3[t])

cooperative activation with GRN model

interpret[{ {A1↦B, type GRN, RGRN R, TGRN C1, hGRN> ... 2754; GRN, TGRN C2},  {A3↦B, type GRN, TGRN C3} }]

{{A1^′[t] 0, A2^′[t] 0, A3^′[t] 0, B^′[t] 1/(1 + ^(-1 - C1 A1[t] - C2 A2[t] - C3 A3[t]))}, {A1, A2, A3, B}}

cooperative activation with primary  NHCA model

TableForm[First[interpret[{ {A1↦B, type  NHCA, mNHCA m1, nNHCA ... B↦A1, type hill, vmax V1, khalf K, nhill n4} }] ]]

A1^′[t]  (V1 B[t]^n4)/(K^n4 + B[t]^n4)
A2^′[t] 0
A3^′[t] 0
B^′[t]  ((1 + TP1 A1[t]^n1)^m1 (1 + TP2 A2[t]^n2)^m1 (1 + TP3 A3[t]^n3)^m1)/(k ... 2)^m1 (1 + TM3 A3[t]^n3)^m1 + (1 + TP1 A1[t]^n1)^m1 (1 + TP2 A2[t]^n2)^m1 (1 + TP3 A3[t]^n3)^m1)

Cooperative Activation with Alternate NHCA model

interpret[{{〈A1, A2, A3〉↦B, type NHCA, TPLUS TP, TMINUS TM, kNHCA k, nNHCA n, mNHCA m}}]

{{A1^′[t] 0, A2^′[t] 0, A3^′[t] 0, B^′[t] &# ...  + TM A2[t]^n + TM A3[t]^n)^m + (1 + TP A1[t]^n + TP A2[t]^n + TP A3[t]^n)^m)}, {A1, A2, A3, B}}


Created by Mathematica  (January 13, 2005) Valid XHTML 1.1!