计算每两条链Biopython之间的氢键

Question

我试图写一个Python函数，它将计算结构中每两个链之间的氢键，并使用Biopython模块从pdb文件返回原子坐标，原子名和残留数表。

限制是：

1. N和O之间的距离<= 3.5。
2. pdb文件中每两个链之间有大约47个Hbonds，因此我们从链A获取ALA / O，从链B获取LEU / N，下一个键将在LEU / O链B和ALA / N之间。链A然后我们再次交换（定期）

这是函数，但它只进行一次迭代（因为中断），结果只写入条件的第二个字符串的输出。如何指定条件或者还有其他我不做的事情？

router.get("/myjobs/:user", function (req, res) { 
   jobs.find({'author.username':user},function (err, myjobs) {
     ...
   }
}

我不是程序员，需要这样才能完成我的论文。非常感谢您的建议和帮助。