如何使用 PCcli 命令生成数百个肽坐标?

How to use PCcli command to generate hundreds of peptide coordinates?

提问人:Yuval 提问时间:10/5/2023 最后编辑:pippo1980Yuval 更新时间:10/6/2023 访问量:25

问:

有一个 PCcli 命令 [https://pypi.org/project/PeptideConstructor/] 可以构建给定肽序列的.pdb文件,但是我不明白如何将其用于数百个肽序列。它对单肽的用法是这样的:

pip install PeptideConstructor

PCcli -s AaDdKSQym -o output.pdb

这里,-s 标志后面的字符是肽序列,-o 标志用于命名输出 .pdb 文件。我在 .csv 文件中有一个肽序列列表,并希望使用 PCcli 命令生成这些序列的 .pdb。

假设我有 100 个序列,并用 阅读“序列”列。pd.read_csv

import pandas as pd

df = pd.read_csv("PeptideSequences.csv")

然后,使用for循环,我尝试将每个序列行放在-s标志之后,并以相同的方式为输出.pdb文件提供序列名称。

for i in range(100):

    PCcli -s str(df.loc[i,'Sequences']) -o str(df.loc[i,'Sequences']).pdb

但是,我收到了无效的语法错误。你能帮我使用这个命令吗?我也希望您能提出其他使用建议。例如,映射序列

pandas 字符串 映射 调用 biopython

评论

0赞 OCa 10/5/2023
欢迎来到 Stack Overflow。你能发布完整的错误吗?(作为文本)。我想知道在您的案例中是否有可能出现最小可重复的例子
0赞 pippo1980 10/6/2023
我仍然不明白单肽折叠是如何折叠的:从头折叠、同源建模、分子动力学 (MD) 模拟或基于深度学习的方法?

答:

0赞 mozway 10/5/2023 #1

您可以使用 subprocess.call

from subprocess import call

for seq in df['Sequences']:
    call(['PCcli', '-s', seq, '-o', f'{seq}.pdb'])