跳转至

第十五章：肌肉读书笔记

作者

本章作者为 James P. Keener 与 James Sneyd。本章把 Vol I Ch 7（钙动力学）和 Ch 12（心肌细胞）建立的"电-化学-力学"工具，扩展到肌肉收缩的横桥动力学。从最经典的 A.V. Hill 1938 力-速度经验公式出发，经 Huxley 1957 离散结合位点模型，到 Peskin 1975/1986 连续结合位点模型的"逆问题"——从实验数据反推横桥速率函数 \(f(x)\)、\(g(x)\)、\(r(x)\)。本章的数学核心是守恒律 PDE (15.21) \(\partial n/\partial t - v(t)\partial n/\partial x = (1-n)f(x) - n g(x)\) 及其在稳态下的力-速度闭式解 (15.37)。

内容概述

本章分为 5 大主题：(1) 15.1 兴奋-收缩耦联——T-tubule + SR + RyR 的 Ca²⁺ 释放链、Ca²⁺ 与肌动蛋白/肌球蛋白相互作用的分子级描述（薄/厚纤维、A 带/I 带/H 带、Z 线、sarcomere ~2.5 μm）、等长张力-长度曲线；(2) 15.2 Hill 模型——力-速度公式 (15.1) \((p+a)v = b(p_0 - p)\)、二元素模型（contractile + series elastic）、Jewell-Wilkie 1958 数据揭示 Hill 模型的失效（步长响应过冲）；(3) 15.3 Huxley 1957 模型——每个横桥只有一个结合位点、\(f(x)\) 和 \(g(x)\) 的线性选择 (15.26-15.27)、稳态解 (15.36) \(n(x; v)\)、力-速度公式 (15.37)、\(\phi = (f_1 + g_1)h/2\) 的特征速度；(4) 15.4 连续结合位点模型（Lacker-Peskin 1986）——结合位点连续分布（"flypaper"模型）、逆问题：从 \(p(v)\) 和 \(\phi(v)\) 实验数据反推 \(F\)、\(g(x)\)、\(r(x)\) 通过 Laplace 变换；(5) 15.5-15.7 钙调控、平滑肌、布朗棘轮（本章未完整展示）。

本章最具数学深度的部分是 15.4.3 逆问题——给定实验测得的 \(p(v)\) 和 \(\phi(v)\)，通过 \(\phi(v) = \rho \epsilon F(1 - N)\) 算出 \(F\)、通过 \(I(v) = \int_h^{-\infty} \exp[\int_x^h g(s)/v \, ds] dx\) 算出 \(I(v)\)、然后用 \(\bar{g}(y) = (1/2\pi i) \int I(1/\sigma) e^{\sigma y} d\sigma\) 的 Laplace 逆变换得到 \(g(x)\)。这是"实验数据→模型参数"的逆向建模范式，是数学生理学的标志性方法。

核心方程与概念

15.1 肌肉结构与兴奋-收缩耦联

sarcomere 结构：~2.5 μm 长；薄（actin + tropomyosin + troponin）+ 厚（myosin）纤维；六角排列（6 薄绕 1 厚）；A 带（含 myosin）、I 带（仅 thin）、H 区（仅 thick）。等长张力-长度曲线（Fig. 15.7，Gordon 1966）：张力在 ~2.0-2.2 μm 长度最大，过短/过长均下降（重叠太多/太少横桥无法形成）。

EC 耦联序列：AP → T-tubule 去极化 → L-type Ca²⁺ 通道（DHPR）→ Ca²⁺ 内流 → RyR 释放 SR Ca²⁺（CICR）→ 细胞质 [Ca²⁺] ↑ → 肌钙蛋白（troponin）结合 Ca²⁺ → 原肌球蛋白（tropomyosin）构象改变 → 暴露肌动蛋白（actin）结合位点 → 横桥结合 → 动力冲程 → ATP 水解 → 重复。

15.2 Hill 力-速度关系

