Pyrosetta 使用指南：蛋白质结构预测和相互作用分析

Pyrosetta 是一款功能强大的蛋白质结构预测软件，主要应用于蛋白质分子建模、蛋白质-蛋白质/配体相互作用预测等生物信息学领域。本指南将带您了解 Pyrosetta 的安装方法、基本用法和示例代码。

安装 Pyrosetta

在安装 Pyrosetta 之前，确保您的电脑已安装 Anaconda 或 Miniconda。在命令行中输入以下命令安装 Pyrosetta：

conda install pyrosetta -c schrodinger

使用 Pyrosetta

使用 Pyrosetta 时，首先需要导入 Pyrosetta 模块。在 Python 脚本中输入以下代码：

import pyrosetta
pyrosetta.init()

这将加载 Pyrosetta 库，使您可以使用 Pyrosetta 提供的各种功能。

蛋白质结构建模

使用 Pyrosetta 进行蛋白质结构建模需要一个蛋白质序列和一个模板蛋白质的结构。以下是一个简单的蛋白质结构建模代码示例：

from pyrosetta import pose_from_sequence, pose_from_pdb, create_score_function
from pyrosetta.rosetta.core.kinematics import FoldTree

# 从序列生成 pose 对象
sequence = 'MKQHKAMIVALIVICITAVVAALVTRKDLCEVHIRTGQTEVAVF'
pose = pose_from_sequence(sequence)

# 从 PDB 文件读取模板蛋白质结构
template_pose = pose_from_pdb('template.pdb')

# 创建得分函数
score_function = create_score_function('ref2015')

# 设置 FoldTree
fold_tree = FoldTree()
fold_tree.add_edge(1, pose.total_residue(), -1)
pose.fold_tree(fold_tree)

# 将模板蛋白质的结构复制到 pose 对象中
for i in range(1, template_pose.total_residue() + 1):
    pose.set_phi(i, template_pose.phi(i))
    pose.set_psi(i, template_pose.psi(i))
    pose.set_omega(i, template_pose.omega(i))

# 运行 relax 程序来优化结构
relax = pyrosetta.rosetta.protocols.relax.FastRelax(score_function)
relax.apply(pose)

蛋白质-蛋白质/配体相互作用预测

使用 Pyrosetta 进行蛋白质-蛋白质/配体相互作用预测需要一个蛋白质-蛋白质/配体复合体的结构。以下是一个简单的蛋白质-蛋白质/配体相互作用预测代码示例：

from pyrosetta import pose_from_pdb, create_score_function
from pyrosetta.rosetta.core.select.residue_selector import ChainSelector
from pyrosetta.rosetta.protocols.docking import setup_foldtree, setup_default_minimizer, \
    setup_packing_rotamers, DockMCMProtocol

# 从 PDB 文件读取蛋白质-蛋白质/配体复合体结构
protein_pose = pose_from_pdb('protein.pdb')

# 创建得分函数
score_function = create_score_function('ref2015')

# 选择蛋白质和配体的链
protein_selector = ChainSelector('A')
ligand_selector = ChainSelector('B')

# 设置 FoldTree
fold_tree = setup_foldtree(protein_pose, 'A_B')
protein_pose.fold_tree(fold_tree)

# 进行蛋白质-蛋白质/配体相互作用预测
dock_protocol = DockMCMProtocol()
dock_protocol.set_scorefxn(score_function)
dock_protocol.set_protein_selector(protein_selector)
dock_protocol.set_ligand_selector(ligand_selector)
setup_default_minimizer(dock_protocol)
setup_packing_rotamers(dock_protocol)
dock_protocol.apply(protein_pose)

以上是 Pyrosetta 的简单使用方法，您可以根据自己的需求进行进一步的学习和探索。

资源和文档

Pyrosetta 官方网站: https://www.pyrosetta.org/
Pyrosetta 文档: https://www.pyrosetta.org/docs/
Pyrosetta 论坛: https://www.pyrosetta.org/forum/

如果您在使用 Pyrosetta 时遇到问题，请随时在论坛中提问。

希望本指南对您有所帮助！