跳转至

第十八章:胃肠道系统读书笔记

书籍信息:Keener & Sneyd, Mathematical Physiology (II: Systems Physiology), Springer 2009

章节信息:Chapter 18: The Gastrointestinal System


一、胃肠道结构概述

胃肠道 tract 从整体上具有共同的基本结构,尽管不同区域存在细节差异。其横截面结构(图18.1)包含以下层次:

外层肌肉层:由大量神经支配的肌肉层组成,排列成环状和纵向方向。这些肌肉层的收缩可以混合消化道内容物,并以可控的方式将食物沿适当方向推进。

黏膜下层:主要由结缔组织组成。

黏膜肌层:一层薄的平滑肌。

固有层:含有毛细血管和多种分泌腺的结缔组织层。

上皮细胞层:其性质因消化道不同区域而异。

神经支配:消化道主要有两层神经元——肠肌间神经丛(位于纵肌和环肌之间)和黏膜下神经丛(位于黏膜下层)。肠肌间神经丛的刺激可增加肠道节律性收缩的频率和强度,并增加兴奋波沿肠壁的传导速度;黏膜下神经丛主要是感觉性的,接收来自肠壁拉伸受体和肠上皮的信息。消化道也受到密集的神经支配,可以控制整个肠道或部分肠道的活动。


二、流体吸收:简单模型

2.1 吸收表面的解剖结构

胃肠道的主要功能是从通过其中的食物和液体混合物中吸收营养物质。肠道的吸收表面由许多褶皱(valvulae conniventes)组成,可将吸收黏膜面积增加约三倍。整个小肠黏膜表面分布着数百万个绒毛(villi),每个绒毛突出黏膜表面约1毫米,使吸收面积再增加十倍。绒毛的吸收表面由具有刷状缘(brush border)的上皮细胞组成,刷状缘包含多达1000个长1微米、直径0.1微米的微绒毛(microvilli),使暴露于肠内容的表面积再增加二十倍。综合所有表面突起,吸收表面积可达约250平方米——相当于一个网球场的面积。

上皮细胞负责吸收营养物质和水分。例如葡萄糖和氨基酸等化学营养物质的吸收与肾脏中相同,通过与 Na⁺ 的共转运体进行。

2.2 水通道蛋白的发现

长期以来,人们认为水是通过膜孔被动渗透而被吸收或转运的。渗透驱动的水流动本质上是细胞体积控制模型的基本假设。然而,1992年发现的水通道蛋白(aquaporins)彻底改变了这一观点。第一个被发现的水通道蛋白是 AQP1,由 Agre 及其同事于1992年发现,Agre 因此获得2003年诺贝尔化学奖。水通道蛋白已在几乎所有生物体中被鉴定,其对水的高度选择性甚至可以排斥 H₃O⁺ 和氢离子。

2.3 简单吸收模型

上皮细胞在其腔侧对水是不渗透的。然而,连接上皮细胞之间的紧密连接处有0.7–1.5纳米的小孔,允许水在腔和间质之间自由扩散。通过这些小孔的水吸收主要由 Na⁺ 在腔和间质之间的梯度驱动。Na⁺ 通过被动转运被输送到上皮细胞内部,然后通过 Na⁺–K⁺ ATP 酶从细胞内部移除到间质。Na⁺ 再由毛细血管血流从间质转运走。

建模时,我们将上皮消化道的一小段考虑为两个充分混合的隔室——腔和间质,由一层膜分开(图18.2)。假设腔内的 Na⁺ 浓度为 \(n_l\),细胞内部为 \(n_i\),间质中所有渗透物的浓度为 \(n\)。Na⁺ 从腔到细胞内部的流动 \(J\) 假定为被动的:

\[J = g(n_i - n_l) \quad (18.1)\]

其中 \(g\) 为某常数。Na⁺ 从细胞内部到间质的流动通过主动的 Na⁺–K⁺ ATP 酶:

\[J = f(n_i) \quad (18.2)\]

其中 \(f\) 为某饱和函数。通过紧密连接的水流 \(q\) 由腔和间质之间的渗透压差驱动:

\[Rq = n - n_l \quad (18.3)\]

