Narzędzie przelicza wzrost populacji w fazie wykładniczej między czterema równoważnymi opisami: liczbą generacji, czasem generacji, właściwą szybkością wzrostu μ i współczynnikiem namnażania. Może również prognozować wielkość populacji albo czas dojścia do zadanego poziomu.

Wynik jest poprawny tylko w granicach modelu. Narzędzie nie rozpoznaje fazy lag, ograniczenia substratu, śmierci, agregacji ani nasycenia pomiaru. Dane z OD, JTK, bezpośredniego liczenia i suchej masy można przeliczać tymi samymi równaniami wyłącznie wtedy, gdy oba punkty pochodzą z tej samej metody i jej zakresu roboczego.

Kalkulatory: parametry wzrostu z dwóch pomiarów, prognoza wzrostu populacji oraz czas wzrostu do wartości docelowej.

Co oblicza tryb „Dwa pomiary”

Dla wartości początkowej N₀, końcowej Nₜ i odstępu Δt kalkulator wyznacza:

  • iloraz Nₜ/N₀, czyli zmianę względną;
  • liczbę generacji n = log₂(Nₜ/N₀);
  • liczbę podwojeń w jednostce czasu k = n/Δt;
  • czas generacji g = Δt/n, gdy populacja wzrosła;
  • właściwą szybkość wzrostu μ = ln(Nₜ/N₀)/Δt;
  • czas podwojenia t_d = ln(2)/μ, równy g dla stałego wzrostu wykładniczego.

Przykład domyślny prowadzi od 10⁶ do 8·10⁶ w trzy godziny. Iloraz wynosi 8, czyli zaszły trzy podwojenia; g = 1 h, k = 1 h⁻¹, a μ ≈ 0,693 h⁻¹.

Jeżeli drugi pomiar jest mniejszy, narzędzie raportuje ujemną zmianę logarytmiczną i czas połowicznego spadku. Nie nazywa ujemnej wartości „ujemnym czasem generacji”, ponieważ byłoby to mylące biologicznie. Spadek może wynikać ze śmierci, sedymentacji, lizy, rozcieńczenia, pobierania próbki albo zmiany odpowiedzi pomiarowej.

Model wykładniczy i jego równoważne zapisy

W idealnym wzroście przez podział binarny:

N(t) = N₀ · 2^(t/g)

Ten sam model można zapisać podstawą liczby Eulera:

N(t) = N₀ · e^(μt)

Z porównania wykładników wynika μ = ln(2)/g. Parametr μ ma wymiar odwrotności czasu. Jeśli czas podano w godzinach, jednostką jest h⁻¹; czas generacji pozostaje w godzinach. Liczba generacji i iloraz populacji są bezwymiarowe.

Niekiedy podręczniki oznaczają liczbę generacji na jednostkę czasu literą k. Wtedy k = 1/g, natomiast μ = k·ln(2). Nie wolno przepisywać wartości między publikacjami bez sprawdzenia definicji symbolu: litera k bywa również używana jako stała śmierci, a r jako tempo wzrostu netto.

Prognoza odejmuje fazę lag, ale nie modeluje jej mechanizmu

Tryb „Prognoza” oblicza liczbę generacji z efektywnego czasu t_eff = max(0, t − lag), a następnie N = N₀·2^(t_eff/g). Jeśli cały zadany czas mieści się w lag, wynik pozostaje równy N₀.

To uproszczenie jest celowo jawne. W prawdziwej hodowli przejście z lag do wzrostu nie musi być nagłe. Poszczególne komórki rozpoczynają podziały w różnym czasie, zmienia się masa na komórkę, a OD może rosnąć przed wzrostem JTK. Parametr lag wpisany do kalkulatora jest założeniem użytkownika, nie wartością wywnioskowaną z danych.

Tabela prognozy pokazuje wartości w równych odstępach całkowitego czasu. Nie jest symulacją zasobu ani krzywej logistycznej. Przy dużej liczbie generacji wynik może być matematycznie ogromny i biologicznie niemożliwy, ponieważ model nie zawiera podłoża, tlenu, przestrzeni ani produktów hamujących.

Czas do celu jest dolnym oszacowaniem w sprzyjających warunkach

Dla Nₜ > N₀ potrzebna liczba generacji wynosi:

n = log₂(Nₜ/N₀)

Czas wzrostu to t = n·g, a po dodaniu jawnie założonej fazy lag: t_całk = lag + n·g. Wynik odpowiada chwili przecięcia idealnej krzywej wykładniczej z poziomem celu.

Wynik nie przewiduje, kiedy rzeczywista hodowla osiągnie zadaną OD, JTK lub masę. Jeżeli cel leży blisko pojemności środowiska, wzrost zwalnia i obliczenie będzie zbyt optymistyczne. Dla planowania doświadczenia wynik powinien służyć jako punkt do ustalenia wcześniejszych i późniejszych pobrań, nie jako dokładny termin końca.

Jak wybrać dwa punkty do obliczenia μ

