Генні регуляторні мережі визначають, коли, де і в якій кількості кожен ген виробляє білок, керуючи всім — від диференціації клітин під час розвитку організму до добових циркадних ритмів. Транскрипційні фактори — це ключові білки, які зв'язуються з ДНК у специфічних ділянках промоторів і або активують, або пригнічують транскрипцію генів-мішеней. Один транскрипційний фактор може регулювати сотні генів, а гени, своєю чергою, регулюють транскрипційні фактори, утворюючи петлі зворотного зв'язку дивовижної складності й точності.
Найпростіші регуляторні мотиви — підграфи мережі з чітко визначеними зв'язками входу й виходу — включають петлю авторегуляції (транскрипційний фактор пригнічує або активує власний ген), тригер (два гени, що взаємно пригнічують один одного і створюють бістабільність — два стійкі стани) та репресилятор (три гени, що по колу пригнічують один одного, породжуючи стійкі коливання, аналогічні біологічному годиннику). Ці мінімальні схеми регулярно трапляються в організмів від бактерій до ссавців, що наводить на думку, що еволюція знайшла оптимальні архітектури мереж.
Цей симулятор моделює динаміку невеликих генних мереж за допомогою звичайних диференціальних рівнянь для концентрацій мРНК і білків, включно з членами кооперативності у формі функцій Хілла. Ви можете побудувати тригери, репресилятори та схеми прямого зв'язку, налаштувати швидкості продукції й деградації та коефіцієнти кооперативності, і поспостерігати, як початкові умови визначають, якого стійкого стану буде досягнуто (бістабільність), або побачити, як у репресилятора виникають коливання, коли кооперативність достатньо висока — це демонструє ключові принципи синтетичної біології та генних схем розвитку.
Що таке транскрипційний фактор і як він контролює експресію генів?
Транскрипційний фактор (ТФ) — це білок, який зв'язується зі специфічними ділянками ДНК (промоторами чи енхансерами) поблизу гена-мішені й впливає на те, чи транскрибуватиме цей ген РНК-полімераза. Активуючі ТФ залучають полімеразу та загальний апарат транскрипції, збільшуючи вироблення мРНК. Пригнічувальні ТФ блокують доступ полімерази або залучають комплекси, що ущільнюють хроматин і вимикають ген. Оскільки самі ТФ кодуються генами, вони можуть регулювати один одного, утворюючи складні регуляторні мережі. Більшість процесів розвитку організму керується комбінаторними кодами ТФ: різні поєднання транскрипційних факторів визначають тип клітини та програму експресії генів.
Що таке бістабільність у генних мережах і чим вона зумовлена?
Бістабільність означає, що генна мережа має два стійкі стаціонарні стани і може перемикатися між ними залежно від початкових умов або короткочасних сигналів. Генетичний тригер — два гени A і B, що взаємно пригнічують один одного — є бістабільним: якщо рівень A високий, він пригнічує B (утримуючи B на низькому рівні), що дає змогу A залишатися високим; і навпаки, якщо високий B, він пригнічує A. Обидва стани самопідтримувані. Бістабільність вимагає позитивного зворотного зв'язку (взаємної репресії або самоактивації) і зазвичай потребує кооперативності зв'язування (коефіцієнт Хілла n > 1). Рішення щодо долі клітини під час розвитку часто використовують бістабільність для надійного й незворотного вибору одного з двох шляхів диференціації.
Що таке репресилятор і як він створює коливання?
Репресилятор, розроблений Елловіцем і Лейблером (2000 рік), — це синтетична мережа з трьох генів, де ген A пригнічує ген B, ген B пригнічує ген C, а ген C пригнічує ген A — кільце репресорів. Коли рівень A високий, він пригнічує B; низький рівень B дає змогу C зростати; високий C пригнічує A; низький A дозволяє B зростати; високий B пригнічує C — і так по колу, породжуючи стійкі коливання концентрацій білків із періодом, який визначається сталими часу продукції й деградації. Це була перша штучно створена синтетична генна схема, що продемонструвала, як проста топологія мережі породжує коливання, і вона суттєво поглибила наше розуміння біологічних годинників.
Функція Хілла H(x) = xⁿ/(Kⁿ + xⁿ) моделює сигмоподібну реакцію транскрипції гена на концентрацію транскрипційного фактора x, де K — концентрація, за якої активація сягає половини максимуму, а n — коефіцієнт Хілла, що вимірює кооперативність. При n=1 (без кооперативності) реакція має вигляд гіперболи — поступова. При n=2 (наприклад, коли ТФ зв'язується у вигляді димера) реакція стає крутішою. При n → ∞ реакція перетворюється на ступінчасту функцію (перемикач). Висока кооперативність (велике n) є необхідною умовою бістабільності й коливань у генних схемах і досягається біологічно через олігомеризацію ТФ або наявність кількох ділянок зв'язування ТФ у промоторі.
Синтетична біологія застосовує інженерні принципи для створення нових біологічних функцій шляхом складання стандартизованих генетичних деталей (промоторів, сайтів зв'язування рибосом, кодувальних послідовностей, термінаторів). Моделі генної регуляції допомагають передбачити, чи працюватиме спроєктована схема так, як задумано, визначити діапазони параметрів для бістабільності чи коливань і оптимізувати часові характеристики та амплітуду. Синтетичні схеми вже розроблено для виявлення онкомаркерів і знищення пухлинних клітин, вироблення цінних біохімічних речовин (інсулін, артемізинін), створення біосенсорів для екологічних токсинів і побудови біологічних комп'ютерів. Математичне моделювання генної регуляції є критично важливим, оскільки біологічні елементи взаємодіють складними нелінійними способами, які неможливо передбачити самою лише інтуїцією.
Зі зростанням n від 1 до 6 функція Хілла H(x) = xⁿ/(Kⁿ+xⁿ) дедалі більше нагадує різкий перемикач замість поступової кривої. На пресеті взаємної репресії «Тригер-перемикач» вище n робить бістабільність надійнішою: два стійкі стани (ген A високий/B низький, або B високий/A низький) стають чіткіше розділеними, і невеликі збурення рідше перемикають систему, тоді як при низькому n мережа може осісти в проміжному, менш чітко визначеному стані.
β (від 1 до 40) задає максимальну швидкість синтезу нового білка, а γ (від 0,2 до 4) — швидкість його деградації чи розведення. За відсутності репресії ген встановлюється на стаціонарну концентрацію приблизно β/γ, тож підвищення β або зниження γ обидва підіймають рівноважний рівень. Відношення β/γ до Порогу K також визначає, чи зможе ген повністю увімкнутись: якщо β/γ набагато більше за K, ген досягає сильної активації після перевищення порогу.
K (від 0,5 до 20) — це концентрація транскрипційного фактора, за якої функція Хілла сягає половини максимального значення; вона зсуває всю сигмоподібну криву реакції ліворуч чи праворуч уздовж осі концентрації, не змінюючи її крутизну (за це відповідає n окремо). Низьке K означає, що ген вмикається при низьких концентраціях ТФ, а високе K вимагає накопичення набагато більшої кількості ТФ, перш ніж ген-мішень відреагує значно.
Ні — Швидкість (від 0,5× до 8×) лише масштабує, наскільки швидко внутрішній годинник симуляції рухається відносно реального часу на екрані, суто для зручності перегляду. Справжній період коливань на пресеті «Репресилятор» визначається сталими часу продукції й деградації (β та γ), що діють на кільце з трьох генів-репресорів, а не повзунком Швидкість, тож зміна Швидкості лише прискорює чи сповільнює відтворення тих самих коливань, не змінюючи їхній базовий період у симуляційних одиницях часу.
У когерентній петлі прямого зв'язку регулятор активує ген-мішень одразу двома шляхами — напряму й опосередковано через проміжний ген, — і мішень вмикається лише коли обидва шляхи узгоджуються. Короткий вхідний імпульс може активувати прямий шлях, але зникає раніше, ніж повільніший проміжний шлях накопичить достатньо продукту, щоб теж перевищити свій поріг, тож ген-мішень відфільтровує короткочасний шум; лише сигнал, що триває достатньо довго для узгодження обох шляхів, вмикає ген-мішень.