蛋白质结构可视化与界面分析
结构图既是展示也是分析：叠合两个结构才能看出构象变化，算出界面残基和埋藏面积才能说明两条链靠什么结合。常用工具是 PyMOL，也有 ChimeraX、VMD 和网页端的 Mol* 等。一张合格的结构图通常有明确的朝向、克制的配色，只保留与结论相关的细节，并标注关键残基和相互作用。界面分析在几何层面给出可量化描述：接触残基列表、氢键与盐桥、埋藏溶剂可及表面积（BSA）、界面残基的极性与疏水组成。典型输入是一个或多个结构文件，输出是图像和界面统计表。
什么时候用
- 论文或汇报中展示蛋白整体折叠、活性口袋或配体结合模式
- 叠合野生型与突变体、apo 与 holo 结构，比较构象变化
- 分析复合物两条链的界面残基、氢键和埋藏面积
- 把保守性、B 因子或 pLDDT 等数值映射到结构上着色
- 批量为多个结构生成统一朝向和配色的图像
不适合的情况
- 根据静态图判断结合强弱：相互作用数量不等于亲和力
- 把预测结构的低置信区域当作真实构象展示
- 用可视化代替动态分析：构象变化幅度与柔性需要轨迹分析
需要准备的数据
- 结构文件
- PDB 或 mmCIF；注意是否含多个生物学组装体、多构象（altloc）和对称拷贝
- 链与残基定义
- 要展示或分析的链、配体残基名、需要标注的残基编号列表
- 数值着色数据（可选）
- 按残基的保守性、pLDDT、突变评分等数值表，需与结构编号一一对应
- 出图要求
- 目标尺寸、分辨率（期刊常要求 300 dpi 以上）、背景透明与否、配色方案
分析步骤
- 1
加载与清理结构
载入后用 remove solvent 去水，按需删除冗余对称拷贝和不相关配体，处理 altloc 只保留一个构象。先看 get_chains 和序列范围，确认链命名和缺失残基位置，缺失片段会在卡通表示中显示为断开。
- 2
建立选择集
用 select 配合 chain、resi、resn、polymer、organic、byres、within 等关键词定义要操作的部分，并给选择集命名。复杂图建议一步步验证每个选择的原子数，避免选到相邻链的同编号残基。
- 3
设置表示与配色
整体用 cartoon，配体和关键残基用 sticks，口袋用 surface 或 mesh。配色宜少：主体一两种浅色，强调对象用高饱和色。show 与 hide 按选择集逐层控制，而不是全局反复切换。
- 4
结构叠合
同源性高用 align（基于序列比对），差异大或序列不同用 super（基于结构）。叠合后记录 RMSD 和参与叠合的原子数，仅用骨架原子（name CA）叠合更稳。局部构象比较时只对保守核心区叠合，再观察目标区域差异。
- 5
界面与相互作用分析
用距离准则选出界面残基，常见做法是取与另一条链任意重原子距离在 4 到 5 Å 以内的残基。再找氢键（供受体距离约 3.5 Å 以内且角度合理）和盐桥（带电基团间约 4 Å 以内），用 distance 对象显示为虚线。
- 6
计算埋藏面积
分别计算各链单独和复合物状态下的溶剂可及表面积，差值的一半通常作为界面埋藏面积。计算前统一 dot_solvent 和 solvent_radius 设置，否则不同次结果不可比。
- 7
出图与复现
固定视角后用 get_view 保存矩阵，设置 ray_trace_mode、抗锯齿和背景，用 ray 指定像素尺寸再 png 导出。把整个流程写成 .pml 脚本，保证图可以重现并批量应用到其他结构。
关键参数与决策
align 与 super 的选择
align 先做序列比对再叠合，适合同一蛋白的不同状态，序列同一性低时比对易出错。super 基于结构特征匹配，对远同源更稳健。两者都可设 cycles 剔除离群原子，报告时注明方法与原子集。
界面距离阈值
界面残基常用 4 Å 或 5 Å 的重原子距离截断。阈值放宽会纳入更多外围残基，界面看起来更大。同一研究内必须统一阈值，并在图注中写明。
表面计算设置
PyMOL 中 dot_solvent 设为 1 计算溶剂可及表面，设为 0 计算分子表面，两者数值差别很大。dot_density 影响精度与速度，比较不同结构时所有参数必须一致。
出图分辨率与配色
单栏图常按宽度约 85 mm、300 dpi 换算像素，再用 ray 指定。配色避免红绿搭配以兼顾色觉障碍读者，背景多用白色或透明，透明表面的 transparency 常设在 0.3 到 0.6。
常见错误
- 直接用不对称单元出图：晶体学不对称单元可能不是生物学组装体，二聚或四聚结构会展示错误，需要用生物学组装文件。
- 不去水、不处理 altloc 就计算界面：水分子和双构象残基会被当作接触原子，界面残基列表虚增。
- 用 align 叠合序列差异大的结构：序列比对错位会让 RMSD 虚高，并把不对应的区域强行重合。
- 用不同表面参数比较埋藏面积：dot_solvent 或 dot_density 不一致时数值不可比，结论会被设置差异左右。
- 视角未固定就分别出图：同一组图朝向不同，读者无法比较，应保存并复用 view 矩阵。
- 展示预测结构时不标注置信度：低 pLDDT 的 loop 在图上与可靠区域看起来一样，容易被过度解读。
结果怎么看
整体结构图
卡通表示的全局视图，配体与关键残基突出显示。看二级结构是否连续、朝向是否让关键区域可见。卡通出现断裂处多为缺失残基，应在图注说明而不是直接连上。
叠合对比图
两个结构以不同颜色叠合，附 RMSD 与叠合原子数。核心区重合而目标 loop 或结构域明显偏移，是有意义的构象差异；若整体处处偏移，先检查叠合方法与所选原子集。
界面相互作用图
界面局部放大，氢键和盐桥以虚线标距离，残基标注编号。关注接触是否形成互补的极性网络；若界面几乎只有疏水接触且距离偏长，结合可能较弱或该构象不是主要结合模式。
界面统计表
列出界面残基数、埋藏面积、氢键与盐桥数量、极性与疏水比例。埋藏面积很小却声称强结合需谨慎；界面集中在柔性末端时，结合模式的可靠性也较低。
和相近方法怎么选
PyMOL 与 ChimeraX
PyMOL 脚本生态成熟、出图风格为许多期刊熟悉，适合批量脚本化出图。ChimeraX 在大体积冷冻电镜密度图、大型复合物和动画上更流畅，内置分析命令更多。两者都支持命令脚本，可按数据类型选择。
命令脚本与界面点选
鼠标点选调整直观，但难以复现，审稿返修时往往无法还原原图。写 .pml 或 Python 脚本可以保存视角、配色和每一步操作，换结构只需改文件名，适合需要统一风格的成组图。
PyMOL 界面分析与专门工具
PyMOL 用距离和表面积即可给出接触残基与埋藏面积，足够一般描述。PDBePISA、PRODIGY 等专门工具还会评估界面是否可能是生物学相关界面或给出结合能估计，适合需要定量比较界面时配合使用。
常见问题
PyMOL 里 align 和 super 有什么区别？
align 先对序列做比对，再用对应残基叠合，适合序列相同或高度相似的结构。super 主要依据结构特征匹配，不要求序列相似，对远同源更稳健。报告 RMSD 时要写清用了哪种方法以及是否只用 CA 原子。
怎么算蛋白界面的埋藏面积？
分别计算两条链单独和形成复合物时的溶剂可及表面积，两者之差即总埋藏面积，习惯上除以二作为界面面积。计算时需固定 dot_solvent、solvent_radius 和 dot_density，并对所有比较对象使用相同设置。
出图要多大分辨率？
多数期刊要求位图不低于 300 dpi，线条图更高。实际做法是按图在版面上的物理宽度乘以 dpi 得到像素数，再用 ray 指定宽高渲染。放大已渲染的小图会模糊，不能替代重新渲染。
怎么按数值给结构着色？
常见做法是把数值写入 B 因子字段，再用 spectrum b 着色并设定数值范围，例如展示 pLDDT 或保守性。注意说明颜色对应的数值区间，否则读者无法判断深浅含义，范围不固定时不同图之间也无法比较。
在 Codar 里怎么做
在 Codar 中提供结构文件并说明要展示或分析什么，它会先确认方案：保留哪些链和配体、叠合方式、界面距离阈值、表面计算参数和出图尺寸配色。随后完成结构清理与选择集构建，进行结构可视化出图，并做结构几何与界面分析，输出接触残基、氢键盐桥与埋藏面积。交付整体结构图、叠合对比图、界面局部图和可复现脚本。
我有一个抗体与抗原的复合物结构，帮我出一张界面放大图标出氢键和盐桥，并给我界面残基列表和埋藏面积。