Narzędzie analizuje serię czasową hodowli okresowej. Z kolumn czasu, biomasy, substratu i produktu wyznacza szybkość wzrostu przez regresję ln(X), czas generacji, bilanse stężeń, uzyski, średnie produkcyjności objętościowe oraz przybliżone szybkości właściwe odniesione do całki biomasy.

Wszystkie wiersze są traktowane jako jeden wybrany odcinek analizy. Przed wklejeniem danych trzeba więc wybrać biologicznie spójny zakres — na przykład fazę wykładniczą dla μ albo pełny przedział produkcyjny dla bilansu. Wysokie R² nie dowodzi, że punkty należą do właściwej fazy.

Przejdź do kalkulatora

Format danych i zakres obliczeń

Każdy wiersz powinien zawierać cztery liczby:

czas; X; S; P

X jest stężeniem biomasy albo spójną miarą proporcjonalną do biomasy, S stężeniem substratu, a P stężeniem produktu. Czas musi rosnąć. Kalkulator odrzuca puste linie, ale nie uzupełnia brakujących wartości. Jeśli produktu nie mierzono, można wpisać zera, pamiętając, że wynik Y_P/S nie ma wtedy wartości eksperymentalnej.

Dane domyślne są przykładem dydaktycznym, nie przebiegiem konkretnego organizmu. Wyniki zależą od wybranego przedziału. Ta sama seria analizowana dla godzin 0–3 może dać inne μ niż dla 2–5 h, ponieważ faza wzrostu się zmienia.

Regresja ln(X) zamiast jednej pary punktów

Wzrost wykładniczy zapisuje się:

X(t) = X₀e^(μt)

Po logarytmowaniu:

ln(X) = ln(X₀) + μt

Kalkulator dopasowuje prostą metodą najmniejszych kwadratów do wszystkich wprowadzonych punktów t, ln(X). Nachylenie jest oszacowaniem μ, a wyraz wolny oszacowaniem ln(X₀) dla tej prostej. Czas generacji wynosi g = ln(2)/μ, jeśli μ > 0.

To podejście wykorzystuje więcej informacji niż obliczenie z początku i końca, ale nie usuwa założeń. Klasyczna regresja zakłada między innymi właściwą postać modelu, niezależność oraz odpowiednią strukturę błędu. OD ma często błąd bardziej zbliżony do addytywnego w skali sygnału, a JTK może mieć zmienność zależną od liczby kolonii. Prosta regresja logarytmiczna jest narzędziem roboczym, nie automatycznie optymalnym modelem statystycznym.

Co mówi R², a czego nie mówi

Współczynnik determinacji pokazuje, jaka część zmienności ln(X) w wybranym zakresie jest wyjaśniona liniową zależnością od czasu. R² bliskie 1 jest zgodne ze wzrostem wykładniczym, ale może wystąpić także dla krótkiego, łagodnie zakrzywionego odcinka.

Wysokie R² nie sprawdza:

  • czy OD pozostaje w zakresie liniowym;
  • czy wszystkie punkty pochodzą z niezależnych próbek;
  • czy faza lag została wykluczona;
  • czy zmieniła się morfologia albo agregacja;
  • czy reszty nie mają systematycznego łuku;
  • czy nachylenie jest powtarzalne biologicznie.

Należy obejrzeć wykres ln(X) i reszty. Kalkulator raportuje liczbę punktów oraz R², ale nie podejmuje za użytkownika decyzji o fazie wzrostu.

Zmiany netto w hodowli okresowej

Dla pierwszego i ostatniego punktu:

ΔX = X_k − X_0

ΔS_użyte = S_0 − S_k

ΔP = P_k − P_0

Znaki są ważne. Zużycie substratu jest dodatnie, gdy końcowe S jest mniejsze. Ujemne ΔS_użyte oznacza wzrost zmierzonego stężenia, błąd, uwolnienie oznaczanego składnika albo niezgodność próbek — nie „ujemne zużycie” do podstawienia w uzysk.

Zmiana produktu może być ujemna, jeśli produkt ulega degradacji lub jest dalej metabolizowany. Wtedy średnia produkcyjność netto nie opisuje całkowitej ilości kiedykolwiek utworzonej.

