⚛️ biophysics

A mechanical model of multicellular remodelling in epithelial monolayers

Dit artikel introduceert een Dynamisch Referentiekader (DRF) binnen een off-lattice celcentrummodel om subcellulaire remodellering in epitheliale monolagen te simuleren, wat het reproduceren van multiscale mechanische responsen en de analyse van hoe de deformatiegeschiedenis de weefselherstel beïnvloedt mogelijk maakt.

Oorspronkelijke auteurs: Osborne, J. M., Van Ammers, R., Davit, Y., Gavaghan, D. J., Byrne, H.

Gepubliceerd 2026-09-14

📖 5 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Osborne, J. M., Van Ammers, R., Davit, Y., Gavaghan, D. J., Byrne, H.

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ⚕️ Dit is een AI-gegenereerde uitleg van een preprint die niet peer-reviewed is. Dit is geen medisch advies. Neem geen gezondheidsbeslissingen op basis van deze inhoud. Lees de volledige disclaimer

De buitenste lagen van ons lichaam en de bekledingen van onze interne organen zijn gemaakt van epitheelweefsels, vellen cellen die fungeren als beschermende barrières en poortwachters voor het lichaam. Deze weefsels zijn niet statisch; ze vormen zichzelf voortdurend opnieuw, rekken uit en ontspannen naarmate organen groeien, genezen of reageren op hun omgeving. Wetenschappers weten al lang dat wanneer deze vellen worden uitgerekt of samengedrukt, ze niet simpelweg terugveren of in één uniforme beweging uitrekken. In plaats daarvan reageren ze in twee verschillende fasen: een snelle, elastische ruk die binnen seconden plaatsvindt, gevolgd door een veel langzamere, kruipende aanpassing die minuten of zelfs uren kan duren. Deze langzame verandering is een vorm van hermodellering, waarbij de interne machinerie van de cellen — specifiek de eiwitsteigers en de kleverige verbindingen die ze bij elkaar houden — zich reorganiseert om de nieuwe vorm te accommoderen. Het begrijpen van dit tweesnelheidsgedrag is cruciaal, omdat het onthult hoe weefsels hun integriteit behouden onder stress en hoe ze kunnen falen bij ziekte, terwijl bestaande computermodellen er moeite mee hebben om deze complexiteit te vangen zonder onmogelijk ingewikkeld te worden.

Een team van onderzoekers heeft nu een nieuwe manier ontwikkeld om dit gedrag te simuleren, door een digitaal model te creëren dat zowel de snelle ruk als de langzame kruip van epitheelvellen vastlegt. In plaats van te proberen elke kleine proteïne binnen elke cel te modelleren, wat enorme rekenkracht zou vereisen, introduceerden de wetenschappers een concept dat ze een "dynamisch referentiekader" noemen. Stel je het weefsel voor als een verzameling cellen die verbonden zijn door veren. In oudere modellen was de lengte van deze veren vast; als je het weefsel uitrekte, rekten de veren uit en behielden ze die nieuwe lengte voor altijd, tenzij de cellen uit elkaar bewogen. Dit nieuwe model voegt een tweede, onzichtbare laag veren toe die de interne geheugen van de cel over zijn vorm vertegenwoordigt. Wanneer het weefsel wordt uitgerekt, bewegen de zichtbare cellen snel, maar deze onzichtbare laag verschuift langzaam zijn eigen doelbreedte om de nieuwe realiteit aan te passen. Deze verschuiving vertegenwoordigt het biologische proces waarbij cellen hun interne structuren reorganiseren om zich aan te passen aan de stress die ze voelen.

