Narzędzie łączy dwa poziomy opisu hodowli ciągłej. Pierwszy jest hydrauliczny: strumień zasilania F, objętość robocza V, szybkość rozcieńczania D = F/V i nominalny czas retencji τ = V/F. Drugi jest biologiczny: model Monoda, w którym szybkość wzrostu zależy od stężenia substratu ograniczającego.

Narzędzie wyznacza stan stacjonarny idealnego chemostatu, przewidywane stężenie substratu i biomasy, produkcyjność biomasy oraz odległość od progu wypłukania. Nie jest sterownikiem procesu ani projektem hodowli. Nie uwzględnia tlenu, maintenance, śmierci, retencji komórek, inhibicji, wielu substratów ani nieidealnego mieszania.

Przejdź do kalkulatora

Najpierw bilans objętości

W idealnym chemostacie objętość robocza jest stała, ponieważ strumień dopływu jest równy strumieniowi odpływu. Szybkość rozcieńczania wynosi:

D = F/V

Jeżeli F podano w L/h, a V w L, wynik ma jednostkę h⁻¹. Nominalny czas retencji jest odwrotnością:

τ = 1/D = V/F

Przy F = 0,25 L/h i V = 1 L, D = 0,25 h⁻¹, a τ = 4 h. Oznacza to, że w ciągu czterech godzin przez reaktor przepływa objętość równa objętości roboczej. Nie oznacza, że każda cząsteczka lub komórka przebywa dokładnie cztery godziny. Idealnie mieszany zbiornik ma rozkład czasów przebywania: część materiału opuszcza go wcześniej, część później.

Jeżeli dopływ i odpływ nie są równe, prosty wzór opisuje chwilowe F/V, ale objętość zmienia się w czasie. Narzędzie zakłada stałe V.

Stan stacjonarny nie oznacza braku aktywności

Stan stacjonarny oznacza stałość obserwowanych stężeń w czasie, mimo ciągłego przepływu i przemian. Biomasa rośnie, a jednocześnie opuszcza reaktor. Substrat wpływa, jest zużywany i wypływa. W najprostszym bilansie biomasy, bez śmierci i retencji:

dX/dt = (μ − D)X

Dla dodatniej, stałej biomasy dX/dt = 0, więc:

μ = D

To klucz chemostatu: poniżej wypłukania szybkość wzrostu w stanie stacjonarnym dopasowuje się do szybkości rozcieńczania. Nie oznacza to, że D jest właściwością organizmu. D ustala operator przez przepływ i objętość, natomiast organizm może albo osiągnąć odpowiednią μ przy pozostałym substracie, albo zostać wypłukany.

Równanie Monoda wiąże wzrost z jednym substratem

Model Monoda zapisuje:

μ(S) = μₘₐₓ · S/(Kₛ + S)

μₘₐₓ jest maksymalną szybkością w warunkach modelu, a Kₛ stężeniem, przy którym μ = μₘₐₓ/2. Kₛ nie jest automatycznie stałą wiązania transportera ani uniwersalnym parametrem gatunku. Zależy od szczepu, substratu, środowiska, historii i sposobu estymacji.

Po połączeniu μ = D z równaniem Monoda otrzymuje się stężenie substratu w stanie nie-wypłukanym:

S* = KₛD/(μₘₐₓ − D)

Wzór wymaga D < μₘₐₓ. Dodatkowo wynik ma sens tylko wtedy, gdy S* < Sᵢₙ. Jeśli wyliczony substrat resztkowy byłby równy lub większy od zasilania, dodatni stan biomasy nie może zostać podtrzymany w tym modelu.

Próg wypłukania zależy od stężenia w zasilaniu

Popularny skrót „wypłukanie przy D = μₘₐₓ” jest granicą dla bardzo dużego Sᵢₙ w prostym modelu. Dokładniej, tuż przy zaniku biomasy substrat w reaktorze zbliża się do Sᵢₙ, więc krytyczna szybkość wynosi:

