A Unified Model for Blood and Lymph Flow with Coupled Nonsmooth Biochemical Dynamics

Titel: Een unificerend model voor bloed- en lymfestroming met gekoppelde niet-gladde biochemische dynamica

Auteurs: Bogna Jaszczak-Dyka en Łukasz Płociniczak
Taal van de samenvatting: Nederlands

1. Probleemstelling

Het lymfestelsel speelt een cruciale rol in het handhaven van weefselhomeostase, het terugvoeren van interstitiële vloeistof en het ondersteunen van het immuunsysteem. Ondanks zijn belang voor de gezondheid (bijv. bij oedeem, kanker en cardiovasculaire ziekten), is het wiskundig modelleren van lymfestroming aanzienlijk minder ontwikkeld dan dat van de hemodynamica (bloedcirculatie).

De belangrijkste uitdagingen zijn:

• Fysische verschillen: Het lymfestelsel is een open, lineair netwerk met lage druk en pulserende stroming, in tegenstelling tot het gesloten, hoge-druk systeem van het hart.
• Actieve pompmechanisme: Lymfe wordt niet door een centraal orgaan gepompt, maar door ritmische samentrekkingen van de vaatwanden (lymfangions), gereguleerd door een complex mechano-chemisch oscillatormechanisme.
• Biochemische koppeling: De opening en sluiting van kleppen wordt gestuurd door de interactie tussen calciumionen ($Ca^{2+}$) en stikstofmonoxide (NO), wat leidt tot niet-gladde (discontinue) dynamiek die moeilijk te modelleren is met standaard differentiaalvergelijkingen.

Het doel van dit onderzoek is een unificerend wiskundig raamwerk te ontwikkelen dat zowel de vloeistofmechanica als de biochemische regulatie van lymfe- en bloedstroom beschrijft, met een specifieke focus op de oscillatoire aard van lymfe-transport.

2. Methodologie

De auteurs hanteren een tweeledige aanpak: een afleiding van de vloeistofmechanica vanuit eerste principes en een analyse van de biochemische kinetiek als een niet-glad systeem.

A. Afleiding van het stromingsmodel (Vloeistofmechanica)

• Basisvergelijkingen: Het model start met de Navier-Stokes-vergelijkingen in cilindrische coördinaten voor een axiaal symmetrische, laminaire stroming.
• Aannames: De vloeistof wordt als Newtoniaans beschouwd (een goede benadering voor lymfe), en de wanden reageren dynamisch op de stroming.
• Schaalanalyse en Lubricatiebenadering: Gezien de kleine aspectverhouding ($\epsilon = R_0/L \ll 1$) van de vaten, wordt een asymptotische expansie toegepast. Dit leidt tot een vereenvoudigd systeem van partiële differentiaalvergelijkingen (PDE's) dat de massabehoud en impulsbehoud beschrijft.
• Constitutieve wet: De relatie tussen transmuraal druk en vaatstraal wordt gemodelleerd met een niet-lineaire functie die de elasticiteit en spanning van het vaatweefsel weergeeft.
• Gekoppelde vergelijkingen: Het resultaat is een gesloten systeem van twee niet-lineaire PDE's voor de gemiddelde snelheid en de vaatstraal (of dwarsdoorsnede), vergelijkbaar met de Saint-Venant-vergelijkingen voor ondiepe waterstroming.

B. Biochemisch Model (Niet-gladde Dynamica)

• Koppeling: De kleppen in de lymfangion worden aangestuurd door de concentraties van $Ca^{2+}$ en NO.
• Vereenvoudiging: Voor de analyse wordt de vaatstraal $R$ als constant beschouwd om de biochemische oscillator te isoleren.
• Niet-gladde ODE's: De dynamiek wordt beschreven door een systeem van twee niet-lineaire gewone differentiaalvergelijkingen (ODE's) met Heaviside-functies (stapfuncties). Deze functies introduceren discontinuïteiten die corresponderen met de drempelwaarden voor klepopeningen en samentrekkingen:
  • $Ca^{2+}$-dynamiek: Bevat afbraak, productie en een term die afhankelijk is van de NO-concentratie en een drempelwaarde voor samentrekking.
  • NO-dynamiek: Bevat afbraak en productie die afhankelijk is van de wandschuifspanning (die op zijn beurt gekoppeld is aan de $Ca^{2+}$-concentratie).
• Filippov-systemen: Het systeem wordt geanalyseerd als een stuksgewijs glad (piecewise-smooth) dynamisch systeem, verdeeld in gebieden ($\Sigma^-$ en $\Sigma^+$) gescheiden door een discontinuïteitsgrens.