Najlepsze są punkty z odcinka, na którym ln(N) jest w przybliżeniu liniowe względem czasu. Wybranie jednego punktu z lag i drugiego z fazy wykładniczej zawyży g. Wybranie punktu z początku i końca hodowli uśredni kilka różnych stanów fizjologicznych.

Praktyczna procedura analizy serii:

  1. odjąć właściwy blank i zachować wartości surowe;
  2. odrzucić punkty poza zweryfikowanym zakresem instrumentu, a nie tylko „ładnie” odstające;
  3. narysować zarówno N(t), jak i ln(N) względem czasu;
  4. wybrać spójny biologicznie odcinek liniowy;
  5. dopasować regresję do wszystkich punktów odcinka zamiast opierać finalny wynik na jednej parze;
  6. sprawdzić reszty, nie tylko R²;
  7. raportować zakres czasu, liczbę punktów, metodę pomiaru i niepewność nachylenia.

Kalkulator dwóch punktów jest dobry do nauki, szybkiej kontroli i odtworzenia prostego rachunku. Regresja jest lepsza do estymacji parametru z eksperymentu.

OD, JTK i liczba komórek nie są wymienne

Gęstość optyczna

OD jest sygnałem rozpraszania zależnym od instrumentu, drogi optycznej, morfologii i agregacji. Jeżeli oba punkty leżą w liniowym zakresie i własności komórek są stabilne, iloraz OD może przybliżać iloraz biomasy. Wartość nasycona przestaje rosnąć proporcjonalnie, więc kalkulator zaniży szybkość.

Rozcieńczenie próbki przywraca zakres tylko wtedy, gdy wynik przeliczy się przez współczynnik rozcieńczenia. Nie wolno porównywać surowego OD próbki nierozcieńczonej z odczytem rozcieńczonym.

JTK

JTK mierzy jednostki zdolne utworzyć kolonię w określonych warunkach. Jedna kolonia może pochodzić z jednej komórki albo agregatu. Zmiana agregacji w czasie zmienia JTK bez proporcjonalnej zmiany biomasy.

Punkty powinny pochodzić z policzalnych płytek i uwzględniać rozcieńczenie oraz objętość posiewu. Zero kolonii nie jest matematycznym zerem populacji i nie może być wprowadzone do logarytmu; oznacza wynik poniżej możliwości wykrycia dla zbadanej porcji.

Liczenie bezpośrednie i cytometria

Licznik komórek lub cytometr może ujmować żywe, martwe i uszkodzone zdarzenia zależnie od barwienia oraz bramkowania. Wzrost liczby całkowitej przy stałym JTK może wskazywać utratę zdolności hodowlanej albo błąd rozdzielenia agregatów.

Sucha masa i produkt

Sucha masa opisuje biomasę, nie liczbę podziałów. Komórki mogą zwiększać rozmiar albo gromadzić polimer bez proporcjonalnej zmiany liczby. Stężenie produktu nie powinno być podstawiane jako N, ponieważ kinetyka produktu może być sprzężona ze wzrostem, mieszana albo niezależna.

Porównanie miar rozwija OD, JTK, liczba komórek i sucha masa — różne miary biomasy.

Rozcieńczenie i pobieranie próbek mogą udawać wzrost lub spadek

W hodowli okresowej pobranie małej próbki zwykle pomija się w prostym modelu. Przy częstym pobieraniu znacznej objętości zmienia się jednak objętość robocza i bilans. Dodanie pożywki rozcieńcza stężenie komórek, choć całkowita liczba może nadal rosnąć.

Jeżeli N oznacza stężenie, po rozcieńczeniu trzeba rozdzielić trzy procesy:

  • namnażanie zmienia liczbę komórek;
  • dodanie lub usunięcie cieczy zmienia objętość;
  • pobranie usuwa część komórek.

Kalkulator nie wykonuje bilansu tych operacji. Do analizy fed-batch lub hodowli ciągłej potrzebny jest bilans masy i objętości.

Dlaczego nie można ekstrapolować bez ograniczeń

Po dziesięciu generacjach iloraz wynosi 2¹⁰ = 1024, po dwudziestu około miliona, a po czterdziestu około biliona. Niewielki błąd g jest wielokrotnie wzmacniany przy długiej prognozie.

Jeśli rzeczywisty czas generacji wynosi 30 minut, lecz przyjęto 27 minut, po dziewięciu godzinach modele różnią się o kilka generacji. To różnica wielokrotna, nie stały błąd addytywny. Dlatego prognoza na wiele generacji powinna pokazywać scenariusze dla granic niepewności.

Drugim ograniczeniem jest pojemność środowiska. Model wykładniczy zakłada stałe warunki na komórkę. W zamkniętej hodowli maleje substrat, zmienia się pH, rośnie stężenie metabolitów, a transfer tlenu może osiągnąć granicę. Krzywa wzrostu jest wynikiem tych procesów, co wyjaśnia Co naprawdę przedstawia krzywa wzrostu populacji.

Niepewność dwóch punktów

Parametr μ jest różnicą logarytmów podzieloną przez czas. Błąd względny małego pomiaru początkowego może silnie wpłynąć na wynik. Dwa punkty nie pozwalają odróżnić szumu od krzywizny i nie dają niezależnego testu liniowości.

