跳转至

第 20 章:内耳(The Inner Ear)

20 章:内耳概述(Chapter Overview)

本节先简述耳蜗的解剖结构(scala vestibuli、scala media、scala tympani 三腔及基底膜)、声音经镫骨传入并转化为电信号的路径,然后交代本章将讨论的三大主题——基底膜的机械调谐、毛细胞的电学谐振、以及非线性耳蜗放大器——并引出"位置理论"作为贯穿全章的主线。耳蜗(cochlea)位于中耳、鼓膜与外耳的深处(图 20.1A),其内部是一个螺旋结构,内含三个分隔的腔:scala vestibuli(前庭阶)、scala media(中阶)和 scala tympani(鼓阶)(图 20.1B、20.2A)。前庭阶与中阶由前庭膜(Reissner's membrane)分隔;中阶与鼓阶则由螺旋板(spiral lamina)和基底膜(basilar membrane)共同分隔。前庭阶和鼓阶内充满外淋巴液(perilymph),其成分与细胞外液相似;中阶则充满内淋巴液(endolymph),具有高 K+、低 Na+ 的特征——这种离子成分差异对毛细胞的转导机制至关重要。声音经中耳传导后被聚焦于镫骨(stapes),镫骨推动卵圆窗(oval window)——卵圆窗是前庭阶的一个开口。所产生的外淋巴波动沿前庭阶传播,并在基底膜和鼓阶中激发互补的波动。在耳蜗顶端,前庭阶与鼓阶之间由蜗孔(helicotrema)连通,使两腔局部压力得以均衡。由于外淋巴本质上不可压缩,鼓阶也必须有一个类似卵圆窗的开口以满足质量守恒——否则镫骨将无法运动;鼓阶上的这一开口称为圆窗(round window)。卵圆窗处镫骨的内向运动,由圆窗处流体的相应外向运动补偿,二者共同构成基底膜位移的流体力学边界条件。声音向电信号的转导由位于基底膜之上的柯蒂氏器(organ of Corti,图 20.2B)完成:柯蒂氏器中的毛细胞(hair cells)顶部有毛束(hair bundles)伸出,毛束与悬于其上的盖膜(tectorial membrane)相接触。基底膜的波动在这些毛上产生剪切力,进而引起毛细胞膜电位变化;该电信号被传递至神经细胞,再由神经细胞送往大脑。整章的展开围绕两条主线:一是声音的物理振动如何被分解为不同频率成分(位置理论),二是耳蜗为何具有极高的灵敏度(主动放大)。

20.1 频率调谐(Frequency Tuning)

本节是 20.1.1 与后续 20.3 节的引子,先把耳蜗所要完成的任务——识别声波的组成频率——梳理清楚,再把不同动物采用的三类主要调谐机制(毛细胞的机械调谐、基底膜的机械调谐、毛细胞的电学调谐)介绍一遍,并指出本章将集中讨论后两类。耳蜗的任务是识别一个声波的组成频率,从而识别声音本身。不同动物完成这一任务的方式可分为三类主要分组:毛细胞的机械调谐(mechanical tuning of the hair cells)、基底膜的机械调谐(mechanical tuning of the basilar membrane)、以及毛细胞的电学调谐(electrical tuning of the hair cells),参 Hudspeth(1985)和 Eatock(2000)的综述。最早的频率调谐理论之一是 Helmholtz(1875)提出的,他假设耳由一组尖锐调谐的元件(sharply tuned elements)组成,每个元件与某一特定频率共振——相当于把耳类比为钢琴弦阵列,每根弦只对一个基频响应。尽管这远非完整的故事,但共振在调谐中确实扮演了重要角色:例如,在许多蜥蜴中,毛细胞上毛束的长度从基部到顶端系统性地增加,正如更长的弦产生更低音调的音符那样,更长的毛束优先响应较低频率的输入,而较短的毛束调谐到较高频率。因此输入频率可由产生最大刺激的位置决定——这是 Helmholtz 位置理论在毛细胞尺度上的具体体现。在哺乳动物中,基底膜本身充当频率分析器,这将在 20.1.1 节讨论。第三种调谐机制源于毛细胞膜上离子通道的性质:每个毛细胞是一个电学谐振器(electrical resonator),具有带通频率响应(band-pass frequency response);产生最大响应的输入频率是毛细胞生物物理性质的函数,而沿耳蜗长度方向这些性质的系统性变化使耳蜗能够基于最大响应的位置区分不同频率。本章不再进一步讨论第一种(蜥蜴式毛束共振)调谐机制,而是集中讨论后两种——基底膜的机械调谐(20.1.1、20.2 节)和毛细胞膜电位的电学谐振(20.3 节)。这两类机制在动物界有清楚的物种差异:低等脊椎动物(龟、牛蛙)以电学谐振为主,哺乳动物以基底膜的机械调谐为主,但即便在哺乳动物中,基底膜的尖锐调谐也离不开毛细胞的主动放大(20.4 节)。

20.1.1 耳蜗宏观力学(Cochlear Macromechanics)

本节先描述哺乳动物中镫骨振动在基底膜上产生的"先增后减"的行波包络、由 von Békésy 测得的频率响应曲线,以及由此归纳出的"位置理论"——再指出基底膜其实不是真正处于张力的膜,而是靠弯曲刚度抵抗形变的结构,最后概述早期和后续两类基底膜模型的发展谱系。在哺乳动物中,镫骨的振动在基底膜上激起一种具有特定形状的波:波幅包络先增大后减小,包络峰值位置依赖于刺激频率(von Békésy, 1960),如图 20.3 所示。图 20.3 中实线表示连续时刻(按时间递增标记为 1、2、3、4)的基底膜偏转;虚线是膜波的包络,包络随时间保持不变——也就是说,虽然瞬时波形在传播,但包络本身是驻定的(stationary)。波速在沿膜传播过程中减小,导致相位持续减小,表观频率则增大。低频刺激的波包络峰值更靠近耳蜗顶端(即蜗孔附近),随着刺激频率增加,包络峰值向耳蜗基部移动(图 20.4)。图 20.4 给出四个不同频率下基底膜波包的相位和振幅随距镫骨距离的变化。包络的幅度是离镫骨距离 x 与刺激频率 ω 的二维函数;图 20.4 中的曲线是固定频率的横截面。另一种呈现方式是取固定距离的横截面,从而给出固定距离下包络幅度作为频率的函数——即基底膜在该固定距离的频率响应。von Békésy 测得的频率响应如图 20.5 所示,可以看到基底膜的每一部分(图中对应不同距离:13、17、20、24、28、31 mm)对某一频率响应最大,且随频率增加,最大响应位置向镫骨方向移动(最近端 13 mm 对高频,最远端 31 mm 对低频)。通过这种方式,耳蜗由基底膜最大幅度位置(即所谓的听觉位置理论,place theory)确定输入信号的频率——这一定律直接对应于 Helmholtz 共振阵列的现代解读:每段基底膜相当于一个针对特定频率尖锐调谐的元件。虽然 von Békésy 的许多结果此后被取代(其实验是在尸体上完成的——虽然有些残忍,但现已知道活人基底膜和毛细胞的性质与之不同,活的基底膜表现出远比尸体尖锐的调谐),他的实验结果及与之相关的位置理论仍是 20 世纪最重要的耳蜗研究之一,对理论研究产生了巨大影响。"基底膜"这一名称具有误导性,因为它并不是真正的膜——这一事实由以下观察证明:将其切开时边缘不会收缩,因此它并不处于张力之下;抵抗运动的力量来自弯曲弹性(bending elasticity)而非面内张力。基底膜的刚度(stiffness)从基部到顶端指数式减小,长度常数约为 7 mm(即按 e^{−x/7mm} 衰减);虽然耳蜗的宽度从基部到顶端减小,但基底膜的宽度沿该方向反而增加,这种几何上的反向变化与刚度变化共同决定了基底膜的振动特性。基底膜上波的模型可按其膜和流体所采用的方程类型加以区分。早期 Ranke(1950)和 Zwislocki(1965)的模型假设外淋巴不可压缩且无粘性,并把基底膜建模为有阻尼的受迫谐振子,沿膜长方向没有弹性耦合;Ranke 使用深水波理论,而 Zwislocki 使用浅水波理论,由此引发了大量关于哪种方法更好的争论。这些模型被许多作者发展,其中最知名的是 Peterson and Bogert(1950)、Fletcher(1951)、Lesser and Berkley(1972)以及 Siebert(1974)。其后 Steele(1974)、Inselberg and Chadwick(1976)、Chadwick et al.(1976)、Chadwick(1980)、Holmes(1980a,b, 1982)等的模型则将基底膜更细致地表示为弹性板(elastic plate),并纳入了流体的粘性和板的几何。实验与理论结果的综述见 Dallos et al.(1990, 1996)。本章概述一些早期和较简单的模型,因为它们能优雅地展示基底膜与外淋巴如何相互作用产生 von Békésy 观测到的波。

20.2.1 不可压缩流体的运动方程(Equations of Motion for an Incompressible Fluid)

本节从流体力学基本方程出发——质量守恒、动量守恒——推导出小振幅下不可压缩无粘流体的线性化方程组,并写出存在速度势 φ 时的特殊形式(不可压缩条件下的 Laplace 方程与 Bernoulli 型压力关系),作为后续基底膜耦合方程的流体侧基础。耳蜗中围绕基底膜的流体不可压缩且假设无粘;这种流体的运动方程是众所周知的,并在许多地方(如 Batchelor, 1967)有推导。设 u = (u1, u2, u3) 为流体速度,p 为压力,ρ 为流体密度。固定体积 V 中流体的质量只能因跨越 V 边界的通量而改变,因此有

\[ \frac{d}{dt} \int_V \rho\, dV = -\int_S \rho (u \cdot n)\, dS \quad (20.1) \]

其中 S 是 V 的表面,n = (n1, n2, n3) 为 V 的向外单位法向。(20.1) 的物理含义是:体积 V 内流体质量的变化率等于穿过边界 S 流出的质量通量的负值——即质量守恒的积分形式。类似地,固定域 V 中流体的动量只能因所施外力或跨越 V 边界的动量通量而改变。因此(对无粘流体,即忽略粘性切应力的情形)动量守恒意味着

\[ \frac{d}{dt} \int_V \rho u_i\, dV = -\int_S [(u \cdot n)\rho u_i + p n_i]\, dS \quad (20.2) \]

面积分中第一项 (u · n) ρ u_i 是对流项,表示动量随流体外流而带走的速率;第二项 p n_i 是压力对边界的法向应力——正是无粘条件下动量通量的来源。利用散度定理把面积分转化为体积分,得到

\[ \int_V \left[ \frac{\partial \rho u_i}{\partial t} + \rho \nabla \cdot (u_i u) + \frac{\partial p}{\partial x_i} \right] dV = 0 \quad (20.3) \]
\[ \int_V \left[ \frac{\partial \rho}{\partial t} + \nabla \cdot u \right] dV = 0 \quad (20.4) \]

由于 V 是任意的,且假设密度 ρ 为常数(不可压缩流体),被积函数必须恒为零,于是

\[ \rho \frac{\partial u}{\partial t} + \rho (\nabla \cdot u) u + \nabla p = 0 \quad (20.5) \]
\[ \nabla \cdot u = 0 \quad (20.6) \]

其中 (20.5) 是 Euler 方程的非线性形式,(20.6) 是不可压缩条件。当流体运动幅度很小(耳蜗中的情形应当满足)时,速度的二阶项 ρ(∇ · u) u 相对于其他项可忽略,于是方程被线性化为

\[ \rho \frac{\partial u}{\partial t} + \nabla p = 0 \quad (20.7) \]
\[ \nabla \cdot u = 0 \quad (20.8) \]

(20.7) 是小振幅下不可压缩无粘流体的线性动量方程。一个重要的特例是 u = ∇φ(无旋流),即速度场可由一个标量势 φ 表示;此时 (20.7) 和 (20.8) 化为

\[ \rho \frac{\partial \phi}{\partial t} + p = 0 \quad (20.9) \]
\[ \nabla^2 \phi = 0 \quad (20.10) \]

(20.10) 表明速度势满足 Laplace 方程——这是不可压缩无旋流动在流体区域内的核心方程;(20.9) 是 Bernoulli 关系的时间依赖形式,给出压力与势的时间导数之间的简单代数关系,其中 p 已归一化使积分常数为零(这相当于在远场或对称位置定义零压)。这一对方程正是 20.2 节基底膜–流体耦合模型中流体一侧的出发点。

20.2.2 基底膜作为谐振子(The Basilar Membrane as a Harmonic Oscillator)

本节描述 Lesser–Berkley 模型的几何抽象(两个矩形流体腔由基底膜分隔),把每点基底膜视为无侧向耦合的有阻尼谐振子,写出其运动方程和流体边界条件,并在"小位移"近似下线性化为驱动项在 y=0 处取值的简化问题。耳蜗最简单的模型之一是把 (20.9) 和 (20.10) 与有阻尼受迫谐振子方程相结合。最清晰的呈现来自 Lesser and Berkley(1972)。在其模型中,耳蜗被假设为图 20.6 所示的构型:耳蜗被简化为两个矩形流体腔,上、下分别填充外淋巴,由基底膜水平分隔;上腔对应 scala vestibuli,下腔对应 scala tympani;为简化起见,scala media(图 20.1 中所示)被省略。下标 1 和 2 分别表示上、下两个腔中的量,因此有两套 (20.9) 和 (20.10):

\[ \rho \frac{\partial \phi_1}{\partial t} + p_1 = \rho \frac{\partial \phi_2}{\partial t} + p_2 = 0 \quad (20.11) \]
\[ \nabla^2 \phi_1 = \nabla^2 \phi_2 = 0 \quad (20.12) \]

其中压力只能确定到一个任意常数(两腔压力以同样任意常数相差)。基底膜的每一点被建模为一个简单的有阻尼谐振子,其质量、阻尼和刚度沿膜长方向变化;由于没有直接侧向耦合,膜任一部分的运动被假设独立于邻近部分。基底膜的位置 y = η(x, t) 由下式给出:

\[ m(x) \frac{\partial^2 \eta}{\partial t^2} + r(x) \frac{\partial \eta}{\partial t} + k(x) \eta = p_2(x, \eta(x,t), t) - p_1(x, \eta(x,t), t) \quad (20.13) \]

其中 m(x) 是基底膜单位面积质量,r(x) 是阻尼系数,k(x) 是单位面积刚度(胡克常数);等式左端是基底膜单位面积上的惯性力 + 阻尼力 + 弹性力,右端是上、下腔压力差提供的驱动力。边界条件要求边界处流体速度与边界速度匹配——即流体不穿透基底膜:

\[ \frac{\partial \eta}{\partial t} = \frac{\partial \phi_1}{\partial y} = \frac{\partial \phi_2}{\partial y}, \quad y = \eta(x,t), \quad 0 < x < L \quad (20.14) \]

再假设顶部没有竖向运动(上壁刚硬):

\[ \frac{\partial \phi_1}{\partial y} = 0, \quad y = l, \quad 0 < x < L \quad (20.15) \]

按原始形式,这是一个非常困难的自由边界问题:膜的位移 η 本身出现在边界条件里,使问题非线性且几何依赖。若假设位移很小,则有显著简化:主导阶的线性问题中,驱动力取 y = 0 处的压力差而非 y = η 处(即把自由边界"冻结"在未扰动位置 y = 0 上):

\[ m(x) \frac{\partial^2 \eta}{\partial t^2} + r(x) \frac{\partial \eta}{\partial t} + k(x) \eta = p_2(x, 0, t) - p_1(x, 0, t) \quad (20.16) \]

边界条件也相应简化为

\[ \frac{\partial \eta}{\partial t} = \frac{\partial \phi_1}{\partial y} = \frac{\partial \phi_2}{\partial y}, \quad y = 0, \quad 0 < x < L \quad (20.17) \]

外部激励的施加有若干方式。Lesser 和 Berkley 的方法是假设与卵圆窗接触的镫骨运动决定卵圆窗的位置。由于 ∂φ/∂x 是流体速度的 x 分量,x = 0 处的边界条件为

\[ \frac{\partial \phi_1}{\partial x} = \frac{\partial F(y, t)}{\partial t}, \quad 0 < y < l \quad (20.18) \]

其中 F(y, t) 是卵圆窗的指定水平位移——也就是说,整个卵圆窗上的水平速度分布由 ∂F/∂t 给出。再假设远端 x = L 处没有水平运动(即耳蜗另一端流体被封闭):

\[ \frac{\partial \phi_1}{\partial x} = 0, \quad 0 < y < l \quad (20.19) \]

这一组方程 (20.11)–(20.19) 构成了 Lesser–Berkley 模型的完整方程组:它把流体 Laplace 方程、基底膜谐振子方程以及镫骨驱动边界条件耦合起来。

20.2.3 解析解(An Analytical Solution)

本节先利用问题关于 y 的奇对称性把求解域减半,然后对单频激励取 Fourier 形式,得到基底膜上的频域代数方程(阻抗 Z = iωm + r + k/(iω)),再对模型方程无量纲化后用 Fourier 级数展开求解,得到一组确定系数 Am 的方程 Amαm = fm,最终展示出与 von Békésy 观测定性一致的"先增后减"行波包络。鉴于问题关于 y 的固有对称性,求解时寻找 y 的奇函数解(Exercise 1),因此只考虑上区域问题,并去掉下标 1。当输入为单频 F(y, t) = F̂(y) e^{iωt} 时,φ(x, y, t) 具有 φ̂(x, y; ω) e^{iωt} 的形式;其他变量类似。对所有变量采用这种形式,得到方程

\[ \nabla^2 \hat{\phi} = 0, \qquad \hat{p} + i \omega \rho \hat{\phi} = 0 \quad (20.20) \]
\[ \frac{\partial \hat{\phi}}{\partial y} = i \omega \hat{\eta}, \qquad i \omega \hat{\eta} Z = -\hat{p}, \quad \text{on } y = 0 \quad (20.21) \]
\[ \frac{\partial \hat{\phi}}{\partial x} = U_0, \quad \text{on } x = 0 \quad (20.22) \]
\[ \frac{\partial \hat{\phi}}{\partial x} = 0, \quad \text{on } x = L \quad (20.23) \]
\[ \frac{\partial \hat{\phi}}{\partial y} = 0, \quad \text{on } y = l \quad (20.24) \]

其中 Z = iωm + r + k/(iω),U_0 = iωF̂。通过在频域中求解,已将基底膜上的微分方程化为代数方程:膜上的偏微分方程原本是时域的二阶 ODE (20.16),但在 Fourier 形式下它变成 iωZ 与 η̂、p̂ 之间的代数关系 (20.21)。iωZ 是有阻尼谐振子的频率响应,阻抗 Z 是 x 的函数(因为 m, r, k 都依赖 x)。注意在 (20.21) 中假设压力是 y 的奇函数——这与 y 方向上基底膜上下的对称性一致(上、下腔压力大小相等、符号相反)。最后,对模型方程做无量纲化:用 L 缩放 x 和 y,用 iωρL 缩放 Z,用 U_0 L 缩放 φ̂;整理后去掉帽子,得

\[ \nabla^2 \phi = 0 \quad (20.25) \]
\[ \frac{\partial \phi}{\partial y} = \frac{2\phi}{Z}, \quad \text{on } y = 0 \quad (20.26) \]
\[ \frac{\partial \phi}{\partial x} = 1, \quad \text{on } x = 0 \quad (20.27) \]
\[ \frac{\partial \phi}{\partial x} = 0, \quad \text{on } x = 1 \quad (20.28) \]
\[ \frac{\partial \phi}{\partial y} = 0, \quad \text{on } y = \sigma \quad (20.29) \]

其中 σ = l/L 是无量纲化的耳蜗高度比。该问题的解析解可通过标准 Fourier 级数求得。寻找如下形式的解

\[ \phi = x \left( 1 - \frac{x}{2} \right) \frac{\sigma - y}{\sigma} + \sum_{n=0}^\infty A_n \cosh[n\pi(\sigma - y)] \cos(n\pi x) \quad (20.30) \]

其中 An 是待定常数;这一构造使得右端第一项自动满足齐次边界条件 (20.27)–(20.29)——可以验证 (1 − x/2)(σ − y)/σ 同时满足 ∂/∂x 在 x = 0 处等于 1、在 x = 1 处等于 0,以及 ∂/∂y 在 y = σ 处等于 0。第二项是 Laplace 方程在矩形域上满足 Neumann 边界条件的标准 Fourier 余弦-双曲余弦展开(其中 cos(nπx) 满足 x = 0、1 处的 Neumann 条件,cosh 在 y = σ 处满足 Neumann 条件)。由于 φ 满足除 (20.26) 之外的所有边界条件,我们用 (20.26) 确定未知系数 An,得到

\[ \sigma + \sum_{n=0}^\infty n\pi A_n \sinh(n\pi\sigma) \cos(n\pi x) - \frac{2}{Z}\left[ x \left(1 - \frac{x}{2}\right) + \sum_{n=0}^\infty A_n \cosh(n\pi\sigma) \cos(n\pi x) \right] = 0 \quad (20.31) \]

两边乘以 cos(mπx) 并从 0 到 1 积分(利用 cos(nπx) 在 [0,1] 上的正交性),得

\[ A_m \alpha_m = f_m \quad (20.32) \]

其中

\[ \alpha_m = \frac{1}{Z} \cosh(m\pi\sigma) - \frac{1}{2 m\pi} \sinh(m\pi\sigma) \quad (20.33) \]
\[ f_m = \sigma \delta_{m0} - \int_0^1 \frac{x(2-x) \cos(m\pi x)}{Z} dx = -\frac{2}{m^2 \pi^2} \quad (20.34) \]

系数 An 现在可以显式求出(Am = fm/αm),并代回 (20.30)。典型结果如图 20.7 所示(参数为 m = 0.05 g/cm²,k = 10^7 e^{−1.5x} dynes/cm³,r = 3000 e^{−1.5x} dynes sec/cm³,ω = 1000/sec,外淋巴密度取 1 g/cm³):波包络与 von Békésy 结果(图 20.3)定性一致,包络峰值随频率增加向耳蜗基部移动。这一 Fourier 解的存在本身就证明了"行波 + 包络"这一图像可由流体–谐振子耦合的线性模型产生。

20.2.4 长波与短波模型(Long-Wave and Short-Wave Models)

本节先介绍两种经典的解析近似——Zwislocki 的长波(浅水)和 Ranke 的短波(深水)近似,及其对耳蜗深度的不同假设,并指出实验倾向于支持短波理论;接着引入 Siebert(1974)的广义模型(在基底膜两端加直接力源),推导无量纲化的 p(x) 主控方程 (20.44);然后在浅水假设下利用 tanh 的 Taylor 展开将问题化为单个二阶常微分方程,再对小参数 ϵ 做多尺度分析,得到振荡解的包络并求出峰值位置 xp;最后在深水假设下利用 tanh ≈ sign(n) 的高频近似,分离正、负空间频率得到 p± 的一阶线性方程组及其解析形式 (20.75),并在生理参数下给出膜位移的简洁表达式 (20.79) 以及特例下的峰值位置 (20.84)。Fourier 解表明 Lesser–Berkley 模型的行为定性正确,但若能对基底膜行为给出更好的解析理解则更理想。模型方程有两种经典近似允许进一步解析研究:Zwislocki 等人研究的长波近似假设波长大于耳蜗深度,而 Ranke 的短波近似则假设相反——耳蜗深度有效无穷。实验表明耳蜗深度对耳蜗波几乎没有影响,从而支持短波理论:即使将耳蜗的一侧完全移除,对波也几乎没有影响。然而两个模型都不能完整描述耳蜗行为(Zwislocki, 1953)。短波和长波模型都可作为 20.2.2 节模型方程的近似特例推得。为此,采用 Siebert(1974)模型(图 20.8),它在前述模型基础上只改动一处:假设在基底膜两端存在直接机械激励。把膜运动方程 (20.16) 修改为包含这种直接激励的形式:

\[ m(x)\frac{\partial^2\eta}{\partial t^2} + r(x)\frac{\partial\eta}{\partial t} + k(x)\eta = p_2(x,0,t) - p_1(x,0,t) + F_0(t)\delta(x) - F_L(t)\delta(x-L) \quad (20.35) \]

与之前一样,假设激励为单频 F_0(t) = F_0 e^{iωt},F_L(t) = F_L e^{iωt}。由 (20.17)、(20.20) 和 (20.35) 得

\[ \nabla^2 p(x,y) = 0 \quad (20.36) \]
\[ -i\omega \eta = \frac{1}{i\omega \rho} \frac{\partial p(x,0)}{\partial y} \quad (20.37) \]
\[ Y(x) p(x,0) = -i\omega \eta(x) + \eta_0 \delta(x) - \eta_L \delta(x-L) \quad (20.38) \]

其中 Y = 2/Z,η_0 = F_0/Z(0),η_L = F_L/Z(L),并已略去 Fourier 变换对应的帽子。由于 Laplace 方程在矩形域上可分离,用 Fourier 级数写出 (20.36) 的解:

\[ p(x,y) = \sum_{n=-\infty}^{\infty} \alpha_n \frac{\cosh[2\pi n(y-l)/L]}{\cosh(2\pi n l/L)} e^{2\pi i n x / L} \quad (20.39) \]

其中用了 y = l 上 ∂p/∂y = 0 的边界条件。由此

\[ -i\omega \eta = \frac{1}{i\omega \rho} \frac{\partial p(x,0)}{\partial y} = -\frac{1}{i\omega \rho} \sum_{n=-\infty}^{\infty} \alpha_n \frac{2\pi n}{L} \tanh(2\pi n l/L) e^{2\pi i n x / L} \quad (20.40) \]

由于我们只关心基底膜行为,从现在起把注意力限于 p(x, 0),记为 p(x)。

浅水近似。 在浅水近似中,假设基底膜上的波长远大于耳蜗深度。因此假设对某整数 N 使得 Nl/L ≪ 1 时,α_n = 0(对所有 n > N)。由于求和只包含 nl/L ≪ 1 的项,每一项中 tanh(2πnl/L) 可由 Taylor 展开的最低阶项近似,即 tanh(2πnl/L) ≈ 2πnl/L,于是 (20.40) 成为

\[ -i\omega \eta \approx -\frac{l}{i\omega \rho} \sum_{n=-\infty}^{\infty} \alpha_n \left(\frac{2\pi n}{L}\right)^2 e^{2\pi i n x / L} \quad (20.41) \]

而对 (20.39) 作类似论证可得

\[ p(x) = \sum_{n=-\infty}^{\infty} \alpha_n e^{2\pi i n x / L} \quad (20.42) \]

结合 (20.41),有

\[ -i\omega \eta \approx \frac{l}{i\omega \rho} \frac{d^2 p}{dx^2} \quad (20.43) \]

结合 (20.37) 和 (20.38),得到 p(x) 的单个方程

\[ Y(x) p(x) = \frac{l}{i\omega \rho} \frac{d^2 p(x)}{dx^2} + \eta_0 \delta_0(x) - \eta_L \delta_L(x) \quad (20.44) \]

要把该方程中的 δ 函数转化为边界条件,从 x = −ϵ 到 x = +ϵ 积分并令 ϵ → 0,可得

\[ \frac{dp}{dx} = -i\omega \rho \eta_0, \quad \text{at } x = 0 \quad (20.45) \]

其中已假设 dp/dx 在 x = 0−(耳蜗边界之外)为零;类似地,从 x = L − ϵ 到 x = L + ϵ 积分并令 ϵ → 0,得

\[ \frac{dp}{dx} = i\omega \rho \eta_L, \quad \text{at } x = L \quad (20.46) \]

注意当 η_0 = 1 且 η_L = 0 时(且由 p + iωρφ = 0),这些边界条件与 Lesser–Berkley 模型中的 (20.27) 和 (20.28) 一致。该方程的解析利用了 Y(x) 是缓慢变化函数的事实:从量纲上看,iωρY(x)/l 具有长度 −2 的单位,它确定了 p(x) 空间振荡的尺度(波长);另一方面,Y(x) 以长度常数 λ−1 ≈ 0.7 cm 指数变化。如果这两个长度常数的比值很小,就称 Y(x) 是缓慢变化函数。此外,可对空间作 x = z/q 重标度,将 (20.44) 化为无量纲形式

\[ \frac{d^2 p}{dz^2} + g^2(\epsilon z) p(z) = 0 \quad (20.47) \]

其中 ϵ 是小正数,g^2(ϵz) = −iωρY(z/q) / (l q^2) 振幅量级为 1 且在 z 上缓慢变化。q 是任意长度尺度,按使 g^2 振幅量级为 1 选取;按假设 λ/q ≪ 1。作为一个具体例子,假设 m = 0、k(x) = k_0 e^{−λx}、r(x) = r_0 e^{−λx},则

\[ \frac{i\omega \rho Y(x)}{l} = -\frac{2\omega^2 \rho}{l k_0} e^{\lambda x} \left(1 + \frac{i\omega r_0}{k_0}\right) \quad (20.48) \]

设 q^2 = 2ω^2 ρ/(l k_0),并定义

\[ \epsilon = \frac{\lambda}{2q} \quad (20.49) \]

如果参数使得 ϵ ≪ 1,则有缓慢变化的振荡。这类问题在振荡理论中是熟知的,可用多尺度分析求解(Kevorkian and Cole, 1996; Keener, 1998)。我们希望求 (20.47) 的近似解。若 g 为常数(ϵ = 0),则 (20.47) 的解简单地为

\[ p(z) = A e^{i g z} + B e^{-i g z} \quad (20.50) \]

但由于 g 假设缓慢变化,我们预期上述基本解在局部(而非全局)是合理的近似。为得到更长程有效的解,引入两个尺度:慢变量 σ = ϵz,以及满足 dτ/dz = f(ϵz) 的快变量 τ,f 是待定函数。于是导数 d/dz 必须替换为偏导数

\[ \frac{d}{dz} = f(\sigma) \frac{\partial}{\partial \tau} + \epsilon \frac{\partial}{\partial \sigma} \quad (20.51) \]

用这两个变量,原常微分方程 (20.47) 变为偏微分方程

\[ f^2(\sigma) \frac{\partial^2 p}{\partial \tau^2} + \epsilon f(\sigma) \frac{\partial^2 p}{\partial \sigma \partial \tau} + \epsilon \frac{\partial}{\partial \sigma}\left[f(\sigma) \frac{\partial p}{\partial \tau}\right] + \epsilon^2 \frac{\partial}{\partial \sigma}\left[f(\sigma) \frac{\partial p}{\partial \sigma}\right] + g^2(\sigma) p = 0 \quad (20.52) \]

显然 f 的最佳选择是 f = g,因为此时关于 ϵ 的最低阶方程为

\[ \frac{\partial^2 p}{\partial \tau^2} + p = 0 \quad (20.53) \]

其通解为

\[ P_0 = A(\sigma) e^{i\tau} + B(\sigma) e^{-i\tau} \quad (20.54) \]

注意 A 和 B 是慢变量 σ 的函数,因为 (20.53) 是偏微分方程。为确定 A 和 B,设 p = P_0 + ϵ P_1 + O(ϵ^2),收集 ϵ 的同次幂,得 P_1 的方程

\[ \left[g^2(\sigma) \frac{\partial^2 P_1}{\partial \tau^2} + P_1\right] = -g(\sigma) \frac{\partial^2 P_0}{\partial \sigma \partial \tau} - \frac{\partial}{\partial \sigma}\left[g(\sigma) \frac{\partial P_0}{\partial \tau}\right] \quad (20.55) \]

由于希望 P_1 是 τ 的周期函数,要求 P_1 "无长期项",即 (20.55) 右端不含与 e^{iτ} 或 e^{−iτ} 成正比的项。于是

\[ \frac{\partial}{\partial \sigma}(g A^2) = 0, \qquad \frac{\partial}{\partial \sigma}(g B^2) = 0 \quad (20.56) \]

\[ A(\sigma) = \frac{A_0}{\sqrt{g(\sigma)}}, \qquad B(\sigma) = \frac{B_0}{\sqrt{g(\sigma)}} \quad (20.57) \]

由此得到关于 ϵ 的最低阶解

\[ p = \frac{1}{\sqrt{g(\sigma)}} \left(A_0 e^{i G(z)} + B_0 e^{-i G(z)}\right) \quad (20.58) \]

其中 G(z) = ∫_0^z g(ϵz) dz。用有量纲变量表示为

\[ p(x) = \phi^{-1/2} \left[ A_1 \exp\left(i \int_0^x \phi(s)\, ds\right) + B_1 \exp\left(-i \int_0^x \phi(s)\, ds\right) \right] \quad (20.59) \]

其中

\[ \phi(x) = \sqrt{-\frac{i\omega \rho Y(x)}{l}} \quad (20.60) \]

常数 A_1 和 B_1 由边界条件 (20.45) 和 (20.46) 确定,再由恒等式 iωη(x) = −Y(x) p(x) 求膜位移。该解的关键特征是它是振荡的且具有包络,包络的最大幅度和位置由频率 ω 决定。在 m = 0、k(x) = k_0 e^{−λx}、r(x) = r_0 e^{−λx} 的特例中,φ(x) = α e^{λx/2},其中 α^2 = 2ω^2 ρ / [l(k_0 + iωr_0)]。令 α = α_r + iα_i 并假设 α_i λ ≫ 1(低频时该假设不成立),则取 η_L = 0,得

\[ \eta(x) = -\frac{1}{i\omega Y(x)} p(x) \approx \hat{A} \exp\left[\frac{3\lambda x}{4} - \frac{2\alpha_i}{\lambda} e^{\lambda x/2} + \frac{2i \alpha_r}{\lambda} e^{\lambda x/2}\right] \quad (20.61) \]

这是一个振幅随相位和幅度都指数增长的振荡:

\[ |\eta| \approx |\hat{A}| \exp\left[\frac{3\lambda x}{4} - \frac{2\alpha_i}{\lambda} e^{\lambda x/2}\right] \quad (20.62) \]

其包络的最大值出现在

\[ x_p = -\frac{2}{\lambda} \ln\left(\frac{4\alpha_i}{3\lambda}\right) \quad (20.63) \]

最大位置依赖于频率,因为当 ω 充分大时,

\[ \alpha_i \approx \sqrt{\frac{\rho \omega}{l r_0}} \quad (20.64) \]

因此对大 ω 有

\[ x_p = -\frac{1}{\lambda} \ln\left(\frac{16 \rho \omega}{9 l \lambda^2 r_0}\right) \quad (20.65) \]

深水近似。 第二种方法——短波或深水近似——假设基底膜波的波长相对于耳蜗深度很短。此时 p(x) 的 Fourier 展开只包含高频,因此存在某个大整数 N 使得 |n| < N 时 α_n = 0;而当 |n| > N 且 l ≫ L 时,tanh(2πnl/L) ≈ sign(n)。于是 (20.40) 变为

\[ -i\omega \eta \approx -\frac{1}{i\omega \rho} \sum_{n=-\infty}^{\infty} \alpha_n \frac{2\pi}{L} |n| e^{2\pi i n x / L} \quad (20.66) \]

现在把求和拆成两部分,定义两个函数

\[ p_+(x) = \sum_{n=0}^{\infty} \alpha_n e^{2\pi i n x / L} \quad (20.67) \]
\[ p_-(x) = \sum_{n=-\infty}^{-1} \alpha_n e^{2\pi i n x / L} \quad (20.68) \]

则 (20.38) 变为

\[ Y p = Y(p_+ + p_-) \approx \frac{1}{\omega \rho}\left(\frac{dp_+}{dx} - \frac{dp_-}{dx}\right) + \eta_0 \delta_0(x) - \eta_L \delta_L(x) \quad (20.69) \]

取其为 p 的控制方程。可通过跨越 x = 0 和 x = L 边界积分并假设耳蜗外 p = 0 来消除 δ 函数影响。这给出

\[ \frac{1}{\omega \rho}[p_+(0) - p_-(0)] = \eta_0 \quad (20.70) \]
\[ \frac{1}{\omega \rho}[p_+(L) - p_-(L)] = \eta_L \quad (20.71) \]

虽然 p_+ 只是正(空间)频率的线性组合,但 Yp_+ 并非如此。然而若假设 Y 是 x 的缓慢变化函数,则 Y 的 Fourier 级数(系数 b_k)满足:对某个相对 p_+ 的主导频率而言较小的某个 k_0,当 |k| > k_0 时 b_k ≈ 0。于是

\[ Y p_+ = \sum_{k=-\infty}^{\infty} c_k e^{2\pi i k x / L} \quad (20.72) \]

其中 c_k = ∑{j=0}^\infty α_j b。如果 p_+ 和 Y 的主导频率按上述方式分离,则对 k ≤ 0 时 c_k 较小,因此可用仅含正频率的 Fourier 级数近似 Yp_+。对 Yp_− 可作类似论证。在这些近似下,(20.69) 分离为关于正、负频率的两对方程:

\[ Y(x) p_\pm(x) \approx \pm \frac{1}{\omega \rho} \frac{d p_\pm}{dx}, \quad 0 < x < L \quad (20.73) \]

这些一阶线性方程可直接积分,得

\[ p_\pm(x) = A_\pm \exp\left(\pm \omega \rho \int_0^x \frac{1}{|Y(\zeta)|} d\zeta\right) \quad (20.74) \]

其中 A_\pm 为常数,所以

\[ p = A_+ \exp\left(\omega \rho \int_0^x \frac{1}{|Y(\zeta)|} d\zeta\right) + A_- \exp\left(-\omega \rho \int_0^x \frac{1}{|Y(\zeta)|} d\zeta\right) \quad (20.75) \]

用 x = 0 和 x = L 处的边界条件确定常数 A_\pm。由 (20.70) 和 (20.71) 知这些常数须满足

\[ A_+ - A_- = \omega \rho \eta_0, \qquad \gamma A_+ - \frac{1}{\gamma} A_- = \omega \rho \eta_L \quad (20.76) \]

其中 γ = exp[ωρ∫_0^L (1/|Y(ζ)|) dζ],由此

\[ A_+ = \frac{\omega \rho}{\gamma^2 - 1}[\gamma \eta_L - \eta_0] \quad (20.77) \]
\[ A_- = \frac{\omega \rho}{\gamma^2 - 1}[\gamma \eta_L - \gamma^2 \eta_0] \quad (20.78) \]

对生理上的 Y 值,除最低频率外 |γ| ≫ 1;例如对图 20.7 中的参数值,ω = 800 时 |γ| = 40,ω = 1500 时 |γ| = 109。由于从物理角度看 η_0 和 η_L 并不大,于是 |A_-| ≈ −ωρη_0 ≫ |A_+|。最后由 (20.38) 和 (20.75),膜位移由下式给出:

\[ \eta = -i\rho \eta_0 Y(x) \exp\left(-\omega \rho \int_0^x \frac{1}{|Y(\zeta)|} d\zeta\right) \quad (20.79) \]

可把 η 的幅度作为 x 的函数画出,给出基底膜上波的包络。固定 x 并把 |η| 画为 ω 的函数即得频率响应。典型结果如图 20.9 所示:与数据的定性吻合良好,频率增加时波包络峰值向镫骨方向移动。在 m = 0、k(x) = k_0 e^{−λx}、r(x) = r_0 e^{−λx} 的特例中,可将波形 (20.79) 计算为

\[ \eta = 2 \eta_0 \xi \exp[\lambda x + \beta(1 - e^{\lambda x})] \quad (20.80) \]

其中

\[ \beta = \frac{2\omega^3 \rho r_0}{\lambda(k_0^2 + \omega^2 r_0^2)} + i \frac{2\omega^2 \rho k_0}{\lambda(k_0^2 + \omega^2 r_0^2)} = \beta_r + i \beta_i \quad (20.81) \]
\[ \xi = \frac{\omega \rho}{i \omega \rho + k_0} \quad (20.82) \]

这里再次看到一个振荡波形,其包络幅度为

\[ |\eta| = 2 \eta_0 |\xi| \exp[\lambda x + \beta_r (1 - e^{\lambda x})] \quad (20.83) \]

最大值出现在

\[ x_p = -\frac{1}{\lambda} \ln\left(\frac{2\omega^3 \rho r_0}{\lambda(k_0^2 + \omega^2 r_0^2)}\right) \quad (20.84) \]

由该式可知 ω 增加时包络峰值向左(即向耳蜗基部)移动。短波模型的主要缺陷是其相位比实验观察增加得更多。Peskin(1976, 1981)也求得了类似的解,他用共形映射和围道积分技术对耳蜗模型的一种特例求出了精确解:在他的模型中,耳蜗膜假设为无限长,r(x) 和 k(x) 选为衰减率为 λ 的指数函数,流体容器高度 λl = π/2。

20.2.5 更复杂的模型(More Complex Models)

本章集中讨论较简单的基底膜模型,这些模型假设耳蜗是二维的、基底膜可用点阻抗函数描述(即膜上每一点都作为有阻尼谐振子,除了通过流体运动间接施加的耦合外,沿膜长方向没有耦合)。尽管基底膜上的波动是听觉过程中的重要组成部分,还涉及许多其他因素(Pickles, 1982; Rhode, 1984; Hudspeth, 1985)。耳蜗响应中的非线性和声发射(acoustic emissions)提示存在主动反馈过程调节波形(20.4 节)。这种反馈可能发生在外毛细胞(outer hair cells)和柯蒂氏器中。简单流体动力学模型不能重现哺乳动物耳蜗中观察到的调谐程度,精确的调谐机制至今仍有争议。许多其他更复杂的模型已被构建(参见例如 Steele, 1974; Steele and Tabor, 1979a,b; Inselberg and Chadwick, 1976; Chadwick et al., 1976; Chadwick, 1980; Holmes, 1980a,b, 1982)。一般而言,这些模型对流体流动使用类似的方程,但将基底膜建模得更细致,包括膜内的空间耦合。得到的膜方程在空间上是四阶的,并在求解中大量使用渐近展开。Steele(1974)做过一项特别详细的研究,他构建了一系列模型——从无限流体中的板,到锥形弹性基底膜,到刚壁耳蜗,再到柔性柯蒂氏拱。最近的模型基于基底膜、柯蒂氏器和外毛细胞的有限元描述,并整合了声学、电学和力学元件(Ramamoorthy et al., 2007),耳蜗的三维螺旋性质则由 Givelberg and Bunn(2003; Givelberg, 2004)建模。

20.3 毛细胞中的电学谐振(Electrical Resonance in Hair Cells)

本节先点出在很多低等脊椎动物中频率分解由毛细胞自身完成而非基底膜上的波,并描述毛细胞顶端的毛束、tip link 以及机械-电转导通道如何把毛束的偏转转化为膜电位变化,再列出毛细胞在阶跃电流下表现出的衰减振荡行为与带通频率响应,最后引出 20.3.1 和 20.3.2 两类模型——Crawford–Fettiplace 和 Ashmore–Attwell 的电路类比模型,以及 Hudspeth–Lewis 的更细致的机制性模型。在许多低等脊椎动物中,频率分解不是由基底膜上的波完成,而是由毛细胞自身完成。龟耳蜗和牛蛙囊(仅举两例)中的毛细胞优先响应某一频率的刺激,这种带通响应由毛细胞膜上的离子通道介导。每个毛细胞顶端是毛束——一组通过纤细纤维(称为 tip link)在尖端彼此相连的静纤毛(stereocilia)。每根静纤毛都是刚性的,尖端受力时绕基部枢转(pivot around its base)而不是弯曲——这种刚性枢转是力学转导的几何基础,力的方向直接决定哪些 tip link 被拉伸。据推测,tip link 像直接连接到离子通道的弹性弹簧:当毛束向一个方向偏转时,tip link 拉开通道(gate open);当毛束向反方向偏转时,tip link 松弛允许通道关闭。机械敏感性离子通道是非选择性的(nonselective),即对 K+、Na+、Ca^{2+} 等阳离子都允许通过;通过这些通道的离子流调制造成毛细胞膜的超极化或去极化。膜电位随后由其他离子通道进一步调制,包括 K+ 通道、Ca^{2+} 敏感的 K+ 通道,以及电压敏感的 Ca^{2+} 通道——这三类通道共同决定了毛细胞的电学响应特性。毛细胞的结构、调谐、敏感性和功能由 Hudspeth(1985, 1989; Hudspeth and Gillespie, 1994)综述,这些论文是该领域可读的综述。在阶跃电流输入下,毛细胞的膜电位表现出衰减振荡(damped oscillations),振荡的周期和幅度依赖于阶跃大小;因此每个细胞有一个固有振荡频率,并对接近该频率的刺激响应最佳。Crawford and Fettiplace(1981)以及 Ashmore and Attwell(1985)已发展出电学谐振的简单模型,虽不基于已知机制的细节,但能很好描述实验结果。Hudspeth and Lewis(1988a,b)之后的工作使用更细致的模型,表明离子电导的实测性质足以解释牛蛙囊毛细胞中的谐振。

20.3.1 电路类比(An Electrical Circuit Analogue)

本节先建立 Crawford–Fettiplace 与 Ashmore–Attwell 模型的基本电路(图 20.10A,含电感 L、电阻 R、电容 Cm 和漏电导 gp),写出其频域方程并展示带通响应 |V_1|,定义品质因子 Q 量化调谐尖锐度;然后把电路与图 20.10B(电压门控 K+ 电导)和图 20.10C(Ca^{2+} 门控 K+ 电导)的生物物理实现等价起来,并指出对生理参数,Ca^{2+} 门控 K+ 电导模型能给出与实验相符的 Q 值,Ashmore–Attwell 因此得出频率调谐更可能源于 Ca^{2+} 敏感 K+ 电导的结论。Crawford and Fettiplace(1981)和 Ashmore and Attwell(1985)的模型基于图 20.10A 所示电路。在电流输入 I 下,电压 V 满足

\[ \frac{d^2 V}{dt^2} + \gamma \frac{dV}{dt} + \omega_0^2 V = f(t) \quad (20.85) \]

其中

\[ \gamma = \frac{g_p}{C_m} + \frac{R}{L} \quad (20.86) \]
\[ \omega_0^2 = \frac{g_p R + 1}{L C_m} \quad (20.87) \]
\[ f(t) = \frac{1}{C_m} \frac{dI}{dt} + \frac{I R}{L C_m} \quad (20.88) \]

取 f(t) = e^{iωt} 时谐振最易显示,此时 V = |V_1(ω)| e^{iωt},其中

\[ V_1(\omega) = \frac{1}{\omega_0^2 - \omega^2 + i\gamma \omega} \quad (20.89) \]

于是 |V_1| 具有带通频率响应,最大值在 ω̂ 处,其中 ω̂^2 = ω_0^2 − γ^2/2。(20.85) 的解形如 exp(−γt/2) exp(±i√(ω_0^2 − γ^2/4)),因此 ω̂ 略小于系统的固有振荡频率。然而若阻尼很小(即 γ 小),频率响应的最大幅度大致出现在系统的固有振荡频率处。|V_1| 峰的尖锐程度是电路调谐程度的一种度量,峰越尖锐表示频率选择性越大。由于

\[ \left.\frac{d^2}{d\omega^2} \frac{1}{(\omega_0^2 - \omega^2)^2 + \gamma^2 \omega^2}\right|_{\omega = \hat{\omega}} = \frac{4}{\gamma^2 (2Q^2 - 1)} \quad (20.90) \]

其中

\[ Q = \frac{\omega_0}{\gamma} \quad (20.91) \]

可见 Q(常称为品质因子)是调谐程度的有用度量。Q 越大,电路的频率选择性越强。现在考虑电路对正弦变化电流的响应。当 I = e^{iωt} 时,

\[ V_1(\omega) = \frac{R}{LC_m} + \frac{i\omega}{C_m} \frac{1}{\omega_0^2 - \omega^2 + i\gamma \omega} \quad (20.92) \]

这仍对应一个带通滤波器,最大响应发生在 ω̂ 处,其中

\[ (\hat{\omega}^2)^2 + 2\left(\frac{R}{L}\right)^2 \hat{\omega}^2 + \left(\frac{R}{L}\right)^2 (\gamma^2 - 2\omega_0^2) = 0 \quad (20.93) \]

Crawford and Fettiplace(1981)用该类型的模型(不含漏电导 g_p)来确定龟耳蜗毛细胞的电学调谐特性。通过把这些调谐曲线与声学刺激毛细胞得到的调谐曲线相比较,他们能够确定电学谐振可以解释毛细胞大部分的频率选择性。虽然上述电路显示出所需的谐振,但若能用电学元件与毛细胞联系更直接的电路来解释电学谐振则更令人满意。这至少有两种方式可实现:图 20.10B 和 C 显示两个电路,分别涉及电压敏感 K+ 电导和 Ca^{2+} 敏感 K+ 电导,至少在形式上与图 20.10A 的电路等价。先考虑图 20.10B 的电路。若漏电导为常数,且 K+ 电导是时间和电压的函数,则

\[ I = C_m \frac{dV}{dt} + g_L(V - V_L) + f g_K(V - V_K) \quad (20.94) \]
\[ \tau \frac{df}{dt} = f_\infty - f \quad (20.95) \]

其中 f_∞ 取线性近似

\[ f_\infty = f_r + \mu(V - V_r) \quad (20.96) \]

这里 V_r 假设为静息膜电位,f_r 是 V = V_r 时 f 的值,μ 是稳态激活曲线的斜率。注意在稳态下,

\[ 0 = g_L(V_r - V_L) - g_K f_r (V_K - V_r) \quad (20.97) \]

因此可以消去 V_L。于是

\[ \frac{dI}{dt} + \frac{I}{\tau} = C_m \frac{d^2 \tilde{V}}{dt^2} + \left(g_L + g_K f_r + \frac{C_m}{\tau}\right) \frac{d\tilde{V}}{dt} + \left(g_L + g_K f_r + \frac{g_K(V_r - V_K)\mu}{\tau}\right) \tilde{V} \quad (20.98) \]

其中 Ṽ = V − V_r,并通过假设 V ≈ V_r 在 V_r 附近线性化。如果

\[ L = \frac{\tau}{g_K(V_r - V_K)\mu} \quad (20.99) \]
\[ R = \frac{1}{g_K(V_r - V_K)\mu} \quad (20.100) \]
\[ g_p = g_L + g_K f_r \quad (20.101) \]

则方程 (20.98) 等价于 (20.85)–(20.88)。对图 20.10C 的电路可遵循类似过程,其中 K+ 电导依赖于 Ca^{2+} 而非电压,但必须对 Ca^{2+} 动力学作额外假设。作为一阶近似,假设 Ca^{2+} 通过通道进入细胞的速率是电压的线性函数(斜率 θ),并按一阶动力学移除(即速率与其浓度成正比)。最后假设开放 K+ 通道的比例与 Ca^{2+} 浓度线性相关。因此

\[ I = C_m \frac{dV}{dt} + g_L(V - V_L) + g_K k c (V - V_K) \quad (20.102) \]
\[ W \frac{dc}{dt} = \frac{I_r + \theta(V - V_r)}{F} - p c \quad (20.103) \]

其中 c 表示 Ca^{2+} 浓度,F 是法拉第常数,W 是细胞体积,p 是 Ca^{2+} 泵送速率,I_r 是 V = V_r 时的稳态 Ca^{2+} 电流。常数 k 是 Ca^{2+} 激活 K+ 电流的速率。同样在稳态附近线性化该系统,可得与 (20.85)–(20.88) 等价的系统,只要

\[ L = \frac{W F}{g_K(V_r - V_K) k \theta} \quad (20.104) \]
\[ R = \frac{p F}{g_K(V_r - V_K) k \theta} \quad (20.105) \]
\[ g_p = g_L + g_K k c_r \quad (20.106) \]

其中 c_r 是静息电位处的稳态 Ca^{2+} 浓度。Ashmore 和 Attwell 表明,虽然带电压敏感 K+ 电导的模型能产生范围广泛的最优频率,但生理参数下所得品质因子 Q 的值比应有的低一个数量级。因此,对合理的参数,模型能区分频率但尖锐度不够。Q 的实验值常常为 5 或以上,而模型中 Q 不超过 0.7。他们认为这是 μ 值太小的结果:μ 的生理值约为 0.33 mV^{−1},但要让模型中 Q 足够大,μ 需要约 3 mV^{−1}。因此 K+ 电流的电压激活斜率似乎不足以解释毛细胞中观察到的谐振。然而在第三个模型中,K+ 电流由 Ca^{2+} 激活可以陡得多。这里 K+ 通道的有效激活斜率为 kθ/(pF),可通过降低泵送速率 p 或提高 K+ 通道对 Ca^{2+} 的敏感性使其变大。Ashmore 和 Attwell 的结论是:毛细胞中的频率调谐更可能是 Ca^{2+} 敏感 K+ 电导的结果,而非电压敏感电导的结果。

20.3.2 频率调谐的机制性模型(A Mechanistic Model of Frequency Tuning)

本节先介绍 Hudspeth–Lewis 的牛蛙囊实验工作,确认 Ca^{2+} 敏感 K+ 电导是频率调谐的定量解释,然后描述其机制性模型:毛束偏转开启转导通道让 K+ 进入 → 去极化激活电压门控 Ca^{2+} 通道 → 细胞内 Ca^{2+} 升高 → Ca^{2+} 敏感 K+ 通道开启 → K+ 外流使细胞复极化 → Ca^{2+} 泵维持平衡,并指出 K+ 既可被动进入又可被动离开细胞这一非典型特征;最后写出三类电流方程 (20.107)(含电压门控 Ca^{2+} 电流、Ca^{2+} 激活 K+ 电流 Ikc、漏电流),指出 Ikc 通道有三闭态两开态的细致模型,参数约 30 个,由电压钳数据约束,最终模型对电流脉冲的响应显示出与实验一致的衰减振荡,转导电流的频率响应为带通,最佳频率 112 Hz、Q ≈ 3,而改变 gkc 可调节谐振频率。这一结论由 Hudspeth and Lewis(1988a,b)的工作进一步支持。基于一系列实验,他们测量了牛蛙囊毛细胞离子电导的动力学性质,构建了这些细胞电学谐振的细致模型。他们的结论是:Ca^{2+} 敏感 K+ 电导的性质(连同电压敏感 Ca^{2+} 电导和漏电导)足以定量解释这些细胞中的频率调谐。其模型的示意图见图 20.11。毛束的机械偏转开启毛束上的转导通道,使正离子(主要是 K+)进入。细胞随之去极化,激活电压门控 Ca^{2+} 通道,细胞内 Ca^{2+} 浓度上升。这又开启 Ca^{2+} 敏感 K+ 通道,K+ 流出细胞,细胞复极化。Ca^{2+} 平衡由把 Ca^{2+} 移出毛细胞的泵维持。该模型的一个关键且相当不寻常的特征是 K+ 既可被动进入也可被动离开细胞。由于毛束伸入 scala media,毛束周围的液体(耳蜗中毛细胞的情形下即内淋巴液)的成分与毛细胞基部周围的液体不同,前者 K+ 高、Na+ 低。这里不展示该模型的所有细节。可以说模型假设存在三类对毛细胞谐振有显著贡献的离子电流:电压门控 Ca^{2+} 电流、Ca^{2+} 激活 K+ 电流和漏电流。因此对施加的电流 I,有

\[ I = C_m \frac{dV}{dt} + I_c + I_{kc} + I_L \quad (20.107) \]

电压门控 Ca^{2+} 电流 I_c 和漏电流由与 Hodgkin–Huxley 模型(第 5 章)类似的方程描述。Ca^{2+} 激活 K+ 通道的模型 I_{kc} 则复杂得多。假设该通道有三个关闭状态和两个开启状态:结合两个 Ca^{2+} 离子把通道转换为可以自发开启的状态,再结合一个 Ca^{2+} 离子可以延长开启时间——也就是说,Ca^{2+} 既是通道开启的触发器,也是维持开启的稳定剂。跃迁速率常数依赖于 Ca^{2+} 和电压。最后,Ca^{2+} 的处理通过假设 Ca^{2+} 通过 Ca^{2+} 通道进入并按一阶过程移除而简单建模。模型的参数(约 30 个)通过让模型与来自单个细胞的电压钳数据吻合来确定,然后考察模型对电流脉冲的响应。结果发现,去极化电流阶跃在模型中诱发衰减膜电位振荡,其频率和幅度依赖于电流阶跃的大小,与实验数据(图 20.12)高度吻合——图 20.12A 显示牛蛙囊毛细胞对不同大小去极化电流阶跃的实际反应,电流越大振荡越明显;图 20.12B 显示模型的对应反应,二者定性一致。为模拟转导电流,在 (20.107) 右端加一项 I_T = g_T(V − V_T)。转导电导假设为毛细胞位移的函数,而毛细胞位移又假设为正弦变化。所得模型频率响应为带通,最大响应在频率 112 Hz 处,品质因子为 3。最大响应处的频率(谐振频率)是模型参数的函数,对模型参数的现实修改可解释牛蛙囊中实验观察到的谐振频率范围。由于谐振频率对 g_{kc} 敏感,模型预测控制 Ca^{2+} 敏感 K+ 通道的数目是细胞调节其频率响应的一种简单方式。

20.4 非线性耳蜗放大器(The Nonlinear Cochlear Amplifier)

本节从第 19 章的视觉系统类比出发,指出耳蜗灵敏度极高(仅受热噪声限制),但耳蜗液本应阻尼低强度信号,由此引出 Gold(1948)的早期"主动机制"假说;接着综述该机制被认为由外毛细胞的主动力学放大提供——外毛细胞在去极化和超极化时分别收缩和伸长,基底膜振动通过侧向运动打开/关闭离子通道,分子机制涉及对膜电位敏感而不需要 ATP 的膜内蛋白——并明确指出该蛋白的身份仍未确定。在第 19 章中我们看到,光感受器能够以可靠且可重复的方式响应单光子;视觉系统因而获得了最大可能的灵敏度。同样,耳蜗最非凡的特征之一是其极高的灵敏度,该灵敏度似乎只受热噪声限制(Manley, 2001)。考虑到耳蜗液本应阻尼低功率信号,耳蜗能检测低幅度信号是令人惊讶的——这里有一个明显的力学悖论:基底膜要振动,但其周围流体应吸收振动能量。正如 Hudspeth(2005)形象地表达的那样,让基底膜在这样的条件下振动无异于让音叉在蜂蜜中振动(akin to asking a tuning fork to vibrate in honey)。早在 1948 年,Gold 就指出了这一困难,并提出存在某种主动机制放大对小声音的响应。如此高的灵敏度(以及尖锐的频率调谐)被认为来自毛细胞的主动力学放大。这种放大的确切分子基础尚不清楚,而且几乎可以肯定在不同物种间存在差异。在哺乳动物中,外毛细胞在去极化和超极化时分别收缩和伸长。基底膜的振动导致外毛细胞顶部毛束的侧向运动,进而开启和关闭静纤毛中的离子通道,引起膜电位变化(Hudspeth, 1997; Nobili et al., 1998; Manley, 2001; Hudspeth et al., 2000)。值得注意的是,毛束对离子通道的开关响应是直接的力学相互作用,因此极快。弹性弹簧(tip link)把一根静纤毛的尖端直接连接到离子通道;当静纤毛运动时,tip link 被拉动,从而开启离子通道。外毛细胞长度的变化似乎来自外膜内数十亿蛋白质的收缩或扩张,并且这些蛋白质对膜电位直接反应而无需 ATP。然而这些蛋白质的精确身份和性质仍不清楚(Ashmore et al., 2000)。

20.4.1 负刚度、适应与振荡(Negative Stiffness, Adaptation, and Oscillations)

本节先描述 Martin 等(2000)的牛蛙囊毛束自发抽动/振荡现象——当被谐振频率的刺激牵引时响应幅度约为初始刺激的两倍——并解释其机制:毛束在小位移下表现出负刚度(与常规弹簧的"扰动→回复力"截然相反,使零位移成为不稳定平衡点,存在两个稳定静止态),且具有位移适应特性(外加位移后力–位移曲线向施加方向平移直至原静止态消失);最后指出这两个特性结合可产生自发振荡(毛束在两个稳定态之间反复翻转伴随曲线持续适应),建模留为 Exercise 8。牛蛙囊的毛束表现出自发抽动或振荡,当被谐振频率的刺激牵引时,响应幅度约为初始刺激的两倍。Martin et al.(2000)提出了一种可能的机制,其基础是毛束在小位移下具有负刚度。这是一个非凡的结果。典型的弹簧(如数学和物理入门课所教)具有正刚度:当从静止位置扰动时,弹簧施加一个试图恢复静止位置的力,即正位移产生正力,负位移产生负力。但牛蛙囊的毛束有更为有趣的响应(图 20.13A)。小的正位移导致负力,倾向于增大位移,从而引起正反馈和静止状态不稳定。类似地,负位移产生正力。因此存在两个稳定静止态:一个为正位移,另一个为负位移。此外,毛束表现出适应,以类似光感受器适应的方式(图 19 章)。这在图 20.13B 中说明。若外加一个固定位移,力–位移曲线同时水平和垂直平移,使施加位移位置处于负刚度区域内。因此,对阶跃性负载增加的响应是初始大幅增加力,随后较慢地弛豫。若维持的位移足够大,力–位移曲线的这种平移导致两个稳态经鞍-结分岔消失。因此自发振荡可以产生。当毛束试图停留在正位移的稳定稳态时,力–位移曲线逐渐适应,向上和向右平移,直到该稳定稳态不再存在。然后毛束被迫翻转到另一个稳定稳态(负位移处),之后过程反向重复。这导致一种自发振荡,其中毛束在两个稳定稳态之间反复翻转,力–位移曲线持续适应。该过程的数学建模留作练习(Exercise 8)。

20.4.2 非线性压缩与 Hopf 分岔(Nonlinear Compression and Hopf Bifurcations)

本节先描述实验观察(活体栗鼠数据):低幅度信号下频率响应曲线尖锐、高幅度下变宽,谐振频率处响应快速饱和("非线性压缩");再指出 Eguíluz 等(2000)认为这两种性质可由接近 Hopf 分岔解释;接着写出 Hopf 标准形 (20.108)–(20.110),加正弦驱动 αe^{iωt} 后求形如 A e^{i(ωt+φ)} 的响应,由模方程 (20.112) 得到 A 与 ω 的关系(图 20.15),低 α 时表现为尖锐调谐与非线性压缩;最后指出在 μ = 0 处 A ∝ α^{1/3},但实验验证仍是开放问题。耳蜗的另一个重要特征是:低幅度信号比高幅度信号调谐更尖锐,因此具有更大的频率选择性。此外,没有任何可听声音下耳蜗响应是线性的——耳蜗响应本质上是非线性的。典型的实验数据(来自活体栗鼠)见图 20.14。在低信号幅度下,频率–响应曲线尖锐成峰,但在高幅度下变得不那么尖锐。此外,在谐振频率处(垂直虚线所示),响应随刺激幅度增加而迅速饱和。这称为非线性响应压缩。换言之,对小幅度刺激的响应速度远大于耳蜗若呈线性行为时所应有的速度。Eguíluz et al.(2000)已表明响应的这两个性质都可能来自接近 Hopf 分岔。由于我们已从生物物理模型看到牛蛙囊毛细胞似乎动力学上接近 Hopf 分岔,这为图 20.13 中看到的响应提供了一个合理且吸引人的解释。在标准形(Wiggins, 2003)中,Hopf 分岔的方程为

\[ \frac{dx}{dt} = \mu x - \omega_0 y - (x^2 + y^2) x \quad (20.108) \]
\[ \frac{dy}{dt} = \omega_0 x + \mu y - (x^2 + y^2) y \quad (20.109) \]

或写成复数形式

\[ \frac{dz}{dt} = (\mu + i \omega_0) z - |z|^2 z \quad (20.110) \]

其中 z = x + iy。该系统在 μ = 0 处有 Hopf 分岔,分岔处出现的振荡周期为 2π/ω_0,形式为 z = √μ e^{iω_0 t}。若用驱动项 α e^{iωt} 驱动该振子,并寻找形如 z = A e^{i(ωt+φ)} 的解,则得

\[ (i A (\omega - \omega_0) - \mu A + A^3) e^{i\varphi} = \alpha \quad (20.111) \]

由于只关心响应的幅度 A,对 (20.111) 两边取模即可,得

\[ \alpha^2 = A^6 - 2\mu A^4 + [\mu^2 + (\omega - \omega_0)^2] A^2 \quad (20.112) \]

对给定的刺激强度(即给定 α),可把 A 画为刺激频率 ω 的函数。结果见图 20.15。当 α 小时,响应有尖锐的频率调谐,且在谐振频率 ω = ω_0 处幅度呈现非线性压缩。这一行为与实验(图 20.14)观察到的定性一致。在谐振频率 ω = ω_0 处,(20.112) 变为

\[ \alpha = A^3 - \mu A \quad (20.113) \]

由此可知,在 Hopf 分岔点处(即 μ = 0),A 与输入幅度的 1/3 次方成正比,因此响应被非线性压缩,正如实验数据所见。任何特定动物的毛细胞是否确实通过处于 Hopf 分岔附近而以这种方式调谐其响应仍不确定,造成非线性压缩、尖锐调谐和放大的精确力学也仍不清楚。尽管如此,该简单理论为这一现象提供了一个引人入胜的解释。

本章个人批注

本章是 Keener 第二卷中"感觉器官"主题的最后一章:第 19 章的视觉和第 20 章的听觉都把"接受器–放大器"机制(视网膜对单光子的捕获 / 耳蜗对小振动的非线性放大)作为贯穿全章的核心框架。这一对照结构在 20.4 节的引子中由 Keener 本人点明——把第 19 章"视网膜能可靠响应单光子"和本章"耳蜗灵敏度仅受热噪声限制"并列——并由此引出"主动放大"这一跨章节的共同主题。在力学层面,本章给出了三个层次的模型:流体层面(不可压缩无粘的 Laplace 方程 + Bernoulli 关系),基底膜层面(点阻抗谐振子或多尺度慢变包络),耦合层面(Siebert 模型的二阶 ODE 或 Ranke 短波近似的一阶 ODE)。这一"流体–膜–驱动"的分层是后续做主动力学有限元模型时(如 Ramamoorthy et al. 2007 的声–电–机耦合模型)必读的基础。值得注意的是 20.2 节虽然篇幅最长,但严格来说它处理的是"被动"耳蜗(即外毛细胞尚未引入放大反馈),而 Keener 在 20.2.5 中坦白:"Simple hydrodynamic models do not reproduce the degree of tuning observed in the mammalian cochlea"——这正是 20.4 节为何必须存在的原因。短波近似下 xp ∝ −(1/λ) ln(ω) 这一对数律与生理观察的"倍频程等距"几乎一致,但它对相位的预测过大——这是被动模型的系统性缺陷,也是非线性 Hopf 解释被提出来的实验动机。Hudspeth–Lewis 模型约 30 个参数、由单细胞电压钳数据约束、得到谐振频率 112 Hz 和 Q ≈ 3——这是一个相对完整的"机制性"模型,其漂亮之处在于把 K+ 既被动进入又被动离开细胞这一解剖学非典型特征直接编码进了方程 (20.107)。本章末尾 Hopf 分岔标准形 (20.108)–(20.110) 是动力系统中与 20.2.4 多尺度分析互补的另一条主线,前者处理时间维度(振幅/压缩)而后者处理空间维度(包络/位置),两者结合才构成"低强度下尖锐调谐且非线性压缩"的完整图像。读完后值得在 questions.md 中追踪:(a) Hudspeth 2005 所说的"tip link 与通道直接的力学耦合"是否已被冷冻电镜结构证实?(b) Hopf 分岔假说在 2009 年之后是否被活体钙成像实验直接支持或反驳?

与上下章的衔接(一段话)

本章位于 Keener 第二卷"细胞生理与器官生理"的收尾,紧接第 19 章视网膜与视觉——Keener 在 20.4 节的引子中明确把第 19 章的"光感受器对单光子的可靠响应"与本章的"耳蜗灵敏度仅受热噪声限制"并列为两类感受器共同面对的物理挑战,并用"honey 中振动的音叉"作类比把主动放大的必要性凸显出来。本章向上承接第 7 章钙动力学的离子通道建模(K+ 通道、Ca^{2+} 通道、Ca^{2+} 激活 K+ 通道在 20.3 节被串联起来)和第 5 章 Hodgkin–Huxley 模型的方程骨架(20.3.2 直接援引 I_c 方程的类似性);向下则与全书的"非线性动力学"主轴呼应——20.2.4 的多尺度分析和 20.4.2 的 Hopf 标准形分别从空间和时间两个维度展示慢变包络与振幅分岔如何共同产生位置理论所要求的尖锐调谐;同时 20.4.1 的鞍-结分岔与第 1 章酶动力学的多稳态结构一脉相承。整章的功能可以概括为:"从流体和膜的宏观方程推到电学谐振的电路类比,再推到主动放大的生物物理机制,最后用 Hopf 分岔标准形给出非线性压缩的极简统一解释"——它是全书中"理论模型如何与具体器官物理耦合"最具代表性的实例。