Uzyski z bilansu końców przedziału

Przy dodatnim zużyciu substratu i dodatnim przyroście biomasy:

Y_X/S = ΔX/ΔS_użyte

Y_P/S = ΔP/ΔS_użyte

Y_P/X = ΔP/ΔX

Jeżeli X, S i P podano jako g/L, jednostkami są odpowiednio gX/gS, gP/gS i gP/gX. Stężenie objętościowe skraca się tylko wtedy, gdy wszystkie wielkości dotyczą tej samej objętości i ta objętość jest stała.

W hodowli okresowej ubytek przez pobieranie próbek, parowanie i dodawanie kwasu lub zasady może zmienić bilans. Dla istotnych zmian objętości należy liczyć całkowite masy m = C·V w każdym punkcie, a nie proste różnice stężeń.

Średnia produkcyjność objętościowa

Kalkulator wyznacza:

r_X,avg = ΔX/Δt

r_P,avg = ΔP/Δt

r_S,avg = ΔS_użyte/Δt

Są to średnie szybkości netto na objętość w całym wybranym przedziale. Nie są chwilowymi pochodnymi. Jeśli produkt powstawał szybko przez godzinę, a następnie pozostawał stały przez cztery godziny, średnia pięciogodzinna rozmyje maksimum.

Do projektowania czasu zbioru użyteczna bywa produkcyjność P(t)/t albo przyrost od zaszczepienia podzielony przez czas. Trzeba jawnie podać, czy uwzględniono fazę przygotowania inokulum i downstream.

Szybkości właściwe z całki biomasy

Zużycie substratu i tworzenie produktu można odnieść do ekspozycji biomasy w czasie:

q_S,avg = ΔS_użyte / ∫X dt

q_P,avg = ΔP / ∫X dt

Kalkulator przybliża całkę regułą trapezów:

∫Xdt ≈ Σ [(X_i + X_(i+1))/2]·(t_(i+1) − t_i)

Jeśli X ma jednostkę gX/L, a S gS/L, q_S ma jednostkę gS/(gX·czas). Wartości są średnimi dla przedziału. Reguła trapezów jest dokładniejsza przy gęstych punktach; przy szybkim wzroście i dużych odstępach liniowa interpolacja X może zawyżać lub zaniżać całkę.

Dla idealnie wykładniczego odcinka między dwoma punktami dokładna średnia logarytmiczna biomasy wynosi (X₂−X₁)/ln(X₂/X₁). Narzędzie używa trapezów, ponieważ seria może obejmować zmienną szybkość i wiele punktów. Przy rzadkich danych warto porównać oba podejścia.

OD jako X: co wolno policzyć

Jeżeli X jest OD, regresja ln(OD) może oszacować μ, o ile OD po blanku jest dodatnie, liniowe względem biomasy i morfologia pozostaje porównywalna. Czas generacji nie wymaga przeliczenia OD na g/L, ponieważ stały współczynnik kalibracji znika w logarytmie ilorazu.

Uzysk Y_X/S z OD ma jednak jednostkę OD na jednostkę substratu i nie jest masowym uzyskiem biomasy. Aby uzyskać gX/gS, trzeba zastosować zweryfikowaną relację OD–sucha masa dla danego organizmu, aparatu, drogi optycznej i zakresu.

Nie należy przeliczać punktów nasyconych przez jedną prostą bez rozcieńczenia. Spektrofotometria i pomiar gęstości optycznej omawia geometrię i zakres.

JTK jako X: interpretacja dotyczy jednostek hodowlanych

JTK/mL można użyć w regresji wzrostu jednostek tworzących kolonie. Wynik μ opisuje zmianę JTK, która łączy podziały, śmierć, agregację, rozpad agregatów i odzysk na podłożu. Nie jest automatycznie szybkością tworzenia wszystkich komórek.

Uzysk JTK/gS jest operacyjnym wynikiem dla metody posiewu, a nie stechiometrycznym uzyskiem biomasy. Jeśli celem jest bilans masy, potrzebna jest sucha masa albo inna ilościowa kalibracja.