其中 \(R\) 为紧密连接的水阻。最后,假设间质有毛细血管血流提供的流入和流出。流体流入为 \(Q\),进入浓度为 \(n_0\),渗透物流出为 \(Q+q\),浓度为 \(n\)。在稳态下,Na⁺ 守恒意味着:

\[g(n_l - n_i) = f(n_i) \quad (18.4)\]
\[(Q + q)n - Qn_0 = f(n_i) \quad (18.5)\]

2.4 无量纲化分析

定义无量纲浓度 \(u_j = n_j/N\),其中 \(y = q/Q\),可得:

\[\rho y^2 + (\rho + u_l)y + \kappa = 0 \quad (18.8)\]

其中 \(\kappa = u_i - u_0 + (1 - \gamma)\beta F(u_i) = 0\)\(\rho = RQ/N\)\(\gamma = g/Q\)\(\beta = gN/Q_f\),以及

\[u_l = u_i + \beta F(u_i) \quad (18.9)\]
\[F(u) = \frac{u^3}{1+u^3}\]

有四个无量纲参数:\(u_0\)(进入间质的渗透物的相对浓度)、\(\rho\)(紧密连接对水的阻力)、\(\gamma\)(腔侧细胞膜对 Na⁺ 的相对渗透性)、\(\beta\)(主动 Na⁺ 转运的最大速度)。

2.5 方程的解与生理意义

方程(18.8)是 \(y\) 的二次多项式,最多有一个正根。当 \(\beta(\gamma - 1)\) 是较大的正数时,\(\kappa\)\(u_i\) 的 N 形函数。相应地,正根也表现出 N 形行为。

重要的生理启示:这意味着可以通过调节腔内水的 Na⁺ 水平来最大化水的吸收。因此,含电解质的液体比纯水更能快速水化。然而,太多 Na⁺ 会产生相反的脱水效果。

当人脱水时,大量醛固酮由肾上腺分泌。醛固酮通过激活通道和泵蛋白的产生来增强上皮细胞对 Na⁺ 的转运,从而增加 Na⁺ 的被动和主动转运。醛固酮的作用可以用增加 \(g\)(腔侧 Na⁺ 转运的电导率)和/或增加 \(Q_f\)(主动 Na⁺ 泵的最大速率)来建模。

2.6 沿着肠道的变化

沿着肠道长度建立水含量和 Na⁺ 浓度的简单模型。假设食糜(从胃进入的食物、水和消化分泌物的混合物)以恒定速度作为塞状流移动。水通过渗透被移除,Na⁺ 被上皮细胞以由局部 Na⁺ 浓度决定的速率移除。在稳态下:

\[\frac{dQ_w}{dx} = -q(n_l) \quad (18.12)\]
\[\frac{d(n_l Q_w)}{dx} = -J(n_l) \quad (18.13)\]

2.7 便秘与腹泻

任何使粪便通过大肠的运动异常缓慢的情况都会导致便秘,使粪便有更多时间失水,从而变硬和干燥。任何阻碍正常运动的肠道病理(包括肿瘤、溃疡或强制抑制正常排便反射)都可引起便秘。

相反的情况——粪便快速通过大肠——称为腹泻。腹泻最常见的原因是感染性腹泻,其中病毒或细菌感染导致黏膜炎症。无论感染在哪里,受感染部位的黏膜分泌率都大大增加,净效应是大量液体被产生以帮助消除感染源。例如,霍乱和其他腹泻细菌的毒素刺激未成熟的上皮细胞释放大量 Na⁺ 和水,可能的目的是通过冲走细菌来对抗疾病。然而,如果这种过量的 Na⁺ 和水分泌不能被成熟、健康细胞的吸收所抵消,结果可能是致命的,因为严重的脱水。


三、驻梯度渗透流

3.1 等渗和逆渗透压运输的问题

为了使水通过渗透转运发生,乍一看细胞外液必须是高渗的,即比细胞质具有更高的溶质浓度。否则,水如何能通过渗透被挤出细胞?然而,多年来已知许多类型的上皮(如肠道和肾脏)能够在没有渗透梯度的情况下运输水,甚至能够抵抗渗透梯度运输。此外,运输的流体本身非常接近等渗。这些观察都与简单的溶质偶联转运理论难以协调。

