Edda Klipp

1. Sommersemester 2010 Systembiologie 3 - Stoichiometry Edda Klipp Humboldt-Universität zu Berlin Institut für Biologie Theoretische Biophysik

2. Stoichiometric Analysis of Cellular Reaction Systems F E 2D 2A B + C v1 v3 v2 G http://www.genome.ad.jp/kegg/pathway/map/map01100.html • Analysis of a system of biochemical reactions • Network properties • Enzyme kinetics not considered

3. Example Catalase Stoichiometric coefficients for Hydrogenperoxid, water, oxygen 1 -2 2 Stoichiometric coefficients can be chosen such that they agree with molecularity, but not necessarily. -1 1 1/2 Their signs depend on the chosen reaction direction. Since reactions are usually reversible, one cannot distinguish between „substrate“ and „product“. -2 -1 2 v- v Stoichiometric Coefficients Stoichiometric coefficients denote the proportions, with which the molecules of substrates and products enter the biochemical reactions.

4. Time Course of Concentrations Usually described by ordinary differential equations (ODE) Example catalase for this choice of stoichiometric coefficienten: 1 -2 2

5. - - æ ö 1 1 0 1 S 1 ç ÷ = - N 0 1 1 0 S ç ÷ Row: Substance 2 ç ÷ 0 0 0 1 S è ø 3 4 1 2 3 Column: reaction Time Course of Concentrations Usually described by ordinary differential equations (ODE) Several reactions at the same time  all rate equations must be considered at the same time. 2 1 3 S S 1 2 4 S 3 External metabolites are not included in N. One can summarize the stoichiometric coefficients in matrix N. The rows refer to the substances, the columns refer to the reactions.

6. Balance equations/Systems equations In general: We consider the substances Siand their stoichiometric coefficients nijin the respective reactionj. r – number of reaction Si – metabolite concentration vj – reaction rate nij– stoichiometric coefficient If the biochemical reactions are the only reason for the change of concentration of metabolites, i.e. if there is no mass flow by convection, diffusion or similar Then one can express the temporal behavior of concentrations by the balance equations.

7. Summary Stoichiometric matrix Vector of metabolite concentrations Vector of reaction rates Parameter vector Metabolite concentrations and reaction rates are dependent on kinetic parameters. WithNcan one write systems equations clearly.

8. The Steady State Reaction systems are frequently considered in steady state, Where metabolite concentrations change do not change with time. This describes an implicite dependency of concentrations and fluxes on the parameters. b.z.w. The flux in steady state is

