Zakrzep rany jako falę podróżna: migracja komórek i front Fisher-KPP

Dwa sposoby zamykania przepustki przez warstwę komórek

Gdy monolawa komórek epiteljowych jest obrażona — w laboratorium często to robiono prosto, krzyżując płytkę z plastiku pipetką, klasyka testu wykrzywiania (scratch assay) — komórki na krawędzi przepustki nie czekają pasywnie, aż nowe komórki dotrą. Zamiast tego aktywnie krążą w głąb, za nimi komórki dzielą się i razem te dwie zachowania zamykają ranę jako jednolite, wspólnie migrujące płytki, a nie jako pojedyncze komórki błądzące losowo.

To dokładnie taka sama sytuacja, jaką analizowali Ronald Fisher oraz niezależnie Andrey Kolmogorov, Ivan Petrovsky i Nikolai Piskunov w 1937 roku dla zupełnie innej sprawy — rozprzestrzeniania się korzystnego gena po populacji. Matematyka jest taka sama, dlatego nazywa się to równaniem Fisher-KPP, niezależnie od tego, czy populacja to geny, słońce lub komórki skóry.

Równanie: difuzja plus rosnąca logistycznie

Ustalmy u(x,t) za lokalną gęstość komórek. Fisher-KPP mówi, że jej prędkość zmiany wynika z sumy terminu rozprzestrzeniającego się (losowe ruchy komórek, modelowane jako difuzja) i terminu wzrostu (proliferacji, modelowane jako rosnące logistycznie w granicach pojemności nosicielnej):

∂u/∂t = Dm * ∂²u/∂x² + k * u * (1 - u/K) ^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^^ ruch/rzucanie się wzrost logistyczny Dm = motylność komórek (koeficjent difuzji) k = prędkość rozrostu K = pojemność nosicielna (gęstość monolitu konfluente) Rozwiązane na prostej z początkowym otworzem, to równanie osiąga stan w postaci fali poruszającej się o stałą prędkość v = 2√(Dm·k). Ta formula jest jednym z najprostszych wyników biologii matematycznej — prędkość uzdrowienia zależy tylko od motylności komórek i szybkości ich rozrodczości, a nie od szerokości pierwotnego otworu.

du/dt = Dm * d2u/dx2   +   k * u * (1 - u/K)
        ^^^^^^^^^^^^^      ^^^^^^^^^^^^^^^^^
        migration/spread    logistic proliferation

Dm = cell motility (diffusion coefficient)
k  = proliferation rate
K  = carrying capacity (confluent monolayer density)

Co się dzieje na krawędzi

Term diffuzji w równaniu Fisher-KPP jest matematycznym zastępcą czegoś bardzo fizycznego: komórki na najbardziej zaawansowanej krawędzi rany rozprostowują się, tworząc płaskie, powłokowate wystania nazywane lamelopodia, które są napędzane przez polymerizację aktynu działającą przeciwko membrannie komórkowej. Monomery aktynu dodają się do rosnących końców łańcuchów wraz z krawędzią liderową, membrana przesuwa się w przód, a nowe adhezje powstają, aby zapewnić przypięcie postępu — cykl ten powtarza się wiele razy na minutę. Komórki kilka rzędów za krawędzią, odporne na to, co dzieje się na krawędzi, głównie rozmnażają się i pasywnie są napięte lub przesuwane przez sąsiednich komórek, co dokładnie reprezentuje podział pracy, który dwie części równania przedstawiają.

Dlaczego fronat nie rozszerza się szybciej niż v

Fisher-KPP pozwala na całą kontynuum prędkości falowych większych lub równych 2√(Dm·k), ale tylko minimalna prędkość jest stabilnie dynamiczna i osiągana z realistycznych, skupionych początkowo danych — ranie o ostrej, skończonej krawędzi zamiast komórkach rozpraszających się do nieskończoności. To dlaczego rzeczywiste eksperymenty z leczeniem ran wiarygodnie pomiarują pojedynczą charakterystyczną prędkość zamknięcia, a nie zakres: niezależnie od dokładnej początkowej postaci szramy, system relaksuje się na tę samą minimalną-prędkościową fronat.

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Wound Healing Cell Migration Simulator i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.