3. Belangrijkste Bijdragen

1. Unificerend Raamwerk: Voor het eerst wordt een systematische afleiding van zowel bloed- als lymfestroming gegeven vanuit dezelfde fundamentele fysische principes, waarbij de verschillen in rheologie en wandmechanica via parameters worden opgelost.
2. Koppeling Mechanica en Biochemie: Het model koppelt de vloeistofmechanica expliciet aan de kinetiek van calcium en NO, waardoor de terugkoppeling tussen mechanische stress (schuifspanning) en chemische signalering wordt vastgelegd.
3. Analyse van Niet-Gladde Systemen: De auteurs bieden een uitgebreide bifurcatie-analyse van het resulterende niet-gladde dynamische systeem. Ze identificeren exacte parameterregimes waarin stabiele limietcycli (oscillaties) ontstaan.
4. Wiskundige Validatie: Het bewijs van het bestaan van een limietcyclus in het vereenvoudigde model biedt theoretische validatie voor de waargenomen ritmische samentrekkingen in het lymfestelsel.

4. Resultaten

• Vloeistofmodel: De afgeleide PDE's beschrijven succesvol de stroming in zowel bloedvaten als lymfevaten. Het model kan zowel passieve stroming als actief gepompte stroming simuleren door de juiste constitutieve wetten en randvoorwaarden toe te passen.
• Oscillaties in het Biochemische Model:
  • Het systeem van twee niet-gladde ODE's vertoont relaxatie-oscillaties onder specifieke voorwaarden.
  • De auteurs bewijzen het bestaan van een stabiele limietcyclus wanneer de evenwichtspunten van beide subsystemen (voor $Ca^{2+}$ en NO) "virtueel" zijn (d.w.z. ze liggen in het gebied waar ze niet fysiek toegankelijk zijn, maar de trajecten worden er wel naartoe getrokken totdat ze de grens kruisen).
  • De voorwaarden voor deze oscillaties zijn: $\frac{\beta}{\alpha} > 1$ en $\frac{\beta + \gamma}{\alpha(\zeta + 1)} < 1$ (waarbij $\alpha, \beta, \gamma, \zeta$ dimensieloze parameters zijn die voortvloeien uit de fysische constanten).
• Bifurcatie-analyse: Er wordt gedetailleerd onderzocht hoe het systeem reageert op veranderingen in parameters (zoals $\beta$ en $\gamma$). Er worden verschillende scenario's geïdentificeerd, waaronder:
  • Grensevenwichtsbifurcaties (Boundary Equilibrium Bifurcations).
  • Overgangen tussen regelmatige evenwichten, virtuele evenwichten en pseudo-evenwichten.
  • Het verschijnen van een "niet-gladde vouw" (non-smooth fold) versus persistentie van de limietcyclus, afhankelijk van de verhouding tussen de parameters.
• Fasevlakken: De analyse toont aan dat de trajecten in het fasevlak $(N, C)$ een gesloten cyclus vormen die de ritmische cyclus van samentrekking (systole) en ontspanning (diastole) van de lymfangion nabootst.

5. Betekenis en Toekomstperspectief

• Fysiologische Validatie: Het model biedt een deterministische verklaring voor de spontane, ritmische pompactie van het lymfestelsel, die eerder vaak als puur empirisch werd beschouwd.
• Medische Toepassingen: Het raamwerk biedt een basis voor het bestuderen van pathologische aandoeningen (zoals lymfoedeem) waarbij de pompmechanica faalt, en voor het ontwerpen van therapeutische interventies.
• Wiskundige Vooruitgang: De toepassing van de theorie van niet-gladde dynamische systemen (Filippov-systemen) op de fysiologie is een innovatieve stap die complex gedrag zoals klepopeningen en -sluitingen nauwkeuriger kan modelleren dan traditionele gladde benaderingen.
• Toekomstig Werk: De auteurs plannen om de volledig gekoppelde PDE-ODE-modellen te analyseren (waarbij de vaatstraal niet langer constant is) en efficiënte numerieke methoden te ontwikkelen om dit niet-lokale, niet-lineaire systeem op te lossen.

Conclusie: Dit artikel levert een robuust en wiskundig onderbouwd model dat de complexe interactie tussen vloeistofmechanica en biochemische regulatie in het lymfestelsel beschrijft. Het bewijst dat de intrinsieke oscillatorische aard van de lymfe-transport een natuurlijk gevolg is van de gekoppelde dynamica van calcium en stikstofmonoxide binnen een niet-glad systeem.