De onderzoekers testten dit idee door computersimulaties van kruip- en spanningsrelaxatie-experimenten uit te voeren, wat standaardmanieren zijn om te testen hoe materialen zich onder kracht gedragen. In een kruipexperiment wordt een constante ruk uitgeoefend op de rand van een digitale weefselplaat, en de onderzoekers observeerden hoe ver het over een bepaalde tijd uitrekte. In een spanningsrelaxatie-experiment wordt het weefsel tot een vaste breedte uitgerekt en daar gehouden, terwijl de kracht die nodig is om het uitgerekt te houden wordt gemeten. Wanneer het model alleen de snelle, elastische respons bevatte, gedroeg het weefsel zich als een eenvoudige rubberband: het rekte snel uit en stopte dan, of de kracht daalde snel en bleef daarna stabiel. Echter, toen de onderzoekers het dynamische referentiekader inschakelden om subcellulaire hermodellering mogelijk te maken, veranderde de simulatie drastisch. Het weefsel reageerde nog steeds eerst snel, maar het bleef vervolgens langzaam uitrekken over de tijd, of de kracht bleef verder dalen, waarmee het tweefasige gedrag dat in echte biologische experimenten wordt gezien, perfect werd nagebootst.

Door de instellingen van dit nieuwe model aan te passen, ontdekte het team dat de snelheid van deze langzame hermodellering afhangt van hoe gevoelig de cellen zijn voor stress en hoe snel ze hun interne "geheugen" bijwerken. Ze ontdekten dat als de cellen traag zijn in het aanpassen, het weefsel zich voornamelijk gedraagt als een eenvoudig elastisch materiaal. Maar als de cellen gevoelig zijn en hun interne structuur in de loop van de tijd bijwerken, vertoont het weefsel de complexe, multisnelheidsrespons die in de natuur wordt waargenomen. Cruciaal is dat het model ook aantoonde dat de geschiedenis van de spanning ertoe doet. Als een weefsel gedurende een lange tijd onder spanning wordt gehouden, past de interne structuur zich aanzienlijk aan. Wanneer de spanning eindelijk wordt losgelaten, veert het weefsel niet onmiddellijk terug; in plaats daarvan herstelt het langzaam, en de snelheid van dit herstel hangt af van hoe lang het onder spanning stond. Dit is een belangrijke bevinding, omdat veel oudere computermodellen voorspelden dat het herstel hetzelfde zou zijn, ongeacht hoe lang het weefsel was uitgerekt, een resultaat dat in strijd is met realiteitsobservaties.

De studie bevestigt dat de langzame, kruipende veranderingen in de vorm van weefsels worden gedreven door de voortdurende reorganisatie van de interne componenten van de cellen, in plaats van door het delen of verplaatsen van cellen naar nieuwe posities. De onderzoekers hebben aangetoond dat dit gedrag gevangen kan worden zonder elke individuele molecuul binnen een cel te hoeven simuleren, wat het mogelijk maakt om grote weefselvellen met duizenden cellen te bestuderen. Hun werk suggereert dat het vermogen van weefsels om zich de mechanische spanningen uit het verleden te herinneren en aan te passen, een fundamentele eigenschap is van hun mechanische ontwerp. Door een kader te bieden dat de microscopische reorganisatie van cellen koppelt aan het macroscopische gedrag van weefsels, biedt dit nieuwe model een duidelijk pad naar het begrijpen van hoe organen zich ontwikkelen, hoe ze genezen en hoe ze kunnen reageren op de mechanische krachten van ziekte. De code en gegevens die zijn gebruikt om deze bevindingen te genereren, zijn nu beschikbaar voor andere wetenschappers om te gebruiken, zodat de gemeenschap meer gedetailleerde en nauwkeurige modellen van levende weefsels kan bouwen.

Technische Samenvatting: Een mechanisch model van multicellulaire remodellering in epitheliale monolagen

Probleemstelling
Epitheliale monolagen, die de basis vormen van veel zoogdierorganen, vertonen complexe mechanische reacties op externe krachten. Experimenteel bewijs (bijv. Harris et al. [6]) toont aan dat deze weefsels reageren op vervorming over meerdere tijdschalen: een snelle elastische respons (seconden) gevolgd door een langzamere relaxatie of kruip (minuten tot uren). Bestaande wiskundige modellen worstelen om dit gedrag te reproduceren zonder celproliferatie of herrangschikking aan te nemen, wat op deze korte tijdschalen verwaarloosbaar is.