D_crit = μ(Sᵢₙ) = μₘₐₓ · Sᵢₙ/(Kₛ + Sᵢₙ)

Jeśli Sᵢₙ nie jest znacznie większe od Kₛ, D_crit może być wyraźnie niższe niż μₘₐₓ. Kalkulator raportuje oba: maksymalną szybkość parametru i próg wynikający z konkretnego zasilania.

W idealnym modelu:

  • D < D_crit pozwala na dodatni stan stacjonarny;
  • D = D_crit jest granicą, przy której X* dąży do zera;
  • D > D_crit prowadzi do wypłukania, X* ≈ 0 i S* ≈ Sᵢₙ.

W rzeczywistym procesie nie planuje się pracy dokładnie na granicy. Dryf pompy, zmiana objętości, temperatury, składu zasilania lub fizjologii może przesunąć próg.

Bilans substratu i przewidywana biomasa

Przy stałym uzysku Yₓ/ₛ, bez maintenance i produktów ubocznych:

X* = Yₓ/ₛ · (Sᵢₙ − S*)

Jeśli Y = 0,5 gX/gS, a zużyto 9 gS/L, model przewiduje 4,5 gX/L. To bilans stechiometryczny, nie bezpośrednia kinetyka. Uzysk może zmieniać się z szybkością wzrostu, dostępnością tlenu i powstawaniem produktów.

Strumień biomasy opuszczającej reaktor na objętość roboczą — produkcyjność objętościowa biomasy — wynosi:

P_X = D · X*

Jednostką jest stężenie biomasy na czas, np. g/L/h. Całkowity strumień masy biomasy to F·X*, np. g/h. Kalkulator podaje oba, zachowując etykiety jednostek użytkownika.

Jak czytać margines do wypłukania

Narzędzie pokazuje stosunek D/D_crit i względny margines (D_crit − D)/D_crit. Są to wskaźniki modelowe. Przykładowo D/D_crit = 0,8 oznacza pracę przy 80% przewidywanego progu, a nie „20% ryzyka”.

Ryzyko zależy od zmienności parametrów. Jeżeli μₘₐₓ oszacowano z niepewnością 15%, margines 5% jest pozorny. Wrażliwość rośnie blisko granicy, ponieważ mianownik μₘₐₓ − D w równaniu S* maleje i niewielka zmiana parametrów powoduje dużą zmianę resztkowego substratu.

Czego nie obejmuje prosty model

Śmierć i maintenance

Jeżeli biomasa traci aktywność z szybkością k_d, bilans przyjmuje postać dX/dt = (μ − k_d − D)X. W stanie dodatnim μ = D + k_d, więc wykorzystanie prostego μ = D zaniża wymaganą szybkość wzrostu.

Maintenance oznacza zużycie substratu na podtrzymanie funkcji bez tworzenia nowej biomasy. Wtedy stały Y przestaje poprawnie opisywać małe szybkości wzrostu. Model Pirt'a rozdziela zapotrzebowanie wzrostowe i podtrzymaniowe.

Retencja komórek

W perfuzji, membranowym bioreaktorze lub układzie z osiadaniem komórki nie opuszczają reaktora z tą samą szybkością co ciecz. Hydrauliczny D = F/V nadal opisuje rozcieńczenie składników rozpuszczonych, ale nie szybkość utraty biomasy. Pojawia się osobny czas retencji komórek.

Wiele substratów i tlen

Wzrost może ograniczać węgiel, azot, fosfor, tlen albo mikroelement. Równanie z jednym S nie rozstrzyga współograniczenia. Jeśli transfer tlenu jest mniejszy od zapotrzebowania, μₘₐₓ z dobrze natlenionej kolby nie opisuje reaktora.

Inhibicja

