[论文解读] On Stretching, Bending, Shearing and Twisting of Actin Filaments II: Multi-Resolution Modelling

Ravinda S. Gunaratne, Carlos Floyd+3 more|arXiv (Cornell University)|Mar 2, 2022

Cellular Mechanics and Interactions被引用 4

一句话总结

该论文提出了一种针对肌动蛋白丝的多分辨率模型，结合了基于单体的椭球体详细表示与粗粒化Cosserat杆模型，以高效模拟拉伸、弯曲、剪切和扭转行为。通过引入源自全原子模拟的曲率依赖刚度和扭转变形-拉伸耦合机制，该方法准确捕捉了扭转行为，尤其在过度扭转时，丝状结构显著收缩——这优于标准Cosserat模型，后者无法再现这种非对称性。

ABSTRACT

We present a multi-resolution methodology for modelling F-actin filaments. It provides detailed microscopic information at the level of individual monomers at a lower computational cost by replacing the monomer-based model in parts of the simulated filament by a rod-based macroscopic model. In the monomer-based description, G-actin is represented by ellipsoids bound at the surface in a double helical configuration to form F-actin. The rod-based model is coarser, in which F-actin is described using a Cosserat model, as seen in the preceding paper [arXiv:2112.01480]. The multi-resolution methodology is illustrated using three case studies, designed to test the properties of F-actin under stretching, bending, shearing and twisting. The methodology is especially suited for situations where filaments are subject to bending deformations. We investigate the limitations of using the standard Cosserat model to capture the complete torsional behaviour of F-actin, presenting its extensions which account for curvature dependent rigidities and a twist-stretch coupling to improve accuracy of the overall multi-resolution scheme.

研究动机与目标

• 开发一种计算高效的F-actin多分辨率模型，结合微观单体细节与宏观杆体力学。
• 解决标准Cosserat模型在模拟F-actin丝在过度扭转与欠扭转状态下非对称扭转响应方面的局限性。
• 基于椭球体模拟数据，推导并实现一种扭转变形-拉伸耦合机制，以提升机械保真度。
• 在保留关键区域机械细节的同时，实现复杂机械载荷下F-actin丝的大规模模拟。
• 通过拉伸、弯曲、剪切和扭转的案例研究验证模型，重点关注扭转行为与曲率依赖刚度。

提出的方法

• 模型采用椭球体表示F-actin的单体结构，沿双螺旋排列，显式表示G-肌动蛋白亚基及其相互作用。
• 在感兴趣区域，采用基于椭球体模型参数导出的Cosserat杆模型，以降低计算成本。
• 引入曲率依赖的扭转刚度函数 ξ(κ̂L,3)，通过根据局部曲率调整单元段的原始长度，模拟由扭转引起的伸长或缩短。
• 使用梯形求积算子 Âih 从离散Voronoi区域曲率估算段曲率，实现机械性能的精确局部适应。
• 基于不同扭转条件下椭球体模拟的数据，参数化扭转变形-拉伸耦合机制，采用六次多项式拟合 p6 到 p0 来定义 ξ(κ̂L,3)。
• 将这些校正项整合到Cosserat杆的正常加速度计算中，实现对扭转变形的动态响应。

实验结果

研究问题

• RQ1 如何设计一种多分辨率模型，以高计算效率模拟F-actin在拉伸、弯曲、剪切和扭转下的行为？
• RQ2 为何标准Cosserat模型无法准确捕捉F-actin丝在过度扭转时的扭转行为？
• RQ3 F-actin中的扭转变形-拉伸耦合机制的本质是什么？如何定量建模以提升机械保真度？
• RQ4 曲率依赖刚度能否提升Cosserat杆模型在模拟肌动蛋白丝力学中的准确性？
• RQ5 在不同扭转状态下，丝的长度如何变化？这种非对称性能否在粗粒化模型中被捕捉？

主要发现

• 标准Cosserat模型在扭转载荷下表现欠佳，无法再现椭球体模型中观察到的非对称长度变化。
• 在过度扭转（如240°）时，丝收缩了原始长度的4.34%；而在欠扭转240°时，收缩量仅为1.26%，表明存在非对称机械响应。
• 通过六次多项式拟合成功参数化了扭转变形-拉伸耦合函数 ξ(κ̂L,3)，其系数为 p6 = -2.6662×10⁶ 至 p0 = -5.0666×10⁻⁴，实现了扭转下长度的精确适应。
• 引入曲率依赖刚度与扭转变形-拉伸耦合显著提升了Cosserat模型在扭转区域的机械保真度。
• 多分辨率模型在远低于全原子或单体基模型的计算成本下，实现了复杂机械行为的精确模拟。
• 案例研究证实，该模型能够准确再现关键机械响应，包括弯曲动力学与非对称扭转变形，验证了其在大规模生物模拟中的适用性。