A multiscale cytoskeletal network model for shear rheological property and its evolutionary mechanism

Technische Samenvatting: Een multiscale cytoskeletnetwerkmodel voor afschuifrheologische eigenschappen en het evolutionaire mechanisme ervan

Probleemstelling
Het cytoskelet is een dynamisch biopolymeernetwerk dat essentieel is voor de cellulaire fysiologie, maar het mechanische gedrag ervan strekt zich uit over meerdere ruimtelijke en temporele schalen. Er bestaat een aanzienlijke uitdaging in het volledig vatten van complexe cellulaire responsen, omdat traditionele modellen moeite hebben met het koppelen van dynamische remodellering (gedreven door cross-linker kinetiek) aan visco-elastische dissipatiemechanismen (gedreven door filamenteigenschappen). Bestaande continuümmodellen slagen er vaak niet in om lokale niet-uniformiteit te vatten, terwijl discrete structurele modellen componenten vaak simplificeren als puur elastische elementen, waarbij ze tijdsafhankelijke eigenschappen en dynamische cross-linker evolutie verwaarlozen. Bovendien zijn veel discrete modellen computationeel beperkt tot korte tijdschalen (≤ 60 s), wat hun vermogen om langetermijn rheologische gedragingen, zoals steady-state kruip, te voorspellen beperkt. Deze studie adresseert de behoefte aan een model dat de bio-chemo-mechanische eigenschappen van lokale eiwitten, de visco-elasticiteit van actinefilamenten en hun vervormingsstaten integreert om de afschuifrheologische eigenschappen over uitgebreide tijdschalen (tot 1.000 seconden) te voorspellen.

Methodologie
De auteurs ontwikkelden een multiscale eindige-elementenmodel (FE-model) gebaseerd op een representatieve eenheid van eindige vierkante cellulaire elementen, afgeleid van elektronenmicroscopische observaties. Het model integreert drie verschillende schalen:

1. Filamentschaal: Actinefilamenten worden gemodelleerd als voorgespannen, driedimensionale visco-elastische staven met drie parameters. Hun constitutieve gedrag wordt beschreven door een tijdsafhankelijke relaxatiemodulus, die de overgang van buig-gedomineerde naar trek-gedomineerde vervorming vastlegt.
2. Noduschaal: Nodussen worden op basis van experimentele statistieken in twee typen onderverdeeld:
   • Fysiek verstrengelde nodussen: Vereenvoudigd als scharnierknopen die niet falen.
   • Chemisch gekruiste nodussen: Gemodelleerd als rotatieveren met kracht-gereguleerde dissociatie- en associatiesnelheden, gestuurd door het Bell-model.
3. Netwerkschaal: Een willekeurig statistisch netwerk wordt gegenereerd door initiële hoeken en willekeurige verplaatsingen aan een uniform rooster toe te voegen. Het model maakt gebruik van een rand-gemodificeerde FE-methode met een incrementeel iteratief algoritme. Dit omvat een geneste Newton-Raphson iteratieprocedure gecombineerd met een gesegmenteerde adaptieve tijdstapstrategie (logaritmische schaal) om het quasi-statisch evenwicht efficiënt op te lossen van milliseconden tot duizenden seconden.

De numerieke oplossing omvat het dynamisch bijwerken van de stijfheidsmatrix om veranderingen in de cross-linker bindingsdichtheid en de visco-elasticiteit van filamenten te reflecteren. Het model berekent de bulk afschuifvervorming, de schijnbare afschuifmodulus en de structurele Shannon-entropie om netwerkremodellering te karakteriseren.

Belangrijkste Bijdragen

• Verlengde Voorspellende Tijdschaal: De studie heeft het voorspellende vermogen van discrete netwerkmodellen succesvol uitgebreid tot ongeveer 1.000 seconden, waarmee de computationele beperkingen van eerdere Brownian/FE-Langevin dynamica modellen zijn overwonnen.
• Geïntegreerde Multiscale Koppeling: Het model koppelt uniek de bio-chemo-mechanische kinetiek van cross-linkers met de visco-elastische vervormingsovergang (van buiging naar trek) van individuele filamenten binnen één enkel numeriek kader.
• Randmodificatie: Een specifieke rand-gemodificeerde FE-methode is ontwikkeld om excentrische randkrachten op randfilamenten te behandelen, wat de nauwkeurigheid van afschuifspanning-voorspellingen verbetert.
• Pathologische en Farmacologische Simulatie: Het kader werd toegepast om cellulaire verweking door kanker en medicamenteuze interventies te simuleren, waarbij microscopische parameterveranderingen werden gekoppeld aan macroscopische rheologische uitkomsten.

Resultaten

• Rheologisch Gedrag: Het model reproduceerde het experimenteel waargenomen afschuifrheologische respons, inclus_ van short- en intermediate-term power-law kruipgedrag en een long-term strain plateau. De power-law exponent bleek af te nemen bij een toenemende filament volumefractie.
• Evolutionair Mechanisme: De simulatie onthulde een vierfasen temporele evolutie:
  1. Snel variërende fase: Snelle groei van de vervorming gedreven door filamentbuiging en hoge netwerkentropie.
  2. Transitie-fase: Vertragende groei naarmate cross-linkers dissociëren, waardoor de vervorming verschuift naar trek-gedomineerde modi en de entropie afneemt.
  3. Intermediate-term: Langzame kruip gedomineerd door de visco-elasticiteit van filamenten naarmate het netwerk stabiliseert.
  4. Long-term plateau: Steady-state gedrag bepaald door de inherente visco-elastische eigenschappen van de filamenten.
• Kankerachtige Transformatie: Simulaties van kankercellen (gekarakteriseerd door verminderde filament voorspanning, viscositeit en dichtheid, en verhoogde cross-linker dissociatie) toonden versnelde netwerkremodellering en een significant verlaagde schijnbare afschuifmodulus.
• Effect van Medicatie: Een single-parameter analyse toonde aan dat:
  • Het verhogen van de filament voorspanning effectief de steady-state afschuifmodulus herstelt.
  • Het verhogen van de filament viscositeit de tijd tot het bereiken van het plateau verlengt.
  • Het verminderen van de cross-linker dissociatiesnelheid de short-term modulus verhoogt en de piek van de structurele entropie vertraagt, wat wijst op een verlengd remodelleringsproces.

Betekenis
De auteurs stellen dat dit werk de kennis over de spatiotemporele evolutie en pathologische mechanismen van cellulaire mechanische responsen doet toenemen. Door te identificeren dat de synergie tussen de cross-linker chemo-mechanische koppeling en de filament vervormingsovergangen de netwerkremodellering en power-law responsen aanstuurt, biedt de studie een mechanistische verklaring voor waargenomen rheologische gedragingen. Het model biedt een theoretische basis voor het reguleren van polymeernetwerkprestaties, en suggereert specifiek dat synergetische interventies die gericht zijn op zowel de filament voorspanning als de cross-linker kinetiek noodzakelijk zijn om de door kanker geïnduceerde verweking effectief te keren. De auteurs merken ook het potentieel op om deze modelleringsbenadering uit te breiden naar andere polymeersystemen, zoals hydrogelen en colloïdale gels, hoewel zij de huidige beperkingen erkennen met betrekking tot het uitsluiten van actieve myosine-eenheden, fysieke verstrengeling mobiliteit en driedimensionale netwerkeffecten.