Как подобрать праймер для ПЦР: расчёт Tm и отжига

Дано: прямой праймер 5'-ATGCCTGACGGTAACCTTAG-3' и обратный 5'-TGAGGCAGCCAATCCTTGAC-3', оба длиной 20 нуклеотидов. Найти: долю ГЦ, температуру плавления $T_m$ каждого праймера и рабочую температуру отжига пары.

Порядок действий короткий: считаем нуклеотидный состав, из него долю ГЦ, затем $T_m$ по правилу Уоллеса и берём температуру отжига на 5 градусов ниже меньшей из двух $T_m$. Ответ: доля ГЦ 50 и 55 %, $T_m$ = 60 и 62 °C, температура отжига $T_a$ = 55 °C; пара согласована, димеров и шпилек по 3'-концам нет. Калькулятор сверху считает те же величины для любой введённой последовательности.

Решение по шагам

Шаг 1. Нуклеотидный состав. Подсчитываем, сколько каких оснований в каждом праймере. Шаг механический, но именно здесь чаще всего теряют балл: одна пропущенная буква сдвигает ответ на 2 или 4 градуса.

Контроль: сумма четырёх чисел в каждой строке обязана равняться длине праймера. У прямого $5+5+5+5=20$, у обратного $5+4+5+6=20$ - состав выписан верно.

Шаг 2. Доля ГЦ. Доля гуанина и цитозина считается от полной длины:

$$\text{GC} = \frac{G + C}{N} \cdot 100\%.$$

Для прямого праймера $\text{GC} = \dfrac{5+5}{20} \cdot 100\% = 50\%$, для обратного $\text{GC} = \dfrac{5+6}{20} \cdot 100\% = 55\%$. Оба значения попадают в рабочее окно 40-60 %, то есть праймеры не переобогащены ни ГЦ-парами, ни АТ-парами.

Шаг 3. Температура плавления по правилу Уоллеса. Для олигонуклеотидов длиной 14-20 нт действует простое эмпирическое правило: каждая пара А-Т даёт 2 градуса, каждая пара Г-Ц - 4 градуса.

$$T_m = 2 \cdot (A + T) + 4 \cdot (G + C).$$

Подставляем числа из таблицы:

$$\begin{aligned} T_m^{\text{пр}} &= 2 \cdot (5 + 5) + 4 \cdot (5 + 5) = 20 + 40 = 60\ ^\circ\text{C}, \\ T_m^{\text{обр}} &= 2 \cdot (5 + 4) + 4 \cdot (5 + 6) = 18 + 44 = 62\ ^\circ\text{C}. \end{aligned}$$

Шаг 4. Согласованность пары. Праймеры работают в одной пробирке при одной температуре, поэтому их $T_m$ должны быть близки:

$$\Delta T_m = |62 - 60| = 2\ ^\circ\text{C} \le 3\ ^\circ\text{C}.$$

Расхождение укладывается в допуск. Будь оно 6-8 градусов, один праймер при рабочей температуре уже отваливался бы, а второй ещё сидел, и реакция пошла бы почти линейно, а не экспоненциально.

Шаг 5. Температура отжига. Рабочую температуру берут ниже меньшей из двух $T_m$, обычно на 3-5 градусов:

$$T_a = \min\left(T_m^{\text{пр}},\, T_m^{\text{обр}}\right) - 5\ ^\circ\text{C} = 60 - 5 = 55\ ^\circ\text{C}.$$

Меньшая, а не средняя: запас считают по слабому праймеру, иначе он просто не сядет на матрицу.

Шаг 6. Проверка 3'-концов. Последний нуклеотид у обоих праймеров - Г или Ц (у прямого G, у обратного C), это и есть ГЦ-замок: пара с тремя водородными связями фиксирует конец, с которого полимераза начинает синтез. В последних пяти звеньях по два ГЦ, то есть замок есть, но конец не переутяжелён.

Кросс-димер проверяем по четырём концевым нуклеотидам: у прямого TTAG, у обратного TGAC. Комплемент к TTAG - это CTAA, совпадения нет. Шпилек тоже нет: ни один четырёхбуквенный фрагмент внутри праймера не комплементарен другому фрагменту того же праймера.