3.2 Diamond-Bossert 模型

为了解释这些观察, Diamond 和 Bossert(1967)提出了驻梯度渗透流模型(standing-gradient osmotic flow model)。该模型分析了具有长而窄通道的上皮细胞,通道一端开放(流体从中流出),另一端封闭。这可以是侧细胞间隙(相邻细胞之间的长而窄的通道),或者是基底膜的折叠。

3.3 模型建立

假设溶质偶联的水转运发生在一个长度为 \(L\)、半径为 \(r\) 的长而窄的圆柱形管中,一端封闭(图18.4)。溶质主动转运过管边界,从细胞质进入管内,水被动跟随,由渗透压差驱动。

\(N(x)\) 为单位面积单位时间内主动转运的速率,\(v(x)\)\(x\) 点处的流体速度,\(c(x,t)\) 为管内溶质的浓度。在圆柱形管中,标准守恒论证给出:

\[\frac{\partial c}{\partial t} = D\frac{\partial^2 c}{\partial x^2} - \frac{\partial (vc)}{\partial x} + \frac{2N(x)}{r} \quad (18.14)\]

类似地,水沿管的速度为 \(v\),单位面积的水流入速率由渗透压 \(P(c(x) - c_0)\) 驱动,其中 \(c_0\) 是圆柱管外各处的溶质浓度。在稳态下:

\[\frac{dv}{dx} = \frac{2P}{r}(c(x) - c_0) \quad (18.16)\]

边界条件为:

\[v(0) = 0, \quad \frac{dc}{dx}(0) = 0 \quad (18.17)\]
\[c(L) = c_0 \quad (18.18)\]

3.4 溶液的主要特征

通过积分可得:

\[v(L)c_0 - D\frac{dc}{dx}(L) = \frac{2}{r}\int_0^L N(x)dx \quad (18.19)\]
\[v(L) = \frac{2LP}{r}(\bar{c} - c_0), \quad \bar{c} = \frac{1}{L}\int_0^c dx \quad (18.20)\]

出口溶质浓度定义为:

\[c_e = \frac{2\pi r \int_0^N(x)dx}{\pi r^2 v(L)} = \frac{2\int_0^N(x)dx}{rv(L)} \quad (18.21)\]

从(18.19)可得:

\[c_e = c_0 - \frac{D\frac{dc}{dx}(L)}{v(L)} \quad (18.22)\]

因此,出口流体速度越高,出口浓度越低,出口流体越接近等渗。

3.5 小扩散极限的奇异扰动分析

对于参数值 \(D = 1000 \mu m^2 s^{-1}\)\(r = 0.05 \mu m\)\(c_0 = 0.3 \mu M\)\(P = 0.2 \mu m s^{-1} \mu M^{-1}\)\(L = 100 \mu m\),有 \(\epsilon = 0.08\)

\(\epsilon = 0\),得到简化系统:

\[\frac{du_o}{dy} = w_o \frac{du_o}{dy} + 2(u_o - 1) = 2n(y) \quad (18.27)\]
\[\frac{dw_o}{dy} = 2(u_o - 1) \quad (18.28)\]

外部解不能满足 \(y = 1\) 处的边界条件。为满足 \(u(1) = 1\),我们在 \(y = 1\) 附近寻求边界层校正。引入缩放变量 \(\xi = (y-1)/\sqrt{\epsilon}\),得到边界层方程。复合解为:

\[u(y) = u_o(y) + (1 - u_o(1))\exp\left(\frac{w_o(1)(y-1)}{\sqrt{\epsilon}}\right) + O(\sqrt{\epsilon}) \quad (18.41)\]

四、逆水运输

4.1 模型扩展

在驻梯度模型中,假设圆柱管外的背景浓度处处相同为 \(c_0\)。该模型表明,溶质浓度的驻梯度可以在具有相同渗透压的区域之间运输水,且运输的流体可以接近该背景渗透压。

