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

Nel mondo microscopico di una cella, il materiale genetico non è sempre una doppia elica statica. A volte, i filamenti di DNA che solitamente si accoppiano tra loro vengono interrotti quando un filamento di RNA interviene e si lega a uno di essi, sporgendo l'altro filamento di DNA. Questa struttura a tre parti, in cui l'RNA e il DNA si tengono per mano mentre un terzo filamento attende nelle retrovie, è nota come R-loop. Queste formazioni sono caratteristiche naturali della biologia, svolgono ruoli nel modo in cui i geni vengono letti e regolati, ma possono anche causare instabilità se compaiono nel posto sbagliato o al momento sbagliato. Comprendere dove e perché queste strutture potrebbero formarsi è una domanda chiave per gli scienziati che studiano i virus, poiché la composizione genetica di un virus può influenzare il modo in cui interagisce con le cellule umane che infetta. La sfida consiste nel distinguere tra una struttura che è semplicemente possibile secondo le leggi della chimica e una che accade realmente all'interno di una cellula vivente.

Un nuovo studio adotta un approccio diverso a questo problema concentrandosi interamente sull'energia chimica necessaria per costruire queste strutture. I ricercatori non hanno guardato all'interno di una cella né hanno misurato direttamente l'attività biologica. Invece, hanno costruito un quadro computazionale per calcolare la stabilità termodinamica di specifiche sequenze virali. In questo contesto, la termodinamica si riferisce ai cambiamenti di energia che si verificano quando le molecole si uniscono; un valore di energia più negativo generalmente significa che un legame è più forte e più propenso a formarsi. Il team si è concentrato sul Virus dell'Immunodeficienza Umana di tipo 1, o HIV-1, utilizzando la sua sequenza genetica di riferimento più comune come mappa. Hanno selezionato sedici regioni distinte lungo il codice genetico del virus, che vanno da segmenti brevi a tratti più lunghi, e hanno sottoposto ciascuna di esse a un rigoroso calcolo in due fasi.

In primo luogo, i ricercatori hanno modellato ciò che accadrebbe se la versione RNA di una specifica regione virale tentasse di legarsi con un filamento di DNA complementare. Hanno utilizzato uno strumento software specializzato progettato per prevedere come si comportano questi ibridi RNA-DNA. Successivamente, hanno modellato ciò che accadrebbe se quella stessa sequenza di DNA tentasse di legarsi con il proprio partner di DNA corrispondente, chiedendosi essenzialmente quanto fosse stabile la naturale doppia elica di DNA in quel punto. Confrontando l'energia dell'accoppiamento RNA-DNA con l'energia dell'accoppiamento DNA-DNA, hanno creato un valore di contrasto specifico per ogni regione. Questo confronto ha permesso loro di vedere se il codice genetico del virus ha una naturale tendenza chimica a formare la struttura a tre filamenti dell'R-loop più facilmente di quanto non mantenga la propria doppia elica.

I risultati hanno rivelato un modello chiaro e coerente attraverso il genoma del virus. Delle sedici regioni testate, quindici hanno mostrato un valore di contrasto negativo, il che significa che l'ibrido RNA-DNA era calcolato essere più stabile della corrispondente struttura DNA-DNA. La differenza media era di otto kilocalorie per mole, un margine significativo nel mondo della chimica molecolare. Solo una regione, un segmento noto come U5, ha mostrato un valore positivo, indicando che in quel punto specifico il DNA preferiva rimanere accoppiato con se stesso. L'analisi statistica ha confermato che questa tendenza non era un evento casuale; la probabilità di vedere così tanti risultati negativi per puro caso era estremamente bassa. I ricercatori hanno anche verificato che il risultato non fosse distorto da sezioni sovrapposte del genoma, e il modello è rimasto costante anche quando tali parti sovrapposte sono state rimosse dal calcolo.

Tuttavia, lo studio è attento a definire i limiti di ciò che significano questi numeri. Le scoperte sono strettamente una previsione basata su modelli chimici e non provano che questi R-loop si formino effettivamente all'interno di un corpo umano o che causino l'integrazione del virus nel DNA umano. I calcoli non hanno tenuto conto del complesso ambiente di una cellula, come la presenza di proteine che potrebbero bloccare queste interazioni, la torsione fisica del DNA o la concentrazione delle molecole. L'autore descrive il proprio lavoro come un modo per dare priorità alle parti del genoma virale che sono chimicamente più probabili nel formare queste strutture, fungendo da guida per esperimenti futuri piuttosto che da risposta definitiva. In sostanza, lo studio fornisce una mappa dettagliata di dove il codice genetico del virus sia più desideroso di creare un tipo specifico di legame, offrendo un'ipotesi testabile che gli scienziati possano verificare in laboratorio.

La ricerca ha anche esplorato se la lunghezza del segmento genetico o la specifica miscela di componenti chimici al suo interno influenzassero i risultati. Sebbene ci siano stati alcuni indizi che regioni con determinate composizioni chimiche mostrassero tendenze più forti, il piccolo numero di regioni testate rendeva queste osservazioni preliminari. La conclusione primaria rimane focalizzata sulla robustezza della differenza di energia stessa. Stabilendo un metodo riproducibile per confrontare la stabilità di RNA-DNA e DNA-DNA attraverso un genoma virale, lo studio offre un nuovo strumento per comprendere il potenziale chimico delle sequenze virali. Suggerisce che per la maggior parte del genoma dell'HIV-1, le condizioni sono chimicamente favorevoli per la formazione di queste strutture a tre filamenti, ma lascia la questione se la natura utilizzi effettivamente questo potenziale alla prossima generazione di esperimenti biologici.