跳转至

第十六章：内分泌系统读书笔记

作者

本章作者为 James P. Keener 与 James Sneyd。本章从 Vol I Ch 1（生化反应）、Ch 4（耦合振荡器）和 Ch 13（造血延迟微分方程）的工具，扩展到内分泌系统——激素脉冲分泌（脉冲周期从毫秒到小时）、下丘脑-垂体-性腺反馈环、月经周期的多室 ODE 模型（Selgrade-Schlosser 1999、Clark 2003）、排卵数控制的全局耦合模型（Lacker 1981, Lacker-Peskin 1981）。本章的数学核心是 16.2.2 Lacker 排卵模型——一个由 \(N\) 个卵泡通过血液循环中雌二醇全局耦合的 ODE 系统，排卵数严格控制在 \(M_1 < M < 2M^*\) 之间，并解释了"为什么大多数卵泡开始发育但最终退化"的生理效率问题。

内容概述

本章分为 3 大主题：(1) 16.1 下丘脑和垂体——下丘脑-垂体的解剖（前叶、后叶、垂体门脉系统）、激素的脉冲分泌（Table 16.1，胰岛素 108-144 脉冲/天、LH 90-121 脉冲/天）、GnRH 脉冲频率对 LH/FSH 释放的调节、Smith 1980 三室 LH 模型、神经脉冲发生器（Brown 1994 用 FHN 方程模拟 LHRH 神经元）；(2) 16.2 哺乳动物排卵——Selgrade-Schlosser 月经周期模型（10 个卵泡发育阶段 + 6 激素反馈）、Lacker 排卵模型（\(N\) 个卵泡通过雌二醇全局耦合）、排卵数的稳定性分析（\(M_1 < M < M_2\) 内可行，\(M^* < M < 2M^*\) 内稳定）；(3) 16.3 胰岛素和葡萄糖（本章未完整展示）。

本章最具洞察力的部分是 Lacker 排卵模型的"大种群-精确时间"效应（16.2.2 末）——卵泡数量从 10 增加到 1000 时，平均排卵时间几乎不变，但排卵时间的标准差显著降低。这解释了为什么绝经期前女性（可用卵母细胞少）通常出现月经不规律——大种群的统计正则性是精确排卵时间的关键。

核心方程与概念

16.1 下丘脑-垂体-性腺轴

解剖：垂体前叶（腺垂体，分泌 LH/FSH/TSH/ACTH/GH/PRL 等）、垂体后叶（神经垂体，储存并释放催产素/加压素）。垂体门脉系统（Fig. 16.1B）：下丘脑神经元分泌释放/抑制激素 → 初级毛细血管丛 → 门静脉 → 次级毛细血管丛 → 垂体前叶内分泌细胞。反馈：垂体后叶分泌的激素可经反向血流回到下丘脑，形成短反馈环。

激素类型：(a) 胺类（肾上腺素、去甲肾上腺素，源自酪氨酸）；(b) 多肽/蛋白质（水溶性，结合膜表面受体）；(c) 固醇类（脂溶性，源自胆固醇，穿过细胞膜结合胞内受体调控基因转录）。信号通路：(a) cAMP/PKA 通路（ACTH、LH、FSH、加压素）；(b) IP3/Ca²⁺ 通路（Ch 7）；(c) 酪氨酸激酶通路（胰岛素）。

Table 16.1 激素脉冲分泌（脉冲数/天）：胰岛素 108-144、胰高血糖素 103-144、生长抑素 72、LH 90-121、TSH 13、ACTH 15、褪黑素 18-24、醛固酮 6-12、皮质醇 15-39。

16.1.1 LH 脉冲分泌（Smith 1980 模型）

三室模型（16.1-16.3）\(\dot{R} = f(T) - b_1 R\)（LHRH），\(\dot{L} = g_1 R - b_2 L\)（LH），\(\dot{T} = g_2 L - b_3 T\)（睾酮），\(f(T)\) 是睾酮对 LHRH 的负反馈函数。Murray 2002 修改：加入延迟 \(g_2 L(t-\tau)\) 反映 LH 通过血液到睾丸的传输时间。

16.1.2 神经脉冲发生器（Brown 1994）

