1. 课程导入与工程视角建立
分子动力学（Molecular Dynamics, MD）模拟在理论上可以给出体系随时间演化的完整轨迹，但在实际工程中，常规 MD 面临一个根本性矛盾：可及时间尺度与生物学/材料学感兴趣事件的时间尺度之间存在数量级鸿沟。一个蛋白质的构象翻转、一次配体解离、一段多肽的折叠，往往发生在微秒到秒量级，而全原子 MD 在通用 GPU 上通常只能稳定覆盖纳秒到微秒量级。增强采样（Enhanced Sampling）正是为跨越这一鸿沟而发展出的一整套方法学，而构象分析工程则是把模拟产生的海量构象转化为可解释、可复现、可交付结论的系统化流程。
本课程把这两部分合并为一个工程视角：增强采样负责“把该看到的状态看到”，构象分析负责“把看到的状态讲清楚”。二者不是先后两个独立阶段，而是共享同一套自由能景观（Free Energy Landscape, FEL）语言、同一套反应坐标（Collective Variable, CV）定义、同一套收敛性判据的闭环系统。理解这一点，是后续 29 章所有方法、参数与案例的共同基础。
核心问题：三个必须同时回答的工程问题
无论使用元动力学（Metadynamics）、伞形采样（Umbrella Sampling）、副本交换（REMD）、自适应偏置（ABF）还是机器学习辅助采样，工程上都要回答三个问题：
- 采样什么？——即反应坐标的选择。CV 必须能区分目标状态、覆盖关键势垒，且维度不宜过高（通常 1–3 维可计算自由能面，4 维以上需特殊处理）。选错 CV，再长的模拟也只是在错误方向上“高效地走错路”。
- 偏置多少？——即偏置势的强度、宽度与更新频率。高斯高度过大导致能量面粗糙、重加权困难；过小则收敛缓慢。工程上常用“势垒高度 / 5～10”作为初始高斯高度的经验起点，再根据重加权质量调整。
- 收敛了吗？——即如何判断自由能面已稳定。常用判据包括：块平均（block averaging）误差、多次独立重复的 FEL 重叠度、重加权后直方图重叠（histogram overlap）、以及关键状态间自由能差的置信区间。
这三个问题构成一个闭环：CV 决定偏置方向，偏置决定采样分布，采样分布决定重加权与收敛判断，而收敛判断又反过来提示是否需要修改 CV 或偏置参数。下图给出这一闭环的工程结构，并标出每一步的典型输入、输出与判断依据。
典型应用场景：从药物设计到材料相变
增强采样与构象分析工程的应用场景高度多样，但共享同一套方法学骨架。以下列举四类典型场景，并说明其关键 CV 与关注量：
| 应用场景 | 典型体系 | 关键 CV | 核心关注量 |
|---|---|---|---|
| 药物—靶点结合/解离 | 激酶 + 小分子抑制剂 | 配体—口袋距离、接触数 | 结合自由能 ΔG、解离路径 |
| 蛋白质构象变化 | GPCR、离子通道 | RMSD、二面角、通道半径 | 开放/关闭态自由能差 |
| 多肽折叠/聚集 | 淀粉样肽段 | 回旋半径、β 折叠含量 | 折叠中间态、聚集倾向 |
| 材料相变与扩散 | 固态电解质、合金 | 配位数、有序参数 | 相变势垒、扩散系数 |
以药物—靶点结合为例，工程上通常先用常规 MD 观察配体在口袋内的稳定性，再选择配体—口袋质心距离与关键氢键数为 CV，采用元动力学逐步填充结合自由能面。若配体在 200 ns 内未发生明显解离，则需检查 CV 是否遗漏了口袋门控残基的构象变化，或考虑加入路径 CV（path CV）以描述更复杂的解离通道。
完整工作流：从结构准备到结论交付
一个可复现的增强采样工程流程通常包含以下阶段，每个阶段都有明确的输入、输出与质量检查点：
- 结构准备与力场选择：从 PDB 获取结构，补全缺失残基与氢原子，选择力场（如 AMBER ff14SB、CHARMM36m）与水模型（TIP3P、TIP4P/2005）。检查质子化态、二硫键、金属配位。输出为可仿真的拓扑与坐标文件。
- 体系平衡：依次进行能量最小化、NVT 升温、NPT 密度平衡，通常各 100–500 ps。检查温度、压力、密度、RMSD 是否稳定。输出为平衡后的构象。
- CV 选择与验证：基于平衡轨迹计算候选 CV 的分布，确认目标态可区分、势垒可覆盖。必要时用短时间试探性偏置模拟验证 CV 有效性。
- 增强采样生产：设置偏置参数（高斯高度、宽度、沉积间隔），运行 100 ns–1 μs 量级模拟，保存偏置势、轨迹与 CV 时间序列。
- 重加权与 FEL 构建：使用 WHAM、MBAR 或直接重加权方法，将偏置采样数据还原为无偏自由能面。检查直方图重叠与统计误差。
- 构象分析与交付：聚类提取代表性构象，分析关键相互作用，绘制 FEL 图，整理收敛性报告与参数表。
课程项目总览：一个贯穿全程的案例
本课程将围绕一个统一的案例展开：某激酶与其小分子抑制剂的结合/解离过程。该案例贯穿 30 章，从结构准备、CV 设计、元动力学采样、重加权、FEL 分析，到最终给出结合自由能与解离路径的工程结论。每一章都会在该案例上增加一个可操作的工程环节，并给出对应的参数、代码片段与结果判读方法。
课程结束时，你将能够独立完成以下任务：给定一个蛋白质—配体复合物结构，选择合适的 CV，设计增强采样方案，运行模拟并判断收敛性，最终交付包含 FEL 图、代表性构象与关键相互作用表的完整分析报告。这一能力可直接迁移到 GPCR、离子通道、多肽聚集、材料相变等不同体系，只需替换 CV 与力场参数即可。