Калькулятор температури плавлення ДНК
Оцініть температуру плавлення (Tm) олігонуклеотиду ДНК за його послідовністю та концентрацією солі. Введіть один ланцюг (A, C, G, T) і концентрацію Na⁺, щоб отримати Tm та вміст GC.
Оцінка температури плавлення ДНК-олігонуклеотиду
Формула температури плавлення ДНК
Де:
- \(A, T, G, C\) — кількість основ аденіну, тиміну, гуаніну та цитозину в послідовності
- \(n\) — загальна кількість основ (довжина послідовності)
- \([\mathrm{Na^+}]\) — сумарна концентрація одновалентних солей у моль/л
Перша формула — правило Воллеса, яке застосовують для коротких олігонуклеотидів до 13 основ; друга — формула з урахуванням концентрації солі (Howley et al., 1979), яку автоматично застосовують від 14 основ. Калькулятор сам обирає відповідну формулу.
Температура плавлення ДНК — приклад обчислення
Приклад: олігонуклеотид із 15 основ із послідовністю ATGCATGCGATTACA в буфері, що містить 50 мМ Na⁺. У послідовності 5 основ A, 4 T, 3 G і 3 C, отже \(G + C = 6\), \(A + T = 9\) і \(n = 15\).
\[ [\mathrm{Na^+}] = 50\ \mathrm{mM} = 0.05\ \mathrm{mol/L}, \quad \log_{10}(0.05) = -1.301 \]
\[ T_m = 100.5 + 41 \cdot \frac{6}{15} - \frac{820}{15} + 16.6 \cdot (-1.301) \]
\[ T_m = 100.5 + 16.4 - 54.67 - 21.60 = 40.6\ \mathrm{°C} \]
Результат: температура плавлення становить \(T_m \approx 40.6\) °C, а вміст GC — 40,0%. Просте правило Воллеса дало б 2 × 9 + 4 × 6 = 42 °C, але для послідовностей від 14 основ і довших його не рекомендують.
Як користуватися калькулятором
- Введіть або вставте послідовність у поле Послідовність ДНК (5′→3′). Приймаються лише літери A, C, G і T; регістр, пробіли та переноси рядків ігноруються.
- Вкажіть сумарну концентрацію одновалентних солей у полі Концентрація солі (Na⁺), у мМ (за замовчуванням 50 мМ). Це значення використовується, лише коли послідовність містить 14 основ або більше.
- Зчитайте температуру плавлення (Tm) у °C та вміст GC у відсотках на панелі результатів.
- Зазирніть у Деталі обчислення, щоб побачити склад основ, застосовану формулу та підставлені числа — так ви зможете перевірити результат вручну.
Що таке температура плавлення?
Температура плавлення, або Tm, — це температура, за якої половина дуплексів ДНК у розчині дисоціює на окремі ланцюги, а половина залишається спареною. Це прямий показник стабільності дуплексу: чим вища Tm, тим більше тепла потрібно, щоб розділити два ланцюги. Tm має значення всюди, де два комплементарні ланцюги ДНК мають надійно знайтися й зв’язатися один з одним — при проєктуванні праймерів для ПЛР, підборі умов гібридизації зондів або плануванні температур відпалу для коротких оліго.
Стабільність дуплексу визначають самі пари основ. Пара G–C тримається трьома водневими зв’язками, тоді як пара A–T — лише двома, тому послідовності, багаті на GC, плавляться за вищих температур. У правилі Воллеса це виражено дуже просто: кожна пара A–T додає до Tm приблизно 2 °C, а кожна пара G–C — приблизно 4 °C. Вміст GC, що виводиться разом з Tm, показує склад послідовності — це зручна перевірка на адекватність, адже вміст GC і Tm завжди змінюються в один бік.
Важливу роль відіграє й сіль. Іони натрію екранують негативні заряди фосфатного остова, тому вища концентрація Na⁺ стабілізує дуплекс і підвищує Tm. Саме тому формула з урахуванням солі містить доданок \(16.6\,\log_{10}[\mathrm{Na^+}]\): кожне десятикратне збільшення концентрації солі додає приблизно 16,6 °C до розрахункової Tm. Порівнюючи оліго, завжди робіть це за однакової концентрації солі й пам’ятайте, що концентрацію слід вводити в мМ — перетворення на моль/л виконується автоматично.
Точність, припущення та практичні поради
- Обидві формули припускають ідеальний дуплекс: послідовність має гібридизуватися з повністю комплементарним ланцюгом без незбігів, петель і прогалин. Навіть одна некомплементарна пара може знизити справжню Tm на кілька градусів.
- Обмеження правила Воллеса: це приблизне емпіричне правило, розраховане на 4–13 основ за концентрації близько 50 мМ Na⁺. Воно не враховує ані концентрацію солі, ані вплив довжини послідовності, тому використовуйте його лише для швидких оцінок дуже коротких оліго.
- Обмеження формули з урахуванням солі: вона враховує довжину й одновалентну сіль, але не концентрацію оліго, кінцеві ефекти, денатуранти (формамід, DMSO) та двовалентні катіони, як-от Mg²⁺, які суттєво стабілізують дуплекси в ПЛР-буферах.
- Для точної роботи використовуйте модель найближчих сусідів: сучасні програми для проєктування праймерів спираються на термодинаміку найближчих сусідів, яка враховує енергію стекінгу кожного динуклеотидного кроку. Обидві формули тут — швидкі наближення з точністю до кількох градусів: вони годяться для планування та перехресних перевірок, але не для остаточного проєктування аналізів.
- Тлумачення результату: як правило, гібридизація або відпал найкраще відбувається приблизно на 5 °C нижче за Tm. Якщо розрахована Tm далеко від цільового значення, змінюйте довжину праймера або вміст GC, а не лише температуру відпалу.
Якщо в тих самих експериментах ви працюєте з пептидами чи білками, вам може знадобитися наш Калькулятор ізоелектричної точки для планування умов електрофорезу, а Калькулятор буферної ємності допоможе перевірити, чи витримує ваш гібридизаційний буфер задане pH під час нагрівання.
Часті запитання
Чому різні калькулятори Tm дають різні значення для однієї й тієї самої послідовності?
Різні калькулятори використовують різні моделі. Прості правила (Воллеса, базові формули з урахуванням солі) ігнорують порядок основ у послідовності, тоді як модель найближчих сусідів враховує енергію стекінгу кожного динуклеотидного кроку, тому зазвичай дає значення, що відрізняються на кілька градусів. На результат також впливають поправки на Mg²⁺, концентрацію оліго та DMSO. Завжди порівнюйте значення Tm, отримані одним методом, — змішувати Tm з різних калькуляторів під час проєктування праймерів є типовою помилкою.
Який вміст GC вважається добрим для праймера або олігонуклеотиду?
Орієнтуйтеся на вміст GC 40–60%. Праймери в цьому діапазоні зв’язуються достатньо міцно, щоб бути специфічними, і водночас надійно денатурують і відпалюються. Ділянки, багаті на GC (понад 65%), можуть утворювати стабільні вторинні структури та потребують вищих температур відпалу, а багаті на AT праймери (менше 40%) можуть зв’язуватися слабо або неспецифічно. Також уникайте довгих серій однієї основи (наприклад, GGGGG) та 3′-кінця з великою кількістю основ G або C, що провокує неспецифічне зв’язування праймера.
Як обрати температуру відпалу для ПЛР на основі Tm?
Типова вихідна точка — температура відпалу приблизно на 3–5 °C нижча за меншу з Tm двох праймерів. Якщо з’являються неспецифічні смуги, підвищуйте температуру відпалу кроками по 1–2 °C; якщо продукту немає — трохи знижте її. Сучасні протоколи часто застосовують градієнтну ПЛР або двоступеневу програму touchdown, щоб знайти оптимум емпірично, адже прості формули Tm лише наближено описують справжню поведінку в ПЛР-буфері.
Чи можна використовувати цей калькулятор для РНК-послідовностей або послідовностей, що містять U?
Ні. Калькулятор приймає лише основи ДНК — A, C, G і T; РНК-послідовності з урацилом (U) відхляються. До того ж гібриди РНК/ДНК і РНК/РНК стабільніші за дуплекси ДНК/ДНК, тому правило Воллеса та формула з урахуванням солі занизили б їхню Tm навіть після заміни U на T. Для РНК-оліго скористайтеся калькулятором на основі моделі найближчих сусідів із РНК-термодинамічними параметрами.
Як Mg²⁺ і DMSO впливають на розраховану Tm?
Обидві формули тут припускають лише одновалентну сіль, тому в типових ПЛР-буферах вони занижують Tm. Mg²⁺ суттєво стабілізує дуплекси, часто підвищуючи ефективну Tm на кілька градусів, а денатуранти, як-от DMSO або формамід, знижують її (приблизно на 0,6 °C на 1% DMSO). Якщо ваш буфер містить Mg²⁺, розглядайте результат як консервативну оцінку та оптимізуйте умови експериментально або скористайтеся інструментом на основі моделі найближчих сусідів, який приймає концентрації двовалентних катіонів.
Що станеться з Tm, якщо мій оліго не повністю комплементарний мішені?
Калькулятор припускає повністю комплементарний дуплекс, тому не може врахувати незбіги. Навіть одна некомплементарна пара, особливо посередині послідовності або біля 3′-кінця праймера, може знизити справжню Tm на 5–15 °C і погіршити ефективність ампліфікації. Для вироджених праймерів або зондів до мінливих ділянок використовуйте інструмент на основі моделі найближчих сусідів, який окремо оцінює кожен тип незбігу, і перевіряйте умови експериментально.
У яких одиницях вводити концентрацію солі — мМ чи моль/л — і як перетворювати?
Вводьте концентрацію Na⁺ у мМ; калькулятор сам перетворює її на моль/л для логарифмічного доданка \(16.6\,\log_{10}[\mathrm{Na^+}]\). Щоб перетворити вручну, поділіть мМ на 1000: 50 мМ = 0,05 моль/л, 500 мМ = 0,5 моль/л. Пам’ятайте, що це сумарна концентрація одновалентних катіонів — у деяких протоколах наводять іонну силу або враховують Tris і K⁺, які варто додати до значення Na⁺ для реалістичної оцінки.