HOOOS

ProteinMPNN设计序列后 怎么用PyMOL分析突变体相互作用力

0 19 分子模拟践行者 PyMOL蛋白质设计
Apple

ProteinMPNN 是目前蛋白质从头设计（De Novo Protein Design）和序列优化领域的主流工具。但它输出的直接结果是 FASTA 格式的氨基酸序列，而非带有三维坐标的 PDB 结构。

要使用 PyMOL 分析突变体与配体或受体之间的相互作用力，必须先经历“序列转结构”的步骤，随后在 PyMOL 中进行结构对齐、微观相互作用表征和定量对比。以下是完整且符合学术研究标准的分析流程。

一、 突变体三维结构的获取路径

直接在 PyMOL 中对野生型（WT）进行单点突变并分析是不够严谨的，因为 ProteinMPNN 往往会引入多点突变，甚至引起骨架的微调。

1. 严谨路径：高精度结构预测（推荐）

将 ProteinMPNN 生成的突变序列输入至 ColabFold / AlphaFold2 或 ESMFold 中进行单体或复合物结构预测。

  • 优势：能够捕捉到因多位点突变导致的骨架微小位移及侧链的最佳构象。
  • 注意：在后续分析前，需检查突变区域的 pLDDT 置信度（通常 > 85 才能进行可靠的侧链相互作用分析）。

2. 快速路径：PyMOL Mutagenesis Wizard（仅适用于单/双点突变）

如果突变位点极少，且确定不会引起骨架变化：

  1. 在 PyMOL 中打开野生型结构。
  2. 点击菜单栏 Wizard -> Mutagenesis。
  3. 在右侧面板选择对应的氨基酸类型，点击要突变的残基。
  4. 利用底部的 Frame 按钮切换不同的侧链旋转异构体（Rotamers），选择排斥力最小（Clash最小，且存在潜在氢键/盐键）的构象。
  5. 点击 Apply 保存突变。

二、 结构对齐与目标区域选定

在 PyMOL 中同时导入野生型（WT）和突变体（MUT）结构，以便进行直观的对比分析。

# 导入结构
load wildtype.pdb, wt
load mutant.pdb, mut

# 将突变体对齐到野生型上
align mut, wt

注：如果序列差异较大导致 align 效果不佳，可使用 super mut, wt 命令，它基于结构特征进行对齐，不受序列特征限制。

为了方便观察，将突变位点及周边 5 Å 范围内的残基单独显示为 Stick 模式：

# 假设突变位点为 A 链的 45 号和 82 号残基
select mut_residues, (mut and chain A and resi 45+82)
select wt_residues, (wt and chain A and resi 45+82)

# 选定突变位点周围 5 埃范围内的环境残基
select env_mut, (byres (mut within 5.0 of mut_residues))
select env_wt, (byres (wt within 5.0 of wt_residues))

# 设定显示样式
preset.technical(selection='all')  # 使用学术预设样式
show sticks, mut_residues or env_mut or wt_residues or env_wt

三、 微观相互作用力分析

蛋白质相互作用的核心在于氢键、盐键、疏水作用、Pi-Pi堆积以及空间位阻（Clash）。我们需要对比 WT 和 MUT 在这些维度上的改变。

1. 氢键与极性相互作用（Hydrogen Bonds）

ProteinMPNN 经常通过改变极性网络来优化结合能。

  • 图形界面操作：
    点击右侧对象列表（Names Panel）中 mut 旁边的 A -> find -> polar contacts -> to any atoms。PyMOL 会自动用黄色虚线标出极性作用。
  • 命令行精确测量（推荐，便于记录测距）：
    # 测量突变残基与周边环境在 3.5 Å 内的氢键
distance mut_hbonds, mut_residues, env_mut, 3.5, mode=2
# 测量野生型对应的氢键
distance wt_hbonds, wt_residues, env_wt, 3.5, mode=2
    mode=2 参数仅显示符合氢键几何学特征的连线，避免虚高。

2. 盐键分析（Salt Bridges）

