Uma olhada prática em como essas bases se comportam no dia a dia
Se você já tentou calcular o conteúdo de GC de um molde de DNA e travou porque esqueceu que adenina e guanina têm tamanhos diferentes, já encontrou um problema real com purinas e pirimidinas. A coisa parece simples na teoria — adenina e guanina são purinas, citosina, timina e uracila são pirimidinas —, mas na prática a diferença estrutural dita comportamento químico, estabilidade térmica e até como seu primer vai falhar.
O que é purinas e pirimidinas, na prática
Purinas possuem um anel duplo (hexagono fundido com pentagono). Pirimidinas têm apenas um anel hexagonal. Isso altera diretamente a largura do par de bases no duplex: A-T e G-C têm praticamente a mesma largura, o que mantém a geometria regular da dupla hélice. Se você inverter a lógica e emparelhar uma purina com outra purina, o espaço fica enorme; pirimidina com pirimidina, fica apertado. O DNA simplesmente não aguenta essa distorção sem colapsar localmente. A parte que os livros pouco enfatizam é a solubilidade. Purinas são menos polares que pirimidinas devido à estrutura biccíclica. Em soluções aquosas com alto teor de sal, a solubilidade diferencial pode gerar precipitação seções seletivas durante a purificação por etanol. Já vi amplicons inteiros precipitarem de forma desigual quando a razão de bases era extremada, e o pellet ficava visivelmente mais amarelado — sinal clássico de excesso de purinas agregadas.
Cálculo de Tm e o erro comum
A equação de Wallace (Tm = 4×(G+C) + 2×(A+T)) funciona bem para oligos curtos, mas ela trata todas as purinas e pirimidinas como se fossero intercambiáveis dentro do grupo. Na verdade, pares G-C contribuem com três pontes de hidrogênio contra duas de A-T, então a presença de homopolímeros de guanina ou citosina eleva o Tm de forma não linear. Quando eu calculava Tm para regiões ricas em G de genes de termófilos, a previsão de Wallace subestimava em cerca de 4 a 6°C na prática. O workaround foi usar o modelo de SantaLucia com ajustes de vizinhança, que leva em conta o efeito do base par adjacente. Um detalhe que quase ninguém lembra: a metilação da citosina (5mC) aumenta o Tm em cerca de 0,5 a 1°C por modificação, e isso pode ser suficiente para que um primer projetado na sequência não metilada falhe no template metilado. Em estudos de metilação com PCR convencional, esse desvio é pequeno, mas em qPCR de alta resolução ele quebra a curva de melting.
👉 Clique no botão abaixo para saber mais sobre o assunto!
Purinas e pirimidinas em sequencing de nova geração
No contexto de NGS, a principal dor de cabeça com essas bases aparece nos homopolímeros. Plataformas de Ion Torrent e algumas variantes do PacBio sofrem diretamente com repetições de uma só purina ou pirimidina, porque a cinética de incorporação fica ambígua. Quando eu tinha um alvo com 7 guaninas seguidas, a chamagem de variante falhava consistentemente entre 6 e 7 cópias, mesmo com cobertura acima de 200×. A solução prática foi usar um polymerase de alta fidelidade com tolerância a repeats e fazer consenso com pelo menos três ciclos de leitura independente. Outro ponto: a oxidação da guanina gera 8-oxoG, que pareia erroneamente com adenina em vez de citosina. Em dados de sequenciamento, isso se manifesta como falsas transições G>T. Se o seu protocolo deLibrary prep não inclui tratamentos com formamida ou UDG, esses artefatos ficam escondidos na chamada. O fix é aplicar um filtro de qualidade por contexto de base e remover variantes que só aparecem em extremos de strand bias.
Design de primer com restrições de purinas e pirimidinas
A regra prática que eu sigo é manter a proporção de purinas e pirimidinas no primer entre 40% e 60%, evitando runs de mais de três bases idênticas. Se o seu gene alvo tem alta riqueza de GC, o primer tende a acumular pirimidinas no extremo 3′, o que reduz a eficiência de extensão porque a DNA polimerase reconhece pior terminações com baixo Tm local. Para contornar isso, eu adicio um tail de 2 a 3 bases ricas em A ou T no 5′ sem alterar a região de annealing, o que equilibra a composição sem comprometer a especificidade. Também costuma ajudar evitar que a última base do 3′ seja uma guanina. A guanina no terminal forma ligações fortes demais e pode induzir primer-dimer se houver complementaridade fraca com outro produto. No meu laboratório, isso resolveu um problema crônico de amplificação inespecífica em templates com alta complexidade secundária.
Quando essa abordagem falha
Não adianta confiar apenas na composição de purinas e pirimidinas para prever comportamento de PCR. Regiões com high GC content ainda exigem aditivos como DMSO, betaina ou formamida, e a concentração ideal varia de protocolo para protocolo. O método também não prediz eficientemente a formação de hairpins em RNA, onde a estabilidade depende de emparelhamentos não canônicos que fogem da lógica simples de par base a base. Nessas situações, simulações termodinâmicas com ferramentas como mfold ou RNAstructure são necessárias, e mesmo assim os resultados precisam de validação empírica.