然而,在许多情况下,上皮细胞层能做更令人惊讶的事情——它能抵抗渗透压差运输水,同时只利用渗透作用。如何实现这一点可以通过对驻梯度模型的简单修改来看。

4.2 分离的细胞内和细胞外区域

代替假设圆柱管外的浓度处处相同,我们将其分为细胞内和细胞外区域,浓度分别为 \(c_0\)\(c_1\)(图18.7),且 \(c_0 > c_1\)。模型方程保持不变,只是 \(x = L\)\(c\) 的边界条件变为 \(c(L) = c_1\)

4.3 净流量方向

如果管内的平均浓度大于细胞内的平均浓度,则开放端有净正流体速度;反之,如果管内的平均浓度小于细胞内的平均浓度,则开放端有净负流体速度,将水吸入细胞。

在小扩散(长管长度)极限中,净流量始终为正,不管 \(c_1\) 有多小。这是因为在小扩散极限中,浓度仅在薄的边界层中低于 \(c_0\),净速度受到的影响很小。

在大扩散(小管长度、小渗透性)极限中,流量可以双向流动:

\[\bar{c} = c_1 + \frac{N_0 L^2}{rD}\left(\frac{1}{3} - \frac{a^2}{3}\right) \quad (18.45)\]

五、胃保护

5.1 黏液和碳酸氢盐的保护作用

消化道内表面是一层柱状上皮细胞,主动分泌黏液和富含碳酸氢盐的液体。黏液高度黏稠,覆盖细胞形成0.5–1.0毫米厚的不溶性层,创造了一个不被其他胃分泌物溶解的润滑边界。该层细胞表面有许多胃小凹(图18.8)。每个胃小凹含有壁细胞,通过主动转运过程分泌盐酸,导致胃腔 pH 约为1。由于供应上皮表面的血液 pH 约为7.4,跨每个上皮细胞的 H⁺ 浓度梯度约为百万倍。显然,上皮细胞必须受到保护免受高腔内酸性的影响。

5.2 稳态模型

为模拟胃保护,假设胃黏膜的腔面是位于 \(x = 0\) 的平面,\(x\) 是垂直于黏膜壁测量的坐标,黏液层厚度均匀为 \(l\)。因此,黏液-腔界面位于 \(x = l\)(图18.9)。

在黏液层内,H⁺ 和 HCO₃⁻ 按以下反应:

\[H^+ + HCO_3^- \rightleftharpoons H_2O + CO_2 \quad (18.46)\]

在稳态下:

\[D_{H^+}\frac{d^2 H^+}{dx^2} = D_{HCO_3^-}\frac{d^2 HCO_3^-}{dx^2} = k_+ H^+ HCO_3^- - k_- [CO_2] \quad (18.48)\]
\[D_{CO_2}\frac{d^2 [CO_2]}{dx^2} = -k_+ H^+ HCO_3^- + k_- [CO_2] \quad (18.49)\]

5.3 边界条件

在腔侧(\(x = l\)),假设 \(H^+ = H_l^+\)\([CO_2] = [CO_2]_l\) 恒定,\(HCO_3^-\) 由平衡关系给出:

\[HCO_3^-]_l = \frac{k_- [CO_2]_l}{k_+ [H^+]_l} \quad (18.50)\]

在上皮侧,假设 \(HCO_3^-\)\(CO_2\) 的通量是已知的:

\[D_{HCO_3^-}\frac{d HCO_3^-}{dx} = -\bar{J} \quad (18.51)\]
\[D_{CO_2}\frac{d [CO_2]}{dx} = -\bar{I} \quad (18.52)\]

\(x = 0\) 处,\(H^+\) 的通量与浓度差成正比:

\[D_{H^+}\frac{d H^+}{dx} = P_{H^+} H^+ \quad (18.54)\]

5.4 无量纲化和奇异扰动

引入无量纲变量 \(y = x/l\)\(u = H^+/H_l^+\)\(v = HCO_3^-/H_l^+\)\(w = [CO_2]/H_l^+\),模型变为:

\[\epsilon \frac{d^2 u}{dy^2} = uv - \zeta w \quad (18.55)\]
\[\epsilon \frac{d^2 u}{dy^2} = \gamma \frac{d^2 v}{dy^2} = -\beta \frac{d^2 w}{dy^2} \quad (18.56)\]

