Biomatematyczny eksperyment nad emergencją życia, symulator ewolucyjny IYEUNI

1. Główna hipoteza

Naszym celem nie jest ręczne zaprogramowanie „żywej istoty”, lecz stworzenie takich warunków, w których coś życiopodobnego stanie się opłacalne samo z siebie.

Hipoteza projektu jest następująca:

Jeśli w świecie 2D damy prostym bytom:

1. lokalne źródła energii,
2. mutujący i dziedziczny zapis struktury,
3. możliwość wiązania się w większe układy,
4. lokalną chemię zależną od dopasowania sekwencji,
5. koszt utrzymania złożoności,

to z samej dynamiki mogą zacząć wyłaniać się struktury przypominające proto-metabolizm, proto-katalizę, proto-organizację i proto-selekcję.

Nie chodzi nam o wierne odtworzenie biologii ziemskiej. Chodzi o znalezienie takiego zestawu reguł minimalnych, przy którym materia zacznie zachowywać się tak, jakby opłacało się jej utrzymywać własną ciągłość organizacyjną.

2. Świat jako chemo-mechaniczna zupa

W naszym modelu przestrzeń nie jest jednorodna. Świat ma wiele nisz środowiskowych i lokalne pola zasobów. Każde miejsce na planszy ma swoją pojemność środowiskową i aktualny poziom dostępnej energii.

To oznacza, że selekcja nie polega tylko na prostym „żyć lub umrzeć”, lecz na pytaniu: czy dany genom umie efektywnie korzystać właśnie z tego lokalnego krajobrazu zasobów.

Środowisko stale się regeneruje, może być okresowo niszczone przez katastrofy, a następnie powoli wracać do równowagi. Dodatkowo dopuszczamy spontaniczne tworzenie się nowych krótkich łańcuchów w pustych obszarach bogatych w energię chemiczną. Dzięki temu świat nie jest statyczną planszą, lecz aktywnym tłem, które samo generuje presję selekcyjną i okazjonalnie otwiera drogę nowym początkom.

3. Genom jako źródło metabolizmu, kształtu i funkcji

Genom jest w naszym podejściu centralny. To sekwencja A/U/G/C, ale nie traktujemy jej jako martwego opisu. Sekwencja bezpośrednio wpływa na:

• profil metaboliczny,
• potencjał katalityczny,
• lokalną komplementarność kontaktów,
• zdolność do kohezji i tworzenia kieszeni reakcyjnych,
• preferencje dotyczące funkcji chemicznych.

Innymi słowy: genotyp nie opisuje jedynie organizmu — genotyp staje się jego zachowaniem.

Profil metaboliczny

Z lokalnych fragmentów sekwencji budujemy profil metaboliczny rozłożony po dostępnych niszach. Jedna z użytych sztuczek matematycznych ma postać:

To celowo sztuczna, ale bogata funkcjonalnie mapa genotyp–środowisko. Nie próbujemy udawać rzeczywistej chemii, tylko tworzymy przestrzeń, w której różne genomy mogą wyspecjalizować się w różnych niszach.

Koszt złożoności

Długość genomu nie jest darmowa. Wprowadzamy nieliniowy koszt utrzymania informacji:

To ważny element całej architektury. Chcemy, aby wzrost złożoności był opłacalny tylko wtedy, gdy daje realną przewagę funkcjonalną. Jeśli długa sekwencja niczego nie poprawia, system sam ją ukarze energetycznie.

4. Ewolucja jako bilans energetyczny

W naszym modelu ewolucja nie jest zewnętrznym algorytmem optymalizacji. Zachodzi jako konsekwencja bilansu przetrwania.

W każdym kroku byt:

1. pobiera energię z otoczenia według swojego profilu metabolicznego,
2. płaci koszt genomu,
3. traci energię przez spontaniczny rozpad,
4. może umrzeć, jeśli jego energia spadnie do zera,
5. jeśli przeżyje, jego genom może mutować.

Dopuszczamy kilka klas zmian:

• mutację punktową,
• wydłużenie sekwencji,
• skrócenie sekwencji,
• duplikację fragmentu,
• abiotyczne wydłużanie w sprzyjającym środowisku.

To oznacza, że selekcja jest u nas emergentna: przetrwa nie to, co my oznaczymy jako dobre, lecz to, co realnie potrafi utrzymać dodatni bilans energii mimo kosztów własnej organizacji.

5. Wiązania jako początek organizmu

Pojedyncza cząstka to za mało, by mówić o proto-życiu. Dlatego wprowadziliśmy lokalne wiązania pomiędzy bytami. Energia wiązania zależy od zgodności lokalnych fragmentów sekwencji oraz ich ekspozycji powierzchniowej. To nie jest binarne „jest kontakt lub nie ma kontaktu”, lecz stopniowalna relacja o sile zależnej od jakości dopasowania.

Uproszczona zależność siły wiązania ma postać:

Im korzystniejsza energia kontaktu, tym mocniejsze połączenie.

To pozwala nam przejść od izolowanych cząstek do struktur, które zaczynają mieć topologię, geometrię i lokalną historię interakcji.

6. Aktywne miejsca i emergencja proto-enzymów

Najważniejsza część całej idei polega na tym, że funkcja chemiczna nie musi należeć do pojedynczego bytu. Może wyłaniać się z odpowiedniego ułożenia kilku bytów połączonych lokalnie.

Dlatego budujemy ActiveSite jako emergentną własność kontaktu, uwzględniając m.in.:

• komplementarność sekwencji,
• kohezję kontaktu,
• lokalne wsparcie sąsiadów,
• wielkość i rozgałęzienie klastra,
• różnorodność sekwencyjną,
• wagę lokalnych połączeń,
• zdolność do tworzenia kieszeni reakcyjnych.

