基于预测性仿真量化脑瘫儿童步态异常中肌肉骨骼损伤的贡献度

本研究针对脑瘫(CP)患儿多水平手术(MLS)疗效差异大、预测难的临床难题，通过个性化肌肉骨骼模型进行预测性仿真，首次在群体层面量化了肌肉无力、挛缩和骨骼畸形对步态运动学的独立与交互贡献。结果显示建模的肌肉骨骼损伤仅能解释平均17%（最高39%）的步态偏差，提示运动控制及未建模损伤（如足踝畸形）的关键作用。该方法为个体化治疗规划提供了新工具。

在儿童康复医学领域，脑性瘫痪(CP)作为最常见的神经发育障碍之一，其治疗始终面临着重大的挑战。这种由围产期脑部损伤引发的疾病，不仅导致运动控制障碍，还会引起肌肉挛缩、无力以及骨骼畸形等一系列肌肉骨骼系统问题，最终表现为异常步态。尽管多水平手术(MLS)是改善患儿行走能力的主要手段，通过延长肌腱、截骨矫形等方式针对性处理这些肌肉骨骼损伤，但其效果却难以预测，个体差异显著。过去二十年间，手术对步态运动学（kinematics）的改善效果似乎进入了平台期。究其根源，在于我们无法清晰解析：在复杂的损伤网络中，每一种特定的肌肉骨骼损伤究竟在多大程度上导致了患儿独特的步态异常？又是如何与运动控制缺陷相互作用的？传统的逆向仿真（inverse simulations）依赖于实验测得的步态数据来估算肌肉力量等参数，却难以回答“如果损伤不存在，步态会如何变化”这类因果性问题。

为了突破这一瓶颈，来自比利时鲁汶大学的研究团队在《Journal of NeuroEngineering and Rehabilitation》上发表了一项创新性研究。他们另辟蹊径，采用了一种名为“预测性仿真”（predictive simulations）的前沿技术。这种方法的核心在于，它不依赖于任何实验步态数据，而是将人体行走抽象为一个最优化控制问题——假设人体会以消耗某种与运动相关的成本（如代谢能量）最小的方式行走。通过求解这个复杂问题，计算机可以“生成”出一个全新的、理论上最优的步态模式。这就好比让一个虚拟的、具有特定损伤的“数字孪生”模型自己学习走路，通过观察其走出的步态与真实患儿步态之间的差异，来反推该损伤的实际贡献。

本研究共纳入了8名计划接受多水平手术的痉挛型脑瘫患儿。研究人员为每位患儿构建了多达8种不同个性化程度的肌肉骨骼模型，形成了一个精巧的“模型组合”。这些模型从仅经过身高体重缩放的基础通用模型（GEN），到逐步引入基于临床检查评估的肌肉无力（GEN_WEAK）、肌肉挛缩（GEN_CTR），以及基于磁共振成像(MRI)精准重建的髋、膝关节骨骼畸形（GEO）等因素，最终到整合所有损伤的完全个性化模型（GEO_FULL）。通过对每一种模型进行预测性仿真，并将产生的虚拟步态（主要是髋、膝、踝关节在矢状面和非矢状面的运动角度）与患儿实际步态分析数据进行比较，研究人员得以像进行一场“控制变量”的科学实验一样，剥离出每一种损伤的独立效应和它们之间的交互作用。

为开展研究，团队运用了几项关键技术：基于个性化肌肉骨骼模型的预测性仿真框架（PredSim），通过求解最优控制问题生成步态；模型个性化方法，即利用临床检查（如徒手肌力测试MMT、被动关节活动度）数据量化肌肉无力和挛缩，并利用MRI影像重建髋膝关节骨骼几何与肌肉路径；以及沙普利值（Shapley value）分析，用于量化不同损伤及其交互作用对步态偏差的贡献度。研究队列来自计划接受多水平手术的脑瘫患儿。

研究结果

