细胞粘附的连续性数学模型及数值模拟.doc

细胞粘附的连续性数学模型 和数值模拟 吕杰 曹金凤 邬黎 摘要 细胞群通过细胞表面的粘附分子不断聚集的过程被称为细胞粘附，是生物学许多研究领域的基础（包括早期的胚胎发育、组织的新陈代谢以及肿瘤生长都涉及到细胞粘附）。本文基于Nicola于2006年提出的连续性细胞粘附模型，通过设计数值计算方法来模拟细胞粘附的整个过程，并在此基础上通过改变细胞粘附力的作用范围、作用效果、以及初始细胞密度等条件来比较模型中各参数对数值模拟结果的影响，最后引入Logistic模型来描述细胞本身的自然增长，构建出一个新的连续性细胞粘附模型，使模型更符合实际。数值模拟结果表明，这一连续性数学模型对不同种类细胞之间粘附的最终形态能做出较准确的预判，能够直接为基础生物学研究和医学临床应用提供参考和帮助。 关键词：细胞粘附、数学模型、数值模拟 绪言 细胞群通过细胞表面的粘附分子（CAMs）不断聚集的过程被称为细胞粘附，是生物学许多研究领域的基础。它与组织的形成（早期的胚胎发育）、稳定（组织间的新陈代谢）、病变（肿瘤细胞的入侵和转移）密切相关。 因此，早在几十年前，国外的许多学者就开始进入这一领域进行研究，通过各种不同的实验来观测、解释和预判细胞粘附这一现象（图1）。其中一个比较著名的实验是由Steinberg在1962年完成的，该实验表明当两种不同种类的细胞混合时，会形成不同的聚集形态，Steinberg还基于热力学模型提出了微分粘附假说，能够较好地预测实验结果，同时揭示了不同种类细胞混合后的最终形态与细胞的自粘附系数和互粘附系数直接相关。 图1 细胞粘附实验 随机地在培养皿中放入细胞，（a）、（b）、（c）分别是5小时后、19小时后以及48小时后细胞粘附的形态 （摘自Armstrong, P.B., 1971） 此后，许多学者在Steinberg实验的基础上，试图通过数学建模，更深入地了解细胞粘附的整个过程。到目前为止，已经形成了几个比较成功的数学模型，但大都采用离散的数学模型，然而，细胞的运动（包括粘附）无论在时间上还是空间上都是一个连续的过程，所以，在2006年，Nicola首先提出了连续性细胞粘附的数学模型。 本文的主要工作是基于Nicola模型，通过设计数值计算方法来模拟细胞粘附的整个过程，并在此基础上通过改变细胞粘附力的作用范围、作用效果、以及初始细胞密度等条件来比较模型中各参数对数值模拟结果的影响，最后引入Logistic模型来描述细胞本身的自然增长，构建出一个新的连续性细胞粘附模型，使模型更符合实际，能更好地模拟细胞粘附所引起的细胞聚集现象，预测最终的细胞群形态。目前，国内对细胞粘附这一问题的研究基本还停留在实验阶段，对此进行数学建模的文献不多见，本文的相关结论亦可为国内基础生物学研究和医学临床应用提供参考和帮助。 本文的内容结构划分如下：第II部分主要介绍Nicola模型；第III部分，通过设计合适的数值计算方法，对Nicola模型进行数值模拟，分别得到一维和二维单组分情况下细胞粘附的最终形态，并通过分析自粘附系数和互粘附系数之间的关系，分别模拟一维和二维情况下不同种类细胞之间粘附的最终形态，并通过与Steinberg和Foty等人的实验结果和基于实验的预测作比较，验证数值计算方法的合理性；第IV部分，通过改变细胞粘附力的作用范围、作用效果、以及初始细胞密度等条件来比较模型中各参数对数值模拟结果的影响，并引入了Logistic模型来描述细胞本身的自然增长，建立新的细胞粘附的连续性数学模型，并在此基础上作相应阐述；第V部分则是对本文的总结与展望。 Nicola模型的简单介绍 Nicola（2006）模型将细胞运动的主要原因归结为细胞的自由扩散和细胞粘附力的作用。以下首先介绍单组分的一维模型。 假设不考虑细胞的增长和死亡，由质量守恒，得到 （1） 其中，表示在横坐标为处，时间为时的细胞密度 表示扩散通量，表示粘附通量 假定系统处于等温等压环境下，则满足定理，有 其中是扩散系数 Nicola假设粘附通量与细胞的密度及细胞间粘附力的大小成正比，与细胞的尺寸成反比，因此可设 其中是一个与粘性线性有关的常数 是施加于细胞的所有外力总和 是细胞的影响半径 在上的细胞受到的外力总和等于所有到的距离为的细胞对它的粘附力的总和。这里的就是影响半径。粘附力的方向取决于处的细胞相对于的位置。基于此，Nicola假设 其中取决于细胞密度，是一个与粘附力大小有关的正的参数，用来确定的方向以及在不同半径处的粘附强度 考虑到细胞聚集的密度总是有限的，Nicola假设 并假设在细胞的影响半径内粘附强度恒为1，则的表达式如下： 因此， （1）式变成 其中 将上