A digital twin for microwave liver treatment replanning

Technische Samenvatting: Een Digitale Tweeling voor Herplanning van Microgolfbehandeling van de Lever

Probleemstelling
Microgolfablatie (MWA) is een toonaangevende minimaal invasieve behandeling voor levertumoren, met voordelen ten opzichte van andere thermische modaliteiten, zoals kortere behandeltijden en een verminderde gevoeligheid voor het "heat-sink"-effect veroorzaakt door bloedvaten. Een succesvolle behandeling is echter sterk afhankelijk van nauwkeurige preoperatieve planning en een precieze intraoperatieve plaatsing van de antenne. In de klinische praktijk maken fysiologische factoren, zoals ademhalingbewegingen, herpositionering van de patiënt en weefseldeformatie, het praktisch onmogelijk om de antenne exact zoals gepland te positioneren.

Wanneer er sprake is van een verkeerde plaatsing, leidt het strikt volgen van het oorspronkelijke plan vaak tot een onvolledige tumorablatie en een potentieel recidief. Omgekeerd is het opnieuw positioneren van de antenne via herinbrenging tijdrovend en brengt het aanzienlijke risico's met zich mee, waaronder tumorzaad (de verspreiding van kankercellen langs de naaldspoor). De huidige klinische planning vertrouwt op door de fabrikant geleverde geometrische modellen die patiëntspecifieke weefseleigenschappen negeren, wat vaak leidt tot een overschatting van de ablatiezones. Er is een kritieke behoefte aan een patiëntspecifiek simulatiehulpmiddel dat capable is van intraoperatieve herplanning om fouten in de plaatsing te corrigeren zonder dat herinbrenging van de antenne nodig is.

Methodologie
De auteurs stellen een "digitale tweeling"-kader voor dat bestaat uit twee hoofdcomponenten: een hoogwaardige MWA-simulator en een optimalisatiealgoritme voor herplanning.

1. Rekenmodel:

   • Fysica: Het systeem lost gekoppelde vergelijkingen voor microgolfvoortplanting en warmteoverdracht op. Warmteoverdracht wordt gemodelleerd met de Pennes Bioheat-vergelijking, terwijl microgolfvoortplanting wordt beheerst door Maxwell-vergelijkingen in het frequentiedomein.
   • Patiëntspecificiteit: Het model integreert patiëntspecifieke diëlektrische en thermische eigenschappen afgeleid van preoperatieve CT- en MRI-beelden. Het houdt rekening met temperatuurafhankelijke weefselveranderingen, waaronder coagulatie, waterverdamping en het afsluiten van bloedperfusie ($\omega_b$).
   • Numerieke Implementatie: Er wordt een Finite Element Method (FEM)-schema gebruikt met de FEniCSx-bibliotheek. Om nauwkeurigheid en rekenkosten in evenwicht te brengen, maken de auteurs gebruik van een zwak gekoppeld benadering waarbij de microgolfvoortplanting één keer wordt opgelost en de externe warmtebron ($Q_{ext}$) wordt benaderd als een temperatuurafhankelijke functie. Dit voorkomt de noodzaak om Maxwell's vergelijkingen bij elke thermische iteratie opnieuw op te lossen.
   • Antennekarakterisering: Omdat interne antennegeometrieën eigendom zijn, schat het model antennekarakteristieken (zoals afmetingen van de choke) door de S11-reflexiecoëfficiënt te optimaliseren tegen bekende externe geometrieën.
   • Validatie: Het model is gevalideerd tegen in vivo varkensexperimenten met drie ablatiescenario's: een baseline ver van grote vaten, één nabij de levervene en één nabij zowel lever- als poortvaten.

2. Herplannings-Optimizer:

   • Doel: Het corrigeren van antenneverplaatsing door drie parameters aan te passen: vermogensinput ($P$), ablatieduur ($T$) en antenneverschuiving langs zijn as ($\vec{t}$).
   • Algoritme: Een Genetisch Algoritme (GA) is geselecteerd vanwege de niet-lineaire aard van het probleem en het ontbreken van gradiëntinformatie. De doelfunctie minimaliseert een gewogen som van de inverse Dice-ähnheidscoëfficiënt en de 95e percentiel Hausdorff-afstand ($dH_{95}$) tussen de gesimuleerde ablatiezone en de doelzone.
   • Beperkingen: De optimalisatie werkt binnen realistische klinische grenzen: vermogen (20–150 W), duur (tot 15 min) en verschuiving ($\pm$5 cm).

Belangrijkste Resultaten

• Simulatie Nauwkeurigheid: De simulator toonde een hoge nauwkeurigheid in het voorspellen van ablatiezones in vergelijking met de grondwaarheid uit in vivo varkensdata. Het behaalde Dice-ähnheidscoëfficiënten van 0,82, 0,81 en 0,79 voor respectievelijk de baseline, levervene en poort/levervene scenario's. Deze resultaten presteerden aanzienlijk beter dan de geometrische modellen van de fabrikant, die Dice-scores opleverden van 0,62, 0,66 en 0,58. De simulator toonde ook robuustheid tegen het heat-sink-effect, met lage Hausdorff-afstanden ($dH_{95} \approx 2$ mm) zelfs in de buurt van grote vaten.
• Herplanningsprestaties: In een herplanningscenario waarbij de antenne kunstmatig 5 mm was verschoven, slaagde de optimizer erin parameters te herstellen om de doelablatiezone te matchen. De herplande simulaties behaalden Dice-scores van 0,83, 0,83 en 0,80, wat een gemiddelde verbetering vertegenwoordigt van respectievelijk 20,3%, 40,7% en 48,1% ten opzichte van de initiële (verkeerd geplaatste) ablatieresultaten. De optimizer herstelde vermogens- en verschuivingsparameters met lage relatieve fouten (bijvoorbeeld <2,3% voor vermogen in baseline gevallen), zelfs in complexe vasculaire omgevingen.

Betekenis en Claims
Het artikel claimt het eerste werk te zijn dat intraoperatieve herplanning van thermische ablatietherapie aanpakt. De primaire betekenis ligt in het aantonen dat een volledige tumorablatie kan worden bereikt zonder herinbrenging van de antenne door specifieke ablatieparameters (vermogen, tijd en axiale positie) te optimaliseren.

De auteurs positioneren dit werk als een proof-of-concept voor een klinisch beslissingsondersteunend hulpmiddel. Hoewel zij erkennen dat de huidige rekentijd (50–95 minuten) real-time gebruik verhindert zonder verdere versnelling (bijvoorbeeld via GPU's of surrogate-modellen), valideert de studie de haalbaarheid van het gebruik van patiëntspecifieke digitale tweelingen om onvermijdelijke intraoperatieve afwijkingen te corrigeren. Het werk suggereert dat een dergelijk hulpmiddel onvolledige ablaties kan voorkomen en de risico's geassocieerd met tumorzaad kan vermijden, met name in complexe gevallen met grote vaten waar standaard planningsmodellen falen.

Beperkingen
De auteurs merken op dat de huidige aanpak een baseline-ablatie vereist om de bloedperfusieparameter ($\omega_0$) te kalibreren, wat in een klinische setting niet haalbaar is zonder eerdere interventie. Zij suggereren dat toekomstig werk methoden kan integreren om perfusie direct te schatten op basis van CT-beelden. Bovendien houdt het huidige model geen rekening met heterogeen tumorweefsel, en is het herplanningshulpmiddel nog niet getest op daadwerkelijke tumormodellen.