Math 如何计算蛋白质db文件中的键角？

如果我有三个蛋白质骨架的x，y，z坐标（N-Ca-C-N-Ca-C…）：

---

如何计算“键角”（Ni-Cai-Ci，Cai-Ci-Ni+1，Ci-Ni+1-Cai+1）？

基本向量几何。两个标准化向量的点积是它们之间角度的余弦

冲洗液：

N -14.152 0.961 4.712
CA -13.296 0.028 3.924
C -11.822 0.338 4.193

（N-Ca）=（-14.152 0.9614.712）-（13.296 0.028 3.924）=（-0.856,0.933,0.778）
标准化：（-0.576,0.628,0.524）

（C-Ca）=（-11.8220.3384.193）-（13.2960.0283.924）=（1.474,0.310,0.269）
归一化：（0.963,0.203,0.176）

点积：（-0.576,0.628,0.524）x（0.963,0.203,0.176）=-0.335

---

angle=acos（-0.335）=109.57度

这可能属于化学部分？如果这是一个编程问题，您可能需要指定您正在使用的语言/库。zero323：谢谢，我没有在数学堆栈交换上看到该网站。Lev Levitsky：谢谢，我正在使用PerlQuick搜索结果：，非常感谢，Beta版。我用PyMol表示了这些点，实际上109.5度是正确的。我最初有类似的点积输出（尽管是负数：0.333）。这给了我70的最终学位，这是错误的！我所做的不同之处在于N-Ca计算；我的x，y，z坐标是你的负坐标（0.856，-0.933，-0.778）。所以我猜向量减法，顺序很重要，是吗？这对我来说仍然很奇怪。。但再次感谢你。

N -14.152 0.961 4.712
CA -13.296 0.028 3.924
C -11.822 0.338 4.193