1. 课程导入与工程全景
在真实药物发现项目中，拿到一个蛋白靶点后，团队最先面对的往往不是“设计分子”，而是一连串更基础的问题：这条序列对应的结构是否可靠？口袋在哪里？这个口袋是否真的可被小分子占据？如果这些问题没有回答清楚，后续的分子对接、虚拟筛选乃至化合物合成都可能建立在错误的地基上。本课程的核心目标，就是建立一套可复现的工程化流程，把“一条蛋白序列”转化为“一组经过验证的可成药位点假设”，并输出可供下游建模使用的结构文件与位点报告。
核心目标与交付物
本课程要交付的不是单一文件，而是一组相互关联的工程产物。最基础的输入通常是靶点的氨基酸序列（FASTA）或一个已知的 PDB 结构编号。经过结构准备后，我们期望得到：
- 清洗后的结构文件：去除非相关水分子、冗余配体、异常残基，补齐缺失原子或缺失环区，并明确质子化状态与二硫键。
- 可成药位点清单：每个位点包含口袋中心坐标、体积、疏水/极性比例、关键残基编号及可成药性评分。
- 位点验证记录：说明该位点来自共晶配体、文献报道还是几何探测，并给出置信等级。
- 可复现的配置与日志：记录所用工具版本、参数、随机种子，使他人能重跑得到一致结果。
这些交付物共同回答一个问题：这个靶点有没有值得投入资源去设计的口袋，以及它在什么位置、由哪些残基构成。
典型工作流
工程全景可以用一条从序列到决策的主线来理解。首先是结构获取与评估：若已有实验结构，检查分辨率、R 因子、缺失残基比例；若无实验结构，则用 AlphaFold2/3 或同源建模生成预测结构，并重点关注 pLDDT 与 PAE 矩阵。第二步是结构准备：加氢、加电荷、修复缺失原子、去除无关配体，必要时保留对口袋有影响的辅因子。第三步是位点探测：既可以从共晶配体直接提取，也可以用 fpocket、SiteMap 等几何方法扫描表面。第四步是可成药性评估：结合口袋体积、疏水分数、封闭程度、氢键供受体分布进行打分。最后是结果整理与交接：输出结构、位点坐标和评估表，供对接或分子生成使用。
常见失败模式与规避思路
从工程角度看，失败往往不是工具不会用，而是关键判断被跳过。第一类问题是结构本身不可用：例如预测结构中口袋区域 pLDDT 低于 70，却仍被当作可靠模板；或者实验结构缺失了口袋周围的环区，导致口袋“漏底”。规避方式是在准备阶段就标记低置信区域，必要时用同源结构补全或缩小分析范围。第二类问题是质子化与电荷处理草率：组氨酸的互变异构、天冬氨酸的质子化状态会直接改变氢键供受体模式，进而影响可成药性判断。第三类问题是把浅槽当口袋：几何探测工具会返回大量候选，但体积小于 100 ų 或疏水比例过低的位点通常不具备实际可成药性。第四类问题是忽略辅因子：激酶结构中的 ATP 结合位点若被当作“空口袋”处理，可能得到与真实结合模式不符的结论。
实战项目设定
为了让课程中的每个环节都有明确对象，本课程设定一个贯穿始终的实战项目：以一个约 300 个残基的示例激酶结构域为靶点，输入为其 FASTA 序列，同时提供一个可选的高分辨率 PDB 模板。项目要求完成以下闭环：
- 获取并评估结构，记录分辨率、缺失残基比例和预测置信度。
- 完成结构准备，输出清洗后的 PDB 或 PDBQT 文件，并说明加氢与质子化策略。
- 运行至少两种位点探测方法，合并候选口袋并去重。
- 对每个候选位点计算体积、疏水比例、封闭度和关键残基，给出可成药性等级。
- 整理为一份位点报告，包含坐标、残基列表、评分和推荐理由。
示例参数仅用于教学演示，例如：口袋体积阈值设为 200 ų，疏水比例阈值设为 0.45，pLDDT 低于 70 的区域不参与口袋中心定位。这些数值不是通用标准，实际项目中应根据靶点家族和文献经验调整。最终判断依据是：位点是否具有足够的空间容纳类药分子，是否由明确的疏水与极性残基共同构成，以及是否与已知配体结合模式一致。