该文档由Tiam Heydari, Maziar Heidari, Omid Mashinchian, Michal Wojcik, Ke Xu, Matthew John Dalby, Morteza Mahmoudi, Mohammad Reza Ejtehadi等作者（分别来自沙里夫理工大学、马克斯·普朗克聚合物研究所、瑞士洛桑联邦理工学院、加州大学伯克利分校、劳伦斯伯克利国家实验室、格拉斯哥大学、德黑兰医科大学、哈佛医学院等机构）完成，以题为《development of a virtual cell model to predict cell response to substrate topography》的研究论文形式发表于《ACS Nano》期刊（2017年，第11卷，第9084-9092页）。这是一项关于细胞响应基质拓扑结构预测模型开发的原创性研究。以下是对该研究的学术报告。
本研究的学术背景属于细胞生物学、生物物理、纳米生物技术和计算生物学交叉领域，特别聚焦于细胞与材料界面相互作用，即细胞如何感知并响应细胞外基质（Extracellular Matrix, ECM）的物理化学特性（如拓扑结构、刚度、化学信息）。已知，工程化的基质（例如具有微/纳米图案或特定机械性能的材料）能够有效引导细胞行为，尤其是间充质干细胞（Mesenchymal Stem Cells, MSCs）的分化方向。例如，纳米级特征可引导MSCs向成骨细胞谱系分化，而材料表面机械性能的变化则能指导干细胞谱系选择。然而，当前领域面临的一个主要挑战是缺乏能够在实验评估之前，就可靠且准确地预测细胞行为的建模方法。这通常意味着需要进行大量耗时的细胞测试迭代才能确定最佳的材料特征。因此，本研究旨在开发一个统一的计算框架，构建一个名为“虚拟细胞模型（virtual cell model）”的多组分细胞模型，以期能够预测细胞在不同基质上培养时，整个细胞和细胞核特征（包括形状、方向乃至染色质构象）的变化，从而为优化面向广泛细胞应用（如干细胞分化）的基质开发提供一个可靠、高效、快速的高通量预测方法。
该研究的工作流程主要包括两大方面：计算模型开发与验证 以及 实验验证。首先，在计算模型开发方面，研究者构建了一个多组分的“虚拟细胞模型”。该模型包含以下关键组成部分：(1) 细胞膜和核膜：使用三角化流体膜模型进行模拟，包含弹性键合势、排斥势、弯曲能及表面积约束能，以反映脂质双层的机械特性（如流动性和弹性）。(2) 细胞质和细胞骨架网络：由一组分布的质量粒子（mass particles）通过粘弹性（Voigt和Maxwell元件）连接形成的三维网络填充在细胞膜与核膜之间，模拟细胞质的粘弹性，其连接弹性参数κ_cyt为k0。(3) 染色质纤维：使用珠-簧链（bead-spring chain）模型模拟，包含键合能、弯曲能和排除体积相互作用，并被限制在核膜内。(4) 细胞外基质模型：采用弹性网络模型模拟，其顶部表面被三角化，基质刚度通过连接键的弹性常数κ_ecm调节。(5) 细胞-基质相互作用：模拟了两种关键相互作用：一是细胞膜对基质的粘附，通过膜单元与基质表面单元之间的类Lennard-Jones势来描述；二是由细胞边缘肌动蛋白聚合产生的有效外向推力（F_active），该力推动膜单元向基质表面靠近，并可设置各向异性以模拟细胞的极性铺展。模型所有参数采用模拟单位，如长度单位a、能量单位k_BT等，并通过C++平台集成实现。模拟过程包括初始化（细胞膜和核膜设为球形）、系统驰豫（约100,000模拟步以达到平衡）、将虚拟细胞放置于不同特征的虚拟基质上、再次驰豫并进行数据采集和分析。
为了验证模型的有效性，研究者设计了与计算模拟相对应的三组实验，对间充质干细胞在具有不同物理特性的基质上的行为进行了验证。第一组实验针对具有沟槽图案的基质：在聚甲基丙烯酸甲酯（PMMA）基板上制造了不同宽度（5 μm和50 μm）和深度（100-400 nm）的微沟槽图案，将MSCs接种于这些基板上培养4天，固定后进行染色（肌动蛋白和DAPI核染色），通过分析显微镜图像计算细胞和细胞核的长宽比（aspect ratio）。第二组实验针对具有不同刚度的基质：使用了不同刚度（1 kPa， 8 kPa， 20 kPa）的聚丙烯酰胺（PA）凝胶作为基底，培养MSCs后，使用随机光学重建显微镜（STORM）对细胞核进行超分辨率成像，并分析其面积、长宽比等形态参数。第三组实验则利用细胞印迹基质：使用原子力显微镜（AFM）获取已发表研究中基于间充质干细胞、软骨细胞和角质形成细胞的形态所制备的细胞印迹基底的拓扑结构，并将其用于指导虚拟细胞模拟。这三组实验分别对应虚拟模型对拓扑引导、刚度响应和复杂形态复制能力的验证。
本研究的主要结果具体体现在以下几个方面。首先，对于沟槽基质的模拟与实验结果显示，虚拟细胞在沟槽基底上能沿沟槽方向伸长，且伸长程度与沟槽深度呈正相关（图1b）。这一定性趋势与实验观察完全一致（图1a）。对实验数据的量化分析（图1c）和虚拟模型的计算分析（图1d）均表明，细胞核的长宽比随着沟槽深度增加而增加。重要的是，虚拟模型的预测显示，细胞的伸长不仅仅是几何覆盖效应（类似于一张墙纸覆盖沟槽后呈现的椭圆形状），而是包含了主动的生物力学响应，其伸长程度超过了单纯的几何覆盖效应。此外，虚拟模型还预测了细胞在沟槽上的取向分布（图1e），显示细胞沿沟槽方向的取向程度也随沟槽深度增加而增强。该模型进一步揭示了基底吸附对细胞核内染色质构象的影响：细胞吸附到基底后，细胞核体积减小，染色质密度增加，导致染色质纤维不同片段间的接触数量增加，并且接触的稳定性增强（即随时间变化的接触百分比降低）（图1f-h）。这些结果暗示，基底拓扑结构可能通过机械传导影响染色质空间组织，从而可能影响基因表达。
其次，在基质刚度的模拟与实验验证方面，STORM实验图像（图2c）显示，基底刚度增加会导致细胞核发生局域形状变形（伸长），其铺展面积增加（图2e）且长宽比增大（图2d）。相应地，虚拟细胞模型在模拟各向异性细胞在不同刚度（κ_ecm从1k0到20k0）的三维弹性网络基底上的响应时，成功预测了类似的趋势：随着基底模型刚度的增加，虚拟细胞核的伸长（长宽比增加，图2g）和铺展面积增加（图2h），同时细胞核体积减小（图2i）。这种体积减小使得核内染色质密度增加，染色质片段间的接触概率空间密度相应增大（图2k）。有趣的是，模型分析表明，在更硬的基底上（对应更高染色质密度），染色质接触随时间变化的百分比反而降低（图2j），意味着这些接触更加稳定。这表明，模型不仅能预测细胞和细胞核的形态变化，还能提供关于核内染色质组织动力学变化的见解，为理解细胞如何感知ECM刚度（机械传导）提供了新的计算视角。
第三，在细胞印迹基质的模拟中，虚拟细胞模型展现了强大的形态预测能力。研究者将实际AFM获取的间充质干细胞和角质形成细胞印迹基底的三维形貌数据（图3a， b上半部分）直接离散化并导入虚拟模型。模拟结果显示，放置在“间充质干细胞”形态印迹上的虚拟细胞，演变成了两极化的、成纤维细胞样的伸长形态和细胞核（图3a下半部分）；而放置在“角质形成细胞”形态印迹上的虚拟细胞，则趋向于形成更加球形的细胞和细胞核形态（图3b下半部分）。这些模拟结果与先前发表的实验观察完全吻合，即干细胞在特定成熟细胞形态的印迹上培养时，会倾向于分化并呈现出该细胞类型的形态特征。这证明了虚拟细胞模型可以复现并解释实验观察到的“形态引导分化”现象，表明细胞能够模拟印迹基底的形状，并在细胞和细胞核层面改变其几何结构，从而可能影响其功能走向。
基于以上结果，本研究得出结论：成功开发了一种能够在计算机（in silico）中探测细胞功能/行为的虚拟细胞模型。通过与实验结果的强相关性验证，该模型展现出作为一种高效、快速的预测工具的潜力，可以为构建用于引导特定细胞功能（如干细胞分化）的最佳化基质提供所需信息。本研究的科学价值在于，首次提供了一个整合了细胞膜、细胞质/骨架、细胞核及染色质、以及可调ECM的统一计算框架，能够定性地预测细胞对基质拓扑、刚度和复杂形态的形态学响应，并深入到染色质构象这一亚细胞层次。其应用价值显著：通过预测性建模，可以减少为寻找最佳材料特性而进行的昂贵且耗时的实验迭代，加速组织工程、再生医学和药物筛选等领域中优化细胞培养基底的设计与开发。该模型为理解物理线索（地形、刚度）如何通过力学信号传导（mechanotransduction）影响细胞命运（特别是干细胞分化）提供了新的机制性见解，例如提示基底物理特性可能通过改变核形状和染色质接触来影响基因表达程序。
本研究的亮点在于：第一，方法论创新：创造性地开发了一个多组分、多尺度的“虚拟细胞”计算模型，将细胞粘附、铺展、细胞核形变乃至内部染色质纤维的构象变化整合在一个统一的力学框架内。第二，预测能力验证全面：模型在三种不同类型的基底（微沟槽、变刚度凝胶、复杂细胞印迹）上均得到实验数据的验证，证明了其广泛的适用性和定性预测的可靠性。第三，深入亚细胞机制：模型不仅预测了细胞和细胞核的宏观形态，还首次在计算模型中探索了基底物理特性对细胞核内染色质空间组织和接触稳定性的潜在影响，将细胞对基质的物理感知延伸到了表观遗传调控的可能层面。第四，潜在的高通量应用前景：研究明确提出了该模型可作为高通量筛选优化细胞培养基底的工具，具有重要的转化医学意义。此外，研究还指出模型的未来发展方向，如通过引入真实细胞参数（杨氏模量、肌动蛋白聚合力等）来实现定量预测，以及考虑培养时间对分化过程的动态影响，从而进一步提升模型的预测能力和生物学真实性。