模型（16.4-16.6）\(\epsilon dv/dt = v(c-v)(v-1) - k_1 g + k_2 a(t)\)（LHRH 神经元，本征双稳态），\(dg/dt = k_3 a(t) + k_4 v - r_g g\)（GABA 神经元，负反馈），\(dz/dt = p(v) - r_z z\)（LH 浓度）。关键观察：这是 FitzHugh-Nagumo 方程的扩展（Ch 5），但兴奋变量 \(v\) 和抑制变量 \(g\) 由同一输入 \(a(t)\) 驱动。频率依赖性：低频输入 → 每个尖峰触发一次脉冲；高频输入 → 抑制占优，\(v\) 不再激活；中间频率 → 不规则脉冲模式。

16.1.3 生长激素（GH）模型

Brown 2004 / MacGregor-Leng 2005 模型：GHRH（兴奋）+ 生长抑素（抑制）共同调节 GH 释放，GH 通过 somatomedins（肝脏分泌）反馈调节下丘脑。实验采样伪随机性：每 10 分钟采样的真实 GH 释放看起来很"随机"，但实际是确定性模型的低频输出。

16.2 月经周期与排卵

16.2.1 Selgrade-Schlosser 月经周期模型

两子模型（pituitary + ovary）：

Pituitary 子模型（16.7-16.10）：\(dLHR/dt = [V_{0L} + V_{1L} E^8/(K_{mL}^8 + E^8)] / (1 + P(t-\tau)/K_{iL}) - k_L (1 + c_{LP} P) LHR / (1 + c_{LE} E)\)。关键特征：(a) 雌二醇 \(E\) 刺激 LH 合成（正反馈）但抑制 LH 释放（负反馈）；(b) 孕酮 \(P\) 抑制 LH 合成但促进 LH 释放；(c) 延迟 \(\tau = 2\) 天模拟血液运输。

Ovary 子模型（16.11-16.19）：10 个卵泡发育阶段（Inactive → Recruited → Secondary → Preovulatory → Ovulatory scars → Corpus Luteum），每个阶段有不同激素分泌率。关键阶段转换：(a) FSH 招募 inactive 卵泡（\(b \cdot FSH\)）；(b) LH 刺激 secondary 卵泡；(c) 5 个时间常数 \(k\) 形成 5 个线性滤波器延迟孕酮分泌（\(S_1, S_2, L_1, L_2, L_3, L_4\)）。

激素输出（16.20-16.22）\(E = e_0 + e_1 SeF + e_2 PrF + e_3 L_4\)，\(P = p(L_3 + L_4)\)，\(I = h_0 + h_1 PrF + h_2 L_3 + h_3 L_4\)。典型周期 28 天，含 LH surge 和多峰雌二醇曲线。

多稳态与 PCOS：模型有两个稳定极限环——正常周期（小振幅、长周期）和 PCOS 周期（小振幅、短周期、低孕酮、高 LH/FSH 比）。外源孕酮足够大可切换——Clark 2003 提议这可能是 PCOS 治疗的基础。

16.2.2 Lacker 排卵模型

核心假设：(1) 卵泡分泌雌二醇的速率是卵泡成熟度的指标；(2) 雌二醇浓度控制 FSH/LH 释放；(3) FSH/LH 控制卵泡成熟率。

关键 ODE（16.23-16.27）\(d\xi_i/dt = f(\xi_i, \xi) = \xi_i \varphi(\xi_i, \xi)\)，\(\varphi(\xi_i, \xi) = 1 - (\xi - M_1 \xi_i)(\xi - M_2 \xi_i)\)，\(\xi = \sum \xi_i\)（总雌二醇浓度，\(i = 1, ..., N\) 个卵泡）。\(M_1, M_2\) 是模型参数，决定排卵数范围。

关键数学性质（16.39）\(\ln(\xi_i - \xi_j)\) 单调性——成熟度的初始顺序永远不变（如果 \(\xi_i(0) > \xi_j(0)\)，则 \(\xi_i(t) > \xi_j(t)\)）。这一单调性使得 Lacker 引入重整化变量 \(\gamma_i = \xi_i / \xi_1\)、\(\tilde{\xi} = \xi / \xi_1\) 和重整化时间 \(\tau = \int_0^t \xi_1^2 ds\)。

