Title: Basics of biochemical systems theory
1Basics of biochemicalsystems theory
- I. Graphical representation of biochemical
networks
2Material flow
3Information flow
FDP
-
HK
4Networks
Glc
G6P
F6P
FDP
-
-
PGI
ATP
ADP
ATP
ADP
PFK
HK
5Fundamental types of variables
- State variables any quantities that may directly
influence a process and whose values might change
(eg. concentrations, temperature) - Flux variables reaction rates
v2
Glc
G6P
F6P
FDP
-
-
v4
v3
v1
PGI
ATP
ADP
ATP
ADP
PFK
HK
v5
6Other types of variables
Independent variables
fixed or isolated variables
Glc
G6P
F6P
FDP
-
-
PGI
ATP
ADP
ATP
ADP
PFK
HK
7Other types of variables
Dependent variables
Dynamic variables
Glc
G6P
F6P
FDP
-
-
PGI
ATP
ADP
ATP
ADP
PFK
HK
8Other types of variables
Aggregated variables
Linear combinations variables
Glc
G6P
F6P
FDP
-
-
PGI
ATP
ADP
ATP
ADP
PFK
HK
ATP ADP ANPtot
v2f - v2r v2
9Other types of variables
Constrained variables
Variables whose values depends on the values of
other variables not through the interactions of
the various processes but through imposed
algebraic relationships (moiety conservation,
equilibria)
Glc
G6P
F6P
FDP
-
-
PGI
ATP
ADP
ATP
ADP
PFK
HK
ATP ANPtot ADP AMP
F6P Keq G6P
10Other types of variables
Implicit variables
Variables that, in a given context, can be
assumed constant or not influencing subsystem
behavior, and are thus neglected
X
Glc
G6P
F6P
FDP
-
-
PGI
ATP
ADP
ATP
ADP
PFK
HK
Y
Z
11Symbology
State variables Xi - Dependent variables -
Independent variables - Aggregated
variables Flux variables vij
v12
X6
X1
X2
X3
X7
v21
v37
-
-
v23
v61
X9
X5
X4
X5
X4
X10
X8
v45
X5 X4 X11
12Strategy for developinga graphical representation
- Sketch graphic in familiar terms
- Potential variables
- Important processes
- Principal interactions
- Make conversion table
- Redraw graph in symbolic terms
- Analyze and refine model
13Common errors and ambiguities
Failure to explicitly account forconsumption/dilu
tion of components
Non-obvious point expansion fluxes dilution
of components owing to volume expansion (e.g.
cell growth).
14Common errors and ambiguities
Confusing convergent processeswith a
multi-substrate reaction
15Common errors and ambiguities
Failure to account for the actual molecularity of
each species in a given reaction
16Common errors and ambiguities
Confusion between material and information flow
-
17System and environment
X7
X4
X5
-
-
X3
X1
X2
X6
18(No Transcript)
19Realism and accuracy vs. feasibility
We should make things as simple as possible, but
not simpler. Albert Einstein
Any intelligent fool can make things bigger,
more complex, and more violent. It takes a touch
of genius and a lot of courage to move in the
opposite direction. Albert Einstein
20Sources of simplification
- Spatial or topological features
(compartmentalization, channeling) - Time-scale separation (evolutionary,
developmental, biochemical, biomolecular) and
quasi-steady-state approximation - Functional simplifications (operating ranges and
linearizations, non-linearity)
21Basics of biochemicalsystems theory
- II. Mathematical representation of processes
22Rate laws
Mathematical functions representing the
instantaneous rate of a process as an explicit
function of all the state variables that have a
direct influence on the process
Mass action rate laws in a well-mixed
homogeneous environment, the rate of an
elementary reaction a1X1 a2X2 anXn ?
is given by
Kinetic orders
Rate constant
Mass action rate laws may also approximate well
the kinetics ofsome more-complex reaction
mechanisms
23Rate laws
Mathematical functions representing the
instantaneous rate of a process as an explicit
function of all the state variables that have a
direct influence on the process
Michaelis-Menten rate laws initial rates of
enzyme-catalyzed reactions in absence of
cooperative effects.
Maximal rate
24Rate laws
Mathematical functions representing the
instantaneous rate of a process as an explicit
function of all the state variables that have a
direct influence on the process
Michaelis-Menten rate laws initial rates of
enzyme-catalyzed reactions in absence of
cooperative effects.
25Rate laws
Mathematical functions representing the
instantaneous rate of a process as an explicit
function of all the state variables that have a
direct influence on the process
Hill rate laws approximation to initial rates of
enzyme-catalyzed reactions with cooperativity.
Maximal rate
Hill number
Km Half-saturation constant
26Rate laws
Mathematical functions representing the
instantaneous rate of a process as an explicit
function of all the state variables that have a
direct influence on the process
Many other rate laws possible
Could we find a convenient (even if only
approximated) generic representation?
27Taylor series
Linear approximation
28First-order (linear) Taylor series approximation
for multiple variables
29Advantages and limitationsof linearized
representation
- Good enough approximation in many cases (small
operating ranges, linear by design) - Strong theory and methods available
- However
- Biochemical processes are strongly nonlinear
30Power-law representation
Taylor series in logarithmic space
First-order truncation
31Power-law representation
Taylor series in logarithmic space First-order
truncation for n variables
32Power-law representation
Meaning of the kinetic orders
gi
With all other conc. constant
gigt(lt)0 ? v increases (decreases) with Xi gi x
? a 1 change in Xi causes a x change in v
33Ranges of values of kinetic orders
Kinetic orders for mass action kinetics
v
Log-log plot
X
v k Xc ? g c
34Ranges of values of kinetic orders
Kinetic orders for substrates under
Michaelis-Menten kinetics
v/Vmax
X/KM
0 g 1
35Ranges of values of kinetic orders
Kinetic orders for inhibitors under
Michaelis-Menten kinetics
v/v(X0)
g-0.5
X/KI
-1 g 0
36Ranges of values of kinetic orders
Kinetic orders for substrates under cooperative
kinetics
v/Vmax
ggt1
X/KM
0 g (4)
37Parameter estimation
- From general considerations
- From experimentally characterized rate expression
- Top-down modeling
38Basics of biochemicalsystems theory
- III. Mathematical representation of networks
39Mass balance
Chemical processes do not destroy or create
matter
Instantaneous rate of X1 accumulation
Instantaneous rate of X4 conversion into X1
(Instantaneous rate of X1 conversion into X2
Instantaneous rate of X1 conversion into X3)
40Steady state
No net accumulation or dilution
O Instantaneous rate of X4 conversion into X1
(Instantaneous rate of X1 conversion into X2
Instantaneous rate of X1 conversion into X3)
41Mass balance equationsfor a network
v51
-
X5
X1
v12
v14
v40
v13
X4
v23
X2
X3
v20
42Mass balance equationsfor a network
General form
vijk unidirectional rate of utilization of Xi
for the production of Xj via the kth
parallel process ?jik number of molecules of Xj
created in each reaction ?jik number of
molecules of Xi consumed in each reaction
43Flux aggregation
Major types
Emphasis
Fluxes
Pools
Strategies
Parallel
Antiparallel
44Flux aggregation
Power-law representation under various types of
aggregation
Emphasis
Fluxes
Pools
Strategies
Parallel
Antiparallel
45Generalized Mass Action systems
Power-law representation on the basis of flux
aggregation emphasizing processes
More simply
Maintains separate identity of fluxes at branch
points
46Flux aggregation
Power-law representation under various types of
aggregation
Emphasis
Fluxes
Pools
Strategies
Parallel
Antiparallel
47Synergistic systems
Power-law representation on the basis of flux
aggregation emphasizing pools
X1, ..., Xn Dependent variables
Xn1, ..., Xnm Independent variables
Kinetic orders
Rate constants
48Some advantages of S systems
- More compact representation
- Fewer parameters
- More mathematically tractable
- Frequently more accurate
49Ranges of values of kinetic orders
Relationship between kinetic orders for the S
systems representation and those for the GMA
representation
V1 ? v211 2v212 v213
v211
X2
v212
X1
v213
X3
For unit enzyme-catalyzed steps, gij is normaly
-4?gij ?4
50Setting up an S system
1. For each dynamic metabolite, aggregate the
production fluxes and the removal fluxes
separately.
51Setting up an S system
2. Identify the variables that influence each
aggregated flux.
52Setting up an S system
3. Represent each aggregated flux as a
power-product of its influencing variables.
53How we got to love the S-system...
54Log gains
Log gains for concentrations
Xi (i1, ..., n) not fixed in this case. The log
gains are systemic properties.
55Log gains
Log gains for aggregated fluxes
Log gains for F factors (turnover numbers,
)
56IX. The distinction between kinetic orders and
log gains of fluxes
57Sensitivities
How parameters influence state variables
For an extensive list of relationships see M. A.
Savageau and A. Sorribas. J.Theor.Biol. 141
(1)93-115, 1989.
58Example
59Stability
Lyapunov stability is derived from the
linearization of the ODEs
For S systems
The analysis proceeds through evaluation of the
Lyapunov exponents as usual
60Recasting
61Further reading
- M. A. Savageau. Biochemical Systems Analysis A
Study of Function and Design in Molecular
Biology, Reading, Mass.Addison-Wesley, 1976. - E. O. Voit. Computational analysis of biochemical
systems. A practical guide for biochemists and
molecular biologists, CambridgeCambridge
University Press, 2000. - M. A. Savageau. Design principles for elementary
gene circuits Elements, methods, and examples.
Chaos 11 (1)142-159, 2001.
62(No Transcript)