Continuümmodellen (bijv. Standard Linear Solid) gaan doorgaans uit van een enkele karakteristieke tijdschaal.
Complexe visco-elastische modellen (bijv. meerdere Maxwell-elementen) kunnen meerdere tijdschalen produceren, maar zijn vaak fenomenologisch, waarbij parameters moeilijk te koppelen zijn aan specifieke subcellulaire mechanismen.
Celgebaseerde modellen (bijv. Vertex Dynamics, Overlapping Spheres) vertrouwen over het algemeen op vaste interactiewetten. Wanneer proliferatie en herrangschikking worden genegeerd, vertonen ze ook vaak een respons met een enkele tijdschaal. Hoewel Subcellular Element (SE) modellen multi-tijdschaalgedrag kunnen vangen via interne reorganisatie, beperkt hun hoge computationele kosten hen tot zeer kleine systemen (<10 cellen), wat hen onpraktisch maakt voor simulaties op weefselniveau.

Methodologie
De auteurs stellen een innovatieve uitbreiding voor van off-lattice celcentrummodellen, specifiek gebruikmakend van een Voronoi-tessellatie (VT) raamwerk. De kerninnovatie is de introductie van een Dynamisch Referentiekader (DRF).

Modelarchitectuur: Het weefsel wordt gerepresenteerd door drie sets celcentra:
1. Fysieke cellen ( $r_i$ ): De werkelijke, vervormende posities van de cellen.
2. Referentiecellen ( $\rho_i$ ): Een dynamische configuratie die over de tijd evolueert en de "geremodelleerde" staat van het weefsel representeert.
3. Statische cellen ( $s_i$ ): Een vaste initiële configuratie die dient als basislijn.
Beherende Vergelijkingen:
- Fysieke Beweging: Fysieke cellen bewegen volgens een krachtbalans tussen viskeuze weerstand, externe krachten en veerkrachten die fysieke cellen verbinden met de referentie-cellen (Gelijk. 2). De veerkracht hangt af van het verschil tussen de fysieke afstand ( $||r_{ij}||$ ) en de referentiedistantie ( $||\rho_{ij}||$ ).
- Referentie Evolutie: De referentiecellen evolueren volgens een aparte bewegingsvergelijking (Gelijk. 3). Ze worden gedreven door twee concurrerende krachten:
  1. Een herstelkracht die hen naar de statische configuratie ( $s_i$ ) trekt.
  2. Een remodelleringskracht die hen naar de tijdsgemiddelde fysieke configuratie ( $\bar{r}_{ij}$ ) over een geheugenvenster $T_M$ trekt.
Dimensieloos maken: Het systeem wordt dimensieloos gemaakt om sleutelparameters te definiëren:
- $\alpha$ (Remodelleringsratio): De relatieve sterkte van de remodelleringskracht versus de herstelkracht.
- $\eta$ (Remodelleringssnelheid): De snelheid waarmee het referentiekader evolueert ten opzichte van het fysieke kader.
- $T_M$ (Remodelleringsgeheugen): Het tijdsvenster waarover de fysieke vervorming wordt gemiddeld.
Simulatieopstelling: Het model simuleert 2D epitheliale monolagen (105 cellen) onder twee experimentele protocollen:
1. Kruiprelaxatie: Constante kracht toegepast; vervorming (strain) wordt geregistreerd.
2. Spanningsrelaxatie: Vervorming wordt tot een maximum gebracht en vastgehouden; spanning (stress) wordt geregistreerd.
Analyse: Tijdreeksen worden gefit met sommen van exponenten. Het aantal tijdschalen wordt bepaald door of een enkele exponent een parsimonieuze fit biedt (fout $E_1 < 0,05$ ) of dat meerdere termen vereist zijn.

Belangrijkste Bijdragen