Hill 公式（15.1）\((p + a)v = b(p_0 - p)\)，\(v = 0\) 时 \(p = p_0\)（等长张力），\(p = 0\) 时 \(v = b p_0/a\)（最大缩短速度）。二元素模型（15.2-15.5）：contractile element（满足 Hill 关系）+ series elastic element（线性，\(P = \alpha(x - x_0)\)），总张力 \(p = P(x)\)，总长度 \(L = l + x\)。ODE（15.5）\(dp/dt = \alpha[dL/dt + b(p_0 - p)/(p + a)]\)。

Jewell-Wilkie 1958 实验揭示 Hill 模型的缺陷：(a) 步长响应后张力恢复太快（模型预测 vs 实测，Fig. 15.11B）；(b) 存在初始振荡。这是Hill 模型的失败——它假设力-速度关系在张力变化后瞬时满足，但实际有弛豫时间。

15.3 Huxley 1957 横桥模型

变量：\(n(x, t)\) = 位移 \(x\) 处的横桥结合分数，\(x > 0\) 产生收缩力，\(x < 0\) 产生阻力，\(x = 0\) 是无应力位置。PDE 守恒律（15.21）\(\partial n/\partial t - v(t) \partial n/\partial x = (1-n) f(x) - n g(x)\)。力（15.23）\(p = \rho \int_{-\infty}^\infty r(x) n(x, t) dx\)（\(r(x) = kx\) 弹簧模型）。能量释放率（15.22）\(\phi = \rho \epsilon \int_{-\infty}^\infty (1-n) f(x) dx\)。

Huxley 速率函数选择（15.26-15.27）： - \(f(x) = 0\)（\(x < 0\)）、\(f_1 x/h\)（\(0 \leq x \leq h\)）、\(0\)（\(x > h\)）——横桥只在 \(0 \leq x \leq h\) 内结合 - \(g(x) = g_2\)（\(x \leq 0\)，快速解离）、\(g_1 x/h\)（\(x > 0\)，与位移正相关）

稳态解（15.36）三段函数（\(x < 0\)、\(0 \leq x \leq h\)、\(x > h\)）。力-速度公式（15.37） $\(p = \frac{\rho k f_1}{f_1 + g_1} \frac{h^2}{2} \left[ 1 - \frac{v}{\phi}(1 - e^{-\phi/v}) + \frac{1}{2G_2^2} \frac{v}{\phi} \right]\)$ 其中 \(\phi = (f_1 + g_1)h/2\) 是特征速度。Huxley 拟合参数 \(F_1 = f_1/(g_1+f_1) = 13/16\)、\(G_2 = g_2/(f_1+g_1) = 3.919\)，对 \(p_0/a \approx 4\) 的 Hill 数据有极好吻合。

15.3.1 等张响应

张力阶跃（15.38-15.41）：肌纤维等长时 \(n_s(x) = f(x)/(f(x)+g(x))\)。张力突然降到 \(p_1\) → 长度阶跃 \(\Delta L\) → 新分布 \(n_s(x + \Delta L)\) → 通过 \(p_1 = \int k(x) n_s(x + \Delta L) dx\) 求 \(\Delta L\)。约束（15.41） $\(-v(t) = \frac{\int k(x) [(1-n) f(x) - n g(x)] dx}{\int k(x) \partial n/\partial x dx}\)$ Civan-Podolsky 1966 实验：相对长度变化 1/300，每个横桥变化 < 10 nm，Huxley 模型无法重现初始振荡——这是 Huxley 模型的"已知缺陷"。

15.4 连续结合位点模型

"flypaper" 模型（Lacker-Peskin 1986）：actin 结合位点连续分布，未结合的横桥无 \(x\) 概念。总结合分数 \(N = \int n dx < 1\)。PDE（15.46）\(\partial n/\partial t - v(t) \partial n/\partial x = (1-N) f(x) - n g(x)\)，关键差异是 结合速率正比于 \((1-N)\) 而不是 \((1-n)\)。

简化形式（15.47-15.48）：所有横桥在 \(x = h\) 处结合，\(f(x) = F \delta(x - h)\)，\(n(h, t) = F(1-N)/v\)。