对称解（16.29-16.30）：\(M\) 个卵泡有相同成熟度时 \(d\xi/dt = \xi + \mu \xi^3\)，\(\mu = -(1 - M_1/M)(1 - M_2/M)\)。关键结果：(a) \(M_1 < M < M_2\) 时 \(\mu > 0\)，\(\xi\) 在有限时间 \(T = \log[(1 + \mu \xi_0^2)/(\mu \xi_0^2)]\) 爆破——排卵发生；(b) \(M < M_1\) 或 \(M > M_2\) 时 \(\mu < 0\)，\(\xi\) 趋于稳态 \(\xi_M = 1/\sqrt{(M - M_1)(M - M_2)}\)——不排卵，卵泡维持中间成熟度。

稳定排卵数（16.48-16.49, 16.16）：变换坐标的特征值 \(\lambda_1 = (M_1 + M_2)M - 2M_1 M_2\)、\(\lambda_2 = -(M_1 + M_2)M + M_1 M_2\)。\(M^* < M < 2M^*\) 时两个特征值均为负——稳定排卵数范围 \([M_1, 2M^*]\)。例如 \(M_1 = 3.85\)、\(M_2 = 15.15\) → \(2M^* = 6.14\) → 实际观测排卵数 4 或 5（与数值模拟 Fig. 16.15 吻合）。

"大种群-精确时间"效应（Fig. 16.17）：种群从 10 → 100 → 1000 时，平均排卵时间几乎不变，但排卵时间标准差显著降低——这是大种群统计正则性的体现。

关键结论

  1. 激素脉冲分泌极其普遍：胰岛素 108-144 脉冲/天、LH 90-121 脉冲/天、生长抑素 72 脉冲/天——脉冲频率是生理调节的关键变量。
  2. GnRH 脉冲频率调节 LH/FSH 比例：恒定 GnRH 释放会减少促性腺激素分泌——这是临床脉冲式 GnRH 给药治疗不孕症的基础。
  3. 神经脉冲发生器（Brown 1994）：LHRH 神经元的本征双稳态 + GABA 抑制反馈 → 模拟 FHN 方程的脉冲动力学。频率依赖性是关键——低频输入触发脉冲，高频输入被抑制。
  4. 月经周期模型的两个稳定极限环：正常周期 vs PCOS 周期，外源孕酮足够大可切换——这是 PCOS 治疗假说的基础。
  5. Lacker 排卵数控制：排卵数严格在 \(M_1 < M < 2M^*\) 之间，其中 \(M^* = M_1 M_2/(M_1 + M_2)\)，且大种群更精确控制排卵时间。
  6. 成熟度顺序不变性（16.39）：\(\xi_i - \xi_j\) 的符号保持——这是 Lacker 重整化变换的数学基础。
  7. 大种群-精确排卵时间（Fig. 16.17）：10 → 1000 个卵泡时，标准差显著降低——解释了绝经期前女性月经不规律。
  8. Selgrade-Schlosser 月经周期通过 6 个 ODE + 5 个延迟链 + 反馈函数重现 28 天周期、LH surge、多峰雌二醇等关键特征。
  9. "大部分卵泡退化"的进化意义：不是浪费，而是"大种群统计正则性"的代价——少数卵泡最终排卵是被"精选"的结果。
  10. 从经验拟合到机理模型：内分泌系统建模大量使用Hill 函数（\(V_1 E^8/(K_m^8 + E^8)\)）和延迟（\(P(t-\tau)\)），但缺乏分子水平的离子通道/受体动力学——这是与神经/心肌建模的本质差异。

