跳转至

第二十章：内耳读书笔记

作者

本章作者为 James P. Keener 与 James Sneyd。本章把 Vol I Ch 4（被动电流、空间常数）、Ch 5（动作电位、Hodgkin-Huxley 模型）和 Ch 7（钙动力学）的工具，扩展到听觉系统——耳蜗宏观力学（基底膜行波）、Lesser-Berkley 1972 解析模型、Shallow-water 和 Deep-water 近似、Ranke vs Zwislocki 之争、毛细胞电共振（Ca²⁺ 激活的 K⁺ 通道）、Hudspeth-Lewis 1988 详细生物物理模型、耳蜗非线性放大器（Gold 1948 假说、外毛细胞 prestin 蛋白）。本章的数学核心是 20.2 耳蜗波动方程 (20.7-20.8, 20.16) 和 20.3 毛细胞电共振电路 (20.85) \(d^2V/dt^2 + \gamma dV/dt + \omega_0^2 V = f(t)\)——前者是反应-扩散 PDE（流体 Laplace 方程 + 基底膜 ODE）的耦合，后者是带通滤波器（RLC 等效电路），频率响应 \(|V_1(\omega)| = 1/|\omega_0^2 - \omega^2 + i\gamma\omega|\)。

内容概述

本章分为 4 大主题：(1) 20.1 频率调谐——Helmholtz 1875 位置理论（place theory）、von Békésy 1960 行波实验（基底膜不同位置对不同频率响应）、三种调谐机制（毛细胞机械调谐、基底膜机械调谐、毛细胞电调谐）；(2) 20.2 耳蜗模型——Ranke 1950 deep-water 近似、Zwislocki 1965 shallow-water 近似、Lesser-Berkley 1972 完整 2D 模型（流体 Laplace + 基底膜 ODE + Fourier 级数解析解）、Siebert 1974 直接力学加载、Steele-Homes 弹性板模型、Peskin 1981 共形映射精确解；(3) 20.3 毛细胞电共振——Ashmore-Attwell 1985 等效电路、电压敏感 K⁺ 通道（Q 因子太低）、Ca²⁺ 激活 K⁺ 通道（Q 因子足够高，Hudspeth-Lewis 1988 验证）、3 闭态 + 2 开态 Markov 模型；(4) 20.4 非线性耳蜗放大器——Gold 1948 假说（耳蜗必须主动放大低强度信号）、外毛细胞 prestin 蛋白的电致伸缩、毛束负刚度导致自发振荡、Martin 2000 慢适应机制。

本章最具洞察力的部分是 20.4.1 毛束的负刚度——毛束在小位移时表现出负刚度（位移增加反而产生回复力减小），这与正常的弹性系统相反。这种负刚度是自发振荡的基础——类比于 Van der Pol 振子，毛束的负阻尼 + 正常流体阻尼 → 极限环振荡。这与 Ch 4 的负阻尼方程、Ch 5 的 Hopf 分岔有共同的数学结构。

核心方程与概念

20.1 频率调谐

三种调谐机制：(a) 毛细胞机械调谐（蜥蜴的毛束长度梯度，基底-顶端逐渐变长，对低频敏感）；(b) 基底膜机械调谐（哺乳动物，宽度 + 刚度梯度，基底-顶端刚度指数递减，长度常数 ~7 mm）；(c) 毛细胞电调谐（龟耳蜗、牛蛙球囊，离子通道的带通响应）。

Helmholtz 1875 位置理论：耳蜗由一组调谐元件（每点谐振于特定频率）构成，输入频率由最大刺激点位置决定。von Békésy 1960 行波实验直接证实——基底膜上形成行波，振幅先增后减，峰值位置依赖频率（高频 → 接近镫骨；低频 → 接近蜗顶）。

基尔霍夫电流守恒（20.1-20.8）\(\rho \partial u/\partial t + \nabla p = 0\)、\(\nabla \cdot u = 0\)（无粘不可压缩流体）；\(\rho \partial \phi/\partial t + p = 0\)、\(\nabla^2 \phi = 0\)（无旋流的势函数 \(\phi\)）。

