这篇研究论文《in silico models for dynamic connected cell cultures mimicking hepatocyte-endothelial cell-adipocyte interaction circle》由意大利比萨大学的研究团队完成，主要作者包括Chiara Andreoni（通讯作者）、Gianni Orsi、Carmelo De Maria、Francesca Montemurro和Giovanni Vozzi。该项研究发表于2014年12月15日的学术期刊《PLOS ONE》上，其DOI号为10.1371/journal.pone.0111946，是一份开放获取（open-access）的学术成果。

这项研究属于生物工程与计算系统生物学交叉领域。研究的核心动机源于对复杂生物系统，特别是涉及多组织器官协同维持体内稳态机制的理解需求。在活体内，维持能量代谢稳态（homeostatic balance）需要肝脏、脂肪组织、内皮等多种细胞类型通过血液流动进行复杂的分子通讯与代谢交互。然而，在实验室中，研究这种涉及三种或更多细胞类型的生化系统，因其相互作用的极度复杂性而变得几乎不可能。传统的“体外”（in vitro）实验虽能提供数据，但往往成本高、耗时长、数据分散，且存在伦理考量。与此同时，许多已有的“硅上”（in silico）计算模型多专注于单一细胞的代谢通路或静态共培养，缺乏模拟动态连接培养系统中多细胞间代谢物互传与稳态平衡的能力。因此，本研究旨在开发一种新颖的计算模型，通过整合已有的肝细胞、内皮细胞和脂肪细胞的单细胞代谢模型，并引入它们之间生理性的互连，来模拟一个动态连接的三培养系统的代谢行为。其最终目标是重现该体外系统在培养基流动条件下，在葡萄糖和脂质代谢方面所观察到的稳态平衡成就。

研究详细工作流程

本研究的工作流程并非线性实验，而是一个迭代的“构建-验证-修正”的计算模型开发过程，可分为以下几个核心阶段：

第一阶段：研究基础与单细胞模型准备 本研究建立在研究团队前期工作的坚实基础上。在此之前，团队已开发并发表了三个独立的虚拟单细胞代谢模型，均使用Simulink（The MathWorks, Inc.）软件平台实现。这些模型分别是：模拟肝细胞代谢的HEMETB模型、模拟内皮细胞代谢的ENMET模型以及模拟脂肪细胞代谢的ADMET模型。这些模型均采用模块化结构，其数学核心基于米氏动力学（Michaelis-Menten kinetics）和非线性微分方程，用以描述酶促反应。它们能够模拟细胞在代谢物浓度变化下的动态行为，但仅限于单培养或静态共培养场景。此外，团队还曾开发过一个名为CREPE的模型，用于模拟静态条件下肝细胞与内皮细胞共培养的代谢串扰。这些前序模型为本研究的集成工作提供了可直接利用的、结构一致的模块化单元。

第二阶段：集成模型构建——从单细胞到三细胞动态系统 本阶段是研究的核心创新环节。研究者并非从零开始，而是将HEMETB、ENMET和ADMET这三个模型作为基础模块进行“连接”。 1. 细胞增殖模型定义：首先，为每个细胞类型定义了更贴近后续验证实验的细胞增殖模型。通过拟合文献中的实验数据，肝细胞增殖采用经典的逻辑斯蒂增长方程（Verhulst方程）进行模拟，而内皮细胞和脂肪细胞（在实验时间尺度内视为成熟且无显著增殖）的细胞数量则被设定为随时间恒定，等于初始接种数。这一步骤确保了模型能反映不同细胞群在培养过程中的数量动态或静态特征。 2. 关键代谢通路扩展与互连机制引入：为了实现多细胞间的代谢交互，研究者在原有单细胞模型的基础上，为每种细胞引入了新的代谢通路和运输机制，这些是模拟稳态平衡所必需的。 * 肝细胞模型（HEMETB扩展）：新增了糖异生（gluconeogenesis）、糖原分解（glycogenolysis）和甘油代谢通路。引入了通过AQP9（水甘油通道蛋白9）介导的甘油摄取（基于菲克扩散定律建模）和通过CD36（脂肪酸转位酶）介导的脂肪酸摄取（基于米氏动力学建模）。同时考虑了葡萄糖和甘油三酯（triglyceride， TG）的释放过程。 * 内皮细胞模型（ENMET扩展）：新增了糖异生和脂肪酸合成通路。最关键的是引入了脂蛋白脂肪酶（lipoprotein lipase， LPL）的活性，该酶位于内皮细胞表面，能够水解循环中的甘油三酯（在本模型中，即来自肝细胞释放的TG），产生游离脂肪酸（free fatty acid， FFA）和甘油，供周围细胞利用。脂肪酸摄取同样通过CD36介导。 * 脂肪细胞模型（ADMET扩展）：引入了甘油三酯/脂肪酸循环中的关键过程，包括甘油三酯水解（脂解）、脂肪酸的摄取（通过CD36）与释放，以及一个特殊的代谢通路——甘油异生（glyceroneogenesis），这是脂肪细胞内从非碳水化合物前体生成甘油-3-磷酸的途径，对于脂肪酸的再酯化至关重要。甘油释放通过AQP7介导。 3. 系统集成与培养基浓度建模：将上述三个扩展后的细胞代谢模块进行物理和生化连接。在Simulink中，它们通过代表“培养基”的共享积分器模块相连。研究人员为葡萄糖、游离脂肪酸和甘油这三种关键代谢物在培养基中的浓度变化建立了核心的状态方程组。这些方程（如文中公式12, 13, 14）本质上是描述每个代谢物在培养基中净变化的物质平衡方程：它们综合了所有细胞类型对该代谢物的摄取（负项）和释放（正项）速率。例如，培养基葡萄糖浓度的变化率等于肝细胞、内皮细胞和脂肪细胞的葡萄糖释放速率之和减去它们的葡萄糖摄取速率之和。这样，模型就能动态模拟三种细胞通过培养基流动（在模型中体现为代谢物的共享和交换）进行的代谢串扰。