挑战和开放性问题

  1. PCOS 的真正机制——Clark 2003 的双稳态是假设性的，临床上 PCOS 的异质性（胰岛素抵抗、高 LH、高雄激素）需要更精细的子模型。
  2. Lacker 模型的局限——(a) 没有显式 FSH/LH 浓度；(b) \(M_1, M_2\) 是经验参数没有分子解释；(c) 仅考虑雌二醇反馈，忽略了抑制素、激活素等其他调节因子。
  3. 胰岛素分泌的钙振荡——16.3 节未在本章展示，但胰岛 β 细胞的钙振荡 → 胰岛素脉冲是 Ch 7 + Ch 16 的交叉点，需要整合两章的数学工具。
  4. 垂体门脉系统的传输延迟——不同激素的延迟不同（LH 延迟 2 天，TSH 延迟分钟），但精确的多延迟耦合模型仍未标准化。
  5. 卵巢衰老与卵母细胞耗竭——Fig. 16.17 提到"为什么绝经期前女性月经不规律"，但 Lacker 模型没有显式追踪卵母细胞数量随年龄减少——需要把种群大小 \(N\) 变成时变量 \(N(t)\)。
  6. 男性 vs 女性内分泌差异——Lacker 模型本质上是对称的，但性别二态性（LH 脉冲频率不同、雌二醇/睾酮主导不同）需要性别特异参数。
  7. 神经-内分泌-免疫交叉——本章只讨论内分泌-内分泌反馈，神经-免疫-内分泌网络（如应激对月经的影响、感染对甲状腺的影响）未涉及。
  8. Lacker 模型在临床辅助生殖（IVF）中的应用——Sarty-Pierson 2005 提到 Lacker 模型用于预测卵巢超刺激反应，但临床决策仍主要基于经验，超刺激失败率高（20-30%）。
  9. 月经周期的种族/个体差异——模型参数 \(V_{0L}, V_{1L}, K_{mL}, K_{iL}\) 等（Table 16.2）来源于特定人群的拟合，参数随种族、年龄、BMI 的变异未量化。
  10. 激素的旁分泌/自分泌作用——本章所有激素作用都是内分泌（通过血液），卵巢内局部旁分泌（如 inhibin 在卵泡液中的局部高浓度）在数学上难以整合。

个人反思与批判性分析

本章的最大优点在于 Lacker 排卵模型——它用一个简洁的 \(N\) 维 ODE 系统（\(N\) 是卵泡数）解释了排卵数的精确控制这一长期困扰生理学家的现象。关键洞察：(a) 排卵数不是被"局部"控制（去除一侧卵巢排卵数不变），而是被全局雌二醇反馈控制；(b) 大种群提供"统计正则性"——排卵时间精度与卵泡数正相关；(c) 数学上的稳定排卵数 \(M_1 < M < 2M^*\) 与不同物种的窝仔数（猪 6-11、羊 2、人 1）有潜在定量对应。

本章的主要缺陷在于 "经验拟合"远多于"机理建模"——Selgrade-Schlosser 模型的 \(V_{0L}, V_{1L}, K_{mL}, K_{iL}\)（Table 16.2）等 30+ 参数都是通过数据拟合得到，缺乏分子水平的受体动力学或信号通路级联。这与 Ch 5（HH 动作电位，3-4 参数）、Ch 11（Guyton 循环，5-6 参数）的"少参数模型"形成鲜明对比。根本原因是内分泌的测量不确定性（激素浓度受采样时间、脉冲相位、个体差异影响）远大于神经电生理（毫秒级电压钳）。

16.1.1 节对 Smith 1980 LH 模型的介绍非常精炼——但临床 LH 脉冲检测（每 5-10 分钟抽血）是侵入性操作，为什么 GnRH 神经元需要脉冲式分泌而非连续分泌至今未完全理解。本章给出 3 种可能解释：(a) 受体脱敏避免；(b) 信号整合效率；(c) 神经元动力学的本征要求——这三种假设都没有被完全验证。

对我自己研究血管生物力学的启发：(1) Lacker 模型的"大种群-精确时间"效应在血管生物学中有直接对应——血管树的分支数（主动脉约 20 级分支、毛细血管约 \(10^9\) 根）提供了巨大的"统计池"，使血压的局部调节（如 myogenic response、metabolic regulation）能在整体血压的恒定上得到平滑；(2) Selgrade-Schlosser 模型的反馈环可以推广到血管-内分泌耦合（如 RAAS 系统：肾素-血管紧张素-醛固酮），用类似 ODE 框架建模；(3) LH 脉冲检测的采样伪随机性（16.1.3 末）的洞察——任何周期性的生理信号在低频采样下都会呈现"随机"特征——这对我们设计血管活性检测的采样频率有直接指导意义。

值得指出，本章没有涵盖：(a) 肾上腺皮质（皮质醇、醛固酮的昼夜节律）；(b) 甲状腺（T3/T4 的反馈调节）；(c) 生长激素-IGF-1 轴；(d) 生殖内分泌的种族差异。这些都是临床内分泌学的核心主题，但本章选择了"激素脉冲 + 排卵"这两个最有数学深度的案例。

重要参考文献

