An Imaging-Informed Reaction-Diffusion Model of Infarct Growth

Wanneer een beroerte toeslaat, wordt de bloedtoevoer naar de hersenen afgesneden, wat een stille, zich uitbreidende crisis in gang zet. Het gebied van het weefsel dat onmiddellijk sterft is de kern, maar daaromheen ligt een kwetsbare zone die de penumbra wordt genoemd. Dit weefsel hunkert naar zuurstof en voedingsstoffen, maar is nog niet dood; het bevindt zich in een precaire staat waarin het ofwel kan herstellen, ofwel kan bezwijken onder de schade die zich vanuit de kern verspreidt. Artsen moeten precies weten hoeveel van deze penumbra uiteindelijk zal veranderen in dood weefsel, aangezien deze voorspelling de levensreddende beslissingen stuurt over welke patiënten onmiddellijke chirurgie nodig hebben en welke gered kunnen worden door medicatie. Momenteel zijn de instrumenten die worden gebruikt om deze voorspellingen te doen vaak een gok. Sommige vertrouwen op complexe computerprogramma's die leren van duizenden eerdere gevallen, maar niet kunnen uitleggen waarom ze een bepaalde voorspelling doen, terwijl andere gedetailleerde biologische modellen gebruiken die te ingewikkeld zijn om uit te voeren op de beelden die artsen daadwerkelijk in het ziekenhuis tot hun beschikking hebben.

Een team van onderzoekers heeft een ander pad gekozen door een model te bouwen dat de zich uitbreidende schade behandelt als een fysieke golf die door een landschap beweegt. In plaats van te gokken op basis van patronen of elke minuscule elektrische signaal binnen een cel te simuleren, gebruikten zij een raamwerk dat beschrijft hoe een verstoring zich door de ruimte en tijd verspreidt. Ze testten dit idee op een kleine groep patiënten, met behulp van standaard magnetische resonantie-imaging scans die werden gemaakt bij hun eerste aankomst in het ziekenhuis. Door de initiële schade en de bloedstroomkaarten in hun vergelijkingen te voeren, simuleerden ze hoe de letsel zich over drie maanden zou uitbreiden. Cruciaal was dat dit een "oracle"-studie was: de onderzoekers gebruikten de uiteindelijke uitkomsten van negentig dagen later om de instellingen van het model af te stemmen, in plaats van een realtime voorspelling te doen voordat de uitkomst bekend was. Deze aanpak stelde hen in staat om de theoretische capaciteiten en structurele beperkingen van het biofysische model te definiëren. De resultaten toonden aan dat deze op fysica gebaseerde aanpak de uiteindelijke grootte en vorm van het beschadigde gebied aanzienlijk beter kon vastleggen dan de standaardmethoden die artsen vandaag de dag gebruiken, wat een nieuwe manier biedt om de toekomst van een beroerte te zien die geworteld is in de werkelijke wetten van hoe weefsel zich gedraagt.

De onderzoekers richtten zich op een specifiek wiskundig concept dat bekend staat als een reactie-diffusieproces. In eenvoudige termen beschrijft dit twee dingen die tegelijkertijd gebeuren: een lokale reactie waarbij een cel besluit te sterven, en een diffusieproces waarbij de stress die die dood veroorzaakt zich naar naburige cellen verspreidt. Stel je een druppel inkt voor die zich in water verspreidt; de inkt beweegt naar buiten, maar in de hersenen is de "inkt" een toxische opbouw van chemicaliën die cellen dwingt op te geven. Het team creëerde een computersimulatie waarbij het hersenweefsel werd verdeeld in kleine kubussen, en ze keken hoe de schade groeide van de initiële kern naar het omliggende, gezond ogende weefsel. Cruciaal was dat ze de schade niet overal gelijkmatig lieten verspreiden. Ze gebruikten de eigen scan van de patiënt om de simulatie te begrenzen, wat betekent dat het model de schade alleen liet uitbreiden naar gebieden waar de bloedstroom al gevaarlijk laag was. Dit zorgde ervoor dat het model naar het juiste gebied keek en gezonde delen van de hersenen die veilig waren, negeerde.