15.4.3 逆问题

思路：实验测得 \(p(v)\)（力-速度）和 \(\phi(v)\)（能量释放率）→ 反推 \(F\)、\(g(x)\)、\(r(x)\)。关键步骤： - 由 \(\phi_0 = \rho \epsilon F(1 - N_0)\) → \(F\) 已知 - \(I(v) = \int_h^{-\infty} \exp[\int_x^h g(s)/v \, ds] dx\) → \(I(v)\) 是 \(1/\bar{g}(y)\) 的 Laplace 变换（15.71） - \(\bar{g}(y) = (1/2\pi i) \int I(1/\sigma) e^{\sigma y} d\sigma\) → 通过 Laplace 逆变换得 \(g\) - 类似地，\(r(x)\) 也从 \(p(v)\) 反推 (15.74)

Hill 数据拟合（15.75-15.78）：用 Hill 的 \(p(v) = (b p_0 - a v)/(v + b)\) 和 \(\phi(v) = \phi_0 + b v(a + p_0)/(v+b)\) 代入，得 \(1/\bar{g}(y)\) 的闭式 Laplace 变换。

15.8 布朗棘轮（部分内容）

闪烁棘轮（flashing ratchet）：非平衡态下，周期性切换的棘轮势 + 扩散产生净粒子流，无需外力梯度。这是热力学第二定律与"分子马达定向输运"的关键交叉点。Eyring 速率理论将势垒高度转换为反应速率。

关键结论

  1. sarcomere 长度 2.5 μm，含 ~\(10^9\) 横桥/cm³ 组织；等长张力-长度曲线峰值 ~2.0-2.2 μm。
  2. Hill 公式 \((p+a)v = b(p_0 - p)\)：1938 年力-速度经验关系，\(a\)、\(b\) 由数据拟合。
  3. Hill 模型在步长响应上过冲——Jewell-Wilkie 1958 实验揭示模型与实测差异，推动横桥动力学发展。
  4. Huxley 1957 模型是第一个从横桥分子机制解释力-速度关系的模型：结合速率 \(f(x)\) 在 \(x > 0\) 时增加，解离速率 \(g(x)\) 在 \(x < 0\) 时增加——这就是横桥的"棘轮"机制。
  5. Huxley 拟合参数：\(F_1 = 13/16\)、\(G_2 = 3.919\)，对应 \(\phi = 400\) nm/s/半 sarcomere，最大缩短速度 1.28 cm/s（与 Hill 1938 4.8 cm/s 吻合）。
  6. Huxley 模型也无法重现步长响应的初始振荡——Podolsky-Nolan 1972-1973 通过修改 \(f\)、\(g\) 修复。
  7. 连续结合位点（flypaper）模型是 Huxley 模型的 \(\Delta x \to 0\) 极限（Peskin 1975）；\(n\) 和 \(f\) 的量纲不同。
  8. 逆问题：从 \(p(v)\) 和 \(\phi(v)\) 实验数据反推 \(F\)、\(g(x)\)、\(r(x)\)，核心是 Laplace 变换 \(I(v) \leftrightarrow 1/\bar{g}(y)\)。
  9. 闪烁棘轮：非平衡态下分子马达可以通过周期势切换 + 扩散产生净输运——这是 ATP 水解驱动肌球蛋白动力冲程的热力学基础。
  10. 力学-化学耦合：每个横桥循环水解一个 ATP，释放 ~\(k_B T\) 自由能，转化为 ~10 pN·nm 机械功（power stroke）。

