想起什么来写什么吧。
整体流程(以Ceas2, TPP, G3P为例):
准备蛋白即配体参数文件(pdb文件需要有header,header的顺序符合cst block的顺序,且residue1和residue2的顺序也必须符合cst block);
设置CST文件;
准备protocol和flag文件;
运行enzyme_design;
结果处理。
CST文件准备:
要想准备好cst文件,需要的先决条件有两个,一个是对CST文件格式的了解,一个是对rosetta参数文件(.params)的了解。
首先,CST格式可见另外一篇博客https://www.cnblogs.com/wq242424/p/9463703.html;
其次,.params文件格式里面,主要注意我们感兴趣的原子,比如HIS的Nhis,确定A1,A2和A3,Nhis作为A1,A2和A3则由.params文件里对应ICOOR行的parent原子决定,A2是A1的parent,A3是A2的parent。如HIS里Nhis对应的原子及A1是ND1,A2为A1的parent为CG,A3为A2的parent为CB,.params文件格式参见.param文件格式,标准氨基酸的.params文件存放在$ROSETTAHOME/main/database/chemical/residue_type_sets/下。
最后,设置好CST文件后,一定要用 CstfileToTheozymePDB.linuxgccrelease 对CST进行验证,因为有时候随着标准氨基酸.params文件的更新,不同原子的定义可能会有改变,命令用法如下:
CstfileToTheozymePDB.linuxgccrelease -extra_res_fa G3P.params -match:geometric_constraint_file G3P.cst
命令会产生出 PDB_Model_G3P.cst_1.pdb 文件,删除pdb文件里的CONNECT行后,使用 pymol PDB_Model_G3P.cst_1.pdb 命令查看motif是否符合设定。
下面列举我遇到过的几种氨基酸的原子信息:
侧链信息:
骨架信息:
上述是列举的官方文档种的残基原子信息。
下面是自己应用过程种用到的:
Ntrp,和Nhis,
OH ,(注意逗号前有空格)
Narg,
结果处理:
对设计出来的众多结果进行筛选,筛选规则的原文如下:
One approach currently used in the Baker group is the following: first, a subset of the - most important criteria is picked, i.e. total_score, ligand binding energy/SR_interface_E_1_2, total constraint score of the catalytic residues (all_cst), packstat, and buried unsatisfied polars of the ligand. Then, for each of these criteria, a minimum value is decided, which all designs considered for expression have to exceed ( i.e. total_score has to be lower than the corresponding Rosetta score of the undesigned scaffold, ligand_binding energy has to be < -10.0, and all_cst has to be < 1.0 ).
大体过程是:
首先,确定4-5个重要的对设计重要的特征,如total_score, ligand binding energy/SR_interface_E_1_2, total constraint score of the catalytic residues (all_cst), packstat, buried unsatisfied polars of the ligand等;
然后,确定每一个特征标准或者称为临界值,例如ligand_binding energy has to be < -10.0,又如 all_cst has to be < 1.0等;
最后,使用 DesignSelect.pl 脚本筛选出符合标准的设计。
DesignSelect.pl脚本的用法:
首先,需要准备一个标准或者临界值文件,如下:
req all_cst value < 1.0
req SR_4_interf_E_1_2 value < -10.0
output sortmin total_score
然后,使用命令:
DesignSelect.pl -d design.out -c <requirements file> -tag_column last > filtered_designs.out
对结果处理,得到的文件filtered_designs.out里包含符合我们条件的entry。
故障报错:
1. ERROR: unknown atom_name: TPP C
准备好flag文件,xml文件等后,运行 rosetta_scripts.linuxgccrelease @flags :
出现错误:ERROR: unknown atom_name: TPP C
没有排查出错误出处,由低版本(2016)换高版本rosetta_scripts(2018)运行后,问题解决。