Wysokie stężenie substratu lub produktu może hamować wzrost. Wtedy μ(S) nie rośnie monotonicznie jak w modelu Monoda, a możliwe są dodatkowe stany stacjonarne. Narzędzie nie stosuje modelu Haldane'a ani inhibicji produktem.

Biofilm i agregaty

Przyczepienie do ścian, nośnika lub sondy zatrzymuje część biomasy. Odpływowa X nie reprezentuje wtedy całkowitej biomasy. Biofilm ma gradienty, których idealnie wymieszany model nie opisuje.

Pompa i objętość są parametrami pomiarowymi

Zadany przepływ pompy nie musi być rzeczywistym F. W pompie perystaltycznej wpływają na niego przewód, docisk, zużycie, ciśnienie wsteczne i lepkość. Kalibracja przez zebranie masy lub objętości w czasie powinna obejmować rzeczywistą konfigurację.

Objętość robocza może różnić się od podziałki naczynia. Piana, zanurzenie przewodów, pobieranie próbek i działanie kontroli poziomu zmieniają V. Jeśli F ma błąd 5% i V błąd 5%, D również ma niepewność, której kalkulator nie pokazuje bez danych wejściowych.

Osiągnięcie stanu stacjonarnego wymaga czasu

Po zmianie D chemostat nie przechodzi natychmiast do nowego stanu. Czas odpowiedzi zależy od rozcieńczania i kinetyki. Jedna objętość wymieniona nie oznacza pełnej stabilizacji. Przy idealnym mieszaniu udział znacznika początkowego maleje wykładniczo jako e^(−Dt); po jednym τ pozostaje około 37%, po trzech około 5%, po pięciu mniej niż 1%.

Biomasa i regulacja metaboliczna mogą stabilizować się wolniej. Kryterium stanu powinno obejmować kilka kolejnych pomiarów biomasy, substratu i produktu, a nie arbitralną liczbę czasów retencji.

Dane potrzebne do estymacji parametrów

μₘₐₓ można oszacować z fazy wykładniczej hodowli okresowej, ale warunki muszą odpowiadać chemostatowi. Kₛ jest trudniejszy: dane tylko przy wysokim S słabo identyfikują parametr, ponieważ μ jest wtedy bliskie μₘₐₓ dla wielu wartości Kₛ.

Lepsza seria obejmuje kilka stanów stacjonarnych D i pomiar S*, z uwzględnieniem błędu obu osi oraz alternatywnych modeli. Linearyzacje odwrotności mogą zniekształcać błędy; dopasowanie nieliniowe do danych surowych jest zwykle właściwsze.

Uzysk wyznacza się z bilansu zużytego substratu i wytworzonej biomasy, w zakresie, gdzie założenie stałości jest uzasadnione. Y z pojedynczej kolby nie musi przenosić się na oddychanie ograniczone tlenem.

Kontrole sensowności

Kontrola Co sprawdzić Typowy błąd
jednostki F i V mają tę samą objętość mL/h podzielone przez L bez konwersji
czas F, μₘₐₓ i wynik używają tej samej jednostki L/min z μₘₐₓ w h⁻¹
zasilanie Sᵢₙ > 0 i odpowiada substratowi modelu całkowity węgiel zamiast składnika ograniczającego
próg porównać D z μ(Sᵢₙ), nie tylko μₘₐₓ zbyt optymistyczny próg przy małym Sᵢₙ
biomasa X* ma zgodną podstawę z Y sucha masa kontra OD bez kalibracji
stan parametry są stabilne w kolejnych próbkach pojedynczy punkt po zmianie przepływu
bilans dopływ i odpływ utrzymują V parowanie albo piana zmienia objętość

Przykład domyślny

Dla F = 0,25 L/h i V = 1 L, D = 0,25 h⁻¹ oraz τ = 4 h. Przy μₘₐₓ = 0,6 h⁻¹, Kₛ = 0,1 g/L i Sᵢₙ = 10 g/L:

