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

Binnen de microscopische wereld van een cel is genetisch materiaal niet altijd een statische dubbele helix. Soms worden de strengen DNA die normaal gesproken met elkaar paren, onderbroken wanneer een RNA-streng langskomt en zich aan een van hen bindt, waardoor de andere DNA-streng opzij wordt geduwd. Deze drieledige structuur, waarbij RNA en DNA elkaars hand vasthouden terwijl een derde streng in de coulissen wacht, staat bekend als een R-loop. Deze formaties zijn natuurlijke verschijnselen in de biologie; ze spelen een rol in hoe genen worden afgelezen en gereguleerd, maar ze kunnen ook instabiliteit veroorzaken als ze op de verkeerde plek of op het verkeerde moment verschijnen. Het begrijpen van waar en waarom deze structuren gevormd kunnen worden, is een cruciale vraag voor wetenschappers die virussen bestudelen, omdat de genetische samenstelling van een virus invloed kan hebben op de manier waarop het interageert met de menselijke cellen die het infecteert. De uitdaging ligt in het onderscheid maken tussen een structuur die louter mogelijk is volgens de wetten van de chemie en een structuur die daadwerkelijk voorkomt in een levende cel.

Een nieuwe studie neemt een andere aanpak van dit probleem door zich volledig te richten op de chemische energie die nodig is om deze structuren op te bouwen. De onderzoekers keken niet in een cel of maten biologische activiteit direct. In plaats daarvan bouwden ze een computationeel kader om de thermodynamische stabiliteit van specifieke virale sequenties te berekenen. In deze context verwijst thermodynamica naar de energieveranderingen die optreden wanneer moleculen samenkomen; een meer negatieve energiewaarde betekent over het algemeen dat een binding sterker is en waarschijnlijker vormt. Het team richtte zich op het Human Immunodeficiency Virus type 1, of HIV-1, en gebruikte de meest voorkomende referentiegenetische sequentie van het virus als kaart. Ze selecteerden zestien verschillende regio's langs de genetische code van het virus, variërend van korte segmenten tot langere stukken, en onderwierpen elk segment aan een rigoureuze tweestapsberekening.

Eerst modelleerden de onderzoekers wat er zou gebeuren als de RNA-versie van een specifieke virale regio probeert te binden met een complementaire DNA-streng. Ze gebruikten een gespecialiseerde softwaretool die ontworig is om te voorspellen hoe deze RNA-DNA-hybriden zich gedragen. Vervolgens modelleerden ze wat er zou gebeuren als diezelfde DNA-sequentie zou proberen te binden met zijn eigen bijbehorende DNA-partner, wat in feite de vraag stelde hoe stabiel de natuurlijke DNA-dubbele helix op diezelfde plek is. Door de energie van de RNA-DNA-paring te vergelijken met de energie van de DNA-DNA-paring, creëerden ze een specifieke contrastwaarde voor elke regio. Deze vergelijking stelde hen in staat om te zien of de genetische code van het virus een natuurlijke chemische neiging heeft om de driestrandige R-loopstructuur gemakkelijker te vormen dan het vast te houden aan zijn eigen dubbele helix.

De resultaten toonden een duidelijk en consistent patroon over het hele genoom van het virus. Van de zestien geteste regio's vertoonden vijftien een negatieve contrastwaarde, wat betekent dat de RNA-DNA-hybride berekend werd stabieler te zijn dan de overeenkomstige DNA-DNA-structuur. Het gemiddelde verschil was acht kilocalorieën per mol, een aanzienlijke marge in de wereld van de moleculaire chemie. Slechts één regio, een segment bekend als U5, vertoonde een positieve waarde, wat aangeeft dat het DNA in dat specifieke punt de voorkeur gaf om met zichzelf gepaard te blijven. De statistische analyse bevestigde dat deze trend geen toevallige gebeurtenis was; de waarschijnlijkheid dat men zoveel negatieve resultaten door toeval zou zien, was extreem laag. De onderzoekers controleerden ook of het resultaat niet werd vertekend door overlappende secties van het genoom, en het patroon bleef standhouden zelfs toen die overlappende delen uit de berekening werden verwijderd.

De studie is echter voorzichtig in het definiëren van de grenzen van wat deze getallen betekenen. De bevindingen zijn strikt een voorspelling gebaseerd op chemische modellen en bewijzen niet dat deze R-lopen daadwerkelijk vormen binnen een menselijk lichaam of dat ze de integratie van het virus in menselijk DNA veroorzaken. De berekeningen hielden geen rekening met de complexe omgeving van een cel, zoals de aanwezigheid van eiwitten die deze interacties kunnen blokkeren, de fysieke draaiing van het DNA, of de concentratie van moleculen. De auteur beschrijft hun werk als een manier om te prioriteren welke delen van het virale genoom chemisch gezien het meest waarschijnlijk deze structuren vormen, wat dient als een gids voor toekomstige experimenten in plaats van een definitief antwoord. In essentie biedt de studie een gedetailleerde kaart van waar de genetische code van het virus het meest begeerig is om een specifiek type binding aan te gaan, wat een toetsbare hypothese biedt voor wetenschappers om in het laboratorium te verifiëren.

Het onderzoek onderzocht ook of de lengte van het genetische segment of de specifieke mix van chemische bouwstenen binnen het segment de resultaten beïnvloedde. Hoewel er enkele aanwijzingen waren dat regio's met bepaalde chemische composities een sterkere neiging vertoonden, waren deze observaties vanwege het kleine aantal geteste regio's voorlopig van aard. De primaire conclusie blijft gericht op de robuustheid van het energieverschil zelf. Door een reproduceerbare methode te vestigen om de stabiliteit van RNA-DNA en DNA-DNA over een viraal genoom te vergelijken, biedt de studie een nieuw instrument voor het begrijpen van het chemische potentieel van virale sequenties. Het suggereert dat voor het grootste deel van het HIV-1-genoom de omstandigheden chemisch gunstig zijn voor de vorming van deze driestrandige structuren, maar het laat de vraag of de natuur dit potentieel daadwerkelijk gebruikt over aan de volgende generatie biologische experimenten.