Om te testen of dit idee werkte, bekeken het team de gegevens van eenentwintig patiënten die een ischemische beroerte hadden gehad. Voor elke persoon namen ze de scan van de dag van de beroerte en draaide hun simulatie vooruit in de tijd, waarbij ze voorspelden hoe de hersenen er negentig dagen later uit zouden zien. Vervolgens vergeleken ze hun voorspelling met de werkelijke scan die negentig dagen na het incident werd gemaakt, die de werkelijke uiteindelijke omvang van het dode weefsel liet zien. De resultaten waren opmerkelijk. De standaardmethode die artsen gebruiken, die simpelweg een lijn trekt rond gebieden met een lage bloedstroom, gokte vaak de verkeerde grootte en voorspelde meestal veel te veel schade. In tegen plaats kwam het nieuwe, op fysica gebaseerde model veel dichter bij de realiteit. Het identificeerde de vorm en omvang van het uiteindelijke letsel op een manier die de eenvoudige bloedstroomkaarten niet konden. Het model was bijzonder goed in het vastleggen van de nuance van hoe de schade zich verspreidt, waarbij het aantoonde dat de snelheid en richting van de verspreiding sterk afhangen van hoe slecht de bloedstroom in elke specifieke plek is.

De studie onthulde ook waarom deze aanpak zo goed werkt. De onderzoekers ontdekten dat de belangrijkste factor niet alleen was hoe snel de cellen stierven, maar hoe gemakkelijk de toxische stress van de ene cel naar de volgende kon bewegen. In gebieden waar de bloedstroom ernstig beperkt was, zwol het weefsel op, waardoor het moeilijker werd voor de stresssignalen om te reizen. Het model hield hier rekening mee door de verspreiding van schade in de meest kritieke gebieden te vertragen, wat voorkwam dat het de uiteindelijke letseling overschatte. Wanneer ze probeerden het model complexer te maken door meer variabelen toe te voegen, verbeterde dit de resultaten niet veel, wat suggereert dat de sleutel tot nauwkeurigheid simpelweg het begrijpen van hoe de fysieke omgeving van de hersenen de verspreiding van schade beïnvloedt was. Deze eenvoud is een groot voordeel, omdat het betekent dat het model gemakkelijker te begrijpen en te vertrouwen is dan de "black box" computerprogramma's die gebruikelijk zijn in de moderne geneeskunde.

De onderzoekers waren echter voorzichtig om de beperkingen van hun werk te vermelden. Dit was een simulatiestudie, wat betekent dat ze de uiteindelijke uitkomst gebruikten om de instellingen van het model af te stemmen, in plaats van te proberen de toekomst in realtime te voorspellen voordat de drie maanden om waren. Ze vonden ook dat het model niet om kon gaan met gevallen waarin de bloedstroom van een patiënt zeer snel succesvol was hersteld, aangezien het model ervan uitgaat dat de schade blijft groeien zolang de bloedtoevoer wordt onderbroken. Bovendien, hoewel het op fysica gebaseerde model de standaardmethoden overtrof, was de verbetering in de primaire nauwkeurigheidsmetriek (Dice-score) marginaal wanneer men rekening houdt met de grote variaties in de resultaten over de kleine groep patiënten. Ondanks deze beperkingen bewijst de studie dat het mogelijk is om de basiswetten van de fysica te gebruiken om de evolutie van een beroerte direct vanuit een ziekenhuisscan te volgen. Het biedt een veelbelovend alternatief voor de huidige methoden, door een instrument te bieden dat zowel wetenschappelijk gefundeerd is als in staat is om zijn eigen redenering uit te leggen, een kwaliteit die essentieel is voor artsen die belangrijke beslissingen nemen voor hun patiënten.