课前准备--氨基酸突变对蛋白质四级结构形成的影响

原创

追风少年i

发布于 2026-09-20 08:22:38

980

作者，Evil Genius

很多学员为文章发愁，没事，不要怕，只要坚持，都会有好的结果。很多人失业，也没事，上天让你好好休息，你就好好休息。

看看我，失业、裁员、破产、负债、网暴、被各种人踩、劳动仲裁、拖欠工资、社保断缴、父母病重等灾难全部经历过了，不也活得好好的，虽然没房没车没存款，但总归还有希望。看见过低三下四的自己，看见过自己低谷别人的态度，自己才有机会真正的站起来。

氨基酸突变对蛋白质四级结构的影响

突变 → 局部结构/理化性质变化 → 单体构象变化 → 蛋白–蛋白界面变化 → 寡聚化状态变化 → 四级结构稳定性/组装能力改变

如果要做科研分析，按照下面这套流程。

一、首先明确要回答的是什么问题

蛋白质四级结构主要涉及：

单体是否还能正确折叠
二聚体/三聚体/四聚体是否还能形成
亚基之间的结合界面是否发生改变
界面结合自由能是否改变
突变是否导致 oligomer → monomer
是否改变二聚体 → 四聚体等平衡
是否产生新的错误聚集界面

例如：

WT：

Chain A
          ↓
     ┌─────────┐
     │         │
     └─────────┘
          ╲
           ╲ interface
            ╲
     ┌─────────┐
     │         │
     └─────────┘
          ↑
        Chain B

          ↓

       Stable dimer

Mutation：

Chain A       Chain B
   │             │
   │  interface  │
   X─────────────X
        ↓
  H-bond / salt bridge
      disrupted

          ↓

   weaker dimerization
          ↓
     monomer ↑

二、最重要的一步：判断突变是否位于蛋白–蛋白界面

这是四级结构分析的核心。

如果突变位点距离另一个亚基很远，例如：

Mutation
   ↓

Chain A
┌───────────────────┐
│                   │
│      ★            │
│                   │
└─────────┬─────────┘
          │
       interface
          │
┌─────────┴─────────┐
│                   │
│      Chain B      │
└───────────────────┘

那么它对四级结构的直接影响通常较弱。

但是不能因此认为“没有影响”。

因为可能存在：

Mutation
   ↓
local structure
   ↓
domain movement
   ↓
interface remodeling
   ↓
oligomerization

也就是说：

远离界面的突变也可能通过变构效应影响四级结构。

三、WT 和 Mutant 结构比较

如果已经有 WT 和突变体结构，首先进行结构比对。

例如：

WT：PDB
Mutant：PDB
AlphaFold/AlphaFold-Multimer预测结构
MD轨迹中的代表构象

重点观察：

1. RMSD

可以计算：

RMSD(WT, Mutant)

但最好分别计算：

分析对象	意义
全蛋白 RMSD	总体结构变化
单个亚基 RMSD	单体折叠变化
界面残基 RMSD	界面重排
界面 backbone RMSD	界面骨架变化
mutation附近 RMSD	局部结构变化

例如：

WT → Mutant

Overall RMSD       1.2 Å
Chain A RMSD       0.8 Å
Chain B RMSD       0.9 Å
Interface RMSD     2.6 Å
Mutation region    3.1 Å

这种情况下就很值得研究：

总体结构变化不大，但是蛋白–蛋白界面发生明显重排。

四、分析界面面积，这是四级结构非常重要的指标

可以计算：

Interface Area / Buried Surface Area

即：

两个亚基结合以后，有多少溶剂可及表面积被埋藏。

例如：

指标	WT	Mutant
Interface area	1850 Å²	1420 Å²
ΔArea	—	−430 Å²

如果突变后界面面积明显下降，说明：

突变可能削弱亚基之间的结合。

反过来，如果：

WT        1850 Å²
Mutant    2100 Å²

也可能意味着突变形成了新的、更紧密的界面。

五、分析突变前后界面的相互作用

这是最直观的一部分。

重点分析：

氢键

WT

ASP ─── H ─── ARG
       ↓
    H-bond

Mutant

ASP ────── ARG
     ×

统计：

H-bond 数量
H-bond occupancy
H-bond distance

盐桥

例如：

GLU−  ···  LYS+

突变：

GLU−  →  ALA

可能导致：

salt bridge lost
      ↓
interface stability ↓

疏水相互作用

例如：

Leu
Val
Ile
Phe
Trp
Met

如果：

Phe → Asp

可能导致：

hydrophobic packing ↓
electrostatic interaction ↑

因此不能只看“氢键”。

六、特别值得分析 ΔΔG

如果目标是判断：

突变到底会不会影响二聚体/多聚体稳定性？

那么建议计算：

例如：

WT dimer
ΔGbind = −45 kcal/mol

Mutant dimer
ΔGbind = −32 kcal/mol

那么：

说明突变以后结合变弱。

需要注意：

ΔΔG 的具体数值高度依赖计算方法，不能把单一软件的结果当成实验真值。

可以结合：

FoldX
Rosetta
MM/GBSA
MM/PBSA
MD
实验数据

进行交叉验证。

七、如果突变位于界面附近，要重点分析“接触残基网络”

例如 WT：

ARG120
               │
               │ salt bridge
               ↓