Zero JTK nie może wejść do ln(X). Wynik poniżej granicy oznaczenia powinien być przechowany jako cenzurowany z limitem, nie zastępowany arbitralnie jedynką.

Substrat ograniczający musi być rzeczywiście zdefiniowany

Y_X/S ma sens, gdy S oznacza składnik, którego ubytek jest mierzony selektywnie. Brix, całkowity węgiel lub ChZT mogą obejmować wiele związków. Spadek jednego cukru może być równoważony pojawieniem się innego związku reagującego w teście.

Przy kilku źródłach węgla komórki mogą zużywać je sekwencyjnie lub równolegle. Jeden bilans całego odcinka ukryje diauksję. Wtedy warto podzielić serię według zmiany substratu i fizjologii, a nie tylko według czasu.

Produkt może być związany ze wzrostem albo nie

W modelu Luedekinga–Pireta chwilowa szybkość tworzenia produktu może mieć część związaną ze wzrostem i część niezależną:

dP/dt = α·dX/dt + β·X

Kalkulator nie estymuje α i β, ponieważ wymagałoby to odpowiedniej liczby punktów, pochodnych i oceny identyfikowalności. Raportuje średnie Y_P/X oraz q_P. Stałość Y_P/X w kilku odcinkach może sugerować sprzężenie ze wzrostem, ale nie jest dowodem mechanizmu.

Produkt wtórny może pojawiać się po spowolnieniu wzrostu. Analiza wyłącznie fazy wykładniczej da wtedy dobry μ i niemal zerowy Y_P/X, choć cały proces jest produkcyjny później.

Korekta objętości i pobierania próbek

Narzędzie zakłada stałą objętość i pracuje na stężeniach. Jeżeli objętość V_i jest znana dla każdego punktu, bilans należy wykonać na masach:

Δm_X = X_kV_k − X_0V_0 + m_X,pobrane

Analogicznie dla substratu i produktu. Suma materiału pobranego we wcześniejszych próbkach nie wraca do reaktora i może być istotna w małej skali.

Dodatek zasady, kwasu, antypianki lub zasilania także zmienia objętość i wnosi składniki. Taki proces nie jest ściśle hodowlą okresową, nawet jeśli główny substrat dodano tylko na początku.

Dane poniżej LoQ i wartości ujemne po blanku

Kalkulator wymaga X > 0, ponieważ oblicza logarytm. Mała dodatnia wartość po odjęciu blanku może jednak mieć ogromny względny błąd. Punktu poniżej LoQ nie należy ratować przez dodanie stałej tylko po to, aby regresja działała.

Dla S i P ujemna wartość po korekcie może sygnalizować szum wokół blanku. Narzędzie przyjmie liczbę do bilansu, ale użytkownik powinien uwzględnić granicę oznaczalności i nie interpretować małego ujemnego stężenia fizycznie.

Kontrole jakości serii

Kontrola Co sprawdza Typowy sygnał błędu
blank X tło biomasy/OD dodatni offset zawyżający wczesne punkty
standard S/P stabilność metody chemicznej dryf całej serii
bilans czasu faktyczny moment pobrania próbki opisane planem, nie zegarem
rozcieńczenie wspólną skalę skok o czynnik 10 lub 100
powtórzenia precyzję pobrania i analizy pojedynczy punkt napędzający nachylenie
wykres reszt postać modelu łuk wskazujący lag lub hamowanie
objętość ważność bilansu stężeń częste pobrania w małym reaktorze

Jak wybierać zakres do regresji

Nie należy wybierać punktów wyłącznie tak, aby maksymalizować R². Najpierw potrzebne są kryteria biologiczne i analityczne: zakończenie lag, dodatni sygnał ponad LoQ, brak nasycenia, brak jawnego załamania przez substrat lub tlen.

Dobrym podejściem jest porównanie kilku sąsiadujących okien, wykres lokalnych nachyleń i reszt, a następnie zapisanie reguły wyboru przed analizą właściwych powtórzeń. Ręczne usunięcie punktu odstającego wymaga przyczyny technicznej albo jawnej analizy wrażliwości.

Jeżeli celem jest maksymalne μ, można raportować najbardziej liniowy odcinek fazy wykładniczej. Jeżeli celem jest średnia szybkość całego procesu, nie wolno nazywać jej μ_max; trzeba użyć bilansu odpowiedniego przedziału.

