Parameteridentifikation
(Bioprozess-Modell)

Wir unterstützen Sie bei der präzisen Parameter-Identifikation am Beispiel eines vereinfachten Upstream-Prozesses zur Herstellung rekombinanten Bioprodukte. Unser Ziel ist es, das wesentliche Systemverhalten Ihres genetisch veränderten Organismus exakt in einem mathematischen Bioprozess-Modell abzubilden. Dieses maßgeschneiderte Modell liefert Ihnen im Anschluss die perfekte Grundlage für Ihre strategischen Entscheidungen: Sie können damit tiefgehende Prozessanalysen wie Sensitivitätsanalysen durchführen oder Ihren bestehenden Bioprozess gezielt optimieren, um beispielsweise Ihre Ausbeute spürbar zu maximieren.

(Stand: Juni 2022)

Anwendungsbeispiel

Die Kultivierung des genetisch veränderten Organismus $\color{#4daf4a}{X}$ (vgl. E. coli) erfolgt in einem ideal durchmischten klassischen Rührkessel-Reaktor. Die Bildung des rekombinanten Bioproduktes $\color{#377eb8}{P}$ wird durch die Zugabe eine Aktivators $a$ (vgl. IPTG) zum Zeit­punkt $t_{\text{a}}$ induziert [1]. Hierbei konkurrieren Bio­mas­se­wachs­tum und Pro­dukt­bil­dung je­weils um dasselbe Substrat $\color{#ff7f00}{S}$ (vgl. Glu­ko­se). Bei Ablauf der Pro­zess­lauf­zeit $t_{\text{harvest}}$ erfolgt die Ernte und Aufarbeitung des gebildeten Bioproduktes. Desweiteren sei bekannt, dass hohe Substrat-Konzentrationen $\color{#ff7f00}{S}$ das Biomassewachstum hemmen. Nach Zugabe des Aktivators wird ein Großteil des Substrates durch Produktbildung verbraucht. Zusätzlich sei eine deutliche Inhibition des Biomassewachstums durch die ansteigende Produktkonzentration gegeben. Nach Auszehrung des Substrates stagniert die Produktbildung. Ebenfalls kann der Zerfall von Biomasse beobachtet werden.

Voraussetzungen zur Durchführung der Parameter-Identifikation sind eine ausreichende experimentelle Datenbasis sowie ein entsprechendes mathematisches Bioprozess-Modell. Die Datenbasis muss das betrachtete System mit mindestens allen essentiellen Eigenschaften abbilden. Hierzu sind eine Vielzahl von Experimenten bzw. Versuchen durchzuführen, um Messdaten für die Parameter-Identifikation zu gewinnen. In den Experimenten werden festgelegte Steuergrößen in geeigneter Kombination variiert. Steuergrößen sind jene Größen, welche einen Eingriff in das jeweilige System erlauben. Im betrachteten Beispiel werden die Start-Biomassekonzentration $X_0$, die Start-Substratkonzentration $S_0$, der Aktivierungszeitpunkt der Produktbildung $t_{\text{a}}$ sowie der Ernte-Zeitpunkt $t_{\text{harvest}}$ als Steuergrößen festgelegt.

Mittels Versuchsplanung inklusive der Aussteuerung von Steuergrößen können relevante Experimente geplant und durchgeführt werden. Im betrachteten Anwendungsbeispiel sei die Datenbasis durch 9 Experimente gegeben, welche in nachfolgender Abbildung visualisiert sind (zum Erweitern anklicken):

Desweiteren sei das notwendige mathematische Bioprozess-Modell bekannt und in Form des nachfolgenden Differentialgleichungssystems gegeben [2,3]:

Die Parameter $k_{\text{d}}$ und $m$ entsprechen hierbei den konstanten Raten des Zerfalls von Biomasse sowie dem Verbrauch von Substrat zur Bedienung des Erhaltungsstoffwechsels der Zelle. Die beiden Ertragskoeffizienten $Y_{\text{X,S}}$ und $Y_{\text{P,S}}$ spiegeln das Verhältnis zwischen verbrauchtem Substrat zu hieraus gebildeter Biomasse bzw. gebildetem Produkt wider. Die nachfolgenden Gleichungen (4) und (5) bilden die Modell-Kinetiken für die Beschreibung von Biomasse-Wachstum bzw. -Zerfall und Produktbildung ab [2]:

Methodischer Hintergrund

Die Parameter-Identifikation, also die Bestimmung der Beträge aller Modell-Parameter, wird als mathematische Minimierung implementiert. Ziel ist es, diese Beträge so zu bestimmen, dass die Abweichung (Fehler) zwischen experimentellen Daten und simulierten Verläufen minimal wird. Hierzu stehen eine Vielzahl von etablierten Verfahren bereit. Im vorliegenden Anwendungsbeispiel wird die Methode der kleinsten Fehlerquadrate (engl.: least squares method) zusammen mit einem direkten heuristischen Suchverfahren (downhill-simplex [4]) eingesetzt. Desweiteren sei gegeben, dass alle Modell-Parameter strukturell sowie praktisch Identifizierbar sind [5].

Das folgende Schema visualisiert den Ablauf der Parameter-Identifikation und ihre wesentlichen Komponenten, deren Zusammenspiel danach beschrieben wird:

Zur Durchführung der Parameter-Identifikation müssen Wertebereiche und Start-Werte $\theta_0$ für die (vektorisierten) Modell-Parameter $\theta$ vorgegeben werden. In den meisten Anwendungsfällen können die Begrenzungen (engl.: constraints oder bounds) aufgrund von Vorwissen sowie formal-kinetischen bzw. mathematischen Gegebenheiten eingeschränkt werden. Nachfolgende Tabelle listet die in diesem Anwendungsbeispiel verwendeten Modell-Parameter sowie deren entsprechende Wertebereiche (Begrenzungen) auf:

Die Festlegung geeigneter Start-Werte stellt im allgemeinen eine große Herausforderung dar. Diese Start-Werte sind im Idealfall durch bestehendes Vorwissen bereits gegeben, oder können bei einfachen Modellen, mit sehr wenigen Parametern, empirisch ermittelt werden. Handelt es sich, wie im vorliegenden Anwendungsbeispiel, um ein dynamisches nicht-lineares Modell mit einer Vielzahl von Modell-Parametern, werden geeignete Start-Werte $\theta_0$ mithilfe verschiedener Methoden ermittelt. Die Generierung dieser Start-Wert-Kombinationen kann beispielsweise (voll-)kombinatorisch oder mittels verschiedener stochastischen bzw. probabilistischen Sampling-Strategien [6–8] erfolgen.

Die optimalen Beträge der Modell-Parameter $\Omega$ werden mittels Minimierung nachfolgender Zielfunktion $J\left(\theta\right)$ ermittelt.

Diese Zielfunktion summiert die quadrierten Differenzen zwischen der beobachteten Größe $y_{\text{data}}\left(t_{j,i}\right)$ und der simulierten Größe $y_\text{model}\left(\theta,t_{j,i}\right)$. Die Indices $i$ ($i=1,\ldots,n$) und $j$ ($j=1,\ldots,k$) referenzieren den diskreten Wert ($i$) innerhalb der Zeitreihe des jeweils betrachteten Experiments ($j$).

Ergebnisse

Für dieses Anwendungsbeispiel wurden 106 Startwert-Kombinationen innerhalb der oben angegebenen Modell-Parameter-Grenzen nach dem Saltelli-Sampling-Schema erzeugt und mittels Zielfunktion $J(\theta)$ bewertet. Die nachfolgende Parameter-Identifikation mittels Minimierung der Zielfunktion ausgehend von geeigneten Startwert-Kombinationen konvergierte abschließend zu den in nachfolgender Tabelle aufgeführten Parameter-Beträgen.

Nachfolgende Abbildung visualisiert, dass die experimentelle Datenbasis durch das parametrierte Bioprozess-Modell in allen wesentlichen Systemeigenschaften abgebildet werden kann (zum Erweitern anklicken):

Das parametrierte Bioprozess-Modell kann im nachfolgenden zur Analyse (z. B. Sensitivitätsanalyse) oder zur Optimierung (z. B. Ausbeute-Maximierung) des bestehenden Bioprozesses eingesetzt werden.