Składanie się białek: Wyjaśnienie modelu siatki HP

Dogmatyka Anfinsena i paradoks Leiwentala

Wczesnym określeńm 1960 roku Christian Anfinsen całkowicie denaturował rybosomę A i obserwował, że spontanicznie się ona odwróciła w aktywną formę po usunięciu denaturanta – ustanawiając dogmatę Anfinsena: struktura domyślna jest stanem termodynamicznie stabilnym, zdeterminowanym wyłącznie sekwencją aminową, bez potrzeby dodatkowego szablonu. Jednak pozostaje zagadka: 100-residołowy protein badający nawet 3 konformacje na resztę ma około 5×10⁴⁷ możliwych struktur – eksploracja wszystkich tych struktur w tempie 10¹³ na sekundę przekraczałaby czas istnienia wszechświata. To paradoks Leiwentala (1969): składanie nie może być losowym poszukiwaniem, więc sekwencja musi zniekształcać drogę do stanu domyślnego.

Model HP: Dwie litery, jedna zasada energii

Proponowany przez Kena Dill w 1985 roku model HP redukuje wszystkie 20 aminokwasów do dwóch typów – H (hydrofobowy, odpychany przez wodę) i P (polarny, komfortowo zachowujący się w wodzie) – umieszczonych na ścieżce unikającej siebie na siatce 2D lub 3D. Energia konformacji pochodzi z jednej zasady: każdy parzyste H-H kontakty na sąsiednich węzłach siatki, ale nie sąsiadujących w sekwencji, przyczynia się do -1; wszystko inne ma wartość zero. Celem jest konformacja o największej liczbie kontaktów H–H. Pomimo ekstremalnego uproszczenia, model ten odzwierciedla rzeczywiste zachowanie – sekwencje z zagęszczonymi aminokwasami H-foldują się kompaktowo, co przypomina prawdziwe hydrofobowe rdzenie białek, ponieważ pogrążanie grup niepolarowych uwalnia uporządkowane cząsteczki wody i zwiększa entropię rozpuszczalnika, wyjaśniając szacunkowo 60–70% rzeczywistej stabilności zgięcia.

E = −1 × (number of non-bonded H–H lattice contacts)

Self-avoiding walks on 2D square lattice (growth ~2.638^n):
  length 20: ~1.7×10^7        length 50: ~2.6×10^19

Energia przepływająca przez kanały i trudność obliczeniowa NP

Nowoczesna teoria składania (Wolynes, Bryngelson, Onuchic) opisuje krajobraz energetyczny: dla losowego heteropolimeru jest niegościnny, pełen pułapek kinetycznych, ale prawdziwe białka – ukształtowane przez ewolucję – mają kanałowy krajobraz, w którym kompaktowość zazwyczaj koreluje z niższą energią potencjalną, rozwiązując paradoks Leiwitala jako kontrolowany spadek zamiast losowego przeszukiwania. Znalezienie dokładnej konfiguracji konformacyjnej białka HP, choćby było możliwe, zostało udowodnione w 1998 roku jako NP-trudne (Crescenzi, Goldman, Papadimitriou, Piccoloboni, Yannakakis): samoizolujące się węzły na 2D siatce rosną w przybliżeniu jak 2,638ⁿ (≈10¹⁹ dla długości 50), więc dokładne przeszukiwanie jest niewykonalne poza krótkimi sekwencjami, a zamiast tego stosuje się heurystyki – Monte Carlo z kryterium Metropolis, algorytmy genetyczne i uczenie ze wzmocnieniem.

Od modeli na podstawie siatki po AlphaFold

Rozszerzenia takie jak HPNX (cztery typy aminowych reszt) i macierz empirycznych kontaktów Miyazawy-Jernigan 20×20 stopniowo zmniejszają dystans między abstrakcją HP a rzeczywistą fizyką białek. Ewolucja od modeli na podstawie siatki przez narzędzia do fragmentowej składarki (Rosetta, I-TASSER) aż po architekturę transformatora równoważnego AlphaFold2, która w 2021 roku osiągnęła dokładność zbliżoną do eksperymentalnej, ilustruje stopniowy postęp w obliczeniach biologicznych. AlphaFold2 nie wykorzystuje siatki, ale ta sama intuicja fizyczna – ukrywanie się hydrofobowych grup jako dominujący czynnik napędowy – jest implikowana w jego danych treningowych zawierających rozstrzygnięte struktury białkowe, co dokładnie robi to, co model HP czyni jawnym i obliczalnym.

Frequently asked questions

Wypróbuj na żywo

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