20.2 耳蜗波动模型

20.2.2 Lesser-Berkley 模型

核心方程（20.11-20.16）：上室（scala vestibuli）\(\rho \partial \phi_1/\partial t + p_1 = 0\)、\(\nabla^2 \phi_1 = 0\)；下室（scala tympani）\(\rho \partial \phi_2/\partial t + p_2 = 0\)、\(\nabla^2 \phi_2 = 0\)；基底膜（20.16）\(m(x) \partial^2\eta/\partial t^2 + r(x) \partial\eta/\partial t + k(x) \eta = p_2 - p_1\)（\(m\) 质量、\(r\) 阻尼、\(k\) 刚度，每点都是独立阻尼振子）；边界（20.14-20.19）\(\partial\eta/\partial t = \partial\phi/\partial y\)（y = 0），\(\partial\phi/\partial y = 0\)（y = l，顶端不动），\(\partial\phi_1/\partial x = \partial F(y, t)/\partial t\)（x = 0，卵圆窗），\(\partial\phi_1/\partial x = 0\)（x = L，远端）。

频率域解（20.20-20.34）：设 \(F(y, t) = \hat{F}(y) e^{i\omega t}\)，\(\phi = \hat{\phi}(x, y; \omega) e^{i\omega t}\) → 偏微分方程化为代数方程，\(\hat{p} + i\omega\rho \hat{\phi} = 0\)。Fourier 解（20.30-20.32）\(\phi = x(1 - x/2) - \sigma y(1 - y/(2\sigma)) + \sum_n A_n \cosh[n\pi(\sigma - y)] \cos(n\pi x)\)，\(A_m \alpha_m = f_m\)。典型结果（Fig. 20.7）：行波包络在 \(x \approx 0.5\) 处最大，随频率增加峰值向基部移动。

20.2.4 Shallow-water 近似（Zwislocki）

核心简化（20.43-20.46）\(-i\omega \eta \approx (l/(i\omega\rho)) d^2p/dx^2\)（短波长极限下 \(p\) 与 \(\eta\) 之间的二阶微分关系）。\(p\) 的 ODE（20.44）\(Y(x) p(x) = (l/(i\omega\rho)) d^2p(x)/dx^2 + \eta_0 \delta(x) - \eta_L \delta(x - L)\)，其中 \(Y = 2/Z\) 是导纳。

多尺度分析（20.47-20.65）：\(d^2p/dz^2 + g^2(\epsilon z) p(z) = 0\)，\(\epsilon = \lambda/(2q) \ll 1\)（\(Y\) 缓慢变化）。WKB 解（20.58-20.65）\(p = (1/\sqrt{g(\sigma)})(A_0 e^{iG(z)} + B_0 e^{-iG(z)})\)，\(G(z) = \int_0^z g(\epsilon z) dz\)。特殊情况（20.65）峰值位置 \(x_p = -(1/\lambda) \ln(16\rho\omega/(9l\lambda^2 r_0))\)——频率每翻倍，峰值位置按对数移动。

20.2.4 Deep-water 近似（Ranke）

简化（20.66-20.78）：设 \(l \gg L\) → \(\tanh(2\pi n l/L) \approx \text{sign}(n)\)。\(p\) 的一阶 ODE（20.73-20.75）\(Y(x) p_\pm(x) \approx \pm (1/\omega\rho) dp_\pm/dx\)，\(p_\pm(x) = A_\pm \exp[\pm\omega\rho \int_0^x |Y(\zeta)| d\zeta]\)。

典型结果（20.83-20.84，Fig. 20.9）\(|\eta| = 2\eta_0 |\xi| \exp[\lambda x + \beta_r (1 - e^{\lambda x})]\)，峰值位置 \(x_p = -(1/\lambda) \ln[2\omega^3 \rho r_0 / (\lambda(k_0^2 + \omega^2 r_0^2))]\)。关键缺陷：模型预测的相位累积远大于实验观测。