其中 \(\beta = D_{CO_2}/D_{H^+}\)\(\gamma = D_{HCO_3^-}/D_{H^+}\)\(\epsilon = D_{H^+}/(k_+ l^2 [H^+]_l)\)\(\zeta = k_-/(k_+ [H^+]_l)\)

参数估计:\(\beta \approx 0.14\)\(\gamma \approx 0.13\)\(k_- = 11 s^{-1}\)\(k_+ = 2.6 \times 10^{10} cm^3 \cdot mol^{-1} \cdot s^{-1}\)\(\bar{J} = 1.4 \times 10^{-10} mol \cdot cm^{-2} \cdot s^{-1}\)\(H_l^+ = 140 mM\)\(l = 0.05 cm\)\(D_{H^+} = 1.75 \times 10^{-5} cm^2 \cdot s^{-1}\)\(P_{H^+} = 1.3 \times 10^{-5} cm \cdot s^{-1}\)。从这些参数值可得 \(\epsilon = O(10^{-7})\)\(\zeta = O(10^{-6})\) 是小参数。

5.5 外部解

\(u = u_0 + \epsilon u_1 + \cdots\),最低阶方程为:

\[0 = f(u_0, y) \quad (18.66)\]
\[\frac{du_0}{dy}(0) = \lambda u_0(0), \quad u_0(1) = 1 \quad (18.67)\]

由于 \(\zeta\) 很小,发现:

\[u_0(y) = 1 - \frac{\lambda + J}{\lambda + 1}(1 - y) + O(\zeta) \quad (18.72)\]
\[w(y) = \alpha + \frac{J+I}{\beta}(1 - y) + O(\zeta) \quad (18.73)\]

因此,最低阶近似下,黏液层中没有 HCO₃⁻,H⁺ 和 CO₂ 沿黏液层线性变化。

5.6 角层校正

外部解不满足 \(y = 0\) 处的边界条件。必须包括角层校正,即在上皮表面附近斜率变化很大的解:

\[\epsilon^{1/2} e^{-\mu y/\epsilon^{1/2}} \quad (18.74)\]

5.7 生理意义

角层的物理意义:最低阶近似下,每个离开上皮的碳酸氢盐分子立即与一个氢离子反应,导致氢离子消失。因此,化学反应的全部发生在角层内。

在上皮表面:

\[H^+_{x=0} = \frac{D_{H^+} H_l^+ - \bar{J}l}{D_{H^+} + P_{H^+}l} \quad (18.79)\]

使用实验参数值,\(H^+_{x=0} = 135 mM\),仅下降3.5%,太小而无法保护上皮细胞免受高腔内酸性的影响。因此,这个黏液层的简单模型不足以解释上皮层是如何受到保护的。

5.8 胃酸分泌与中和

上述胃保护模型的主要困难是碳酸氢盐通量 \(J\) 太小,无法充分减少上皮细胞表面的氢离子。一个通过检查盐酸分泌与碳酸氢盐释放之间关系的模型由 Lacker 和他的同事构建。

盐酸由泌酸腺的壁细胞分泌。首先,细胞内的水分离解为氢离子和羟基离子。氢离子通过 H⁺–K⁺ ATP 酶主动分泌。同时,氯离子主动分泌,Na⁺ 主动吸收,分别通过不同的 ATP 酶。结果是腔内高浓度盐酸。同时,二氧化碳与羟基离子结合(在碳酸酐酶催化下)形成碳酸,进而形成碳酸氢盐。碳酸氢盐从细胞扩散到细胞外介质,并通过毛细血管血流运输。

方向:毛细血管血流从胃小凹中的泌酸细胞流向腔的上皮 lining。由于上皮细胞在胃小凹的下游,它们从血液中吸收碳酸氢盐,然后分泌到黏液中。因此,随着酸产量的增加,上皮细胞分泌到腔内的碳酸氢盐速率也增加。