Domyślny przykład krok po kroku

Seria rośnie od 0,10 do 2,00 gX/L w 5 h. Iloraz końców wynosi 20, ale μ kalkulatora pochodzi ze wszystkich sześciu punktów. Substrat maleje z 20,0 do 4,8 gS/L, więc zużycie netto wynosi 15,2 gS/L. Produkt rośnie z 0 do 3,2 gP/L.

Z bilansu końców:

Y_X/S = (2,00−0,10)/(20,0−4,8) = 1,90/15,2 = 0,125 gX/gS

Y_P/S = 3,2/15,2 ≈ 0,211 gP/gS

Y_P/X = 3,2/1,9 ≈ 1,684 gP/gX

Te liczby są tylko arytmetyką przykładu. Tak wysoki Y_P/X może być możliwy dla produktu wydzielanego, ale powinien skłonić do sprawdzenia podstawy jednostek i bilansu węgla.

Dlaczego bilans węgla jest mocniejszy niż pojedynczy uzysk

Substrat nie znika wyłącznie do biomasy i produktu docelowego. Węgiel może przejść do CO2, produktów ubocznych, EPS i związków pozostających w supernatancie. Tlen wpływa na stopień utlenienia produktów. Uzysk masowy bez korekty składu pierwiastkowego nie zamyka bilansu węgla.

Jeśli suma przypisanych produktów przekracza zużyty substrat na tej samej podstawie masowej, możliwe są inne źródła w pożywce, woda w produkcie, niezgodne jednostki albo błąd metody. Jeśli suma jest bardzo mała, brakuje głównego strumienia, zwykle gazowego lub nieoznaczonego metabolitu.

Powtórzenia i agregowanie wyników

Każdą niezależną hodowlę najlepiej dopasować osobno, a następnie porównać μ, uzyski i produkcyjności. Sklejenie wszystkich punktów trzech hodowli w jedną regresję sztucznie zwiększa liczbę obserwacji i ignoruje strukturę powtórzeń.

Średnia OD w każdym czasie może być użyteczna do wykresu, ale regresja średniej nie zachowuje pełnej zmienności. Raport powinien podawać liczbę hodowli biologicznych, sposób obchodzenia się z powtórzeniami technicznymi i miarę rozrzutu.

Powtórzenie techniczne a biologiczne rozwija różnicę, a Karty kontrolne, badania biegłości i efekt partii pokazuje, jak wykrywać zmiany systemowe.

Jak raportować analizę

Minimalny raport zawiera tabelę surowych danych, jednostki, korektę blanku, współczynniki rozcieńczenia, wybrany zakres, liczbę punktów, równanie dopasowania, R², μ, g, definicje X, S i P, zmiany objętości oraz wzory uzysków.

Przykład:

Dla punktów 1–5 h dopasowano ln(X) do czasu metodą najmniejszych kwadratów. X oznaczało suchą masę w g/L; każdy punkt pochodził z osobnej próbki tej samej hodowli. Raportowano μ, g i R²; uzyski obliczono z różnicy pierwszego i ostatniego stężenia po korekcie o łączną objętość próbek. Szybkości właściwe oszacowano z całki X(t) regułą trapezów.

Powiązane narzędzia

Kalkulator wzrostu populacji i czasu generacji służy do szybkiego rachunku z dwóch punktów i prognozy. Kalkulator chemostatu i progu wypłukania stosuje bilans ciągły. Osobne narzędzia dla JTK i rozcieńczeń oraz dla pełnego bilansu substratu rozwijają etapy, których ten analizator nie automatyzuje.

Źródła i dalsza lektura

Zakres obliczeń oparto na ćwiczeniach 2, 3 i 7 książki Mikrobiologia ogólna i przemysłowa — ćwiczenia laboratoryjne oraz częściach D1–D3 Mikrobiologii — krótkich wykładów, zakupionych z licencją do wykorzystania wiedzy. Zapis regresji, uzysków i szybkości właściwych skonfrontowano z otwartymi materiałami inżynierii bioprocesowej.bioprocess,nist-regression