9. v1 v2 v3 d S1 1 -1 0 dt S2 0 1 -1 = dS1 / dt = v1-v2 dS2 / dt = v2-v3 System equations . . = S N v Matrix formalism dSi /dt = 0 Steady state Assumption: Linear kinetics Nv = 0is usually a non-linear equation system, which cannot be solved analytically (necessitates knowledge of kinetic(). Example Unbranched pathway variabel

10. The Stoichiometric Matrix N - Characterizes the network of all reactions in the system - Contains information about possible pathways

11. The Kernel Matrix K In steady state holds Non-trivial solutions exist only if the columns ofNare linearly dependent. Mathematically, the linear dependencies can be expressed by a matrixK with the columnskwhich each solve K – null space (Kernel) ofN The number of basis vectors of the kernel ofNis

12. Calculation of the Kernel Matrix The Kernel matrixK can be calculated with the Gaussian Elimination Algorith for the solution of homogeneous linear equation systems. Example Alternative: calculate with computer programmes Such as „NullSpace[matrix]“ in Mathematica.

13. Representation of Kernel Matrix The Kernel matrixKis not uniquely determined. Every linear combination of columns is also a Possible solution. Matrix multiplication with a regular MatrixQ „from right“ gives another Kernel matrix. For some applications one needs a simple ("kanonical") representation of the Kernel matrix. A possible and appropriate choice is Kcontains many zeros. I – Identity matrix

14. Informations from Kernel Matrix K • Admissible fluxes in steady state • Equilibrium reactions • Unbranched reaction sequences • Elementary modes

15. Admissible Fluxes in Steady State Examples v1 v2 v3 S0 S1 S2 S3 v2 v1 v4 S0 S1 S2 S3 v3 v1 v3 S v2

16. Admissible Fluxes in Steady State With the vectorski (k1, k2,…) is also every linear combination A possible columns ofK. for example: insteadand also In steady state holds All admissible fluxes in steady state can be written as linear combinations of vectorski : v1 v3 für S v2 The coefficientsihave the respective units, eg. or .

17. Equilibrium Reactions Case: all elements of a row inKare 0 Then: the respective reaction is in every steady state in equilibrium. Example

18. Unbranched Reaction Steps v1 v2 v3 S1 S2 S3 S0 v4 S4 The basis vectors of nullspace have the same entries for unbranched reaction sequences. Unbranched reaction sequences can be lumped for further analysis.

19. Kernel Matrix –Dead Ends v1 v2 v3 S1, S2, S3 intern, S0, S4 extern S1 S2 S4 S0 v4 S3 Necessary and sufficient condition for a „Dead end“: One metabolite has only one entry in the stoichiometric matrix (is only once Substrate or product). Flux in steady state through this reaction must vanish in steady state (J4 = 0). Model reduction: one can neglect those reactants for steady state analyses.

20. Kernel Matrix – Irreversibility S1, S2 internal, S0 , S3, S4 external v1 v2 v3 S1 S2 S3 S0 v4 S4 Other choice of basis vectors Mathematically possible, biologically not feasible The basis vectors of a null space are not unique. The direction of fluxes (signs) do not necessarily agree with the direction of irreversible reactions. (Irreversibility limits the space of possible steady state fluxes.)

21. Michaelis-Menten kinetics Isolated reaction: Pyruvatkinase, Na/K-ATPase Conservation Relations: Matrix G If compounds or groups are not added to or deprived of a Reaction system, then must their total amount remain constant. Examples

22. Conservation Relations - Calculation If there exist linear dependencies between the rows of the stoichiometric matrix, then one can find a matrixGsuch as N – stoichiometric matrix holds Due to The integration of this equation yields the conservation relations.

23. Conservation Relations – Properties of G The number of independent row vectorsg (= number of Independent conservation relations) is given by (n = number of rows of the stoichiometric matrix = number of metabolites) GTis theKernel matrixof NT, and can be found in the same way asK. (Gaussian elimination algorithm) The matrixGis not unique, with Pregular quadratic matrix is again conservation matrix. Separated conservation conditions:

24. Conservation of atoms or atom groups, e.g. Pyruvatdecarboxylase (EC 4.1.1.1) CH3CO-group Protons Carboxyl group carbon oxygen hydrogen Elektric charge Conservation Relations – Examples

25. RearrangeN, L – Linkmatrix (independent upper rows, dependent lower rows) RearrangeS respectively (indep upper rows, dep lower rows) For dependent concentrations hold Reduced ODE system Conservation relations – Simplification of the ODE system If conservation relations hold for a reaction system, then the ODE system can be reduced, since some equations are linearly dependent.

26. v2 v1 S1 v3 S2 S4 v4 S3 v5 Basic Elements of Biochemical Networks Systems equations r – number of reactions Si – metabolite concentrations vj – reaction rates nij – stoichiometric coefficients Individual reaction properties Network properties Matrix representation

27. Nachtrag vom 10. Mai 2010

28. Non-negative Flux Vectors In many biologically relevant situations have fluxes fixed signs. We can define their direction such that Example: opposite uni-directional rates instead of net rates, - Description of tracer kinetics or dynamics of NMR labels - different isoenzymes for different directions of reactions - for (quasi) irreversible reactions Sometimes is the value of individual rates fixed. Both conditions restrict the freedom for the choice of Basis vectors for K.

29. Elementary Flux Modes Situation: some fluxes have fixed signes, others can operate in both directions. Which (simple) pathes connect external substrats? v2 v2 P2 P2 v1 v1 P1 P1 S S v3 v3 P3 P3

30. Elementary Flux Modus • An elementary flux mode comprises all reaction steps, • Leading from a substrate S to a product P. • Each of these steps in necessary to maintain a steady state. • The directions of fluxes in elementary modes fulfill • the demands for irreversibility

31. Number of Elementary Flux Modes v1 v2 v3 S1 S2 S3 S0 v4 S4 v1 v2 v3 S1 S2 S3 S0 v4 S4 The number of elementary modes is at least as high as the number of basis vectors of the null space.