D_crit = 0,6·10/(0,1+10) ≈ 0,594 h⁻¹

S* = 0,1·0,25/(0,6−0,25) ≈ 0,0714 g/L

Przy Yₓ/ₛ = 0,5 gX/gS:

X* = 0,5·(10−0,0714) ≈ 4,964 g/L

P_X = 0,25·4,964 ≈ 1,241 g/L/h

Model przewiduje dodatni stan, bo D/D_crit ≈ 0,421. Wynik nie jest gwarancją, że rzeczywisty organizm osiągnie prawie 5 g/L; wymaga poprawnych parametrów i spełnienia założeń.

Przykład blisko progu

Przy tych samych parametrach i D = 0,58 h⁻¹, S* = 2,9 g/L, a X* = 3,55 g/L. Niewielkie zwiększenie D do D_crit ≈ 0,594 h⁻¹ prowadzi w modelu do S* → 10 g/L i X* → 0.

To pokazuje nieliniowość granicy. Biomasa nie maleje proporcjonalnie do różnicy μₘₐₓ−D; blisko progu resztkowy substrat szybko rośnie.

Jak raportować obliczenie

Raport powinien podać rzeczywisty przepływ z kalibracji, objętość roboczą, D, τ, organizm, temperaturę, skład zasilania, definicję substratu ograniczającego, źródło μₘₐₓ, Kₛ i Y, zastosowany model oraz obserwowane kryterium stanu.

Przykładowy zapis:

Chemostat 1,00 L zasilano z rzeczywistym przepływem 0,248 L/h, co dawało D = 0,248 h⁻¹ i τ = 4,03 h. Parametry Monoda dopasowane w tych samych warunkach wynosiły μₘₐₓ = 0,60 h⁻¹ i Kₛ = 0,10 g/L. Dla Sᵢₙ = 10,0 g/L model przewidywał D_crit = 0,594 h⁻¹; pracowano przy 42% progu. Stan uznano za ustalony po trzech kolejnych bilansach biomasy i substratu bez trendu przekraczającego kryterium metody.

Dlaczego wynik kalkulatora trzeba porównać z bilansem rzeczywistym

W stanie ustalonym strumień substratu dopływającego to F·Sᵢₙ, a odpływającego F·S*. Różnica powinna być zgodna z ilością przypisaną biomasie, produktom i oddychaniu. Jeśli obliczone zużycie znika bez odpowiednika, przyczyną może być niepełny zestaw produktów, sorpcja, wytrącanie, błąd analityczny albo niewłaściwa definicja substratu.

Analogicznie strumień biomasy F·X* warto porównać z szybkością zużycia substratu oraz bilansem węgla. Zgodność pojedynczego OD nie wystarcza, ponieważ współczynnik OD–sucha masa może zmieniać się z morfologią. Zamknięcie bilansu nie dowodzi poprawności modelu, ale duże niezamknięcie jest sygnałem, że co najmniej jedno założenie, oznaczenie lub strumień wymaga ponownej kontroli.

Powiązane narzędzia i artykuły

Kalkulator wzrostu populacji i czasu generacji wyznacza μ z pomiarów okresowych. Hodowla okresowa, zasilana i ciągła porównuje architekturę procesów. Co kontroluje bioreaktor wyjaśnia pętle pomiarowe, a Mieszanie, napowietrzanie i transfer tlenu pokazuje, dlaczego kinetyka biologiczna nie wystarcza bez zdolności aparatury.

Źródła i dalsza lektura

Równania opracowano na podstawie części D3 książki Mikrobiologia — krótkie wykłady oraz rozdziałów o hodowli ciągłej w Mikrobiologii i Mikrobiologii ogólnej i przemysłowej — ćwiczenia laboratoryjne, zakupionych z prawem do wykorzystania wiedzy. Terminologię i bilanse skonfrontowano również z materiałem NCBI.ncbi-continuous