Ответ. Прямой праймер: 20 нт, ГЦ 50 %, $T_m = 60\ ^\circ\text{C}$. Обратный: 20 нт, ГЦ 55 %, $T_m = 62\ ^\circ\text{C}$. Пара согласована ($\Delta T_m = 2\ ^\circ\text{C}$), рабочая температура отжига $T_a = 55\ ^\circ\text{C}$.

Формулы Tm: правило Уоллеса и поправка на соль

Правило Уоллеса - грубая сумма вкладов водородных связей. Пара А-Т держится двумя связями, пара Г-Ц - тремя, поэтому ГЦ-богатый дуплекс расходится при более высокой температуре. Коэффициенты 2 и 4 подобраны эмпирически для гибридизации коротких зондов в солевом буфере, и область применимости узкая: 14-20 нуклеотидов. За её пределами правило быстро начинает завышать ответ.

Вторая формула учитывает концентрацию одновалентных катионов и длину:

$$T_m = 81{,}5 + 16{,}6 \lg\left[\mathrm{Na^+}\right] + 0{,}41 \cdot \text{GC}\% - \frac{600}{N}.$$

При стандартных 50 мМ одновалентной соли $\lg 0{,}05 = -1{,}30$, и для наших праймеров получается:

$$\begin{aligned} T_m^{\text{пр}} &= 81{,}5 - 21{,}6 + 0{,}41 \cdot 50 - 30 = 50{,}4\ ^\circ\text{C}, \\ T_m^{\text{обр}} &= 81{,}5 - 21{,}6 + 0{,}41 \cdot 55 - 30 = 52{,}5\ ^\circ\text{C}. \end{aligned}$$

Разрыв почти в 10 градусов пугает, но он закономерен. Правило Уоллеса выведено для высокой ионной силы и для 20-мера слегка завышает; формула с логарифмом соли выведена для длинных зондов (больше 50 нт) и, наоборот, занижает - вклад члена $600/N$ при малых $N$ становится непропорционально большим. Точный расчёт делают по модели ближайших соседей, которая учитывает не отдельные основания, а соседние пары; для нашего прямого праймера она даёт около 57 градусов, то есть ровно между двумя оценками.

Практический вывод: формула на бумаге - стартовая точка, а не приговор. Реальную $T_a$ подбирают градиентной ПЦР, прогоняя одну смесь в диапазоне 52-60 градусов и выбирая ту температуру, где в геле одна чёткая полоса. Как выглядит правильный и неправильный результат в геле, разобрано в статье про электрофорез ДНК.

Проверка на димеры, шпильки и ГЦ-замок

Температура - не единственное требование к праймеру. Три структурных дефекта губят реакцию даже при идеальной $T_m$.

Самодимер и кросс-димер. Если 3'-конец праймера комплементарен участку того же или второго праймера, две молекулы слипаются и полимераза достраивает их друг на друге. Короткий побочный продукт размножается быстрее целевого и виден в геле как жирная полоса ниже 50 пар нуклеотидов. Проверка вручную: выписать последние 4-5 нуклеотидов, построить к ним комплемент в обратном порядке и поискать этот фрагмент в обеих последовательностях.

Шпилька. Праймер сворачивается сам на себя, если внутри него есть инвертированный повтор: фрагмент и его обратный комплемент, разделённые петлёй. Стебель из 4 пар с петлёй в 8-10 нуклеотидов ещё терпим, стебель из 5-6 пар уже опасен. Ищут так же: перебирают четвёрки букв и смотрят, не встречается ли дальше их обратный комплемент.

ГЦ-замок. Последний нуклеотид лучше делать Г или Ц: именно с 3'-конца полимераза начинает синтез, и прочная тройная связь снижает долю ложных стартов. Но больше трёх ГЦ в последних пяти звеньях брать не стоит - слишком липкий конец садится и на неполностью комплементарные места.

Отдельно следят за гомополимерными участками: четыре и больше одинаковых букв подряд (например, GGGG) дают проскальзывание полимеразы. В наших праймерах максимальная серия - две одинаковые буквы.

Специфичность: почему одной Tm мало

Праймер не ищет ген как объект - он садится на любую комплементарную ему последовательность, а в геноме человека таких мест может оказаться несколько. Поэтому после расчёта $T_m$ пару прогоняют через поиск по геному (BLAST или Primer-BLAST) и смотрят, есть ли вторая площадка посадки на разумном расстоянии. Если её нет, продукт будет один.

