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 del mundo microscópico de una célula, el material genético no siempre es una doble hélice estática. A veces, las hebras de ADN que normalmente se aparean entre sí se ven interrumpidas cuando una hebra de ARN aparece y se une a una de ellas, desplazando la otra hebra de ADN. Esta estructura de tres partes, donde el ARN y el ADN se toman de la mano mientras una tercera hebra espera en la sombra, se conoce como bucle R (R-loop). Estas formaciones son rasgos naturales de la biología, desempeñando papeles en cómo se leen y regulan los genes, pero también pueden causar inestabilidad si aparecen en el lugar equivíndido o en el momento equivocado. Comprender dónde y por qué se forman estas estructuras es una pregunta clave para los científicos que estudian los virus, porque la composición genética de un virus puede influir en cómo interactúa con las células humanas que infecta. El desafío radica en distinguir entre una estructura que es simplemente posible según las leyes de la química y una que realmente sucede dentro de una célula viva.

Un nuevo estudio adopta un enfoque diferente para este problema al centrarse enteramente en la energía química necesaria para construir estas estructuras. Los investigadores no miraron dentro de una célula ni midieron la actividad biológica directamente. En su lugar, construyeron un marco computacional para calcular la estabilidad termodinámica de secuencias virales específicas. En este contexto, la termodinámica se refiere a los cambios de energía que ocurren cuando las moléculas se unen; un valor de energía más negativo generalmente significa que un enlace es más fuerte y es más probable que se forme. El equipo se centró en el Virus de la Inmunodeficiencia Humana tipo 1, o VIH-1, utilizando su secuencia genética de referencia más común como mapa. Seleccionaron dieciséis regiones distintas a lo largo del código genético del virus, que variaban desde segmentos cortos hasta tramos más largos, y sometieron cada una a un riguroso cálculo de dos pasos.

Primero, los investigadores modelaron lo que sucedería si la versión de ARN de una región viral específica intentara unirse con una hebra de ADN complementaria. Utilizaron una herramienta de software especializada diseñada para predecir cómo se comportan estos híbridos de ARN-ADN. Después, modelaron lo que sucedería si esa misma secuencia de ADN intentara unirse con su propia pareja de ADN correspondiente, preguntando esencialmente qué tan estable es la doble hélice de ADN natural en ese mismo punto. Al comparar la energía del apareamiento ARN-ADN frente a la energía del apareamiento ADN-ADN, crearon un valor de contraste específico para cada región. Esta comparación les permitió ver si el código genético del virus tiene una tendencia química natural a formar la estructura de bucle R de tres hebras con mayor facilidad de la que mantiene su propia doble hélice.

Los resultados revelaron un patrón claro y consistente a través del genoma del virus. De las dieciséis regiones probadas, quince mostraron un valor de contraste negativo, lo que significa que el híbrido de ARN-ADN fue calculado como más estable que la estructura de ADN-ADN correspondiente. La diferencia promedio fue de ocho kilocalorías por mol, un margen significativo en el mundo de la química molecular. Solo una región, un segmento conocido como U5, mostró un valor positivo, indicando que en ese punto específico, el ADN prefería permanecer apareado consigo mismo. El análisis estadístico confirmó que esta tendencia no fue una ocurrencia aleatoria; la probabilidad de ver tantos resultados negativos por azar era extremadamente baja. Los investigadores también verificaron para asegurar que el resultado no estuviera sesgado por secciones superpuestas del genoma, y el patrón se mantuvo constante incluso cuando esas partes superpuestas fueron eliminadas del cálculo.

Sin embargo, el estudio es cuidadoso al definir los límites de lo que significan estos números. Los hallazgos son estrictamente una predicción basada en modelos químicos y no prueban que estos bucles R se formen realmente dentro de un cuerpo humano o que causen que el virus se integre en el ADN humano. Los cálculos no tuvieron en cuenta el complejo entorno de una célula, como la presencia de proteínas que podrían bloquear estas interacciones, el retorcimiento físico del ADN o la concentración de moléculas. El autor describe su trabajo como una forma de priorizar qué partes del genoma viral son químicamente más propensas a formar estas estructuras, sirviendo como una guía para futuros experimentos en lugar de una respuesta final. En esencia, el estudio proporciona un mapa detallado de dónde el código genético del virus está más ansioso por realizar un tipo específico de enlace, ofreciendo una hipótesis comprobable para que los científicos la verifiquen en el laboratorio.

La investigación también exploró si la longitud del segmento genético o la mezcla específica de componentes químicos dentro de este influía en los resultados. Aunque hubo algunos indicios de que las regiones con ciertas composiciones químicas mostraban tendencias más fuertes, el pequeño número de regiones probadas significó que estas observaciones fueron preliminares. La conclusión principal sigue centrada en la robustez de la diferencia de energía en sí misma. Al establecer un método reproducible para comparar la estabilidad de ARN-ADN y ADN-DNA a través de un genoma viral, el estudio ofrece una nueva herramienta para comprender el potencial químico de las secuencias virales. Sugiere que, para la mayor parte del genoma del VIH-1, las condiciones son químicamente favorables para la formación de estas estructuras de tres hebras, pero deja la pregunta de si la naturaleza realmente utiliza este potencial a la siguiente generación de experimentos biológicos.