20.2.5 Peskin 1981 精确解

特殊条件：\(r(x), k(x) \propto e^{-\lambda x}\)，流体容器高度 \(\lambda l = \pi/2\)，无限长基底膜。方法：共形映射（conformal mapping）+ 围道积分——获得了真正的精确闭式解，是后续 WKB 近似（Lesser-Berkley）的对照基准。

20.3 毛细胞电共振

20.3.1 简化电路模型

核心方程（20.85-20.89）\(d^2V/dt^2 + \gamma dV/dt + \omega_0^2 V = f(t)\)，\(\gamma = g_p/C_m + R/L\)，\(\omega_0^2 = g_p R + 1/(LC_m)\)。频率响应（20.89）\(|V_1(\omega)| = 1/|\omega_0^2 - \omega^2 + i\gamma\omega|\)，带通滤波器。

关键观察（20.90-20.91）\(Q = \omega_0/\gamma\) 是品质因子——越大越尖锐。实验观测：\(Q\) 通常 ≥ 5。

20.3.2 电压敏感 K⁺ 通道（Fig. 20.10B）

电路（20.94-20.101）\(I = C_m dV/dt + g_L(V - V_L) + f g_K(V - V_K)\)，\(\tau df/dt = f_\infty - f\)，\(f_\infty = f_r + \mu(V - V_r)\)。线性化（20.98）等价于 (20.85)，但 \(Q\) 仅约 0.7（实验要求 ≥ 5）——生理参数不够尖锐。

20.3.3 Ca²⁺ 激活 K⁺ 通道（Fig. 20.10C）

电路（20.102-20.106）\(I = C_m dV/dt + g_L(V - V_L) + g_K k c(V - V_K)\)，\(W dc/dt = I_r/F + \theta(V - V_r)/F - p c\)。关键改进（20.104-20.105）\(L = WF/(g_K(V_r - V_K) k\theta)\)、\(R = pF/(g_K(V_r - V_K) k\theta)\)——\(Q\) 可达实验要求（减小 \(p\) 或增大 \(k\theta\) 即可）。Ashmore-Attwell 1985 结论：Ca²⁺ 激活 K⁺ 通道是毛细胞频率调谐的主要机制。

20.3.4 Hudspeth-Lewis 1988 详细模型

3 闭态 + 2 开态 Markov 模型（20.107）\(I = C_m dV/dt + I_c + I_{kc} + I_L\)。关键参数（约 30 个）通过单细胞电压钳数据拟合。典型结果（Fig. 20.12）：对阶跃电流的阻尼振荡，最大响应 ~112 Hz，\(Q \approx 3\)。参数敏感：共振频率敏感于 \(g_{kc}\)，毛细胞可以通过调节 Ca²⁺ 激活 K⁺ 通道的数量来调谐自己的频率响应。

20.4 非线性耳蜗放大器

Gold 1948 假说：耳蜗对低强度信号的极高灵敏度不可能来自纯被动机制——流体阻尼会"在蜂蜜中抑制音叉"——必须存在主动机械放大。现代证实（Hudspeth 1997, Nobili 1998, Manley 2001）：外毛细胞通过 prestin 蛋白的电致伸缩主动收缩/伸长（不需要 ATP，直接由膜电位驱动），放大基底膜振动 ~100 倍，提高频率选择性 ~10 倍。

毛束负刚度（Martin 2000）：毛束在小位移时负刚度（\(F = -kx\) 变 \(F = +kx\) 在负刚度情况下）——这是自发振荡的基础。这与 Ch 4 的负阻尼方程、Ch 5 的 Hopf 分岔有共同的数学结构。