碱性潮汐(alkaline tide):随着 H⁺ 分泌到腔内的速率增加,下游 H⁺(即胃静脉 H⁺)减少。这强化了该模型的主要思想——上皮细胞分泌 HCO₃⁻ 是由泌酸细胞分泌 H⁺ 驱动的,因此当腔内 H⁺ 增加时,胃保护自动增加。


六、小肠中的耦合振荡器

6.1 收缩的时间控制

整个消化道的平滑肌细胞在其膜电位中表现出振荡,周期从2到40次/分钟不等。这种电控制活动(ECA)的典型例子如图18.11所示。虽然细胞膜电位的去极化可引起肌肉收缩,但这仅在膜电位去极化超过阈值时发生,在这种情况下,电位开始以更高频率振荡或爆发。是否发生爆发取决于神经元或化学刺激的水平。电爆发称为电反应活动(ERA)。肌肉收缩的频率不能高于 ECA 的频率,因此局部 ECA 的特性限制了可能的肌肉收缩类型。

6.2 电活动的波动

如果每个局部振荡器与其邻居不耦合,预期不会有沿肠道传播的有组织收缩波。然而,本节的主要观点是:振荡器之间的弱耦合导致 ECA 和 ERA 沿肠道的波动传播。

频率 plateau:在完整的肠道中,ECA 频率在幽门附近整个区域是恒定的,尽管内在频率沿该区域稳步下降。在约60厘米处,ECA 频率开始下降。在频率 plateau(恒定频率区域),每个振荡器与其邻居锁相,导致沿肠道远离幽门传播的有组织波。

相位锁定:振荡以相同频率振荡不一定意味着波动行为。相位波发生在沿肠道长度从一个点到下一个点有恒定的相位提前(或延迟)时。

当各段解耦时,每个都表现出振荡 ECA,但内在频率沿远离幽门的方向下降;频率 plateau 在分离的段中消失。在完整的肠道中,最靠近幽门的最高频率段使附近的振荡器同步,这些振荡器具有相似但较低的频率。然而,当内在频率差异太大时,同步是不可能的:频率 plateau 崩溃,波失去规律性。

6.3 耦合振荡器模型

数值模拟:许多研究者使用数值模拟来研究小肠中耦合振荡器链的行为。Diamant 等(1970)将5到25个 van der Pol 振荡器耦合,频率沿链下降。每个振荡器与其频率较低的最近邻居耦合,耦合是阻性的。数值模拟表明,振荡器组织成频率 plateau,其长度随耦合强度增加而增加。

相位方程:对电活动的数学研究从假设有 \(n+1\) 个耦合振荡器开始:

\[\frac{du_i}{dt} = F_i(u_i) + \epsilon \sum_{j=1}^{n+1} a_{ij} H(u_j) \quad (18.80)\]

使用多尺度或平均技术,这些方程可以简化为描述各振荡器相位偏移的方程:

\[\frac{d}{d\tau}\delta\theta_i = \xi_i - \frac{1}{2\pi}\left(1 + \sum_{j=1}^{n+1} a_{ij} h(\delta\theta_j - \delta\theta_i)\right) + O(\epsilon^2) \quad (18.81)\]

其中 \(\tau = \epsilon t\) 是慢时间。

当每个振荡器仅在线性链中与其最近邻居耦合时:

\[\frac{du_i}{dt} = F_i(u_i) + \epsilon(u_{i+1} - u_i) + \epsilon(u_{i-1} - u_i) \quad (18.82)\]

相邻相位差 \(\phi_i = \delta\theta_{i+1} - \delta\theta_i\) 的方程为:

\[\frac{d\phi_i}{d\tau} = \Delta_i + [h(\phi_{i+1}) + h(-\phi_i) - h(\phi_i) - h(-\phi_{i-1})] + O(\epsilon^2) \quad (18.83)\]

最后,相位差方程变为:

\[\dot{\phi} = \beta\Delta + KH(\phi) \quad (18.84)\]

6.4 频率 Plateau 的形成

Ermentrout 和 Kopell(1984)的研究表明:在每个 plateau 上,相位差不完全恒定,而是做小的振荡,仅在平均意义上锁定,这种现象有时称为相位捕获(phase trapping)。

