在中枢神经系统疾病治疗领域，血脑屏障如同一道坚固的防线，阻碍了大多数药物通过血液循环进入大脑。为了突破这一屏障，医学界将目光投向了直接通过脑脊液给药的创新途径。这种给药方式能够绕过血脑屏障的限制，将药物直接送达病灶区域，为阿尔茨海默病、帕金森病等神经系统疾病的治疗带来了新的希望。然而，这条看似捷径的道路却充满了未知与挑战——药物在复杂脑脊液系统中的分布规律难以预测，个体解剖结构差异、脑脊液流动动力学、颅脊顺应性等多重因素相互交织，使得给药方案的优化如同在迷宫中寻找出口。传统的动物实验虽然能够提供部分答案，但成本高昂且难以实现多参数系统优化。特别是在非人灵长类动物模型中开展研究，更是面临伦理和资源的双重压力。正是基于这样的背景，Mohammadreza Khani等研究人员在《Fluids and Barriers of the CNS》上发表了他们的最新研究成果，致力于通过数字医学技术破解这一难题。研究团队创新性地构建了首个包含全脑脊液系统的非人灵长类动物数字模型，该模型独特之处在于整合了真实的生理特征：不仅包含完整的脑室系统、基底池、皮质蛛网膜下腔等关键解剖结构，还精细模拟了31对脊髓神经根丝的复杂几何形态，以及心脏搏动和呼吸运动共同驱动的脑脊液振荡现象。更重要的是，研究人员通过动态网格技术成功实现了颅脊顺应性的模拟，使模型能够再现脑脊液搏动沿神经轴的非均匀衰减这一关键生理现象。在技术方法上，研究团队基于4岁雄性食蟹猴的活体MRI数据，通过T2加权成像获取了0.38毫米各向同性的高分辨率解剖图像，并利用相位对比MRI在六个椎体水平量化了脑脊液流动波形。采用OpenFOAM软件平台中的多相混合物模型进行计算流体动力学模拟，同时通过3D打印制作了几何结构完全相同的实验验证模型，以荧光素钠作为示踪剂进行比对验证。
模型几何结构特征研究构建的模型总脑脊液体积为16.2毫升，其中脊髓部分占7.8毫升，颅内容积为8.3毫升。模型精细再现了非人灵长类动物特有的第四脑室五个开口的解剖特征，与人类三个开口的结构形成明显差异。脊髓神经根丝直径范围在0.25-0.55毫米之间，平均每个椎体水平每侧有2.25个根丝，总表面积达22.7平方厘米，为药物与组织相互作用提供了重要界面。刚性数字模型的实验验证通过对比刚性数字模型与3D打印实验模型的结果，研究发现两者在时空示踪剂分布上表现出高度一致性（R2=0.88）。注射后60分钟，两个模型中示踪剂的空间分布模式和区域药物百分比均高度吻合，验证了数字模型预测的可靠性。这种强相关性为后续顺应性模型的扩展应用奠定了坚实基础。顺应性对脑脊液动力学的影响引入顺应性特征后，模型表现出显著的生理学变化：硬膜径向位移在颈椎区域最大达到120微米，腰椎区域约为20微米。与刚性模型相比，顺应性模型中的脑脊液峰值流速在腰椎区域降低了约10倍（0.065 vs 0.65毫升/分钟），导致整个神经轴的平均雷诺数从250降至100以下。这种流体动力学特性的改变直接影响了药物的分散行为。
顺应性对药物分散的影响研究发现顺应性是决定药物分布的关键因素。注射后1小时，顺应性模型中91.9%的给药剂量保留在腰椎区域，而刚性模型仅为72.9%。在刚性模型中，药物能够在1小时内扩散至颅底，而顺应性模型中的药物甚至未能超过T11椎体水平。药物在注射部位的曲线下面积值也呈现显著差异（顺应性模型2.8%-小时 vs 刚性模型1.2%-小时），表明顺应性显著延缓了药物的神经轴分散。