关键结论

  1. von Békésy 1960 行波实验直接证实 Helmholtz 位置理论——基底膜不同位置对不同频率响应。
  2. Lesser-Berkley 1972 完整模型用 Fourier 级数解析解再现行波包络（Fig. 20.7）——包络峰值随频率增加向基部移动。
  3. Shallow-water（Zwislocki）与 Deep-water（Ranke）近似通过 WKB 多尺度分析给出峰值位置公式 \(x_p \propto \ln\omega\)。
  4. Peskin 1981 共形映射精确解是无穷长基底膜 + 特殊参数的闭式解，是 WKB 近似的对照。
  5. 毛细胞电共振（20.85）是带通滤波器，\(Q = \omega_0/\gamma\) 衡量调谐尖锐度。
  6. Ca²⁺ 激活 K⁺ 通道（不是电压敏感 K⁺ 通道）是毛细胞调谐的主要机制——Ashmore-Attwell 1985 比较两种模型后确定。
  7. Hudspeth-Lewis 1988 3 闭态 + 2 开态 Markov 模型重现毛细胞阻尼振荡和 ~112 Hz 共振频率。
  8. Gold 1948 假说 + 外毛细胞 prestin 蛋白的主动放大是耳蜗极高灵敏度（受限于热噪声）的关键。
  9. 毛束负刚度（Martin 2000）是自发振荡的物理基础——与 Van der Pol 振子的负阻尼同构。
  10. 位置理论 + 电调谐 + 主动放大三机制协同：基底膜粗调（位置）、毛细胞电共振细调、prestin 蛋白放大。

挑战和开放性问题

  1. 基底膜的精确弹性模型——20.2 节展示的 2D 流体 + 1D 基底膜简化是教科书级别，真实的基底膜是 3D 弹性板 + organ of Corti + 内外毛细胞耦合，完整 3D 有限元模型（Steele 1974, Ramamoorthy 2007）才接近真实。
  2. 毛细胞的非线性响应——20.4 节只展示了主动放大的定性描述，非线性 ODE/PDE 模型（如 fast-slow 双时间尺度）尚未标准化。
  3. prestin 蛋白的分子动力学——20.4.1 说"exact identity and properties of these proteins remains elusive (Ashmore et al., 2000)"，20 年后仍未完全解决。
  4. 毛束的负刚度机制——Martin 2000 假说需要更精细的肌动蛋白-肌球蛋白动力学模型支撑。
  5. 耳蜗的反馈-前馈耦合——外毛细胞放大的正反馈如何在避免自激振荡（耳声发射）的同时保持稳定？需要主动控制理论视角。
  6. 听力损失——老年性聋（presbycusis）、噪声性聋、药物性聋（庆大霉素）的数学模型——毛细胞、神经节细胞、神经纤维的退化动力学。
  7. 前庭系统——20 节开头提到前庭系统（半规管、椭圆囊、球囊）检测运动和加速度，但其数学模型（壶腹嵴、椭圆囊斑的流体力学 + 毛细胞转导）未在本章展示。
  8. 听觉神经编码——听觉神经纤维的频率-强度-时间响应（频率调谐曲线、锁相、适应性）的概率模型。
  9. 中枢听觉处理——耳蜗核、上橄榄复合体、外侧丘系、下丘、内侧膝状体、初级听皮层的多级信息处理。
  10. 声源定位——双耳时间差（ITD）、双耳强度差（ILD）的 Jeffress 模型及其在鸟类和哺乳动物中的实现。

个人反思与批判性分析

本章的最大优点在于 20.3.2 Ca²⁺ 激活 K⁺ 通道是电调谐的主要机制——这是一个反直觉但通过数学分析建立的结论。直觉上"电压敏感 K⁺ 通道"似乎更直接负责电共振，但计算表明其 Q 因子仅 0.7（远低于实验要求的 ≥5），而Ca²⁺ 激活 K⁺ 通道通过整合 Ca²⁺ 动力学（缓慢累积 + 快速释放）能产生足够尖锐的频率响应。这是数学生理学"否定直觉"的经典案例——数学分析给出的结论（"voltage alone is not enough, you need Ca²⁺ dynamics"）比"细胞生物学家头脑中"先验认为的机制（"电压敏感通道最直接"）更有力。