Przykładowo stabilność aktywnego miejsca składamy jako ważoną kombinację kilku cech:

A jego finalną moc opisujemy syntetycznie przez:

To jest sedno naszej hipotezy: coś na kształt enzymu może pojawić się jako relacja geometryczna i sekwencyjna, a nie jako pojedynczy z góry zdefiniowany obiekt.

7. Chemia: polimeryzacja, ligacja, cięcie, restrykcja

Gdy aktywne miejsce osiągnie odpowiednią jakość, może uruchomić reakcję. W naszym modelu aktywne miejsca mogą prowadzić do czterech typów operacji:

• polimeryzacji,
• ligacji,
• cięcia,
• restrykcji.

Prawdopodobieństwo reakcji zależy od wielu lokalnych czynników: siły wiązania, mocy aktywnego miejsca, zgodności motywu, stabilności kontaktu i dostępnej energii wykonawcy.

Dzięki temu chemia nie jest tłem dekoracyjnym. Chemia zmienia materiał dziedziczny, długości sekwencji i zdolność układów do dalszego organizowania się.

8. Mechanika jako warunek powstania ciała

Sama chemia nie wystarczy. Jeśli struktura nie potrafi utrzymać własnej topologii podczas ruchu, nie stanie się zaczątkiem organizmu.

Dlatego wprowadziliśmy:

• ruch ciągły z tłumieniem,
• pamięć poprzedniego przemieszczenia,
• lokalny consensus ruchu struktur związanych,
• projekcję wiązów,
• procedury rescue,
• procedury spill,
• naprawę topologii po rasteryzacji,
• lokalne wyrównanie kątów.

To wszystko służy jednemu celowi: chcemy, aby struktura wielocząstkowa nie rozpadała się tylko dlatego, że świat jest liczony na siatce.

Właśnie tu znajduje się dziś jeden z najważniejszych problemów badawczych. Granica między „naturalnym rozpadem” a „artefaktem numerycznym” jest bardzo cienka. Jeśli organizmy rozpadają się przez błędy mapowania ruchu na raster, nie uczymy się niczego o emergencji życia — uczymy się tylko ograniczeń własnej implementacji. Dlatego tyle uwagi poświęcamy naprawie wiązań, wyrównaniu geometrii i stabilności struktur po ruchu.

9. Energetyka współpracy

Złożone struktury nie mogą istnieć za darmo. Jeśli kilka bytów tworzy wspólną organizację, ta organizacja musi na siebie zarobić.

Dlatego wprowadziliśmy:

• koszt utrzymania wiązań,
• koszt naprężenia dalekich wiązań,
• ryzyko zerwania przeciążonych relacji,
• lokalne współdzielenie energii po wiązaniach,
• dyfuzję energii pomiędzy sąsiadami.

Dyfuzja energii ma uproszczoną postać:

Koszt utrzymania połączenia opisujemy z kolei jako:

W ten sposób próbujemy sprawdzić, czy układy związane mogą uzyskać przewagę netto dzięki współpracy, czy też będą systematycznie przegrywać z prostszymi, tańszymi formami.

10. Jakie życie próbujemy zobaczyć

Nie oczekujemy, że z tego modelu od razu wyłoni się pełny organizm w biologicznym sensie. Szukamy wcześniejszego progu.

Interesuje nas pojawienie się układów, które:

• dłużej utrzymują własną strukturę,
• stabilizują przepływ energii,
• tworzą funkcjonalne lokalne kontakty,
• potrafią chemicznie modyfikować siebie lub sąsiadów,
• dziedziczą korzystne cechy i podlegają selekcji.

Jeśli zobaczymy trwałe struktury, które utrzymują się nie dlatego, że są ręcznie wymuszone, lecz dlatego, że opłaca się im istnieć, to właśnie tam zaczniemy mówić o proto-życiu.

11. Najbliższe cele

Z obecnej postaci kodu wynikają dla nas bardzo konkretne najbliższe zadania.

1. Dalsza stabilizacja topologii struktur

Chcemy ograniczyć przypadki, w których organizacja wielocząstkowa rozpada się wskutek konfliktów rasteryzacji, złego rozwiązywania ruchu lub nadmiernego rozciągania wiązań.

2. Lepszy balans chemii

Musimy nadal stroić parametry tak, aby reakcje chemiczne były produktywne i selekcyjnie istotne, ale nie dominowały świata czystym chaosem.

3. Realna opłacalność współpracy

Zależy nam na tym, by układy związane uzyskiwały przewagę tylko wtedy, gdy ich organizacja naprawdę poprawia metabolizm, stabilność lub katalizę.

4. Długie obserwacje emergencji

Symulacja jest zbudowana tak, by można ją było streamować i oglądać przez długi czas. To oznacza, że chcemy nie tylko uruchamiać krótkie eksperymenty, ale śledzić wielogodzinne przebiegi i patrzeć, czy pojawiają się trwałe reżimy organizacji.

12. Synteza końcowa

Nasz projekt jest próbą odpowiedzi na pytanie, czy w świecie 2D można stworzyć warunki, w których materia zacznie organizować się tak, jakby odkryła własny interes w przetrwaniu.

Nie chcemy pisać gotowego życia. Chcemy pisać warunki, z których życie mogłoby zacząć się wyłaniać.

Dlatego łączymy genom, energetykę, chemię kontaktową, mechanikę wiązań i ewolucję w jednym systemie. Jeśli z takiego układu zaczną pojawiać się trwałe, selekcyjnie uprzywilejowane struktury zdolne do podtrzymywania własnej organizacji, będzie to dla nas sygnał, że zbliżyliśmy się do uchwycenia czegoś bardzo istotnego: nie biologii w pełnym sensie, lecz momentu poprzedzającego biologię.