Chain A ─── MUTATION ─── Chain A
                         │
                         │
                      PHE145
                         │
                    hydrophobic
                         │
                         ↓
                    Chain B

突变以后可能发生：

Mutation
   ↓
ARG120 position changes
   ↓
PHE145 rotates
   ↓
new contact formed
   ↓
interface rearrangement

因此建议建立：

Residue–Residue Contact Map

比较：

WT contact map

        Chain B
       120 130 140
100     █
110       ███
120         ██
130       ███

和：

Mutant contact map

        Chain B
       120 130 140
100
110      █
120       █
130         ████

这个对于解释四级结构改变非常有用。

八、如果没有突变体实验结构，可以用 MD

这是比较推荐的路线：

WT
               │
               ↓
          Protein prep
               │
               ↓
        WT dimer/tetramer
               │
               ↓
          MD simulation
               │
               ↓
     ┌─────────┴─────────┐
     ↓                   ↓
interface              dynamics
analysis                analysis


             Mutant
               │
               ↓
        Mutation modeling
               │
               ↓
       Mutant dimer/tetramer
               │
               ↓
          MD simulation
               │
               ↓
     ┌─────────┴─────────┐
     ↓                   ↓
interface              dynamics
analysis                analysis

然后比较 WT vs Mutant。

九、MD中建议重点计算这些指标

指标	判断什么
RMSD	整体结构稳定性
RMSF	哪些残基变得更灵活
Rg	蛋白整体紧致程度
SASA	暴露程度变化
Interface area	亚基接触面积
H-bond	氢键网络变化
Salt bridge	盐桥变化
Contact number	亚基接触数量
Contact occupancy	接触持续时间
Interface distance	亚基间距离
PCA	构象运动模式
MM/PBSA	结合能变化
MM/GBSA	结合能变化

尤其推荐：

RMSF + Interface contact + H-bond + salt bridge + MM/PBSA

联合分析。

十、如果是四聚体，分析方法要进一步升级

例如：

Chain A ─── Chain B
           │           │
           │           │
        Chain D ─── Chain C

这时候不能只分析：

A-B

而应该分析：

A-B
B-C
C-D
D-A

甚至：

A-C
B-D

根据实际四级结构判断。

可以建立一个：

Interface matrix

Interface	WT	Mutant	Δ
A-B	−35.2	−27.1	+8.1
B-C	−34.8	−34.0	+0.8
C-D	−35.1	−27.5	+7.6
D-A	−34.9	−33.8	+1.1

这样马上就能发现：

Mutation主要破坏 A-B / C-D 界面。

十一、还要考虑“寡聚化平衡”

这是很多结构分析容易忽略的。

假设：

WT

Monomer ⇌ Dimer ⇌ Tetramer
              ↑
           dominant

突变以后：

Mutant

Monomer ⇌ Dimer ⇌ Tetramer
   ↑
  increased

那么即使突变体的单体结构：

RMSD ≈ WT

也不能说明没有影响。

因为真正变化的可能是：

association equilibrium

也就是：

因此如果研究问题是“突变影响四级结构形成”，最好最终结合实验：

SEC
SEC-MALS
AUC
native PAGE
crosslinking
analytical ultracentrifugation

验证实际 oligomeric state。

十二、推荐的一套完整分析路线

如果有一个：

WT PDB + 突变位点

可以直接按照：

WT PDB
  │
  ├──→ 确定四级结构
  │       ↓
  │   Chain A/B/C/D
  │
  ├──→ 定位 mutation
  │
  ├──→ 判断是否位于 interface
  │
  ├──→ 计算 interface area
  │
  ├──→ 分析 H-bond
  │
  ├──→ 分析 salt bridge
  │
  ├──→ 分析 hydrophobic contacts
  │
  ├──→ 分析 residue contact network
  │
  └──→ 计算 ΔΔG
          ↓
      WT vs Mutant
          ↓
       MD simulation
          ↓
 ┌────────┼────────┐
 ↓        ↓        ↓
RMSF   Interface  H-bond
       contacts
 ↓        ↓        ↓
PCA     SASA     Salt bridge
          │
          ↓
      MM/PBSA
          │
          ↓
  四级结构稳定性改变
          │
          ↓
  oligomerization改变？
          │
          ↓
     实验验证

如果是做的“蛋白突变 + 分子对接/MD”方向

尤其建议把分析分成 三个层次：

Level 1：结构层面

mutation → backbone/side-chain → interface

Level 2：能量层面

mutation → H-bond/salt bridge/hydrophobic → ΔΔG

Level 3：动态层面

mutation → MD dynamics → interface fluctuation → oligomerization stability

最终可以形成一个非常清晰的科研结论链：

Mutation → local structural perturbation → interface remodeling → binding-energy change → altered inter-subunit dynamics → quaternary-structure stability/oligomerization change

完整的流程：从 PDB → PyMOL界面定位 → FoldX/Rosetta ΔΔG → GROMACS WT/Mutant MD → RMSD/RMSF/界面接触/H-bond/MM-PBSA → 四级结构结论

生活很好，有你更好。

原创声明：本文系作者授权腾讯云开发者社区发表，未经许可，不得转载。

如有侵权，请联系 cloudcommunity@tencent.com 删除。

数据分析

#分子对接

#分子动力学

#突变