肌肉骨骼损伤对步态缺陷的贡献

整体而言，引入损伤模型改善了仿真步态与实测步态在髋屈曲和膝屈曲角度上的一致性，但对髋内收、髋旋转和踝背屈的模拟改善有限。在单独建模的损伤中，肌肉挛缩对改善仿真吻合度的贡献最大。与基础模型（GEN）相比，单独建模挛缩（GEN_CTR）使髋屈曲和膝屈曲的均方根误差（RMSD）分别平均降低了3.9°和2.8°。单独建模无力（GEN_WEAK）主要影响髋部运动学，使髋屈曲的RMSD降低了1.9°。而单独建模骨骼畸形（GEO）也使髋屈曲的RMSD降低了3.9°。当合并建模挛缩和无力（GEN_FULL）时，改善效果更为明显，髋屈曲和膝屈曲的RMSD分别降低了5.7°和4.0°。

损伤间的交互作用

沙普利值分析揭示，损伤之间存在非线性的交互作用。肌肉挛缩是减少矢状面步态偏差（降低RMSD）最主要的因素（中位数贡献为降低6.4°），而单纯的无力或骨骼畸形贡献很小。然而，当挛缩与无力或骨骼畸形共同建模时，会产生额外的协同效应，导致RMSD的降低幅度大于各损伤独立贡献之和（额外降低约1.9°）。这表明，在存在挛缩的情况下，无力或骨骼畸形会“放大”挛缩对步态的不良影响。

个体间差异性

建模的肌肉骨骼损伤对步态异常的解释程度存在显著的个体差异。对于基线时仿真与实测步态差异较大（即步态异常更严重）的患儿，模型个性化能明显改善仿真效果（RMSD降低更多）。例如，对于步态呈蹲踞模式（crouch gait）的患儿CP6，完全个性化模型（GEO_FULL）能更好地模拟其增加的髋膝屈曲。反之，对于步态偏差较小的患儿CP3，引入损伤模型并未能一致性地改善仿真效果。此外，引入骨骼畸形模型（GEO）降低了不同患儿间仿真误差的变异性，提示骨骼畸形解释了一部分个体间步态差异的来源。

结论与意义

本研究通过系统的预测性仿真分析，得出核心结论：在计划接受多水平手术的脑瘫患儿中，当前方法所能建模的肌肉骨骼损伤（肌肉无力、挛缩及髋膝骨骼畸形）仅能解释一小部分（平均17%，最高不超过39%）的步态运动学偏差。这一有限且多变的贡献度，与多水平手术疗效的有限性和可变性相互印证。研究发现，肌肉挛缩是影响矢状面步态的最主要损伤因素，且其效应会被无力或骨骼畸形放大。然而，仍有大部分步态偏差无法通过肌肉骨骼损伤模型解释。

这强烈提示，运动控制障碍（如痉挛、选择性运动控制减少）以及未建模的骨骼畸形（如胫骨扭转、足部畸形） 在脑瘫步态异常中扮演着至关重要的角色。例如，足踝的复杂畸形在本研究的模型中没有得到充分体现，而这在脑瘫患儿中极为普遍。此外，当前基于临床检查数据的模型个性化方法可能存在精度限制。

尽管如此，这项研究标志着在理解脑瘫步态异常病因方面迈出了重要一步。它首次在群体层面量化了多种肌肉骨骼损伤的独立与交互贡献，为建立更全面的病因学模型奠定了基础。研究所展示的预测性仿真结合个性化建模的方法，是一个极具前景的工具平台。未来，通过整合更精准的损伤模型（尤其是运动控制障碍和足踝畸形），以及利用更定量化的仪器评估数据（如测力计替代徒手肌力测试），该平台有望更准确地模拟个体患者的步态，从而为临床医生评估不同手术方案潜在的生物力学后果、实现真正的个体化治疗规划提供强大的决策支持。最终，帮助避免无效手术，将治疗资源精准地导向最可能受益的患儿。