第三阶段：模型验证与参数评估 构建好的集成模型需要通过实验数据进行验证。研究者选择了三篇逻辑连贯的文献（Vinci等人2010, 2012; Iori等人2012）作为验证基准，这些文献报道了使用模块化多室生物反应器（MCMB）进行的相应体外实验数据。 1. 单培养验证（基线验证）：首先，将单个细胞模型（肝细胞、内皮细胞、脂肪细胞模型）的模拟结果与文献中对应细胞在静态和动态（有培养基流动）条件下的单培养实验数据进行比较。实验数据包括培养48小时内，培养基中葡萄糖、游离脂肪酸和甘油浓度的变化。为了模拟“流动”对细胞代谢的影响（实验观察到动态条件下代谢活性通常更高），研究者并非修改模型结构，而是引入了一套可启发式调整的“校正因子”。通过在模型界面中选择“静态”或“动态”模式，系统会自动对一系列关键酶的动力学参数（如Vmax）应用不同的比例系数，以此来重现不同培养条件下观察到的代谢表型差异。通过调整这些校正因子，模型成功再现了大部分单培养实验趋势，例如肝细胞在动态条件下更强的葡萄糖和脂肪酸摄取、脂肪细胞持续的脂肪酸和甘油释放等。 2. 三细胞连接系统验证（核心目标验证）：随后，将集成的三细胞模型（使用从单培养动态条件优化得到的参数）的模拟结果，与文献中动态连接的三培养（肝-内皮-脂肪）系统的实验数据进行比较。初始模拟结果并不理想：模型预测出培养基中葡萄糖被持续消耗、脂肪酸被完全清除、甘油浓度下降，这与实验观察到的稳态平衡（三种代谢物浓度在48小时内均无明显净变化）大相径庭。这表明，简单地将单培养模型参数相加，无法重现多细胞系统涌现出的稳态特性。 3. 模型修正与关键参数识别：基于生物化学知识和对模型行为的分析，研究者推断，在连接系统中，占主导地位的肝细胞可能调整了其代谢策略以支持其他细胞。因此，他们有选择地降低了肝细胞的几个关键摄取速率：将肝细胞葡萄糖摄取率降低了100倍，CD36介导的脂肪酸摄取Vmax降低了5倍，AQP9介导的甘油摄取率降低了1000倍。同时也略微降低了脂肪细胞的脂肪酸摄取能力。经过这些针对性的参数调整后，修改后的模型模拟结果与实验数据取得了良好的一致，能够重现葡萄糖和游离脂肪酸浓度在培养基中基本保持稳定的稳态现象。

第四阶段：数据分析方法 研究中的数据“分析”主要体现在模型仿真与实验数据的对比上。实验数据从文献图表中提取，表示为均值±标准差。模拟采用Simulink中的ODE 23s（刚性/修正罗森布洛克）数值积分方法进行。对比重点关注0小时和48小时两个时间点代谢物浓度的净变化，与实验文献的分析方法保持一致。模型的内在数据（如细胞内代谢物浓度）由于缺乏实验测量值，其初始条件通常设为零，其动态变化完全由模型方程驱动。

研究主要结果

