[论文解读] Mathematical models in biology

Barbara Bori C. Mazzag|eScholarship (California Digital Library)|Jun 16, 2003

Slime Mold and Myxomycetes Research被引用 29

一句话总结

本论文建立了并分析了三种关键生物过程的数学模型：*Azospirillum brasilense* 的趋氧性、cAMP 和钙信号介导的生长锥导向，以及流体剪切力作用下内皮细胞的形变。通过反应-扩散方程、随机模拟和黏弹性网络模型，研究揭示了分子适应性与机械特性如何实现稳健的梯度感应与定向运动；关键结果表明，在趋氧性中形成条带结构，而在振荡流作用下内皮细胞表现出频率依赖的形变特征。

ABSTRACT

Aerotaxis is the particular form of chemotaxis in which oxygen plays the role of both the attractant and the repellent. Aerotaxis occurs without methylation adaptation, and it leads to fast and complete aggregation toward the most favorable oxygen concentration. Biochemical pathways of aerotaxis remain largely elusive, however, aerotactic pattern formation is well documented. This allows mathematical modeling to test plausible hypotheses about the biochemical mechanisms. Our model demonstrates that assuming fast, non-methylation adaptation produces theoretical results that are consistent with experimental observations. We obtain analytical estimates for parameter values that are difficult to obtain experimentally. Chemotaxis in growth cones differs from gradient sensing in other animal cells, because growth cones can change their attractive or repulsive response to the same chemical gradient based on their internal calcium or cAMP levels. We create two models describing different aspects of growth cone guidance. One model describes the internal switch that determines the direction of movement. However, this model allows chemotaxis under certain conditions only, so a second model is created to propose a mechanism that allows growth cone guidance in any environment. Endothelial cells go through extensive morphological changes when exposed to shear stress due to blood flow. These morphological changes are thought to be at least partially the result of mechanical signals, such as deformations, transmitted to the cell structures. Our model describes an endothelial cell as a network of viscoelastic Kelvin bodies with experimentally obtained parameters. Qualitative predictions of the model agree with experiments.

研究动机与目标

• 通过基于氧气依赖性转向频率的细菌转向行为随机模型，理解 *Azospirillum brasilense* 趋氧性的机制。
• 研究生长锥如何通过钙信号与cAMP信号动力学检测并响应导向信号分子（如 netrin-1）。
• 利用 Kelvin 体的串联与并联排列构建黏弹性网络模型，模拟内皮细胞细胞骨架的力学响应。
• 分析参数变化对稳态与振荡流条件下内皮细胞形变与受力分布的影响。
• 将分子尺度的信号动力学与宏观细胞行为（如条带形成与定向迁移）相联系。

提出的方法

• 基于氧气依赖性转向频率开发了细菌趋氧性的随机模型，并对系统进行无量纲化以利于分析。
• 构建了生长锥内 cAMP 与钙信号动力学的反应-扩散模型，引入 Goldbeter-Koshland 开关以实现适应与放大。
• 通过 Kelvin 体的串联与并联组合构建黏弹性网络模型，以表征内皮细胞骨架的力学行为。
• 采用四阶龙格-库塔法对稳态与振荡流条件下的常微分方程组进行数值模拟。
• 开展参数敏感性分析，评估黏弹性参数（μ, η）对形变与受力分配的影响。
• 使用 MATLAB 代码模拟并可视化内皮细胞模型中随时间演化的形变与受力分布。

实验结果

研究问题

• RQ1 氧气梯度如何导致 *Azospirillum brasilense* 群体在趋氧行为下形成稳定的条带结构？
• RQ2 钙浓度在调节生长锥对 netrin-1 梯度响应时的转向方向中起什么作用？
• RQ3 内皮细胞的黏弹性特性如何影响其在稳态与振荡流体剪切力下的形变行为？
• RQ4 在振荡流作用下，内皮细胞的峰值形变与稳态形变对频率的依赖关系如何？
• RQ5 当材料参数变化时，黏弹性网络内部的机械力分布如何改变？

主要发现

• 数值模拟证实，趋氧性条带结构源于氧气梯度与氧气依赖性转向频率，条带在空间与时间尺度上均保持稳定。
• 模型预测，高细胞质钙浓度可导致生长锥暴露于更高配体浓度的一侧产生更高活性物质（A）浓度，而低钙浓度则反转该梯度，从而解释转向行为。
• 内皮细胞形变对黏度参数 η₁₂ 极为敏感，在振荡流作用下，随着 η₁₂ 增大，峰值形变亦随之增加。
• 在双体并联 Kelvin 模型中，稳态形变与弹簧常数 μ₀₂ 呈非线性依赖关系，最大形变出现在中间值附近。
• 振荡流下的峰值形变表现出类似共振的行为，特定频率下响应最大，表明存在频率选择性的机械敏感性。
• 网络模型中的受力分配系数显著依赖于 μ₀₂ 与 η₁₂，表明机械载荷分布高度依赖于黏弹性参数。