ProteinMPNN 极易在蛋白质表面或界面设计出带电荷残基。强烈的静电吸引（盐键）对稳定结构至关重要。

  • 原理：带正电残基（Lys, Arg, His）的侧链末端氮原子，与带负电残基（Asp, Glu）侧链末端氧原子，距离在 4.0 Å 以内。
  • PyMOL 识别命令：
    # 选定突变体中的酸碱残基
select mut_positive, mut and (resn LYS+ARG+HIS)
select mut_negative, mut and (resn ASP+GLU)
# 寻找它们之间的盐键
distance mut_salt, mut_positive, mut_negative, 4.0, mode=2
    对比 WT 与 MUT 之间的盐键网络数量。若突变体新增了盐键，说明局域稳定性或结合亲和力得到了提升。

3. 疏水堆积与疏水斑块（Hydrophobic Interactions）

当 ProteinMPNN 将表面的亲水残基突变为疏水残基时，可能造成蛋白质聚集；但若在核心区（Core）或紧密界面优化了疏水堆积，则会显著提高热稳定性。

  • 可视化疏水表面：
    通过将疏水氨基酸高亮，观察其空间排布。
    # 定义疏水残基
select hydrophobic, resn ALA+VAL+LEU+ILE+PHE+TRP+MET+PRO+GLY
color gray80, all
color orange, hydrophobic
  • 观察堆积密度：
    将突变区域显示为圆盘空间填充模型（Spheres）或表面（Surface）：
    show spheres, mut_residues or (env_mut and hydrophobic)
set sphere_scale, 0.5
    检查突变是否填补了原有的“空腔”（Pocket/Cavity），或者是否造成了疏水残基非正常暴露。

4. 空间位阻与撞击分析（Steric Hindrance / Clashes）

若 ProteinMPNN 在空间受限区域设计了较大侧链的氨基酸（如 Ala -> Phe），可能会引入致命的位阻。

  • 排斥力检测：
    # 显示范德华半径重叠
show vdw, mut_residues
    如果侧链的球体出现大面积重叠，或者使用 PyMOL 的 Measurement -> Clashes 插件检测到红色报警盘，则需要对突变体进行局部能量最小化（如使用 Rosetta Minimize 或 Amber 进行简短弛豫）。

5. Pi-Pi 堆积与阳离子-Pi 作用（Pi-Pi & Cation-Pi）

这在蛋白质相互作用界面（PPI）中非常常见，涉及芳香环残基（Phe, Tyr, Trp）。

  • 观察构象：
    检查突变残基（如突变为了 Phe）的芳香环与周围芳香环的距离。
    • Face-to-face 堆积：两环中心距离 3.5 - 5.0 Å，且环面接近平行。
    • Edge-to-face 堆积（T-shaped）：一个环的边缘指向另一个环的中心，距离在 5.0 - 7.0 Å 之间。

四、 结果对比与学术作图输出

在完成相互作用分析后，需要输出能够清晰展示突变效果的对比图。

1. 双窗口对比（Split-screen）

为了避免 WT 和 MUT 的线条重叠在一起导致画面混乱，可以使用 PyMOL 的多网格模式：

set grid_mode, 1

此时 WT 和 MUT 会分左右（或上下）两个窗口同步联动旋转，极大地方便了突变前后氢键/盐键改变的直观对比。

2. 静电势表面分析（APBS）

如果 ProteinMPNN 改变了局部的电荷分布，建议使用 APBS 插件生成静电势表面（Electrostatic Potential Surface）：

  1. 安装并打开 APBS Electrostatics 插件。
  2. 对 WT 和 MUT 分别运行静电势计算。
  3. 观察结合界面是从“电荷相斥”变为了“电荷互补”，这也是解释 ProteinMPNN 设计成功率的强力证据。

3. 高清渲染输出

在导出论文插图前，进行如下参数优化以去除锯齿、增加阴影立体感：

bg_color white                  # 设置背景为白色
set ray_shadows, 1              # 开启阴影
set cartoon_fancy_helices, 1    # 优化螺旋视觉效果
set ambient, 0.35               # 增加环境光亮度
ray 2400, 1800                  # 高清渲染（分辨率可根据需要调整）
png mutant_analysis.png         # 保存图片

点评评价

captcha
健康