¶ 生物物理学
生物物理学是运用物理学的理论、方法和技术研究生物系统的结构与功能的交叉学科。它从分子、细胞到生物体多个尺度上揭示生命过程的物理机制，是现代生命科学的重要组成部分。
¶ 学科概述
生物物理学（Biophysics）位于物理学与生物学的交叉地带，其核心目标是：
- 定量理解生物系统的行为规律
- 用物理原理解释生命现象的本质机制
- 发展实验与理论工具来测量和模拟生物过程
¶ 研究尺度与子领域
| 尺度 | 典型对象 | 对应子领域 | 特征尺寸 |
|---|---|---|---|
| 分子级 | DNA、蛋白质、脂质 | 分子生物物理 | 1–10 nm |
| 亚细胞级 | 细胞膜、细胞骨架、细胞器 | 膜生物物理 | 10–1000 nm |
| 细胞级 | 神经元、肌肉细胞 | 细胞生物物理 | 1–100 μm |
| 组织/器官级 | 神经网络、骨骼 | 神经物理、生物力学 | 1 mm–1 m |
¶ 一、分子生物物理
¶ 1.1 蛋白质折叠与能量景观
蛋白质从一条无序的多肽链折叠为特定的三维结构，是生命中最精妙的物理过程之一。
能量景观（Energy Landscape）模型：蛋白质折叠可以理解为一个高维能量曲面上的搜索问题。
假设蛋白质的构象空间有 个自由度，其势能函数为 ，其中 表示所有原子的坐标。天然态（折叠态）对应全局能量极小值。
能量景观示意：
能量 E
↑
| ╱╲ ╱╲
| ╱ ╲╱╲ ╱ ╲ ← 局部极小（折叠中间态）
| ╱ ╲╱╲ ╲
|╱ ╲ ╲
| ╲ ╲
| ╲ ╲
| ╲ ╲___← 漏斗底部（天然态，全局极小）
└──────────────────────────────→ 构象坐标
具体数值例子：一个由 100 个氨基酸组成的蛋白质，若每个氨基酸有 3 种主链构象，则总构象数为 。如果以每秒 个构象的速度搜索，需要 秒——远超宇宙年龄（ 秒）。然而实际折叠时间只有毫秒到秒级别，这就是 Levinthal 悖论。
实际折叠数据（以核糖核酸酶 A 为例）：
| 参数 | 数值 | 说明 |
|---|---|---|
| 氨基酸数 | 124 | 4 个二硫键 |
| 折叠时间（25°C） | ~1.5 s | 从完全变性状态 |
| 天然态半衰期 | ~8 h（37°C） | 热稳定性 |
| 折叠能垒 | ~20–50 kJ/mol | 由过渡态决定 |
| 天然态自由能 | 至 kJ/mol | 相对于未折叠态 |
¶ 1.2 DNA 与 RNA 的物理性质
双螺旋结构的力学参数：DNA 可以看作一根弹性杆，其构象由三个力学常数决定：
| 力学参数 | 符号 | 数值 | 物理意义 |
|---|---|---|---|
| 持续长度 | ~50 nm (≈150 bp) | 超过此长度 DNA 开始弯曲 | |
| 扭转持久长度 | ~100 nm | 扭转刚度 | |
| 拉伸模量 | ~1000 pN | 拉伸刚度 | |
| 螺旋螺距 | 3.4 nm (10.5 bp) | B-DNA 螺距 | |
| 直径 | 2.0 nm | 双螺旋直径 |
单分子力谱实验：用光镊或磁镊可以测量 DNA 在外力下的拉伸行为。
其中 为施加力， 为拉伸量， 为松弛长度， 为玻尔兹曼常数， 为绝对温度。
实验数据对比（lambda噬菌体 DNA，48,502 bp）：
| 外力 (pN) | 相对伸长 () | 观察到的现象 |
|---|---|---|
| 0.1 | 1.01 | 热波动主导，弹性类似蠕虫状链 |
| 1.0 | 1.10 | 熵弹性占主导 |
| 10 | 1.30 | 焓弹性开始显现 |
| 65 | 1.70–1.80 | B-S 转变，螺旋结构被强制拉伸 |
| >100 | 2.00+ | DNA 双链开始解离（熔解） |
¶ 1.3 分子马达
分子马达是将化学能（ATP 水解）转化为机械能的蛋白质机器，著名例子包括肌球蛋白、驱动蛋白和动力蛋白。
驱动蛋白（Kinesin）的步进过程：
步进循环：
ATP ADP+Pi
│ │
┌──▼───────▼──┐
│ │ │ ← 头1结合微管, 头2空闲
│ ──┴── │
└─────────────┘
↓ ATP 结合头2
┌─────────────┐
│ ││ │ ← 头2向微管旋转
│ ──┴▼─── │
└─────────────┘
↓ 头2结合，头1释放 ADP
┌─────────────┐
│ ││ │ ← 前进约 8 nm
│ ──┴── │
└─────────────┘
↓ 头1释放，进入下一周期
关键性能参数：
| 参数 | 数值 | 驱动蛋白 | 肌球蛋白V |
|---|---|---|---|
| 步长 | nm | 8.3 | 36 |
| 最大速度 | nm/s | 800 | 450 |
| ATPase 速率 | s⁻¹ | ~100 | ~15 |
| 单步耗能 | ATP/步 | 1 | 1 |
| 做功 | k_B T/步 | ~100 | ~100 |
| 效率 | % | ~50 | ~40 |
| 持续步数 | 步 | >100 | ~10 |
¶ 二、膜生物物理
¶ 2.1 脂质双层的物理模型
细胞膜的基本结构是脂质双层，厚度约 4–5 nm。其力学性质可以用连续弹性理论描述。
弯曲能量密度（Helfrich 模型）：
其中 为弯曲模量， 为高斯弯曲模量，, 为主曲率， 为自发曲率。
典型数值：
| 脂质类型 | (× ) | 厚度 (nm) | 相转变温度 (°C) |
|---|---|---|---|
| DOPC | 20–25 | 3.7 | -17 |
| DPPC | 30–40 | 4.1 | 41 |
| 细胞膜（混合） | 15–30 | 4–5 | 可变 |
| 胆固醇-MPC | 25–35 | 4.5 | 取决于比例 |
¶ 2.2 离子通道与膜电位
Nernst 方程：离子在膜两侧的平衡电位由浓度差决定。
其中 为气体常数（8.314 J/(mol·K)）， 为绝对温度， 为离子价数， 为法拉第常数（96,485 C/mol）。
具体数值计算：对于哺乳动物细胞在 37°C（310 K）时的钾离子：
典型离子浓度与平衡电位：
| 离子 | 细胞内 (mM) | 细胞外 (mM) | 平衡电位 (mV) |
|---|---|---|---|
| K⁺ | 140 | 5 | -89 |
| Na⁺ | 15 | 150 | +62 |
| Ca²⁺ | 10⁻⁴ (游离) | 2 | +123 |
| Cl⁻ | 10 | 120 | -65 |
实际膜电位：静息膜电位通常约为 mV，靠近 K⁺ 的平衡电位，因为静息时膜对 K⁺ 的渗透性远大于其他离子。
Goldman-Hodgkin-Katz 方程（考虑多种离子的渗透性）：
使用典型相对渗透度 ，计算得 mV。
¶ 三、神经生物物理
¶ 3.1 动作电位的 Hodgkin-Huxley 模型
Hodgkin 和 Huxley（1952 年诺贝尔奖）用鱿鱼巨轴突建立了动作电位的定量模型。
等效电路模型：
细胞外
───────┬────────┬────────┬────────┬──────
║ ║ ║ ║
║ C_m ║ ║ ║
║ ║ ║ ║
───────┴────────┴────────┴────────┴──────
[G_Na] [G_K] [G_L] 电池
细胞内
膜电流由三部分构成：
门控变量动力学：
其中 和 是电压依赖的速率常数。
具体参数（鱿鱼巨轴突，6.3°C）：
| 参数 | 符号 | 数值 |
|---|---|---|
| 膜电容 | 1 F/cm² | |
| 最大钾电导 | 36 mS/cm² | |
| 最大钠电导 | 120 mS/cm² | |
| 漏电导 | 0.3 mS/cm² | |
| 钾平衡电位 | -77 mV | |
| 钠平衡电位 | +50 mV | |
| 漏电平衡电位 | -54.4 mV |
动作电位各阶段的时间过程：
| 时间 (ms) | 膜电位 (mV) | 钠活化 | 钾活化 | 钠失活 | 事件 |
|---|---|---|---|---|---|
| 0 | -65 | 0.05 | 0.32 | 0.60 | 静息状态 |
| 0.5 | -50 | 0.20 | 0.32 | 0.55 | 阈值达到，钠通道开始打开 |
| 1.0 | +10 | 0.72 | 0.35 | 0.28 | 钠电导达到峰值 |
| 1.5 | +40 | 0.95 | 0.42 | 0.12 | 动作电位峰值 |
| 2.0 | +20 | 0.90 | 0.58 | 0.04 | 钠通道失活，钾通道打开 |
| 3.0 | -20 | 0.60 | 0.75 | 0.02 | 钾外流占主导 |
| 4.0 | -60 | 0.20 | 0.60 | 0.08 | 超极化（低于静息） |
| 5.0 | -65 | 0.10 | 0.40 | 0.30 | 恢复到静息状态 |
注意：表中数值是近似值，用于展示各变量在动作电位过程中的相对变化趋势。
¶ 3.2 信息编码：发放率 vs 时间编码
神经元如何编码信息是一个核心问题。两种主要理论：
发放率编码：信息由平均发放率 携带
举例：一个运动神经元在自主运动时的发放率从 10 Hz 增加到 50 Hz，编码了肌肉收缩力的增加。
时间编码：信息也可以由精确的脉冲时间携带
| 编码方式 | 精度要求 | 信息传输率 | 典型场景 |
|---|---|---|---|
| 发放率编码 | 10–100 ms | ~30 bits/s | 运动控制、肌肉力量 |
| 时间编码 | 1 ms | >200 bits/s | 听觉定位、触觉辨位 |
| 相位编码（theta 节律） | 5–10 ms | ~100 bits/s | 海马体位置细胞 |
实际信息量计算（以飞蛾听觉神经元为例）：
- 时间精度： ms
- 每秒独立时间窗：1000
- 每个时间窗平均发放率：0.1 spikes
- 信息量： bits/s
¶ 四、生物力学
¶ 4.1 细胞力学
弹性模量：不同类型细胞的刚性差异巨大。
| 细胞类型 | 弹性模量 (Pa) | 方法 | 典型应用场景 |
|---|---|---|---|
| 红细胞 | 1–5 × 10³ | 光镊拉伸 | 变形通过毛细血管 |
| 成纤维细胞 | 0.5–3 × 10³ | 原子力显微镜 | 组织修复 |
| 软骨细胞 | 2–10 × 10³ | 微压痕 | 关节力学 |
| 骨细胞 | 10–50 × 10³ | 微吸管 | 骨骼力学感知 |
| 癌细胞 | 0.1–1 × 10³ | 原子力显微镜 | 转移能力（更软） |
癌细胞与正常细胞的力学区别：研究表明，癌细胞比正常细胞更软。例如：
- MCF-7 乳腺癌细胞： kPa
- MCF-10A 正常乳腺细胞： kPa
- 这一差异可达到 2–3 倍，被认为是癌细胞易于侵入和转移的力学基础之一。
¶ 4.2 肌肉收缩的滑动丝模型
Hill 模型——肌肉收缩力与速度的经验关系：
其中 为肌肉力， 为收缩速度， 为最大等长收缩力， 和 为 Hill 常数。
实验数据（蛙缝匠肌，0°C）：
| 相对荷载 | 收缩速度 (mm/s) | 功率 (mW/g) |
|---|---|---|
| 0.0 | 4.5 | 0 |
| 0.2 | 3.2 | 0.64 |
| 0.4 | 1.9 | 0.76 |
| 0.6 | 0.9 | 0.54 |
| 0.8 | 0.3 | 0.24 |
| 1.0 | 0 | 0 |
最大功率出现在约 处，这是肌肉收缩速度与力量的最优平衡点。
¶ 4.3 骨骼力学
骨骼是一种复合生物材料，包含羟基磷灰石（无机相，65%）和胶原蛋白（有机相，25%）。
力学性能：
| 骨类型 | 抗压强度 (MPa) | 抗拉强度 (MPa) | 弹性模量 (GPa) | 断裂韧性 () |
|---|---|---|---|---|
| 皮质骨（纵向） | 170–193 | 50–151 | 17–20 | 2–7 |
| 皮质骨（横向） | 133 | 27–42 | 6–13 | — |
| 松质骨 | 2–12 | 1–5 | 0.5–5 | 0.01–0.1 |
¶ 五、生物物理技术
¶ 5.1 结构生物学三大技术
| 技术 | 分辨率 | 样品要求 | 优点 | 局限 |
|---|---|---|---|---|
| X射线晶体学 | 0.5–3 Å | 需生长晶体（数周–数月） | 原子级分辨率 | 需结晶，非生理态 |
| 核磁共振 (NMR) | 1–5 Å | 需高浓度（mM级），~30 kDa | 溶液态，动态信息 | 分子量受限 |
| 冷冻电镜 (Cryo-EM) | 2–4 Å | 快速冷冻，~100 kDa+ | 无需结晶，大分子复合物 | 对小分子分辨率低 |
实际数据处理举例：Cryo-EM 工作流程
样品快速冷冻 ─→ 电子束照射，收集2D投影 ─→
│
2D图像分类（数万–数十万张）
│
3D重构（傅里叶变换）
│
原子模型搭建
分辨率与信息获取（以核糖体研究为例）：
| 分辨率 | 可识别的结构特征 | 典型应用 |
|---|---|---|
| >10 Å | 亚基轮廓，二级结构大致位置 | 低分辨拓扑 |
| 5–10 Å | -螺旋（管状密度），-折叠 | 骨架路径追踪 |
| 3–5 Å | 氨基酸侧链密度，部分水分子 | 药物结合位点 |
 Å |