临界点方程的解存在当且仅当 \(K^{-1}(-\beta\Delta)\) 的每个分量介于 \(m\)\(M\) 之间。当 \(\beta\) 越过临界值 \(\beta_0\) 时,一对 plateau 出现;随着 \(\beta\) 进一步增加,多个 plateau 出现。

6.5 相位波与实际波的差异

重要的是要注意,\(k \neq j\) 的相位差 \(\phi_k\)\(\phi_{k-}\) 附近做小振荡,但不恒等于 \(\phi_{k-}\)。因此,在每个频率 plateau 上,相位不是锁定的,而只是在每个周期中平均"捕获"。因此,与一些实验数据不同,不应期望在每个 plateau 上看到完全规则的传播波,因为这样的波要求在每个瞬间相位锁定并沿 plateau 恒定相位差。

6.6 Cajal 间质细胞

近年来,胃肠道的电活动控制涉及的因素远比耦合振荡器复杂。特别是,Cajal 间质细胞(ICC)在 ECA 调节中起极其重要的作用。

ICC 由 Santiago Ramon y Cajal 于1911年首次发现,他称之为"原始神经元",并理论认为它们是平滑肌收缩的调节器,受神经系统调节。ICC 可分为两组:肠肌间 ICC(在大多数胃肠道区域,在纵肌和环肌层之间形成薄层)和肌内 ICC(分布在平滑肌细胞本身中)。

现在已知 ICC 是驱动胃肠道平滑肌电慢波的起搏细胞,分离的平滑肌细胞不能独立振荡。ICC 还受肠神经系统调节;看起来胃肠道区域的神经支配更多指向 ICC 而非平滑肌本身。

ICC 与心脏之间的重要相似性:两者中,神经元都支配起搏细胞,起搏细胞以某种方式协调其活动以控制肌肉运动。在起搏细胞和平滑肌细胞中,内部 Ca²⁺ 的控制(通过 IP3 受体、Ryanodine 受体和 Ca²⁺ 敏感离子通道)起关键作用。

6.7 生物物理和解剖学模型

Miftakhov 等的模型:迄今为止最详细的胃肠道电活动生物物理模型。该模型包括肠神经系统的许多不同神经元(每个都有树突、胞体和轴突成分)以及突触中神经递质和 Ca²⁺ 的动力学。

Pullan 等的模型:采用不同方法,构造了复杂的解剖学胃模型。使用来自可见人项目的数据,构造了几何上准确的胃,在此之上可以叠加有限元网格并求解模型方程。然而,他们仅使用高度简化的生物物理模型( FitzHugh-Nagumo 方程)来模拟电活动。

胃电图(EGG):胃的电活动可以通过胃电图测量,就像心电图测量心脏的电活动一样。然而,即使在今天,EGG 记录,特别是旨在测量肠道电活动的记录,也难以进行和解释。每个心跳是一个单一的、局部的、相对精确的、连续重复多次的事件。胃肠道电活动没有这些属性,它是弥散的、沿胃肠道长度性质和频率变化的、持续发生的。


七、公式汇总表