本章的主要缺陷在于 20.2 节基底膜模型的过度简化——教科书级别的 2D 流体 + 1D 基底膜 ODE + 解析 Fourier 解，距离真实的耳蜗（3D 螺旋、内外毛细胞、organ of Corti、神经纤维）至少 3 层抽象。没有这些简化模型的预测能与临床听力图（audiogram）精确对应——临床听力学的核心是每个频率的阈值，但本章的预测精度仅能解释"行波包络峰值移动"这一粗略特征。

20.4.1 的 毛束负刚度 是本章最具洞察力的概念——它意味着毛束本身是一个非保守系统，类似于 Van der Pol 振子的负阻尼。负刚度 + 正阻尼 → 极限环振荡 → 自发振荡（毛束的 twitches）。这种"生物材料的反常力学行为" 挑战了"生物组织是粘弹性固体"的传统观念，需要主动的、非平衡的模型来解释。

对我自己研究血管生物力学的启发：(1) 基底膜的行波机制与血管壁的脉搏波传播在数学上完全相似——血管壁的弹性梯度（从主动脉到外周动脉）类似于基底膜的刚度梯度，都支持"行波 + 局部最大振幅"模式。脉搏波速度（PWV）的频率依赖性（PWV 随频率增加而增加）与耳蜗波速度随频率变化在形式上同构；(2) Ca²⁺ 激活 K⁺ 通道作为频率调谐机制的发现，提示血管平滑肌的 Ca²⁺-K⁺ 通道（如 BK 通道）可能在血管张力振荡（vasomotion, Ch 19）中起类似作用——BK 通道激活 → K⁺ 外流 → 超极化 → Ca²⁺ 通道关闭 → 钙瞬变衰减 → 收缩结束——形成自持振荡；(3) Gold 假说 + 主动放大机制与血管的肌源性反应（myogenic response, Bayliss 效应）的"主动调节"本质相同——血压升高 → 血管平滑肌收缩 → 阻力增加 → 血流恢复。耳蜗的 prestin 蛋白与血管平滑肌的肌动蛋白-肌球蛋白虽然分子不同，但主动放大的数学框架（带通滤波 + 负阻尼振荡器）相同。

值得指出，本章对前庭系统（半规管、椭圆囊、球囊）几乎没有涉及——这些是平衡觉和本体感觉的关键器官，涉及流体动力学（内淋巴在壶腹嵴上的流动）、惯性（耳石对加速度的反应）、毛细胞转导（与耳蜗毛细胞同源）。前庭系统的数学模型与 Ch 19 的血管振荡、Ch 18 的小肠蠕动、Ch 20 的耳蜗调谐共同构成"生物系统的流-固-电耦合"的精彩案例。

重要参考文献

