# RDCM 分析

## 概述

残基距离接触矩阵（RDCM, Residue Distance Contact Matrix）计算残基对之间的平均距离，分析残基间空间关系和长程接触。适用于分析接触和距离、识别长程相互作用、研究折叠/展开转变、验证实验接触图。

## 工作流程

### 输入文件
- `-s`：结构文件（.tpr）
- `-f`：轨迹文件（.xtc/.trr）

### 命令

```bash
echo -e "C-alpha\n" | gmx mdmat -s md.tpr -f md.xtc -mean rdcm.xpm
```
- **目的**：计算残基对之间的平均距离矩阵
- **参数说明**：
  - `-s md.tpr`：输入拓扑文件
  - `-f md.xtc`：输入轨迹文件
  - `-mean rdcm.xpm`：输出平均距离矩阵（XPM 格式）
- **管道输入**：`"C-alpha"` 选择计算距离的原子组
- **选择原因**：C-alpha 代表残基位置，计算残基间距离
- **原理**：对轨迹中每一帧计算残基对距离，然后对所有帧取平均

### 输出
- **rdcm.xpm**：平均距离矩阵（单位：nm）
  - 矩阵元素 (i,j) = 残基 i 和残基 j 之间的平均距离
  - 对角线附近为近距离（序列相邻残基）

### 可视化
```bash
dit xpm_show -f rdcm.xpm -o rdcm.png
```

## 结果解释

### 距离标度（颜色表示）

| 距离 (nm) | 含义 |
|----------|------|
| < 0.5 | 紧密接触，空间相邻 |
| 0.5-1.0 | 可能短程相互作用 |
| 1.0-2.0 | 无直接接触 |
| > 2.0 | 远距离，无相互作用 |

### 结构特征
- **对角线**：短距离，序列相邻残基
- **非对角线模式**：长程接触表示三级结构
- **结构域接触**：非对角线块表示结构域间相互作用

### 二级结构模式
- **α-螺旋**：特征对角线模式，残基 i 与 i+3, i+4 距离近
- **β-折叠**：β-链对在非对角线区域显示近距离

## 常见问题

**Q: 距离矩阵显示均匀高距离？**  
A: 可能蛋白质展开或原子选择不正确。

**Q: 距离矩阵显示噪声？**  
A: 增加模拟时间获得更好的平均，或检查轨迹质量。

**Q: 意外的距离模式？**  
A: 验证蛋白质序列和结构，检查模拟稳定性。