挑战和开放性问题

  1. Huxley 模型的初始振荡失效——Podolsky-Nolan 通过修改 \(f\)、\(g\) 函数修复，但修改后函数缺乏分子基础。"应该用什么样的 \(f\)、\(g\)" 至今无共识。
  2. 心肌的 Frank-Starling 机制——15.1 末尾明确说"the same is not true for cardiac muscle, which considerably complicates theoretical studies"。心肌的力-长度依赖是心输出量调节的核心，但 Huxley 框架对心肌的扩展不完整。
  3. 平滑肌的 latch state——平滑肌可在低 ATP 消耗下维持张力（latch state），但 15.1 节说"the physiology is considerably different"。门控横桥 + latch 横桥的双态模型（Khaki & Murphy 1991 等）未在本章展示。
  4. 钙调控的非平稳性——15.3 节假设 f 和 g 与 [Ca²⁺] 无关（tetanus 状态），但生理状态下 [Ca²⁺] 动态变化（Ch 7 的钙瞬变）。将 Ca²⁺ 动力学耦合到横桥动力学是开放问题。
  5. 肌节的异质性——肌节长度不均匀性导致弱肌节先断裂（sarcomere length non-uniformity, SLNU），这可能是肌源性力-速度曲线下凹的物理原因。Campbell 2015 等人的工作未在本章展示。
  6. 横桥的分子动力学——15.4 节给出"过氧化 ATP → ADP + Pi" 的化学描述，但Pi 释放顺序、ADP 释放顺序、ATP 结合顺序的精确时间常数仍未在数学模型中完全整合。
  7. 横桥协同性——多横桥之间的力学耦合（一个横桥的力影响相邻横桥的结合概率）通常被忽略。真实肌肉中横桥是协同的，相邻结合位点距离 ~38 nm 横桥间相互作用不可忽略。
  8. 完整肌肉的 3D 力学——sarcomere 模型被简化为 1D 沿纤维方向，忽略了肌原纤维的横向力学和Z 盘 / 肌膜的横向耦合。整个肌纤维的 3D 有限元模型与 1D 横桥模型的耦合是开放问题。

个人反思与批判性分析

本章的最大优点在于 15.4.3 逆问题——这是数学生理学的标志性方法：从可测量的整体行为（\(p(v)\)、\(\phi(v)\)）反推不可测量的微观动力学（\(f\)、\(g\)、\(r\)）。这一方法在实验生物物理中价值巨大：它把"显微观察"和"宏观测量"通过 Laplace 变换联系起来，避免了对单个横桥的实时追踪（目前还不可能）。

本章的主要缺陷在于 15.3.1 步长响应的"修补"——Podolsky-Nolan 通过修改 \(f\)、\(g\) 函数修复 Huxley 模型的振荡，但这种"为拟合而调整"缺乏分子基础。这与 Ch 12 心肌离子通道模型的发展路径类似：Noble 1962 → MNT 1975 → Beeler-Reuter 1977 → Luo-Rudy 1994，每一代都"为拟合而调整"，最终积累了"参数化经验"。何时是"机理模型"、何时是"经验拟合"是科学哲学问题，本章未触及。

15.2 节的 Hill 模型和 15.3 节的 Huxley 模型是肌力学的"牛顿力学 vs 量子力学"对比——Hill 模型是现象学（黑箱二元素 + 经验公式），Huxley 模型是机理（横桥分子机制）。两者的关系本质上是宏观与微观的桥梁——Hill 的 \(a\)、\(b\) 可以从 Huxley 的 \(f_1\)、\(g_1\)、\(h\) 推出（\(\phi \approx b\)，\(p_0/a \approx 4\)）。

对我自己研究血管生物力学的启发：(1) 横桥力-速度曲线与血管平滑肌的力-速度曲线在形式上完全相同（Hill 公式适用）——血管壁在正常生理范围内主要通过横桥动力学调节张力，钙瞬变通过调节横桥结合概率来调节收缩力；(2) 闪烁棘轮机制可能解释血管壁的分子马达（如 kinesin、myosin V 在细胞骨架上的输运）如何产生长程定向输运而无需外力梯度；(3) Peskin 连续结合位点模型的逆问题方法可以推广到血管壁的弹性反问题——从压力-直径数据反推血管壁的多层结构（内膜、中膜、外膜）的弹性模量，与临床血管成像结合可建立个体化的血管力学模型。

