Згортання білків: Модель HP-латиної решітки

Магнітне поле симуляції - Закон Гаусса

На початку 1960-х років Хр. Анфензен повністю денатурував рибонуклеїнову кислоту А та спостерігав, як вона самовільно згортається у свою активну форму після видалення денатуруючого агента — встановлюючи догму Анфензена: природна структура є термодинамічно стабільною станом, що визначається лише амінокислотним послідовством, без додаткової шаблону. Але залишається загадком: 100-residюбна білок з вибіркою навіть 3 конформацій на залишку має приблизно $5 imes 10^{47}$ можливих структур — дослідження всіх цих структур при швидкості $10^{13}$ за секунду перетнуло б існування Всесвіту. Це парадокс Левінштейна (1969): згортання не може бути випадковим пошуком, тому послідовність повинна схиляти шлях до природної форми.

Закон электромагнітної індукції Фарадея та закон Ленца

Пропонований Кеном Діллом у 1985 році, модель HP (H-P) зменшує всі 20 амінокислот до двох типів — H (гідрофобна, відштовхується водою) та P (полярна, комфортна у воді), розміщених на самонесучому шляху через 2D або 3D сітку. Енергія конформації походить від одного правила: кожний пар гідрофобних (H) атомів, що знаходяться сусідніми на сітці, але не є послідовними, сприяє -1; все інше сприяє нулю. Мета — конформація з найбільшою кількістю H-H контактів. Незважаючи на надмірне спрощення, це відтворює реальну поведінку — послідовності з групуванням H атомів складаються компактно, що відображає реальні гідрофобні ядра білків, оскільки захоплення неполярних груп вивільняє впорядковані молекули води та збільшує сольовий ентропій, що пояснює приблизно 60–70% реальної стабільності згортання.

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

Енергетичні каналі та NP-складність

Сучасна теорія згортки (Воленес, Брінгєльсон, Онухік) описує енергетичний ландшафт: для випадкового гетерополімеру він буває нерівним та сповнений кінетичних пасток, але справжні білки — сформовані еволюцією — мають «трубкоподібний» ландшафт, де компактність зазвичай корелює з нижчою вільною енергією, вирішуючи парадокс Левінталя як спрямоване спускання, а не випадковий пошук. Однак, точне визначення конформації основної станів білка було доведено NP-складним у 1998 році (Кресцензі, Голдман, Пападімітроу, Піколбоні, Янакакіс): самовільні переміщення на 2D решітці зростають приблизно як 2.638ⁿ (≈10¹⁹ для довжини 50), тому точний пошук неможливий за межами коротких послідовностей, і використовуються евристичні методи — Монте-Карло з критерієм Метаполіса, генетичні алгоритми та навчання з підкріпленням.

Магнітне поле симуляції - Закон Гаусса

Розширення, такі як HPNX (чотири типи амінокислот) та матриця емпіричних контактів Miyazawa-Jernigan 20x20, поступово зменшують відстань між абстракцією HP і фізикою реальних білків. Шлях від моделей решіток через інструменти з'єднання фрагментів (Rosetta, I-TASSER) до архітектуриequivariant трансформаторів AlphaFold2, яка досягла майже експериментальної точності у 2021 році, демонструє багатошаровий прогрес у обчислювальній біології. AlphaFold2 не використовує решітку, але та сама фізична інтуїція — гідрофобне занурення як основна рушійна сила — явно закодована в його навчальних даних розв'язаних структур, що є тим, що робить модель HP вираженою та обчислюваною.

Frequently asked questions

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте the simulation і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.