Dynamisch Referentiekader (DRF): Het artikel introduceert een mechanisme om de instantane vervorming van een cel te ontkoppelen van zijn mechanische referentieconfiguratie. Dit stelt het model in staat om subcellulaire remodellering (cytoskelet reorganisatie, myosine-omzet, remodellering van adhesiebindingen) te simuleren zonder expliciet de subcellulaire componenten te modelleren, waardoor de computationele kosten van SE-modellen worden vermeden.
Emergentie van Multi-tijdschalen: De auteurs demonstreren dat het DRF-raamwerk van nature multi-tijdschaalresponsen genereert (snelle elastische + trage relaxatie/kruip) in celgebaseerde modellen zonder proliferatie.
Parametersensitiviteit: De studie kwantificeert hoe de remodelleringsparameters ( $\alpha, \eta, T_M$ ) de overgang tussen single- en multi-tijdschaalgedrag controleren.
Geschiedenisafhankelijkheid: Het model reproduceert succesvol de experimentele observatie dat weefselherstel bij het vrijlaten afhankelijk is van de duur van de vervorming, een fenomeen dat standaard elastische modellen niet kunnen vatten.

Resultaten

Kruipexperimenten:
- Zonder remodellering ( $\alpha=0$ ) vertoont het weefsel een single-tijdschaal respons (snelle elastische vervorming gevolgd door een plateau).
- Met remodellering ( $\alpha > 0$ ) ontstaat een tweede, tragere tijdschaal naarmate het referentiekader evolueert, wat verdere vervorming mogelijk maakt.
- De emergentie van meerdere tijdschalen is robuust over verschillende toegepaste krachtmagnitudes, maar is zeer gevoelig voor $\alpha$ (remodelleringssterkte) en $\eta$ (remodelleringssnelheid). Het vergroten van het geheugenvenster ( $T_M$ ) bevordert eveneens multi-tijdschaalgedrag.
Spanningsrelaxatieexperimenten:
- Vergelijkbaar met kruip introduceert remodellering een secundaire relaxatiefase.
- De ramp-tijd ( $T_R$ ) fungeert als controleparameter; tussenliggende ramp-tijden leggen onderscheidende post-ramp relaxatietijdschalen bloot, terwijl zeer korte of zeer lange rampen deze kunnen maskeren.
- Het vergroten van $T_M$ breidt de parameterruimte uit waarin multi-tijdschaalresponsen worden waargenomen.
Herstel van Vervorming:
- In modellen zonder remodellering is de snelheid van herstel na het vrijlaten van een vastgehouden vervorming onafhankelijk van de vasthoudduur ( $T_H$ ).
- In het DRF-model is de herstelsnelheid omgekeerd evenredig met $T_H$ . Het model voorspelt dat langere vasthoudtijden het referentiekader verder laten vervormen (wat de geaccumuleerde subcellulaire remodellering representeert), wat resulteert in een lagere spanning bij vrijlating en een trager daaropvolgend herstel. Dit gedrag komt overeen met de algemene experimentele bevinding van Harris et al. [6] dat herstel afhankelijk is van de duur van de vervorming, terwijl de specifieke $1/T_H$ relatie een kwantitatief resultaat is van het voorgestelde model.

Betekenis en Claims
Het artikel claimt een computationeel efficiënt raamwerk te bieden voor het bestuderen van subcellulaire remodellering op weefselniveau. Door standaard celcentrummodellen uit te breiden met een DRF, kunnen de auteurs de complexe, multi-tijdschaal mechanische gedragingen reproduceren die worden waargenomen in epitheliale monolagen, zonder terug te vallen op prohibitief dure subcellulaire simulaties of fenomenologische continuümmodellen met onmeetbare parameters.

De auteurs stellen dat dit raamwerk "gedetailleerdere beschrijvingen en analyses van de ontwikkeling en functie van biologische weefsels" mogelijk maakt. Zij benadrukken dat het model modulair is: het DRF-concept kan worden toegepast op andere off-lattice modellen (bijv. Vertex Dynamics) en kan worden uitgebreid naar 3D of kan processen zoals celgroei en differentiatie bevatten door tussenliggende referentietoestanden toe te voegen. Het werk suggereert dat de geschiedenis van de vervorming, gemedieerd door subcellulaire omzet, een cruciale factor is in weefselmechanica die nu systematisch kan worden onderzocht met deze aanpak.

Meer zoals dit