值得指出，本章对 ATP 能量学的讨论非常简略。单个 ATP 水解释放约 80 kJ/mol = 130 zJ/molecule，但单个横桥动力冲程做功约 10 pN·nm = 10 zJ——效率约 8%，其余以热释放。这与肌肉热力学效率（~25%）的差异是"内部摩擦"（同一时刻有多个横桥处于不同结合状态，互相干扰）。这种集体效应（不是单个横桥的效率）才是真实肌肉效率的来源。

重要参考文献

[X1] Hill, A. V. (1938). The heat of shortening and the dynamic constants of muscle. Proceedings of the Royal Society B, 126, 136-195. — 15.2 节力-速度公式 (15.1) 的原始论文；肌肉力学的奠基之作。 [X2] Huxley, A. F. (1957). Muscle structure and theories of contraction. Progress in Biophysics and Biophysical Chemistry, 7, 255-318. — 15.3 节 Huxley 1957 模型的原始论文；横桥动力学的开创性 PDE 模型。 [X3] Jewell, B. R., & Wilkie, D. R. (1958). An analysis of the mechanical components in frog's striated muscle. Journal of Physiology, 143, 515-540. — 15.2.2 节 Hill 模型失败的关键实验；推动从 Hill 模型到 Huxley 模型的转变。 [X4] Peskin, C. S. (1975). Mathematical aspects of heart physiology. PhD thesis, Albert Einstein College of Medicine. — 15.4 节连续结合位点模型和逆问题的原始论文。 [X5] Lacker, H. M., & Peskin, C. S. (1986). A mathematical model of the volution of myosin cross-bridge states. In Some Mathematical Questions in Biology — Muscle Physiology (American Mathematical Society). — 15.4.1 节连续结合位点模型的另一独立提出。 [X6] Gordon, A. M., Huxley, A. F., & Julian, F. J. (1966). The variation in isometric tension with sarcomere length in vertebrate muscle fibres. Journal of Physiology, 184, 170-192. — 15.1 节 Fig. 15.7 等长张力-长度曲线的经典实验数据。 [X7] Podolsky, R. J., & Nolan, A. C. (1972, 1973). Series of papers on crossbridge modification to fit step response data. Journal of Physiology and Proceedings of the National Academy of Sciences USA. — 15.3.1 节末提及的 Huxley 模型修补。 [X8] Civan, M. M., & Podolsky, R. J. (1966). Contraction kinetics of striated muscle fibres following quick changes in load. Journal of Physiology, 184, 511-534. — 15.3.1 节 Huxley 模型失效的实验。 [X9] Brokaw, C. J. (1976). Computer simulation of movement-generating cross-bridges. Biophysical Journal, 16, 1013-1027. — 15.3 末 Huxley 模型参数选取的来源。 [X10] White, D. C. S., & Thorson, J. (1975). The kinetics of muscle contraction. Progress in Biophysics and Molecular Biology, 27, 175-255. — 15.1-15.2 节 Fig. 15.4, 15.8, 15.10, 15.11 数据的综合出处。 [X11] Berne, R. M., & Levy, M. N. (1993). Physiology (3rd ed.). Mosby. — 15.1 节 Fig. 15.1, 15.2 肌细胞和 sarcomere 图示的来源。 [X12] Hill, T. L. (1974, 1975). Free Energy Transduction in Biology. Academic Press. — 15.4 节末提到的横桥生物化学热力学的经典专著。 [X13] Pate, E., & Cooke, R. (1989). A model of crossbridge action: the effects of ATP, ADP and Pi. Journal of Muscle Research and Cell Motility, 10, 181-196. — 15.4 节末横桥能量学的现代模型。 [X14] Eisenberg, E., & Hill, T. L. (1978). A cross-bridge model of muscle contraction consistent with the law of the heart. Progress in Biophysics and Molecular Biology, 33, 55-82. — Hill 的热力学横桥模型与心肌 Frank-Starling 定律的整合。 [X15] Bers, D. M. (2002). Cardiac excitation-contraction coupling. Nature, 415, 198-205. — 15.1 节 Fig. 15.3 AP + Ca²⁺ 瞬变 + 收缩三阶段图的来源。