Rosetta Relax教程：快速优化蛋白质结构

本教程将指导您使用Rosetta软件中的relax工具对蛋白质结构进行最简单的优化。

步骤

安装Rosetta软件并设置环境变量。
准备蛋白质结构PDB文件，保存为example1.pdb。
打开终端并输入以下命令：

rosetta_scripts.default.linuxgccrelease -s example1.pdb -parser:protocol <(echo '<ROSETTASCRIPTS>
    <SCOREFXNS>
        <ScoreFunction name='score3' weights='score3'/>
    </SCOREFXNS>
    <RESIDUE_SELECTORS>
        <Chain name='chainA' chains='A'/>
    </RESIDUE_SELECTORS>
    <TASKOPERATIONS>
        <RestrictToRepacking name='repack_only' />
    </TASKOPERATIONS>
    <MOVERS>
        <FastRelax name='relax' scorefxn='score3' task_operations='repack_only' />
    </MOVERS>
    <PROTOCOLS>
        <Add mover_name='relax'/>
    </PROTOCOLS>
</ROSETTASCRIPTS>') -out:pdb -out:prefix example1_relaxed

参数说明：

-s: 指定输入的PDB文件。
-out:pdb: 指定输出PDB文件。
-out:prefix: 指定输出文件名前缀。

等待程序运行完成。 程序运行结束后，您将在当前目录下找到名为 example1_relaxed.pdb 的文件，这就是经过relax后的蛋白质结构。

代码解析

以上命令使用Rosetta脚本语言定义了一个简单的relax流程，主要包括以下几个部分：

SCOREFXNS: 定义了使用的评分函数，这里使用默认的score3函数。
RESIDUE_SELECTORS: 定义了需要进行操作的残基，这里选择了所有属于A链的残基。
TASKOPERATIONS: 定义了对选中残基的操作，这里使用了RestrictToRepacking，表示只进行侧链堆积优化。
MOVERS: 定义了进行结构优化的mover，这里使用了FastRelax，并指定了评分函数和task operation。
PROTOCOLS: 定义了执行mover的顺序，这里只执行一次FastRelax。

总结

本教程介绍了使用Rosetta进行蛋白质结构relax的基本流程。您可以根据自己的需求修改脚本中的参数，例如选择不同的评分函数、task operation和mover等，以实现不同的优化目标。