Заодно фиксируют длину ампликона - расстояние между внешними границами праймеров. Пусть она равна 480 парам нуклеотидов. Тогда время элонгации при 72 градусах берут по скорости Taq-полимеразы около 1000 нуклеотидов в минуту:

$$t_{\text{эл}} = \frac{480}{1000} \cdot 60 \approx 30\ \text{с}.$$

Длина ампликона служит и контролем: полоса в геле обязана встать на отметке 480 пар нуклеотидов. Сам цикл реакции с тремя температурными ступенями разобран в статье ПЦР: принцип метода, здесь нас интересовал только подбор затравок. Если работаешь с РНК через обратную транскрипцию, праймеры дополнительно сажают на стык экзонов - тогда примесь геномной ДНК не даст продукта. А нуклеотидный состав матрицы считается правилом Чаргаффа: этот расчёт разобран в задаче как определить состав ДНК.

Частые ошибки

• Берут среднее $T_m$ пары вместо минимального. При $T_m$ = 60 и 68 среднее 64 градуса означает, что слабый праймер не сядет вовсе. Запас считают от меньшей температуры плавления.
• Применяют правило Уоллеса к праймеру в 28 нуклеотидов. Формула $2(A+T) + 4(G+C)$ даёт для него 80-90 градусов, чего не бывает: за пределами 14-20 нт нужна поправка на длину и соль.
• Путают долю ГЦ праймера и долю ГЦ матрицы. В расчёт $T_m$ идёт состав самого праймера; состав матрицы влияет лишь на выбор участка.
• Считают $T_a$ равной $T_m$. При $T_a = T_m$ половина дуплексов уже разошлась, выход падает в разы. Рабочий отступ - 3-5 градусов вниз.
• Забывают проверить 3'-концы на комплементарность. Праймеры с идеальными $T_m$ и долей ГЦ дают в геле одни димеры, если их концы слипаются.
• Оставляют разницу в длине. Праймеры 18 и 26 нуклеотидов разойдутся по $T_m$ на 10 и более градусов, даже если доли ГЦ у них одинаковы.

FAQ

Какой длины должен быть праймер? Рабочий диапазон - 18-25 нуклеотидов. Короче 18 нт последовательность перестаёт быть уникальной: 4 в степени 17 уже сопоставимо с размером генома человека. Длиннее 25 нт праймер дорожает, хуже диффундирует и чаще сворачивается в шпильку.

Почему доля ГЦ должна быть 40-60 процентов? Это компромисс. Ниже 40 % праймер слабо держится и требует низкой $T_a$, на которой растёт неспецифический отжиг. Выше 60 % он садится слишком прочно, в том числе на частично комплементарные места, и в геле появляются лишние полосы.

Что делать, если температуры плавления пары разошлись на 8 градусов? Не усреднять, а править последовательности: укоротить или удлинить один праймер на 1-2 нуклеотида либо сдвинуть его на матрице. Сдвиг на три-четыре позиции обычно меняет $T_m$ на 2-6 градусов и сводит пару.

Обязательно ли считать $T_m$ вручную, если есть онлайн-калькуляторы? На контрольной - да, там спрашивают именно правило Уоллеса. В работе формулу тоже нужно понимать: калькуляторы используют разные модели и разные концентрации соли, и без этого выбрать между их ответами, расходящимися на 5-10 градусов, невозможно.

Коротко

1. Выписать нуклеотидный состав праймера и проверить, что $A + T + G + C = N$.
2. Найти долю ГЦ по формуле $(G+C)/N \cdot 100\%$ и убедиться, что она попала в окно 40-60 %.
3. Посчитать $T_m = 2(A+T) + 4(G+C)$ - правило применимо для 14-20 нуклеотидов; для длинных праймеров брать формулу с $\lg[\mathrm{Na^+}]$ и членом $600/N$.
4. Свести пару: расхождение $\Delta T_m$ не больше 3 градусов, иначе править длину одного из праймеров.
5. Взять $T_a$ на 5 градусов ниже меньшей $T_m$, проверить 3'-концы на димеры, шпильки и ГЦ-замок. Для примера из условия: ГЦ 50 и 55 %, $T_m$ = 60 и 62 °C, $T_a = 55\ ^\circ\text{C}$.