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

Resumo Técnico: Um Modelo de Rede Citoesquelética Multiescala para Propriedades Reológicas de Cisalhamento e Seu Mecanismo Evolutivo

Definição do Problema
O citoesqueleto é uma rede de biopolímeros dinâmica, essencial para a fisiologia celular, mas seu comportamento mecânico abrange múltiplas escalas espaciotemporais. Existe um desafio significativo em capturar de forma abrangente as respostas celulares complexas, pois os modelos tradicionais têm dificuldade em acoplar a remodelação dinâmica (impulsionada pela cinética dos ligantes cruzados/cross-linkers) com os mecanismos de dissipação viscoelástica (impulsionados pelas propriedades dos filamentos). Os modelos contínuos existentes frequentemente falham em capturar a não uniformidade local, enquanto modelos estruturais discretos frequentemente simplificam os componentes como elementos puramente elásticos, negligenciando propriedades dependentes do tempo e a evolução dinâmica dos ligantes cruzados. Além disso, muitos modelos discretos são computacionalmente limitados a escalas de tempo curtas (≤ 60 s), o que limita sua capacidade de prever comportamentos reológicos de longo prazo, como o creep (fluência) em estado estacionário. Este estudo aborda a necessidade de um modelo que integre as propriedades bio-quimico-mecânicas de proteínas locais, a viscoelasticidade dos filamentos de actina e seus estados de deformação para prever propriedades reológicas de cisalhamento ao longo de escalas de tempo estendidas (até 1.000 segundos).

Metodologia
Os autores desenvolveram um modelo de elementos finitos (FE) multiescala baseado em uma unidade representativa de elementos celulares quadrilaterais finitos, derivada de observações de microscopia eletrônica. O modelo integra três escalas distintas:

1. Escala do Filamento: Os filamentos de actina são modelados como hastes sólidas viscoelásticas de três parâmetros e pré-tensionadas. Seu comportamento constitutivo é descrito por um módulo de relaxação dependente do tempo, capturando a transição da deformação dominada pela flexão para a dominada pela tensão.
2. Escala do Nó: Os nós são categorizados em dois tipos baseados em estatísticas experimentais:
   • Nós fisicamente emaranhados: Simplificados como nós de articulação (hinge nodes) que não falham.
   • Nós quimicamente ligados (cross-linked): Modelados como molas rotacionais com taxas de dissociação e associação reguladas pela força, governadas pelo modelo de Bell.
3. Escala da Rede: Uma rede estatística aleatória é gerada introduzindo ângulos iniciais e deslocamentos aleatórios a uma grade uniforme. O modelo utiliza um método de FE com modificação de contorno por meio de um algoritmo iterativo incremental. Isso inclui um esquema de iteração Newton-Raphson aninhado combinado com uma estratégia de passo de tempo adaptativo segmentado (escala logarítmica) para resolver eficientemente o equilíbrio quase-estático de milissegundos a milhares de segundos.

A solução numérica envolve a atualização dinâmica da matriz de rigidez para refletir as mudanças na densidade de ligação dos ligantes cruzados e na viscoelasticidade do filamento. O modelo calcula a deformação de cisalhamento volumétrica, o módulo de cisalhamento aparente e a entropia de Shannon estrutural para caracterizar a remodelação da rede.

Principais Contribuições

• Escala de Tempo Preditiva Estendida: O estudo expandiu com sucesso a capacidade preditiva dos modelos de rede discretos para aproximadamente 1.000 segundos, superando as limitações computacionais de modelos anteriores de dinâmica Browniana/FE-Langevin.
• Acoplamento Multiescala Integrado: O modelo acopla de forma única a cinética bio-quimico-mecânica dos ligantes cruzados com a transição de deformação viscoelástica (flexão-para-tensão) de filamentos individuais dentro de um único arcabouço numérico.
• Modificação de Contorno: Um método específico de FE com modificação de contorno foi desenvolvido para lidar com forças de superfície excêntricas em filamentos de contorno, aumentando a precisão das previsões de tensão de cisalhamento.
• Simulação Patológica e Farmacológica: O arcabouço foi aplicado para simular o amolecimento canceroso e intervenções farmacológicas, vinculando mudanças de parâmetros microscópicos a resultados reológicos macroscópicos.

Resultados

• Comportamento Reológico: O modelo reproduziu a resposta reológica de cisalhamento observada experimentalmente, incluindo o comportamento de creep (fluência) de lei de potência de curto e médio prazo e um platô de deformação de longo prazo. O expoente da lei de potência mostrou-se decrescente com o aumento da fração volumétrica do filamento.
• Mecanismo Evolutivo: A simulação revelou uma evolução temporal de quatro estágios:
  1. Estágio de variação rápida: Crescimento rápido da deformação impulsionado pela flexão do filamento e alta entropia da rede.
  2. Estágio de transição: Crescimento desacelerado à medida que os ligantes cruzados se dissociam, deslocando a deformação para modos dominados por tensão e reduzindo a entropia.
  3. Médio prazo: Creep lento dominado pela viscoelasticidade do filamento conforme a rede se estabiliza.
  4. Platô de longo prazo: Comportamento de estado estacionário determinado pelas propriedades viscoelásticas inerentes dos filamentos.
• Transformação Cancerosa: Simulações de células cancerosas (caracterizadas por redução de pré-tensão, viscosidade e densidade do filamento, e aumento da dissociação do ligante cruzado) demonstraram uma remodelação acelerada da rede e uma diminuição significativa do módulo de cisalhamento aparente.
• Efeitos de Drogas: Uma análise de parâmetro único mostrou que:
  • O aumento da pré-tensão do filamento restaura efetivamente o módulo de cisalhamento em estado estacionário.
  • O aumento da viscosidade do filamento prolonga o tempo para atingir o platô.
  • A redução da taxa de dissociação do ligante cruzado aumenta o módulo de curto prazo e retarda o pico da entropia estrutural, indicando um processo de remodelação prolongado.

Significância
Os autores afirmam que este trabalho avança a compreensão da evolução espaciotemporal e dos mecanismos patológicos das respostas mecânicas celulares. Ao identificar que a sinergia entre o acoplamento quimico-mecânico do ligante cruzado e as transições de deformação do filamento governa a remodelação da rede e as respostas de lei de potência, o estudo fornece uma explicação mecanística para os comportamentos reológicos observados. O modelo oferece uma base teórica para regular o desempenho de redes poliméricas, sugerindo especificamente que intervenções sinérgicas visando tanto a pré-tensão do filamento quanto a cinética do ligante cruzado são necessárias para reverter efetivamente o amolecimento induzido pelo câncer. Os autores também observam o potencial de estender esta abordagem de modelagem para outros sistemas poliméricos, como hidrogéis e géis coloidais, embora reconheçam as limitações atuais relativas à exclusão de unidades de miosina ativas, mobilidade de emaranhamento físico e efeitos de rede tridimensional.