Kalkulator temperatury topnienia DNA
Oszacuj temperaturę topnienia (Tm) oligonukleotydu DNA na podstawie jego sekwencji i stężenia soli. Wpisz sekwencję jednej nici (A, C, G, T) oraz stężenie Na⁺, aby otrzymać Tm i zawartość GC.
Szacowanie temperatury topnienia oligonukleotydu DNA
Wzór na temperaturę topnienia DNA
Gdzie:
- \(A, T, G, C\) — liczba zasad adeniny, tyminy, guaniny i cytozyny w sekwencji
- \(n\) — łączna liczba zasad (długość sekwencji)
- \([\mathrm{Na^+}]\) — całkowite stężenie soli jednowartościowych w mol/l
Pierwszy wzór to reguła Wallace'a, stosowana dla krótkich oligonukleotydów do 13 zasad; drugi to wzór z korektą na sól (Howley i wsp., 1979), stosowany automatycznie od 14 zasad. Kalkulator sam wybiera odpowiedni wzór.
Temperatura topnienia DNA – przykład obliczenia
Przykład: 15-zasadowy oligonukleotyd o sekwencji ATGCATGCGATTACA w buforze zawierającym 50 mM Na⁺. Sekwencja zawiera 5 zasad A, 4 T, 3 G i 3 C, więc \(G + C = 6\), \(A + T = 9\) oraz \(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} \]
Wynik: temperatura topnienia wynosi \(T_m \approx 40.6\) °C, a zawartość GC — 40,0%. Prosta reguła Wallace'a dałaby 2 × 9 + 4 × 6 = 42 °C, ale nie zaleca się jej stosowania dla sekwencji liczących 14 zasad i więcej.
Jak korzystać z kalkulatora
- Wpisz lub wklej sekwencję w pole Sekwencja DNA (5′→3′). Akceptowane są wyłącznie litery A, C, G i T; wielkość liter, spacje i podziały wierszy są ignorowane.
- Podaj całkowite stężenie soli jednowartościowych w polu Stężenie soli (Na⁺), w mM (domyślnie 50 mM). Wartość ta jest wykorzystywana tylko wtedy, gdy sekwencja ma co najmniej 14 zasad.
- Odczytaj temperaturę topnienia (Tm) w °C oraz zawartość GC w procentach w panelu wyników.
- Zajrzyj do sekcji Szczegóły obliczeń, aby sprawdzić skład zasad, zastosowany wzór i podstawione liczby — dzięki temu możesz zweryfikować wynik ręcznie.
Czym jest temperatura topnienia?
Temperatura topnienia, oznaczana jako Tm, to temperatura, w której połowa dupleksów DNA w roztworze rozpada się na pojedyncze nici, a połowa pozostaje sparowana. To bezpośrednia miara stabilności dupleksu: im wyższa Tm, tym więcej ciepła potrzeba do rozdzielenia obu nici. Tm ma znaczenie wszędzie tam, gdzie dwie komplementarne nici DNA muszą niezawodnie się odnaleźć i połączyć — przy projektowaniu starterów PCR, doborze warunków hybrydyzacji sond czy wyznaczaniu temperatur przyłączania krótkich oligonukleotydów.
Stabilność dupleksu wynika z samych par zasad. Parę G–C utrzymują trzy wiązania wodorowe, a para A–T ma ich tylko dwa, dlatego sekwencje bogate w GC topią się w wyższych temperaturach. W regule Wallace'a wyraża się to bardzo wprost: każda para A–T dodaje do Tm około 2 °C, a każda para G–C — około 4 °C. Podawana obok Tm zawartość GC pokazuje, jak zbudowana jest sekwencja — to szybki test poprawności, ponieważ zawartość GC i Tm zawsze rosną i maleją razem.
Sól również odgrywa ważną rolę. Jony sodu ekranują ładunki ujemne szkieletu fosforanowego, więc wyższe stężenie Na⁺ stabilizuje dupleks i podwyższa Tm. Właśnie dlatego wzór z korektą na sól zawiera człon \(16.6\,\log_{10}[\mathrm{Na^+}]\): każdy dziesięciokrotny wzrost stężenia soli dodaje do przewidywanej Tm około 16,6 °C. Porównując oligonukleotydy, zawsze porównuj je przy tym samym stężeniu soli i pamiętaj, że stężenie podaje się w mM — przeliczenie na mol/l odbywa się automatycznie.
Dokładność, założenia i praktyczne wskazówki
- Oba wzory zakładają idealny dupleks: sekwencja musi hybrydyzować z w pełni komplementarną nicią, bez niedopasowań, pętli ani przerw. Nawet pojedyncze niedopasowanie może obniżyć rzeczywistą Tm o kilka stopni.
- Ograniczenia reguły Wallace'a: to zgrubna reguła przeznaczona dla sekwencji liczących około 4–13 zasad przy stężeniu Na⁺ zbliżonym do 50 mM. Nie uwzględnia stężenia soli ani efektów długości sekwencji, więc stosuj ją wyłącznie do szybkich szacunków dla bardzo krótkich oligonukleotydów.
- Ograniczenia wzoru z korektą na sól: uwzględnia on długość i sole jednowartościowe, ale nie stężenie oligonukleotydu, efekty końcowe, czynniki denaturujące (formamid, DMSO) ani kationów dwuwartościowych, takich jak Mg²⁺, które w buforach PCR istotnie stabilizują dupleksy.
- Do dokładnych prac używaj modelu najbliższych sąsiadów: nowoczesne programy do projektowania starterów opierają się na termodynamice najbliższych sąsiadów, uwzględniającej energię stapiania każdego kroku dinukleotydowego. Oba wzory dostępne tutaj to szybkie przybliżenia o dokładności rzędu kilku stopni — idealne do planowania i sprawdzania krzyżowego, ale nie do ostatecznego projektowania testów.
- Interpretacja wyniku: jako regułę można przyjąć, że hybrydyzacja lub przyłączanie przebiega najlepiej około 5 °C poniżej Tm. Jeśli obliczona Tm daleko odbiega od wartości docelowej, zmodyfikuj długość startera lub zawartość GC, a nie samą temperaturę przyłączania.
Jeśli w tych samych doświadczeniach pracujesz z peptydami lub białkami, przydatny może się okazać również nasz kalkulator punktu izoelektrycznego przy planowaniu warunków elektroforezy, a kalkulator pojemności buforowej pomoże sprawdzić, czy bufor hybrydyzacyjny utrzyma pH podczas podgrzewania.
Najczęściej zadawane pytania
Dlaczego różne kalkulatory Tm podają różne wartości dla tej samej sekwencji?
Różne kalkulatory opierają się na różnych modelach. Proste reguły (Wallace'a, podstawowe wzory z korektą na sól) pomijają kolejność zasad w sekwencji, natomiast model najbliższych sąsiadów uwzględnia energię stapiania każdego kroku dinukleotydowego, więc jego wyniki zwykle różnią się o kilka stopni. Korekty na Mg²⁺, stężenie oligonukleotydu i DMSO również przesuwają wynik. Zawsze porównuj wartości Tm otrzymane tą samą metodą — mieszanie Tm z różnych kalkulatorów przy projektowaniu starterów to częsty błąd.
Jaka zawartość GC jest dobra dla startera lub oligonukleotydu?
Celuj w zawartość GC 40–60%. Startery w tym zakresie wiążą się na tyle silnie, by być specyficzne, a jednocześnie niezawodnie się denaturują i przyłączają. Regiony bogate w GC, powyżej 65%, mogą tworzyć stabilne struktury drugorzędowe i wymagają wyższych temperatur przyłączania, natomiast startery bogate w AT, poniżej 40%, mogą wiązać się słabo lub niespecyficznie. Unikaj też długich ciągów tej samej zasady (np. GGGGG) oraz końca 3′ z wieloma zasadami G lub C, który sprzyja niespecyficznym inicjacjom.
Jak wybrać temperaturę przyłączania w PCR na podstawie Tm?
Częstym punktem wyjścia jest ustawienie temperatury przyłączania około 3–5 °C poniżej niższej Tm z pary starterów. Jeśli pojawiają się niespecyficzne prążki, podnoś temperaturę przyłączania krokami co 1–2 °C; jeśli nie ma produktu, nieco ją obniż. Nowoczesne protokoły często stosują PCR gradientowy lub dwuetapowy program touchdown, aby empirycznie znaleźć optimum, ponieważ proste wzory na Tm są jedynie przybliżeniem rzeczywistego zachowania w buforze PCR.
Czy mogę użyć tego kalkulatora do sekwencji RNA lub sekwencji zawierających U?
Nie. Kalkulator akceptuje wyłącznie zasady DNA: A, C, G i T — sekwencje RNA z uracylem (U) są odrzucane. Hybrydy RNA/DNA i RNA/RNA są też bardziej stabilne niż dupleksy DNA/DNA, więc reguła Wallace'a i wzór z korektą na sól zaniżałyby ich Tm nawet po podstawieniu T w miejsce U. Dla oligonukleotydów RNA użyj zamiast tego kalkulatora opartego na modelu najbliższych sąsiadów z termodynamicznymi parametrami dla RNA.
Jak Mg²⁺ i DMSO wpływają na obliczoną Tm?
Oba wzory zakładają wyłącznie sole jednowartościowe, więc w typowych buforach PCR zaniżają Tm. Mg²⁺ silnie stabilizuje dupleksy, często podnosząc efektywną Tm o kilka stopni, natomiast czynniki denaturujące, takie jak DMSO czy formamid, ją obniżają (mniej więcej 0,6 °C na każdy 1% DMSO). Jeśli bufor zawiera Mg²⁺, traktuj wynik jako ostrożne oszacowanie i optymalizuj eksperymentalnie albo skorzystaj z narzędzia opartego na modelu najbliższych sąsiadów, które uwzględnia stężenia kationów dwuwartościowych.
Co się dzieje z Tm, gdy mój oligonukleotyd ma niedopasowanie do sekwencji docelowej?
Kalkulator zakłada idealnie komplementarny dupleks, więc nie uwzględnia niedopasowań. Nawet pojedyncze niedopasowanie, zwłaszcza w środku sekwencji lub przy końcu 3′ startera, może obniżyć rzeczywistą Tm o 5–15 °C i osłabić wydajność amplifikacji. Dla starterów zdegenerowanych lub sond celujących w zmienne regiony użyj narzędzia najbliższych sąsiadów, które ocenia każdy typ niedopasowania osobno, i przetestuj warunki eksperymentalnie.
Stężenie soli podaję w mM czy w mol/l i jak je przeliczyć?
Podaj stężenie Na⁺ w mM; kalkulator przelicza je wewnętrznie na mol/l do członu logarytmicznego \(16.6\,\log_{10}[\mathrm{Na^+}]\). Aby przeliczyć ręcznie, podziel mM przez 1000: 50 mM = 0,05 mol/l, 500 mM = 0,5 mol/l. Pamiętaj, że chodzi o całkowite stężenie kationów jednowartościowych — niektóre protokoły podają siłę jonową albo uwzględniają Tris i K⁺, które należy doliczyć do wartości Na⁺, aby uzyskać realistyczne oszacowanie.