[X1] von Békésy, G. (1960). Experiments in Hearing. McGraw-Hill. — 20.1 节 Fig. 20.3-20.5 行波包络实验的原始数据；1961 年获 Nobel 奖。 [X2] Helmholtz, H. (1875). On the Sensations of Tone. Longmans, Green. — 20.1 节位置理论的原始论文。 [X3] Lesser, M. B., & Berkley, D. A. (1972). Fluid mechanics of the cochlea. Journal of Fluid Mechanics, 51, 497-512. — 20.2.2-20.2.3 节完整 2D 模型的原始论文。 [X4] Zwislocki, J. (1965). Analysis of some auditory characteristics. In Handbook of Mathematical Biology. — 20.2.4 节 shallow-water 模型的代表。 [X5] Ranke, O. F. (1950). Theory of operation of the cochlea: A contribution to the hydrodynamics of the cochlea. Journal of the Acoustical Society of America, 22, 772-777. — deep-water 模型的代表。 [X6] Siebert, W. M. (1974). Ranke revisited — a simple short-wave cochlear model. In Psychophysics and Physiology of Hearing. Academic Press. — 20.2.4 节直接力学加载的改进。 [X7] Peskin, C. S. (1981). A particular solution of the cochlear problem. In Mathematical Sciences Department. New York University. — 20.2.5 节共形映射精确解。 [X8] Steele, C. R. (1974). Behavior of the basilar membrane with pure-tone excitation. Journal of the Acoustical Society of America, 55, 148-162. — 20.2.5 节 3D 弹性板模型的代表。 [X9] Holmes, M. H. (1980a,b, 1982). A nonlinear oscillator model of the cochlea. Journal of Sound and Vibration. — 非线性基底膜模型的代表。 [X10] Crawford, A. C., & Fettiplace, R. (1981). An electrical tuning mechanism in turtle cochlear hair cells. Journal of Physiology, 312, 377-412. — 20.3.1 节电共振电路模型的原型。 [X11] Ashmore, J. F., & Attwell, D. (1985). Models for electrical tuning in hair cells. Proceedings of the Royal Society B, 226, 325-344. — 20.3.2-20.3.3 节两种 K⁺ 通道模型的比较，证明 Ca²⁺ 激活机制更优。 [X12] Hudspeth, A. J., & Lewis, R. S. (1988a,b). A model for electrical resonance and frequency tuning in saccular hair cells of the bullfrog. Journal of Physiology, 400, 275-297. — 20.3.4 节 3 闭态 + 2 开态 Markov 模型的原型。 [X13] Hudspeth, A. J. (1985). The cellular basis of hearing: the biophysics of hair cells. Science, 230, 745-752. — 毛细胞综述的经典。 [X14] Hudspeth, A. J. (1997). Mechanical amplification of stimuli by hair cells. Current Opinion in Neurobiology, 7, 480-486. — 主动放大综述。 [X15] Hudspeth, A. J. (2005). How the ear's works work. Nature, 431, 893-896. — "在蜂蜜中抑制音叉" 比喻的出处。 [X16] Gold, T. (1948). Hearing. II. The physical basis of the action of the cochlea. Proceedings of the Royal Society B, 135, 492-498. — 耳蜗主动放大假说的原始论文。 [X17] Nobili, R., Mammano, F., & Ashmore, J. F. (1998). How well do we understand the cochlea? Trends in Neurosciences, 21, 159-167. — 耳蜗综述。 [X18] Manley, G. A. (2001). Evidence for an active process and a cochlear amplifier in nonmammals. Journal of Neurophysiology, 86, 541-549. — 主动放大的非哺乳动物证据。 [X19] Martin, P., Mehta, A. D., & Hudspeth, A. J. (2000). Negative hair-bundle stiffness betrays a mechanism for mechanical amplification by the hair cell. Proceedings of the National Academy of Sciences USA, 97, 12026-12031. — 20.4.1 节负刚度假说的原型。 [X20] Ashmore, J. F., Géléoc, G. S. G., & Harbott, L. (2000). Molecular mechanisms of sound amplification in hair cells. Philosophical Transactions of the Royal Society B, 355, 839-846. — prestin 蛋白分子机制的综述。 [X21] Hudspeth, A. J., Choe, Y., Mehta, A. D., & Martin, P. (2000). Putting ion channels to work: Mechanoelectrical transduction, adaptation, and amplification by hair cells. Proceedings of the National Academy of Sciences USA, 97, 11765-11772. — 毛细胞综述。 [X22] Dallos, P., Popper, A. N., & Fay, R. R. (1990, 1996). The Cochlea. Springer. — 耳蜗解剖和生理的综合专著。 [X23] Rhode, W. S. (1984). Cochlear mechanics. Annual Review of Physiology, 46, 231-246. — 耳蜗力学综述。 [X24] Berne, R. M., & Levy, M. N. (1993). Physiology (3rd ed.). Mosby. — 20.1-20.2 节 Fig. 20.1, 20.2 耳蜗解剖图的来源。 [X25] Eatock, R. A. (2000). Adaptation in hair cells. Annual Review of Neuroscience, 23, 285-314. — 毛细胞适应综述。