Калькулятор температуры плавления ДНК
Оцените температуру плавления (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^+}]\): каждое десятикратное увеличение соли добавляет к расчётной Tm около 16,6 °C. Сравнивая олигонуклеотиды, всегда делайте это при одинаковой концентрации соли и помните, что концентрацию нужно вводить в мМ — пересчёт в моль/л выполняется автоматически.
Точность, допущения и практические советы
- Обе формулы предполагают идеальный дуплекс: последовательность должна гибридизоваться с полностью комплементарной цепью без несовпадений, петель и разрывов. Даже одно несовпадение может снизить истинную Tm на несколько градусов.
- Ограничения правила Уоллеса: это грубое правило для последовательностей примерно из 4–13 оснований при концентрации около 50 мМ Na⁺. Оно не учитывает концентрацию соли и влияние длины, поэтому годится лишь для быстрой оценки очень коротких олигонуклеотидов.
- Ограничения формулы с поправкой на соль: она учитывает длину и одновалентную соль, но не концентрацию олигонуклеотида, концевые эффекты, денатуранты (формамид, ДМСО) и двухвалентные катионы вроде Mg²⁺, которые заметно стабилизируют дуплексы в ПЦР-буферах.
- Для точных расчётов используйте модель ближайших соседей: современные программы для дизайна праймеров применяют термодинамику ближайших соседей, учитывающую энергию стэкинга каждого динуклеотидного шага. Приведённые здесь две формулы — быстрые приближения с точностью до нескольких градусов; они хороши для планирования и перекрёстной проверки, но не для окончательной разработки методики.
- Интерпретация результата: как показывает практика, гибридизация и отжиг лучше всего идут примерно на 5 °C ниже Tm. Если расчётная Tm сильно отличается от нужной, меняйте длину праймера или содержание GC, а не только температуру отжига.
Если в тех же экспериментах вы работаете с пептидами или белками, возможно, вам пригодится наш калькулятор изоэлектрической точки для планирования условий электрофореза, а калькулятор буферной ёмкости поможет проверить, сохранит ли гибридизационный буфер свой pH при нагревании.
Часто задаваемые вопросы
Почему разные калькуляторы Tm дают разные значения для одной и той же последовательности?
Разные калькуляторы используют разные модели. Простые правила (Уоллеса, базовые формулы с поправкой на соль) не учитывают порядок оснований, тогда как модель ближайших соседей учитывает энергию стэкинга каждого динуклеотидного шага, поэтому её результаты обычно отличаются на несколько градусов. Поправки на Mg²⁺, концентрацию олигонуклеотида и ДМСО тоже сдвигают результат. Всегда сравнивайте значения Tm, полученные одним методом, — смешивать Tm из разных калькуляторов при дизайне праймеров является распространённой ошибкой.
Каким должно быть содержание GC в праймере или олигонуклеотиде?
Стремитесь к содержанию GC 40–60 %. Праймеры в этом диапазоне связываются достаточно прочно, чтобы быть специфичными, но при этом надёжно денатурируют и отжигаются. Участки, богатые GC (более 65 %), могут образовывать устойчивые вторичные структуры и требуют более высоких температур отжига, а AT-богатые праймеры с содержанием GC ниже 40 % связываются слабо или неспецифично. Также избегайте длинных повторов одного основания (например, GGGGG) и 3′-конца с большим числом G или C — это провоцирует неспецифичный прайминг.
Как выбрать температуру отжига для ПЦР на основе Tm?
Обычно начинают с температуры отжига примерно на 3–5 °C ниже меньшей Tm пары праймеров. При появлении неспецифичных полос повышайте температуру отжига шагами по 1–2 °C; если продукта нет — немного снизьте её. Современные протоколы часто используют градиентную ПЦР или двухступенчатую программу touchdown, чтобы эмпирически найти оптимум, поскольку простые формулы Tm лишь приближённо описывают поведение в ПЦР-буфере.
Можно ли использовать этот калькулятор для последовательностей РНК или последовательностей с U?
Нет. Калькулятор принимает только основания ДНК — A, C, G и T; последовательности РНК с урацилом (U) отклоняются. Кроме того, гибриды РНК/ДНК и РНК/РНК устойчивее дуплексов ДНК/ДНК, поэтому даже при замене U на T правила Уоллеса и формула с поправкой на соль занижали бы их Tm. Для олигонуклеотидов РНК используйте калькулятор на основе модели ближайших соседей с термодинамическими параметрами для РНК.
Как Mg²⁺ и ДМСО влияют на расчётную Tm?
Обе формулы здесь предполагают только одновалентную соль, поэтому в типичных ПЦР-буферах они занижают Tm. Mg²⁺ сильно стабилизирует дуплексы, часто повышая эффективную Tm на несколько градусов, а денатуранты вроде ДМСО или формамида её снижают (примерно на 0,6 °C на каждый 1 % ДМСО). Если в буфере есть Mg²⁺, считайте результат консервативной оценкой и оптимизируйте условия экспериментально либо используйте инструмент на основе модели ближайших соседей, допускающий ввод концентраций двухвалентных катионов.
Что происходит с Tm, если олигонуклеотид имеет несовпадение с мишенью?
Калькулятор предполагает полностью комплементарный дуплекс и не может учесть несовпадения. Даже одно несовпадение, особенно в середине последовательности или вблизи 3′-конца праймера, способно снизить истинную Tm на 5–15 °C и ухудшить эффективность амплификации. Для вырожденных праймеров или зондов к вариабельным участкам используйте инструмент на основе модели ближайших соседей, который отдельно оценивает каждый тип несовпадения, и проверяйте условия экспериментально.
В каких единицах вводить концентрацию соли — в мМ или моль/л, и как её перевести?
Введите концентрацию Na⁺ в мМ; калькулятор сам переведёт её в моль/л для члена \(16.6\,\log_{10}[\mathrm{Na^+}]\). Для пересчёта вручную разделите мМ на 1000: 50 мМ = 0,05 моль/л, 500 мМ = 0,5 моль/л. Учтите, что это суммарная концентрация одновалентных катионов: в некоторых протоколах указывается ионная сила или включаются Tris и K⁺, которые стоит прибавить к значению Na⁺ для реалистичной оценки.