1. 成功开发了一个模块化的三细胞动态连接代谢模型：该模型整合了肝细胞、内皮细胞和脂肪细胞的代谢网络，能够模拟它们通过培养基流动进行的代谢物交换。
2. 单培养模型验证成功：通过引入基于培养条件（静/动态）的启发式校正因子，单个细胞模型的模拟结果能够较好地重现对应体外单培养实验中葡萄糖和脂代谢的主要趋势。这证明了基础单细胞模型的可靠性以及通过参数调整模拟环境刺激影响的可行性。
3. 揭示了简单加和模型的局限性：首次尝试运行三细胞连接模型时，直接使用单培养（动态）参数得到的结果与三培养实验数据严重不符。这一关键结果表明，动态连接的三培养系统并非其组成细胞单培养行为的简单加和，而是一个具有涌现属性的复杂系统。肝细胞的代谢主导性在简单加和模型中过度凸显，破坏了稳态。
4. 通过针对性参数调整实现了稳态平衡的模拟：通过有目的地降低肝细胞对葡萄糖、脂肪酸和甘油的竞争性摄取能力，修改后的模型成功模拟出了实验中观察到的代谢稳态。这一修正过程本身具有科学发现价值，它提示在真实的连接系统中，可能存在类似的代谢适应机制，例如肝细胞在正常血糖条件下降低自身葡萄糖摄取，以优先满足内皮和脂肪等葡萄糖依赖型细胞的需求。
5. 模型性能评估：最终修正后的模型对葡萄糖和游离脂肪酸浓度动态的模拟与实验数据吻合良好。对甘油浓度的模拟虽有改善，但匹配度相对较低，研究者分析这可能与模型中脂肪细胞和内皮细胞的甘油释放量不足，以及细胞内代谢物初始条件设置为零的局限性有关。

研究结论与意义

本研究成功构建并初步验证了一个能够模拟肝细胞-内皮细胞-脂肪细胞动态连接培养系统代谢行为的计算模型。研究表明，通过集成已有的单细胞模型并精心设计细胞间的代谢互连，可以重现体外实验中观察到的多细胞系统层面的特性——稳态平衡。更重要的是，研究过程揭示了从单细胞行为预测多细胞系统行为的复杂性，强调了系统涌现属性的存在。为实现稳态，模型需要对关键参数进行特定调整，这反过来为理解体内可能的代谢调节机制（如器官间的代谢优先级调整）提供了计算视角和假设。

该研究的科学价值在于：1）方法学创新：提供了一种从模块化单细胞模型构建多细胞、多组织相互作用模型的可扩展框架。2）研究工具：所开发的模型是一个可用的“硅上”研究平台，能帮助其他研究者以低成本、无伦理顾虑的方式，进一步探究整合代谢、代谢串扰及其调控机制。3）加深理解：通过模型构建和验证过程，加深了对肝脏、脂肪和内皮组织在能量代谢稳态中协同作用机制的理解。

应用价值主要体现在生物医学研究和药物开发领域：该模型可用于模拟代谢性疾病（如肥胖、糖尿病、代谢综合征）状态下多组织互作的紊乱，辅助病理机制研究和药物疗效的初步预测。

研究亮点

1. 研究目标新颖：首次尝试构建一个模拟三种不同细胞类型（肝、内皮、脂肪）在动态流动培养条件下代谢交互的计算模型，直接针对能量代谢稳态这一核心生理/病理过程。
2. 方法具有继承性与创新性：基于团队成熟的前期单细胞模型进行集成，保证了模型基础的可靠性；通过引入细胞间代谢物交换的数学模型和针对培养条件的参数校正策略，实现了从静态单细胞到动态多细胞系统的功能跨越。
3. 揭示了“涌现属性”：研究过程清晰展示了复杂系统研究的典型挑战——整体不等于部分之和。模型初期失败与后期成功的对比，强有力地证明了多细胞系统稳态是一种涌现属性，必须通过系统层面的互连和调节才能实现。
4. 良好的验证与迭代过程：研究严格遵循计算模型研究的范式，使用独立的体外实验数据进行多层级验证（单培养基线验证、多培养核心验证），并根据 discrepency 进行合理的模型修正和假设生成，体现了研究的严谨性。

其他有价值内容

研究在讨论部分展望了模型的未来发展方向，包括：扩展模型以纳入其他关键代谢器官（如胰腺）、增加新的代谢方面（如乳酸代谢、激素调节）、引入基因表达调控模型以使系统能更动态地适应刺激，以及尝试结合优化算法进行参数估计等。这些展望为后续研究指明了潜在路径。文中还详细列出了所使用的所有酶动力学参数表、状态方程和化学计量方程（作为补充材料），增强了研究的可重复性和透明度。