解析TMalign文本文件中的转换矩阵
TM-align
将两个蛋白质结构通过旋转和位移对齐后:
TMalign test1.pdb test2.pdb -m mtx.txt
输出转换矩阵,文件内容为:
------ The rotation matrix to rotate Chain_1 to Chain_2 ------
m t[m] u[m][0] u[m][1] u[m][2]
0 0.7438770778 0.9928096071 0.1017845140 0.0629999746
1 3.6134235067 0.0559138503 0.0710345055 -0.9959054877
2 18.3578484222 -0.1058429281 0.9922671070 0.0648325754
Code for rotating Structure A from (x,y,z) to (X,Y,Z):
for(i=0; i<L; i++)
{
X[i] = t[0] + u[0][0]*x[i] + u[0][1]*y[i] + u[0][2]*z[i];
Y[i] = t[1] + u[1][0]*x[i] + u[1][1]*y[i] + u[1][2]*z[i];
Z[i] = t[2] + u[2][0]*x[i] + u[2][1]*y[i] + u[2][2]*z[i];
}
- 旋转矩阵
u
是一个 3x3 的矩阵,用于定义结构 A 到结构 B 的旋转。 - 平移向量
t
是一个长度为 3 的向量,