🔬 Візуалізація згортання білка

HP-модель (гідрофобно-полярна) — це спрощена ґраткова модель, введена Лау і Діллом у 1989 році. Вона зводить кожну амінокислоту до одного з двох типів: гідрофобний (H) або полярний (P), розміщує їх на квадратній ґратці як самоуникний шлях і визначає енергію як від'ємну кількість незв'язаних H-H контактів. Попри свою простоту, вона правильно передбачає, що гідрофобні залишки збираються у прихованому ядрі — та сама рушійна сила, що керує реальним згортанням білків.

Білки згортаються тому, що приховування гідрофобних залишків від води термодинамічно вигідне. Молекули води утворюють впорядковані водневозв'язані клітки навколо відкритих неполярних груп, що коштує ентропії. Приховування цих груп у компактному ядрі вивільняє молекули води для більш невпорядкованих конфігурацій, знижуючи загальну вільну енергію системи. Цей гідрофобний ефект є домінантною фізичною силою, що керує згортанням, доповнюваною водневими зв'язками, електростатикою та ван-дер-ваальсовими взаємодіями.

Кожна пара гідрофобних залишків, що є сусідами на ґратці, але не послідовними в ланцюзі, дає внесок мінус одну одиницю енергії. Загальна енергія тому є від'ємним цілим числом; чим вона від'ємніша, тим стабільніша конформація. Симуляція прагне мінімізувати цю енергію, що відповідає максимізації кількості H-H контактів і, отже, розміру гідрофобного ядра. Графік енергії в реальному часі дозволяє спостерігати за прогресом пошуку.

На кожному кроці пропонується випадковий рух: або кінцевий рух (один кінець переходить на нову сусідню позицію ґратки), або поворотна ротація (один сегмент ланцюга повертається на 90 градусів навколо випадково обраного залишку-осі). Якщо результуюча конформація перетинає сама себе, вона одразу відхиляється. Інакше обчислюється зміна енергії ΔE. Якщо ΔE від'ємна (нижча енергія), рух приймається безумовно; якщо додатна — приймається з імовірністю e^(−ΔE/kT). Цей критерій Метрополіса гарантує, що ланцюг вибирає конформації відповідно до розподілу Больцмана при температурі kT.

Безрозмірна температура kT задає масштаб теплових флуктуацій. При високій kT (до 5) фактор Больцмана e^(−ΔE/kT) залишається близьким до 1 навіть при великому зростанні енергії, тож ланцюг широко досліджує конформаційний простір і рідко застрягає. При низькій kT (до 0.1) практично приймаються лише рухи, що знижують енергію, тож ланцюг різко спускається у найближчий локальний мінімум. Кнопка «Знайти оптимум» примусово встановлює kT = 0.1 на 8000 кроків, щоб провести пошук глибоко в найдоступнішому басейні.

Парадокс Левінталя, сформульований Сайрусом Левінталем у 1969 році, вказує на те, що білок з N залишків має порядку 3^N можливих конформацій. Для скромного білка зі 100 залишків це приблизно 10^47 станів. Якби ланцюг перебирав їх зі швидкістю 10^6 за секунду, пошук мінімально-енергетичного згортання методом грубої сили зайняв би значно довше за вік Всесвіту. Проте реальні білки надійно згортаються за мікросекунди–мілісекунди. Розгадка в тому, що згортання — не випадковий пошук: енергетичний ландшафт має форму лійки з багатьма спадними шляхами, що ведуть ланцюг до нативного стану.

Обчислювальна задача пошуку конформації з мінімальною енергією в HP-моделі на двовимірній квадратній ґратці була доведена NP-повною Хартом і Істрейлом у 1996 році. Це означає, що не існує відомого алгоритму, який гарантує глобальний оптимум за поліноміальний час зі зростанням довжини послідовності. На практиці для пошуку хороших, хоча не обов'язково оптимальних, згортань використовуються евристики: Монте-Карло, імітаційний відпал, генетичні алгоритми. Кнопка «Знайти оптимум» у цій симуляції запускає жадібний пошук при низькій температурі, а не вичерпний.

Теорія енергетичного ландшафту, розроблена насамперед Волінесом, Онучичем і Діллом у 1990-х, описує згортання як рух по багатовимірній енергетичній поверхні над конформаційним простором. Для послідовностей, що надійно згортаються, ландшафт має форму лійки: конформації на краю лійки мають високу енергію та невпорядковані, конформації біля дна — низьку енергію та компактні. Лійка забезпечує термодинамічне зміщення, що послідовно скеровує ланцюг вниз, пояснюючи, чому згортання швидке та відтворюване, попри величезну кількість можливих станів.

Модель навмисно спрощена. Реальні білки використовують 20 типів амінокислот, а не 2. Вони згортаються у тривимірному просторі, а не на плоскій ґратці. Важливі взаємодії, такі як водневі зв'язки, жорсткість кістяка, геометрія бічних ланцюгів, електростатика та сольватація, повністю відсутні. HP-модель також не може розрізняти різні гідрофобні залишки між собою. Тим не менш, вона передає основну фізику гідрофобного колапсу та лійкоподібного ландшафту, роблячи її цінним концептуальним інструментом, хоча і непридатним для кількісних прогнозів щодо реальних білків.

N позначає N-кінець — кінець поліпептидного ланцюга, що починається з аміногрупи (-NH2), а C позначає C-кінець з карбоксильною групою (-COOH). У реальних білках послідовність амінокислот завжди зчитується від N до C, що також є напрямком біосинтезу на рибосомі. У HP-моделі ці позначки просто відзначають два кінці самоуникного шляху, щоб можна було відстежувати орієнтацію ланцюга під час зміни конформації в процесі згортання.