Powtórzenia techniczne opisują zmienność pobrania i pomiaru tej samej próbki. Powtórzenia biologiczne obejmują niezależne hodowle i pokazują zmienność procesu. Uśrednienie wszystkich odczytów przed logarytmowaniem nie zawsze jest równoważne uśrednieniu indywidualnych szybkości.

Jeżeli dostępne są niepewności N₀, Nₜ i czasu, należy propagować je przez logarytm albo użyć bootstrapu/regresji. Sam kalkulator celowo nie tworzy pozornie dokładnego przedziału bez tych danych.

Kontrole sensowności wyniku

Przed zaakceptowaniem rachunku warto sprawdzić:

Kontrola Pytanie Sygnał problemu
faza czy oba punkty są w tym samym odcinku wzrostu? krzywizna ln(N)
zakres czy instrument odpowiada liniowo? wynik blisko nasycenia lub blanku
jednostka czy oba pomiary mają tę samą skalę? OD surowe kontra OD po rozcieńczeniu
objętość czy nie zmieniła się przez zasilanie/pobranie? skok stężenia bez zmiany liczby całkowitej
populacja czy miara obejmuje tę samą frakcję? zmiana agregacji albo żywotności
czas czy znaczniki czasu dotyczą faktycznego pobrania? opóźnienie między pobraniem i odczytem
biologia czy g jest możliwe przy danych warunkach? wynik wymagający nierealnego bilansu

Kalkulator ostrzega o liczbach skrajnych i spadku, ale nie ma informacji o organizmie, temperaturze ani podłożu. Kontrola biologiczna należy do użytkownika.

Jak raportować wynik

Minimalny zapis powinien brzmieć na przykład:

W odcinku 2–5 h OD po korekcie blanku wzrosło z 0,08 do 0,61. Z dwóch punktów otrzymano μ = 0,677 h⁻¹ i g = 1,02 h; oba odczyty przeliczono do wspólnej drogi optycznej. Wartość służy jako szybkie oszacowanie, a parametr końcowy wyznaczono regresją siedmiu punktów z trzech niezależnych hodowli.

Należy podać organizm lub układ, podłoże, temperaturę, atmosferę, metodę, blank, zakres, przedział czasu, liczbę punktów i sposób dopasowania. Sam zapis „czas generacji 30 min” nie pozwala odtworzyć wyniku.

Przykłady kontrolne

Cztery podwojenia

Jeżeli populacja rośnie z 2,5·10⁵ do 4,0·10⁶, iloraz wynosi 16, a n = 4. Przy odstępie 120 min czas generacji wynosi 30 min, k = 0,0333 min⁻¹, a μ ≈ 0,0231 min⁻¹.

Niepełna liczba generacji

Wzrost z 10⁶ do 5·10⁶ daje n = log₂(5) ≈ 2,322. Pokolenia populacji nie muszą być liczbą całkowitą, ponieważ komórki dzielą się asynchronicznie, a wynik jest opisem zbiorczym.

Spadek

Zmiana z 8·10⁶ do 2·10⁶ to iloraz 0,25 i n = −2. W ciągu czterech godzin odpowiada to dwóm połowicznym spadkom, czyli czasowi połowicznego spadku 2 h przy stałym tempie netto. Nie dowodzi to, że każda komórka umiera z jedną stałą szybkością.

Prognoza przekraczająca pojemność

Start od 10⁶ i czas generacji 20 min przez 10 h daje 30 generacji oraz około 1,07·10¹⁵ w jednostce wejściowej. Matematyka jest poprawna dla modelu, ale wynik biologiczny powinien wywołać ostrzeżenie: zamknięta hodowla niemal na pewno przestanie być wykładnicza wcześniej.

Granice zastosowania

Narzędzie nie służy do prognozowania zakażenia u człowieka, dawki patogenu ani dynamiki epidemii. „Czas generacji” w mikrobiologicznej hodowli oznacza czas podwojenia populacji w danych warunkach. W epidemiologii generation interval oznacza odstęp między zakażeniem źródła i osoby wtórnej; to inna wielkość, opisana w Reprodukcja zakażenia, wzrost epidemii i czas generacji.

Kalkulator nie zastępuje również modelu chemostatu, Monoda, Gompertza, Baranyiego ani modelu logistycznego. Każdy dodaje inne mechanizmy i parametry. Prostota modelu wykładniczego jest zaletą, gdy analizuje się właściwy odcinek, i wadą, gdy próbuje się nim opisać całą hodowlę.

Źródła i dalsza lektura

Równania i interpretację wzrostu opracowano na podstawie rozdziałów D1–D3 książki Mikrobiologia — krótkie wykłady, rozdziałów o hodowli i pomiarze w Mikrobiologii oraz Mikrobiologii ogólnej i przemysłowej — ćwiczenia laboratoryjne, zakupionych z licencją do republikacji. Definicje skonfrontowano z podręcznikiem OpenStax i materiałem NCBI o matematyce wzrostu bakterii.openstax,ncbi-growth