《PLOS Computational Biology》：Cross-bridge model for predicting muscle short-range stiffness during movement
编辑推荐：
肌肉骨骼模拟可深入揭示肌肉骨骼系统特性如何影响日常运动的生物力学。其中一项特性是肌肉对拉伸的初始抵抗，即短程刚度（short-range stiffness, SRS），它是响应外部扰动稳定运动的关键。现有肌肉骨骼模拟采用缺乏激活依赖刚度特性的现象学Hill型
肌肉骨骼模拟可深入揭示肌肉骨骼系统特性如何影响日常运动的生物力学。其中一项特性是肌肉对拉伸的初始抵抗，即短程刚度（short-range stiffness, SRS），它是响应外部扰动稳定运动的关键。现有肌肉骨骼模拟采用缺乏激活依赖刚度特性的现象学Hill型模型，难以捕捉短程刚度；现有模拟亦未捕捉肌肉缩短后短程刚度的历史依赖性降低，即肌肉触变性（muscle thixotropy）。交叉桥（cross-bridge, XB）模型虽能复现短程刚度，但何种模型特性为捕捉其历史依赖性所必需仍不清楚。研究人员利用已有11条透化大鼠比目鱼肌纤维数据集，测试了各类交叉桥模型在广泛行为相关长度变化与激活水平下复现经验短程刚度及其历史依赖变化的能力。研究人员以缩短前后观测短程刚度之比量化肌肉触变性，计算各模型预测短程刚度比与实测刚度比的根均方偏差（root-mean-square deviation, RMSD）。研究发现，仅当纳入细肌丝与粗肌丝的协同激活时，交叉桥模型方能在小幅与大幅历史依赖刚度降低条件下捕捉短程刚度变化（RMSD ≤ 0.1）；反之，Hill型模型及无协同肌丝激活的交叉桥模型低估短程刚度，且在大幅历史依赖刚度降低条件下无法捕捉其变化（RMSD > 0.2）。采用高斯近似求解交叉桥分布可得相似结果，而计算成本约降低八倍。因此研究人员提议将带协同肌丝激活的高斯近似交叉桥模型植入肌肉骨骼模拟，以改进运动中短程刚度的预测。
研究背景与问题提出
骨骼肌既是肌肉骨骼系统的动力源也是制动器，其对拉伸的初始快速抵抗即短程刚度（short-range stiffness, SRS）在扰动稳定中发挥关键作用，而缩短后的刚度降低称为肌肉触变性（muscle thixotropy）。现有肌肉骨骼模拟普遍采用现象学Hill型模型，该类模型缺乏激活依赖刚度且无法描述触变性的历史依赖；即便加入经验扩展也仅区分有无先前运动，不能刻画激活水平、恢复时间与振幅的连续影响。早期两态交叉桥（cross-bridge, XB）模型虽能定性复现短程刚度，但在次最大激活下无法重现触变性随激活增强而增大的经验规律。因此，研究人员旨在从透化大鼠比目鱼肌纤维（n=11）的拉伸-缩短-拉伸数据中，甄别出能跨激活水平与恢复时间捕捉短程刚度历史依赖的交叉桥模型，并将其工程化以便植入全身运动模拟。该文发表于《PLOS Computational Biology》。
主要关键技术方法
研究人员以已有11条透化成年雌性Sprague-Dawley大鼠比目鱼肌慢肌纤维的体外拉伸-缩短-拉伸数据集为基准，在Matlab中拟合与验证四类交叉桥模型：两态无协同、两态细肌丝协同、三态细+粗肌丝协同（含super-relaxed态）、四态再加forcibly-detached态；同时以带/不带串联弹性（series compliance）的Hill型模型作基准。交叉桥应变分布分别用离散分箱法与高斯近似法求解Huxley偏微分方程。参数分固定文献值与在两组拟合试验（大振幅短恢复、无预处理参考）上拟合两部分；以首10 ms力-长斜率定短程刚度，以测试拉伸比参考拉伸的刚度比定量触变性，用根均方偏差（RMSD）与AIC比较。
研究结果
交叉桥动力学单独不足以解释次最大激活拉伸-缩短力轨迹：Hill无串联弹性模型拟合差；Hill带串联弹性与两态无协同XB模型在最大激活可追随拉伸力上升，但在次最大激活下既低估条件拉伸初始力率，又无法复现测试拉伸相对条件拉伸的短程刚度下降。
加入协同肌丝动力学改善次最大激活初始力率：两态XB coop与三态XB coop在最大激活下改善缩短后等长力发展；在次最大激活下提升条件拉伸力率并复现测试拉伸刚度下降，但三态模型等长恢复偏慢、末段力略高。
建模协同肌丝动力学提升拉伸-缩短-拉伸力轨迹精度，而额外forcibly-detached态影响甚微：纤维平均RMSD随模型复杂度下降，三态/四态XB coop最低；Hill类与两态XB在激活、振幅、恢复时间扫描中RMSD较高。
建模协同肌丝激活对跨条件捕捉短程刚度变化至关重要：仅三态/四态XB coop复现相对短程刚度随pCa呈U形、随振幅增大而降、随恢复时间延长而升（SRS RMSD ≤0.1）；Hill与两态XB误预测为单调下降或恒定。
带协同肌丝激活的交叉桥模型预测行为相关条件下的肌肉触变性：在大幅刚度降低试验子集，四态XB coop SRS RMSD=0.101最低；全84试集合三态/四态XB coop力RMSD≈0.103且AIC最优；模型还正确预测正常姿势摆动（振幅>3%、恢复<0.5 s）下刚度降>30%，及长恢复时触变性消失。
高斯近似交叉桥分布对力轨迹影响小但计算成本降一个数量级：高斯近似三态XB coop单步耗时约0.000017 s，变步长积分较特征线法同精度快约8.6倍，力误差<0.01% F0。
讨论部分总结
研究人员指出，常用Hill型模型即使带串联弹性也只能由指数力-长关系产生伪历史依赖，无法解释次最大激活触变性最强这一经验事实；而单纯 actin-myosin 交互的两态XB同样失败，关键在细肌丝协同（thin filament cooperativity）与粗肌丝协同（thick filament cooperativity，力依赖的super-relaxed?disordered-relaxed转换）。前者在低钙下因交叉桥脱附削弱细肌丝活化而延长刚度恢复，后者经缩短致力下降使肌球蛋白头进入super-relaxed态进一步拖慢再附着。四态forcibly-detached态在当前慢速短拉伸数据集未显著改善力平台，AIC不支持其引入。高斯近似因近期模型采用高斯附着函数与指数脱附函数而比Zahalak当初假设更自洽，精度损失可忽略却换得约8倍加速。作者明示模型未含肌联蛋白（titin）动态、肌节非均质与丝间距变化，故仍存力轨迹偏差；此类机制可能是残余力增强/抑制与rheopexy的补充来源。将协同交叉桥模型植入全身模拟可弥补Hill型在偏心收缩与姿势扰动下的短板。
结论部分翻译
本研究旨在依据骨骼肌短程刚度历史依赖降低（即肌肉触变性）的捕捉能力甄别肌肉模型。研究发现，纳入肌丝协同动力学可改善历史依赖短程刚度降低的预测；采用交叉桥分布的高斯近似在模型精度与计算需求间取得良好折衷，便于植入肌肉骨骼模拟。基于带协同肌丝动力学的交叉桥肌肉模型的肌肉骨骼模拟，可阐明历史依赖短程刚度在全身运动中的作用。
