PL 
EnglishDomyślny
Українська🧬 Zegar molekularny — oscylator PER/CRY/BMAL1
Symuluj pętlę ujemnego sprzężenia zwrotnego transkrypcja-translacja PER/CRY/BMAL1 równaniami oscylatora Goodwina. Zobacz oscylacje stężeń białek z okresem 24 h.
O tej symulacji
Ta symulacja animuje oscylator Goodwina, klasyczny minimalny model molekularnego zegara ssaków. Trzy sprzężone zmienne — mRNA (M), cytoplazmatyczne białko PER/CRY (Pc) i jądrowe PER/CRY (Pn) — są całkowane za pomocą solvera Rungego-Kutty czwartego rzędu. Jądrowe PER/CRY hamuje własny gen poprzez funkcję Hilla z wykładnikiem n, i to właśnie ta opóźniona pętla ujemnego sprzężenia zwrotnego, a nie zewnętrzny zegar, tworzy samopodtrzymujący się rytm ~24-godzinny widoczny na wykresie szeregu czasowego po prawej stronie.
🔬 Co pokazuje
Lewy panel to schematyczna komórka: mRNA opuszcza jądro, jest tłumaczone na cytoplazmatyczne PER/CRY, a to białko przemieszcza się do jądra (tempa k₁, k₂), gdzie wyłącza własną transkrypcję. Prawy panel przedstawia M, Pc i Pn w ciągu ostatnich 24 symulowanych godzin, z pasmami dnia/nocy i bieżącym oszacowaniem okresu na podstawie kolejnych przejść M przez próg.
🎮 Jak korzystać
Przeciągnij współczynnik Hilla n, degradację mRNA k_d, degradację białka v_d oraz tempa transportu cytoplazma↔jądro k₁/k₂, aby zmienić kształt cyklu, lub przejdź od razu do gotowego scenariusza: Normalny zegar, Mutacja Tau (szybki, ~20 godz.), Opóźniona faza (wolny) lub Arytmiczny (n=4). Prędkość kontroluje tempo odtwarzania; Odtwórz/Pauza i Reset pozwalają zatrzymać lub zrestartować całkowanie w dowolnym momencie.
💡 Czy wiesz, że...
Poniżej współczynnika Hilla wynoszącego około n = 8 pętla sprzężenia zwrotnego jest zbyt "miękka", by podtrzymać oscylacje, a zegar zanika do płaskiego stanu ustalonego — to scenariusz Arytmiczny. Prawdziwe geny zegarowe osiągają efektywne n znacznie powyżej tego progu dzięki wielu kolejnym etapom fosforylacji i tworzenia kompleksów — dokładnie to zaburza klasyczna mutacja tau kinazy CK1ε, skracając cykl do około 20 godzin.
Najczęściej zadawane pytania
Co dokładnie modeluje tu oscylator Goodwina?
To trójzmienny układ równań różniczkowych reprezentujący własną pętlę sprzężenia zwrotnego transkrypcji-translacji genu: mRNA (M) jest tłumaczone na białko cytoplazmatyczne (Pc), które przemieszcza się do jądra (Pn) i hamuje dalszą transkrypcję własnego genu poprzez wyraz represji typu Hilla. To podręcznikowy mechanizm proponowany dla pętli PER/CRY/BMAL1, która napędza dobowy zegar ssaków.
Dlaczego układ oscyluje, zamiast się ustabilizować?
Hamowanie transkrypcji przez białko jądrowe jest opóźnione (musi zostać najpierw przetłumaczone w cytoplazmie i przetransportowane do jądra) i silnie nieliniowe, gdy współczynnik Hilla n jest wystarczająco duży. Właśnie to połączenie opóźnienia i stromej nieliniowości odpycha układ od stabilnej równowagi w kierunku powtarzającego się cyklu granicznego, zamiast prostego zaniku do stanu ustalonego.
Co dokładnie kontroluje współczynnik Hilla n?
Określa, jak gwałtownie włącza się represja w miarę gromadzenia się jądrowego PER/CRY. Niskie n daje stopniową, łagodną odpowiedź, która zanika; symulacja wymaga w przybliżeniu n większego lub równego 8, aby sprzężenie zwrotne było wystarczająco strome do podtrzymania rytmicznych oscylacji, a przy n = 4 (scenariusz Arytmiczny) zegar całkowicie przestaje oscylować.
Co powoduje, że scenariusz Mutacja Tau przyspiesza zegar?
Scenariusz Mutacja Tau podnosi tempo degradacji białka v_d oraz tempa transportu k1/k2, naśladując efekt mutacji kinazy CK1(tau), która sprawia, że białko PER rozkłada się szybciej u prawdziwych zwierząt. Szybszy obrót białka skraca opóźnienie w pętli sprzężenia zwrotnego, co skraca okres do około 20 godzin zamiast normalnych ~24.
Jak mierzony jest na ekranie okres oscylacji?
Symulacja obserwuje ślad mRNA M, wyczekując momentów, w których przekracza on rosnąco ustalony próg (0,5), i zapisuje symulowany czas między kolejnymi przejściami jako szacowany okres, pokazywany na pasku statystyk. Ponieważ zależy to od dwóch kolejnych przejść, odczyt pozostaje pusty przez kilka sekund po Reset lub po zmianie parametrów, aż zostanie uchwycony nowy pełny cykl.