A General Computational Framework for Sequence-Level Thermodynamic Analysis of Modeled RNA–DNA and DNA–DNA Duplex States: A Proof-of-Concept Study in HIV-1

Dentro do mundo microscópico de uma célula, o material genético nem sempre é uma dupla hélice estática. Às vezes, as fitas de DNA que costumam se emparelhar uma com a outra são interrompidas quando uma fita de RNA surge e se liga a uma delas, afastando a outra fita de DNA. Essa estrutura de três partes, onde o RNA e o DNA dão as mãos enquanto uma terceira fita espera nos bastidores, é conhecida como um laço R (R-loop). Essas formações são características naturais da biologia, desempenhando papais em como os genes são lidos e regulados, mas também podem causar instabilidade se aparecerem no lugar errado ou na hora errada. Compreender onde e por que essas estruturas podem se formar é uma questão fundamental para cientistas que estudam vírus, pois a composição genética de um vírus pode influenciar como ele interage com as células humanas que infecta. O desafio reside em distinguir entre uma estrutura que é meramente possível de acordo com as leis da química e uma que realmente acontece dentro de uma célula viva.

Um novo estudo adota uma abordagem diferente para este problema, focando inteiramente na energia química necessária para construir essas estruturas. Os pesquisadores não olharam para dentro de uma célula ou mediram a atividade biológica diretamente. Em vez disso, construíram uma estrutura computacional para calcular a estabilidade termodinâmica de sequências virais específicas. Neste contexto, a termodinâmica refere-se às mudanças de energia que ocorrem quando as moléculas se unem; um valor de energia mais negativo geralmente significa que uma ligação é mais forte e tem maior probabilidade de se formar. A equipe focou no Vírus da Imunodeficiência Humana tipo 1, ou HIV-1, usando sua sequência genética de referência mais comum como um mapa. Eles selecionaram dezesseis regiões distintas ao longo do código genético do vírus, variando de segmentos curtos a trechos mais longos, e submeteram cada uma a um rigoroso cálculo de duas etapas.

Primeiro, os pesquisadores modelaram o que aconteceria se a versão de RNA de uma região viral específica tentasse se ligar com uma fita de DNA complementar. Eles usaram uma ferramenta de software especializada, projetada para prever como esses híbridos de RNA-DNA se comportam. Em seguida, modelaram o que aconteceria se essa mesma sequência de DNA tentasse se ligar com seu próprio parceiro de DNA correspondente, essencialmente perguntando quão estável é a dupla hélice de DNA natural naquele mesmo ponto. Ao comparar a energia do emparelhamento RNA-DNA com a energia do emparelhamento DNA-DNA, eles criaram um valor de contraste específico para cada região. Essa comparação permitiu ver se o código genético do vírus tem uma tendência química natural de formar a estrutura de três fitas do laço R mais facilmente do que mantém sua própria dupla hélice.

Os resultados revelaram um padrão claro e consistente através do genoma do vírus. Das dezesseis regiões testadas, quinze mostraram um valor de contraste negativo, significando que o híbrido RNA-DNA foi calculado como sendo mais estável do que a estrutura DNA-DNA correspondente. A diferença média foi de oito quilocalorias por mol, uma margem significativa no mundo da química molecular. Apenas uma região, um segmento conhecido como U5, mostrou um valor positivo, indicando que, naquele ponto específico, o DNA preferia permanecer pareado consigo mesmo. A análise estatística confirmou que essa tendência não foi uma ocorrência aleatória; a probabilidade de ver tantos resultados negativos por acaso era extremamente baixa. Os pesquisadores também verificaram se o resultado não era enviesado por seções sobrepostas do genoma, e o padrão manteve-se constante mesmo quando essas partes sobrepostas foram removidas do cálculo.

No entanto, o estudo é cuidadoso ao definir os limites do que esses números significam. As descobertas são estritamente uma previsão baseada em modelos químicos e não provam que esses laços R realmente se formam dentro de um corpo humano ou que eles causem a integração do vírus no DNA humano. Os cálculos não levaram em conta o ambiente complexo de uma célula, como a presença de proteínas que podem bloquear essas interações, o torcimento físico do DNA ou a concentração de moléculas. O autor descreve seu trabalho como uma forma de priorizar quais partes do genoma viral são quimicamente mais propensas a formar essas estruturas, servindo como um guia para experimentos futuros, em vez de uma resposta final. Em essência, o estudo fornece um mapa detalhado de onde o código genético do vírus está mais ansioso para fazer um tipo específico de ligação, oferecendo uma hipótese testável para cientistas verificarem em laboratório.

A pesquisa também explorou se o comprimento do segmento genético ou a mistura específica de blocos de construção químicos dentro dele influenciava os resultados. Embora houvesse alguns indícios de que regiões com certas composições químicas mostravam tendências mais fortes, o pequeno número de regiões testadas significava que essas observações eram preliminares. A conclusão primária permanece focada na robustez da própria diferença de energia. Ao estabelecer um método reproduzível para comparar a estabilidade de RNA-DNA e DNA-DNA através de um genoma viral, o estudo oferece uma nova ferramenta para compreender o potencial químico de sequências virais. Sugere que, para a maior parte do genoma do HIV-1, as condições são quimicamente favoráveis para a formação dessas estruturas de três fitas, mas deixa a questão de se a natureza realmente utiliza esse potencial para a próxima geração de experimentos biológicos.