[X1] Smith, W. R. (1980). Hypothalamic regulation of pituitary secretion of luteinizing hormone. Bulletin of Mathematical Biology, 42, 35-78. — 16.1.1 节三室 LH 模型的原型；反馈压抑延迟微分方程。 [X2] Cartwright, M., & Husain, M. (1986). A model for the effect of estrogen on the pulsatile release of luteinizing hormone. Journal of Theoretical Biology, 123, 239-250. — LH 雌激素反馈模型的扩展。 [X3] Murray, J. D. (2002). Mathematical Biology II: Spatial Models and Biomedical Applications (3rd ed.). Springer. — 16.1.1 节末提及的延迟 LH 模型。 [X4] Brown, D., Herbison, A. E., Robinson, J. E., Marrs, R. W., & Leng, G. (1994). Modelling the luteinizing hormone-releasing hormone pulse generator. Journal of Neuroendocrinology, 6, 55-63. — 16.1.2 节神经脉冲发生器 FHN 模型的原型。 [X5] Knobil, E. (1981). Patterns of hypophysiotropic signals and gonadotropin secretion in the rhesus monkey. Biology of Reproduction, 24, 44-49. — GnRH 脉冲频率对 LH/FSH 调节的实验证据。 [X6] Wildt, L., Marshall, G., & Knobil, E. (1981). Experimental induction of puberty in the infantile rhesus monkey. Science, 207, 1373-1375. — GnRH 脉冲频率与排卵恢复的临床实验。 [X7] Lacker, H. M. (1981). Regulation of ovulation number in mammals: A follicle interaction hypothesis. Journal of Mathematical Biology, 13, 233-248. — 16.2.2 节 Lacker 排卵模型的开创性论文。 [X8] Lacker, H. M., & Peskin, C. S. (1981). A model of the hypothalamic pulse generator. Journal of Mathematical Biology, 13, 153-161. — Lacker 模型与下丘脑脉冲发生器的整合。 [X9] Selgrade, J. F., & Schlosser, P. M. (1999). A model for the hypothalamic-pituitary-ovarian axis with a detailed description of feedback regulation. Mathematical Biosciences, 158, 1-31. — 16.2.1 节 Selgrade-Schlosser 月经周期模型的原型。 [X10] Clark, L. H., Schlosser, P. M., & Selgrade, J. F. (2003). Multiple stable periodic solutions in a model for the hypothalamic-pituitary-ovarian axis. Journal of Mathematical Biology, 46, 218-252. — 16.2.1 节 PCOS 双稳态极限环的发现。 [X11] Bertram, R., Egli, M., Toporikova, N., & Freeman, M. E. (2006). A mathematical model for the mating-induced prolactin rhythm. Biological Reproduction, 74, 50-57. — 16.1.3 末提及的催乳素节律模型。 [X12] Brown, D., MacGregor, R., & Leng, G. (2004). Data-driven modeling of the growth hormone pulse generator. Journal of Theoretical Biology, 228, 179-189. — 16.1.3 节 GH 脉冲发生器模型。 [X13] MacGregor, D. J., & Leng, G. (2005). Modelling the hypothalamic control of growth hormone secretion. Journal of Neuroendocrinology, 17, 788-803. — GH 模型的进一步发展。 [X14] Brabant, G., Prank, K., & Schöfl, C. (1992). Pulsatile hormone secretion: Generation, gating and modulation. Trends in Endocrinology and Metabolism, 3, 101-108. — Table 16.1 激素脉冲分泌频率的综述。 [X15] Sarty, G. E., & Pierson, R. A. (2005). An application of a dynamic follicle model to the prediction of ovulation in the mare. Animal Reproduction Science, 86, 213-225. — Lacker 模型在辅助生殖中的应用。 [X16] MacLachlan, E. A., Santoro, N., Land, H., et al. (1990). The measurement of gonadotropin and estradiol levels over the menstrual cycle. Fertility and Sterility, 53, 769-774. — Fig. 16.9 月经周期激素数据来源。 [X17] Knobil, E. (1980). The neuroendocrine control of the menstrual cycle. Recent Progress in Hormone Research, 36, 53-88. — 月经周期的下丘脑-垂体-卵巢轴控制的综述。 [X18] Liu, B., & Deng, D. (1991). A mathematical model for the hypothalamic-pituitary-ovarian axis in the rat. Journal of Biological Systems, 1, 5-20. — 第一个与实验数据拟合的内分泌模型。