Stöchiometrische Matrix
Die stöchiometrische Matrix ist eine Matrix, welche die Stöchiometrie eines Reaktionsnetzwerkes in kompakter Form repräsentiert.
Sie wird meist mit <math>N</math> abgekürzt. In der Regel entsprechen die Spalten von <math>N</math> den Reaktionen des Systems während die Zeilen den chemischen Spezies entsprechen. Spezies einer Reaktion, welche in Summe konsumiert werden, erhalten einen Eintrag mit negativen Vorzeichen und Spezies, welche in Summe produziert werden, erhalten einen positiven Eintrag an der Position von <math>N</math>, welche der gegebenen Reaktion und Spezies entspricht. Die Änderung der Spezies nach der Zeit ist dann gegeben durch <math>\text{d}c/\text{d}t = Nv</math> wobei <math>v</math> den Vektor der Reaktionsraten (auch „Flussvektor“ genannt) repräsentiert. In einem stationären Zustand gilt demnach <math>Nv = \mathbf{0}</math>, d. h. das System befindet sich in einem Fließgleichgewicht.
Die stöchiometrische Matrix lässt Rückschlüsse auf die Reaktionsraten der stationären Zustände zu. Im stationären Zustand muss dann mindestens <math>v\in\text{kern}(N)</math> gelten, d. h. die Menge aller solcher <math>v</math> liegen im Nullraum von <math>N</math>. Diese Bedingung ist unabhängig von der Kinetik, welcher das gegebene chemische System zugrunde liegt.
Der Rang der stöchiometrischen Matrix gibt die Zahl linear unabhängiger Reaktionen an.
Beispiele
Die vier Reaktionen
<math> \begin{array}{crcl} R_1 := & \text{A} & \rightarrow & \text{B},\\ R_2 := & 4 \text{C} + 2 \text{B} & \rightarrow & 5 \text{C} + \text{B},\\ R_3 := & \text{C} & \rightarrow & \text{D},\\ R_4 := & \text{D} & \rightarrow & \text{A}, \end{array} </math>
lassen sich wie folgt als Matrix kodieren:
<math> N = \left[\begin{array}{rrrr} -1 & 0 & 0 & 1\\
1 &-1 & 0 & 0\\ 0 & 1 &-1 & 0\\ 0 & 0 & 1 &-1
\end{array}\right], </math>
wobei die Spalten der Reihe nach den Reaktionen <math>R_1, R_2, R_3, R_4</math> und die Zeilen den Spezies <math>\text{A}, \text{B}, \text{C}, \text{D}</math> entsprechen. Reaktion <math>R_1</math> konsumiert eine Einheit von Spezies <math>\text{A}</math> und produziert eine Einheit von Spezies <math>\text{B}</math> (erste Spalte). Gleichfalls wird in Reaktion <math>R_2</math> netto eine Einheit von Spezies <math>\text{C}</math> produziert während netto eine Einheit von Spezies <math>\text{B}</math> konsumiert wird (zweite Spalte). Ist im Beispiel <math>v = [1,1,1,1]^T</math>, dann ergibt sich <math>Nv=\mathbf{0}</math> und das System befindet sich in einem stationären Zustand.
Wie man an diesem Beispiel sieht kann das originale chemische System nicht allein durch Kenntnis der stöchiometrischen Matrix rekonstruiert werden. Für Reaktion <math>R_2</math> existieren unendliche viele Möglichkeiten, welche die gleiche Spalte in <math>N</math> erzeugen:
<math> \begin{array}{rcl} 3 \text{C} + \text{B} & \rightarrow & 4 \text{C},\\ 4 \text{C} + 2 \text{B} & \rightarrow & 5 \text{C} + \text{B},\\ 5 \text{C} + 3 \text{B} & \rightarrow & 6 \text{C} + 2 \text{B},\\
& \vdots &
\end{array} </math>
Demnach enthält die stöchiometrische Matrix weniger Information als der ursprüngliche Satz an Reaktionen.
Reduzierung einer stöchiometrischen Matrix auf die sogenannten "unabhängigen" Reaktionen
Beispielsweise sollen 6 bekannte Reaktionen simultan ablaufen<ref>Autorenkollektiv (K. Budde federführend): Reaktionstechnik III (Reihe Verfahrenstechnik), VEB Verlag für Grundstoffindustrie, Leipzig 1977, Beispiel 3 S. 23, Berechnung eines komplexen Reaktionssystems nicht-volumenkonstanter Gasreaktionen mit der stöchiometrischen Matrix und den Umsatzvariablen der unabhängigen Reaktionen</ref>:
1) <math> 1 \mathrm A + 2\mathrm B \; \longrightarrow \; 1\mathrm C </math>
2) <math> 1 \mathrm C \; \longrightarrow \; 1\mathrm A + 2\mathrm B </math>
3) <math> 2 \mathrm A + 2\mathrm B \; \longrightarrow \; 1\mathrm D +1\mathrm E </math>
4) <math> 1 \mathrm D + 1\mathrm E \; \longrightarrow \; 2\mathrm A +2\mathrm B </math>
5) <math> 1 \mathrm A + 3\mathrm B\; \longrightarrow \; 1\mathrm D + 1\mathrm F</math>
6) <math> 1 \mathrm D + 1\mathrm F \; \longrightarrow \; 1\mathrm A + 3\mathrm B </math>
Wie man an den obigen chemischen Reaktionsgleichungen sieht, sind jeweils zwei Gleichungen als Gleichgewichtsreaktion zusammengehörig:
- 1) und 2)
- 3) und 4)
- 5) und 6)
Sie können daher nicht alle voneinander "unabhängig" sein. Eine Reduzierung der Anzahl der Gleichungen ist daher sinnvoll und möglich. Welche der genannten Reaktionsgleichungen letztendlich als unabhängige verbleiben ist faktisch eine zufällige oder willkürliche Entscheidung, die vom Ergebnis des Gaußschen Verfahrens abhängt.
Nimmt ein Stoff/Reaktand i nicht an der Reaktion j teil, so ist sein stöchiometrischer Koeffizient für diese Reaktion j in der Koeffizientenmatrix mit Null anzusetzen. Dies gilt generell für Inertstoffe, falls diese mit bilanziert werden müssen.
| Reaktand i | ||||||
|---|---|---|---|---|---|---|
| Reaktion j | A | B | C | D | E | F |
| 1) | −1 | −2 | +1 | 0 | 0 | 0 |
| 2) | +1 | +2 | −1 | 0 | 0 | 0 |
| 3) | −2 | −2 | 0 | +1 | +1 | 0 |
| 4) | +2 | +2 | 0 | −1 | −1 | 0 |
| 5) | −1 | −3 | 0 | +1 | 0 | +1 |
| 6) | +1 | +3 | 0 | −1 | 0 | −1 |
Nun erfolgt die lineare Vereinfachung mit dem Gaußschen Algorithmus. Im Beispiel durch Addition von Zeilen (Reaktionen bzw. von deren stöchiometrischen Koeffizienten).
| Reaktand i | |||||||
|---|---|---|---|---|---|---|---|
| Reaktion j | A | B | C | D | E | F | Hinweise |
| 1) | −1 | −2 | +1 | 0 | 0 | 0 | bleibt (unabhängige Reaktion) |
| 2+1) | 0 | 0 | 0 | 0 | 0 | 0 | fällt weg (Reaktion 2) |
| 3) | −2 | −2 | 0 | +1 | +1 | 0 | bleibt (unabhängige Reaktion) |
| 4+3) | 0 | 0 | 0 | 0 | 0 | 0 | fällt weg (Reaktion 4) |
| 5) | −1 | −3 | 0 | +1 | 0 | +1 | bleibt (unabhängige Reaktion) |
| 6+5) | 0 | 0 | 0 | 0 | 0 | 0 | fällt weg (Reaktion 6) |
Ergebnis der Bearbeitung:
| Reaktand i | ||||||
|---|---|---|---|---|---|---|
| Reaktion j | A | B | C | D | E | F |
| 1) | −1 | −2 | +1 | 0 | 0 | 0 |
| 3) | −2 | −2 | 0 | +1 | +1 | 0 |
| 5) | −1 | −3 | 0 | +1 | 0 | +1 |
Es verbleiben also im Beispiel nur drei Reaktionen als linear unabhängige Reaktionen. Der Rang der Koeffizientenmatrix ist daher r=3.
Weitere Eigenschaften
Seien <math>\mathbb{P} := \{x\in\mathbb{R}\ \vert\ x>0\}</math> die Menge aller reellen Zahlen größer Null und <math>\overline{\mathbb{P}} := \mathbb{P}\cup\{0\}</math> die Menge aller reellen Zahlen größer oder gleich Null<ref name=Gunawardena2003>J. Gunawardena: Chemical reaction network theory for in-silico biologists. 2003</ref>. Weiterhin sei <math>r</math> die Anzahl der Reaktionen im gegebenen chemischen System. Die folgenden Mengen an Vektoren <math>v</math> repräsentieren stationäre Zustände chemischer Systeme mit unterschiedlichen Randbedingungen<ref name=Gagneur2004>J. Gagneur and S. Klamt: Computation of elementary modes: a unifying framework and the new binary approach. BMC Bioinformatics 2004</ref>:
- nur reversible Reaktionen, d. h. Reaktionen, welche in Vorwärts- und Rückwärtsrichtung ablaufen können: <math>\text{ker}(N)</math>;
- nur irreversible Reaktionen, d. h. Reaktionen, welche nur in Vorwärtsrichtung ablaufen können: <math>\text{ker}(N)\cap\overline{\mathbb{P}}^r</math>;
- eine Mischung aus reversiblen und irreversiblen Reaktionen: <math>\{v\in\mathbb{R}^r\text{ }\vert\text{ }Nv=\mathbf{0}\text{ and }v_i \geq 0 \text{ } \forall i\in Irr\}</math>, wobei <math>Irr</math> die Menge der Indices der irreversiblen Reaktionen bezeichnet.
Anwendung
Die stöchiometrische Matrix ist ein zentrales Werkzeug der Systembiologie. Sie ermöglicht eine systematische Analyse der Flussvektoren von stationären Zuständen eines chemischen oder biologischen Systems. Im Allgemeinen ist hierbei der Größe des zu analysierenden Systems fast keine Grenze gesetzt, da ihre Verwendung nur Methoden aus der linearen Algebra erfordert. Methoden, welche auf die stöchiometrische Matrix zurückgreifen, sind z. B.: FBA (Flux Balance Analysis)<ref name=Orth2010>J. D. Orth and I. Thiele and B. O. Palsson: What is flux balance analysis? Nat Biotechnol 2010</ref>, FCA (Flux Coupling Analysis)<ref name=Burgard2004>A. P. Burgard and E. V. Nikolaev and C. H. Schilling and C. D. Maranas: Flux coupling analysis of genome-scale metabolic network reconstructions. Genome Res 2004</ref>, FVA (Flux Variability Analysis)<ref name=Mahadevan2003>R. Mahadevan and C. Schilling: The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metabolic engineering 2003</ref>, das Konzept der EFMs (Elementary Flux Modes)<ref name=Schuster2000>S. Schuster and D. A. Fell and T. Dandekar: A general definition of metabolic pathways useful for systematic organization and analysis of complex metabolic networks. Nat Biotechnol 2000</ref><ref name=Larhlimi2008>A. Larhlimi and A. Bockmayr: On inner and outer descriptions of the steady-state flux cone of a metabolic network. Computational Methods in Systems Biology 2008</ref> und ähnliche Methoden wie Extreme Currents<ref name=Larhlimi2008></ref><ref name=Clarke1988>B. L. Clarke: Stoichiometric network analysis. Cell Biochem Biophys 1988</ref> und Extreme Pathways<ref name=Larhlimi2008></ref><ref name=Schilling2000>C. H. Schilling and D. Letscher and B. O. Palsson: Theory for the systemic definition of metabolic pathways and their use in interpreting metabolic function from a pathway-oriented perspective. J Theor Biol 2000</ref>, DFBA (Dynamic FBA)<ref name=Mahadevan2002>R. Mahadevan and J. S. Edwards and F. J. Doyle: Dynamic flux balance analysis of diauxic growth in Escherichia coli. Biophys J 2002</ref> und CRNT (Chemical Reaction Network Theory)<ref name=Gunawardena2003></ref><ref name=Horn1972>F. Horn and R. Jackson: General mass action kinetics. Arch Rational Mech Anal 1972</ref><ref name=Feinberg1979>M. Feinberg: Lectures on chemical reaction networks. 1979</ref><ref name=Feinberg1995>M. Feinberg: The existence and uniqueness of steady states for a class of chemical reaction networks. Arch Rational Mech Anal 1995</ref>.
In der Verfahrenstechnik werden komplexe Reaktionssysteme, in denen also mehrere "unabhängige" chemische Reaktionen gleichzeitig im Reaktionsvolumen ablaufen, mithilfe der Stöchiometrischen Matrix der unabhängigen Reaktionen bilanziert. Dazu benutzt man Umsatzvariablen und modifizierte Umsatzvariablen dieser Reaktionen.
Literatur
- Autorenkollektiv (federführend K. Budde): Reaktionstechnik III (Reihe Verfahrenstechnik), VEB Deutscher Verlag für Grundstoffindustrie, Leipzig 1977, S. 21–27 (Bilanzierungsbeispiel Nr. 3 einer komplexen chemischen Gasreaktion. Aufstellung und Vereinfachung der Koeffizienten-Matrix. Bestimmung des Ranges der Matrix (Anzahl linear unabhängiger Reaktionen). Vereinfachung der Matrix mit dem Gaußschen Algorithmus. Bilanzierung der nicht-volumenkonstanten Gasreaktionen als Stoffmengenbilanzen mittels normaler Umsatzvariable für einen kontinuierlichen Prozess.); sowie Budde: Reaktionstechnik I, Leipzig 1974, Kapitel "Stöchiometrische Grundlagen komplexer Reaktionen" S. 35
- Anwendung des Gibbsschen Phasengesetzes auf chemische Reaktionen, S. 105–109 in Karl-Heinz Näser, Dieter Lempe, Otfried Regen: Physikalische Chemie für Techniker und Ingenieure, VEB Deutscher Verlag für Grundstoffindustrie, Leipzig 1990, 19. unveränderte Auflage, ISBN 3-342-00545-9
Siehe auch
- Freiheitsgrad (die Anzahl der unabhängigen Reaktionsgleichungen hängt von der Anzahl der Reaktanden und Phasen etc. ab)
- Gibbssches Phasengesetz
Einzelnachweise
<references/>
Weblinks
- Stöchiometrische Analyse zellulärer Reaktionssysteme (PDF; 168 kB)