原子水平结构，水分子，金属离子 | 催化机理研究 |
¶ 5.2 单分子技术
| 技术 | 力范围 | 位移精度 | 力精度 | 典型应用 |
|---|---|---|---|---|
| 光镊 (Optical Trap) | 0.1–100 pN | 0.1–10 nm | 0.1 pN | 分子马达步进、DNA 拉伸 |
| 磁镊 (Magnetic Tweezers) | 0.01–100 pN | 1–100 nm | 0.01 pN | DNA 扭转、拓扑异构酶 |
| 原子力显微镜 (AFM) | 10–10⁴ pN | 0.5–0.1 nm | 1 pN | 蛋白质机械展开 |
| 荧光共振能量转移 (FRET) | — | 1–10 nm | — | 构象变化动态测量 |
AFM 单分子力谱应用实例——蛋白质机械展开：
以**肌联蛋白（Titin）**为基础的单分子力拉伸实验：
- 肌联蛋白的 Ig 结构域在 ~150–250 pN 时解折叠
- 每个解折叠事件对应长度增加 ~25–30 nm
- 典型的力-伸长曲线呈现锯齿状，每个"牙齿"代表一个结构域的展开
力 (pN)
↑
250┤ ╱╲
│ ╱ ╲╱╲
200┤ ╱ ╲╱╲ ← 每个峰对应一个Ig结构域展开
│ ╱ ╲╱╲
150┤ ╱ ╲╱╲
│ ╱ ╲╱╲
100┤╱ ╲╱╲
└──────────────────────────────→ 伸长 (nm)
0 100 200 300
¶ 六、生物物理的前沿方向
¶ 6.1 系统生物物理
整合多种尺度的建模方法：
| 尺度 | 建模方法 | 计算复杂度 | 时间尺度 |
|---|---|---|---|
| 原子级 | 分子动力学 (MD) | 数百万原子·ns | ns–s |
| 粗粒化 | CG-MD | 快 100–1000 倍 | s–ms |
| 介观 | 布朗动力学 | 进一步加速 | ms–s |
| 连续/宏观 | 有限元分析 | 数万元素 | s–hr |
多尺度建模具体案例——离子通道：
- 全原子 MD：计算 K⁺ 选择性过滤器的离子通过能垒（2 Å 精度，~50 ns 模拟）
- 粗粒化模型：模拟整个通道蛋白的构象变化（~1 s）
- 连续模型：用 Poisson-Boltzmann 和 Nernst-Planck 方程模拟离子通量
¶ 6.2 AI 与生物物理的融合
AlphaFold 为生物物理带来了革命性变化：
蛋白质结构预测精度演进（CASP 竞赛数据）：
| 年份 | CASP 轮次 | 最佳方法 | 全局距离测试评分（GDT） |
|---|---|---|---|
| 2018 | CASP13 | AlphaFold 1 | ~60 |
| 2020 | CASP14 | AlphaFold 2 | ~90 |
| 2022 | CASP15 | AlphaFold 2 + RoseTTAFold | ~90+ |
| 2024 | — | AlphaFold 3 (多模态) | 结合小分子、核酸 |
¶ 参考文献
- Phillips, R., Kondev, J., Theriot, J., & Garcia, H. G. (2012). Physical Biology of the Cell (2nd ed.). Garland Science.
- Howard, J. (2001). Mechanics of Motor Proteins and the Cytoskeleton. Sinauer Associates.
- Hille, B. (2001). Ion Channels of Excitable Membranes (3rd ed.). Sinauer Associates.
- Hodgkin, A. L., & Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. The Journal of Physiology, 117(4), 500–544.
- Jumper, J., et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nature, 596, 583–589.
- Alberts, B., et al. (2014). Molecular Biology of the Cell (6th ed.). Garland Science.
- Nelson, P. (2014). Biological Physics: Energy, Information, Life (Updated ed.). W. H. Freeman.
- Bustamante, C., et al. (2000). Single-molecule studies of DNA mechanics. Current Opinion in Structural Biology, 10(3), 279–285.