编号 公式名称 公式 描述
(18.1) Na⁺ 被动通量 \(J = g(n_i - n_l)\) Na⁺ 从腔到细胞内部的被动流动
(18.2) Na⁺ 主动通量 \(J = f(n_i)\) Na⁺–K⁺ ATP 酶主动转运
(18.3) 渗透水流 \(Rq = n - n_l\) 水通过紧密连接的流动
(18.4) 稳态 Na⁺ 守恒 \(g(n_l - n_i) = f(n_i)\) 被动和主动 Na⁺ 转运平衡
(18.8) 无量纲水流方程 \(\rho y^2 + (\rho + u_l)y + \kappa = 0\) 水流问题的无量纲形式
(18.9) 无量纲浓度关系 \(u_l = u_i + \beta F(u_i)\) 腔浓度与细胞浓度的关系
(18.15) 溶质守恒(驻梯度) \(D\frac{d^2 c}{dx^2} - \frac{d}{dx}(vc) + \frac{2N(x)}{r} = 0\) 圆柱管内溶质稳态方程
(18.16) 水流速度方程 \(\frac{dv}{dx} = \frac{2P}{r}(c(x) - c_0)\) 渗透驱动水流
(18.20) 出口速度 \(v(L) = \frac{2LP}{r}(\bar{c} - c_0)\) 管端流体速度
(18.21) 出口浓度 \(c_e = \frac{2\int_0^N(x)dx}{rv(L)}\) 运输流体的溶质浓度
(18.46) 碳酸氢盐缓冲反应 \(H^+ + HCO_3^- \rightleftharpoons H_2O + CO_2\) 胃黏液中的酸碱反应
(18.55) 无量纲 H⁺ 方程 \(\epsilon \frac{d^2 u}{dy^2} = uv - \zeta w\) 胃保护的偏微分方程
(18.72) 外部解(H⁺) \(u_0(y) = 1 - \frac{\lambda + J}{\lambda + 1}(1 - y)\) 黏液层 H⁺ 分布
(18.80) 耦合振荡器方程 \(\frac{du_i}{dt} = F_i(u_i) + \epsilon \sum a_{ij} H(u_j)\) 振荡器链的动力学
(18.84) 相位差方程 \(\dot{\phi} = \beta\Delta + KH(\phi)\) 相位锁定分析的核心方程
(18.101) 平均频率 \(\omega(\epsilon) + \lim_{T\to\infty}\frac{1}{T}\int_0^T \delta\theta_k'(t)dt\) 第 k 个振荡器的平均频率
(18.45) 逆水运输平均浓度 \(\bar{c} = c_1 + \frac{N_0 L^2}{rD}(\frac{1}{3} - \frac{a^2}{3})\) 短管极限的平均溶质浓度

八、关键概念与生理意义总结

8.1 流体吸收

水通道蛋白(aquaporins)的发现彻底改变了我们对水跨膜转运的理解。水通道蛋白允许水特异性通过,而大多数其他溶质和离子被排斥在外。这表明水通过细胞膜的转运不仅仅是简单的渗透驱动,而是受到高度调控的过程。驻梯度渗透流模型解释了上皮细胞如何能够运输等渗流体甚至抵抗渗透梯度运输水——关键在于溶质主动泵入狭窄的侧细胞间隙,在该空间中维持高于细胞质的溶质浓度,从而将水吸入该空间并沿流体流出。

8.2 胃保护

胃上皮面临一个显著挑战:腔内 pH 约为1,而血液 pH 约为7.4,跨上皮细胞的 H⁺ 浓度梯度约为一百万倍。黏液-碳酸氢盐屏障提供化学保护,黏液层作为扩散屏障,碳酸氢盐与渗透过来的 H⁺ 中和。然而,简单的扩散-反应模型表明,仅靠这一机制不足以解释观察到的上皮保护程度。胃酸分泌与碳酸氢盐释放之间的耦合——"碱性潮汐"现象——提供了一个自动调节的保护机制,当酸分泌增加时,碳酸氢盐分泌也增加。

8.3 肠动力学的耦合振荡器

消化道运动依赖于平滑肌的协调收缩,而这种协调部分由称为电控制活动(ECA)的内在电振荡控制。重要的是,这些局部振荡器通过 gap junctions 弱耦合,形成空间扩展系统。这种耦合导致有趣的空间组织现象——频率 plateau,其中相邻振荡器锁定到公共频率,产生沿肠道传播的规则波。Ermentrout-Kopell 的数学分析揭示了这种 plateau 形成的机制,涉及相位方程的不变流形和"相位捕获"现象。Cajal 间质细胞(ICC)的发现进一步揭示了胃肠道电活动的复杂性——ICC 是真正的起搏细胞,类似于心脏的窦房结。

8.4 建模方法的启示

本章展示了生理学建模的几个重要方面:简化模型的价值(即使失去定量准确性,简单模型也能揭示基本机制)、奇异扰动技术在处理多尺度问题中的威力、耦合振荡器理论在理解空间组织模式中的应用,以及复杂生物系统建模中 identifiability 和模型复杂性的权衡问题。


参考书籍:Keener & Sneyd, Mathematical Physiology (II: Systems Physiology), Springer 2009, Chapter 18

笔记整理日期:2026年5月11日