时间不充分，简单记录下自己实践过程中的做法:

1. 首先，非标准残基都需要转换成.params文件，使用 <path-to-Rosetta>/main/source/scripts/python/public/molfile_to_params.py -n TPP TPP.mol2 --chain=F --clobber 转换，此命令会产生TPP_0001.pdb和TPP.params两个文件；

2. 其次，如果配体需要进行构象搜索，可以使用obabel，这个免费且简单，见博客http://www.cnblogs.com/wq242424/p/8231600.html，需要做限制性构象搜索的话，就比较难了，我用的是schrodinger的conformational search里的高级搜索；

3. 如果经过上一步产生了ligand_conformers.sdf，就需要为 <path-to-Rosetta>/main/source/scripts/python/public/molfile_to_params.py -n GAD -p GAD --chain=X --center=-31,-35.5,8 --clobber --conformers-in-one-file GAD_confs.sdf 命令多加一个参数 --conformers-in-one-file GAD_confs.sdf ，这个参数会在产生的.params文件末尾加一句 PDB_ROTAMERS GAD_conformers.pdb 来表明你的多重构象的文件，此命令会产生GAD.pdb，GAD.params和GAD_conformers.pdb三个文件，其中GAD.pdb是以conformers中第一个构象为结构；

4. 多底物对接的问题，经过查阅手册和论坛，得出的结论是通过控制需要dock的链("X/Y/F/。。。")来指定需要对接的小分子，比如我要对接TPP和GAD，就需要在single movers后面及前面LIGAND_AREAS和INTERFACE_BUILDERS那里指定分别对“F”和“X”链进行dock，并在MOVERS处指定两底物分子分别对应的操作，options不必特意根据双底物改变，下面给出了我做的两个XML文件的示例；

5. StartFrom,Transform,Translate和Rotate的差别和比较。StartFrom是在底物不在蛋白口袋内部，需要程序先将底物以Coordinates为基准放到口袋中的程序；Transform是使用蒙特卡洛的方法，对对接过程的底物和蛋白小分子随机取样的过程；Translate和Rotate则是对处于口袋中的小分子进行平移旋转的操作。底物在蛋白口袋中时，就可以不使用StartFrom；使用Transform时，则需要指定GRID。

6. 至于下面的high_resolution_dock XML文件，用途是在你已经知道你的底物的位置的情况下，对底物及周边残基进行高分辨率的对接，与其说对接，可能说packing和minimizing更恰当，因为它是对周围的残基优化的同时变换底物的构象，以求找到最优结果。

分享一下自己的两个xml文件及对应的flags：

low_resolution_dock+high_resolution_dock

This protocol will simply do low-resolution followed by high-resolution docking.

It will also report the binding energy (ddg) and buried-surface area (sasa) in the score file.

<ROSETTASCRIPTS>

	<SCOREFXNS>

		<ligand_soft_rep weights="ligand_soft_rep">

		</ligand_soft_rep>

		<hard_rep weights="ligand">

		</hard_rep>

	</SCOREFXNS>

	<LIGAND_AREAS>

		<docking_sidechain_X chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true" minimize_ligand="10"/>

		<final_sidechain_X chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true"/>

		<final_backbone_X chain="X" cutoff="7.0" add_nbr_radius="false" all_atom_mode="true" Calpha_restraints="0.3"/>

	</LIGAND_AREAS>

	<INTERFACE_BUILDERS>

		<side_chain_for_docking ligand_areas="docking_sidechain_X"/>

		<side_chain_for_final ligand_areas="final_sidechain_X"/>

		<backbone ligand_areas="final_backbone_X" extension_window="3"/>

	</INTERFACE_BUILDERS>

	<MOVEMAP_BUILDERS>

		<docking sc_interface="side_chain_for_docking" minimize_water="true"/>

		<final sc_interface="side_chain_for_final" bb_interface="backbone" minimize_water="true"/>

	</MOVEMAP_BUILDERS>

	<MOVERS>

	single movers

		<StartFrom name="start_from_X" chain="X">

			<Coordinates x="-31" y="-35.5" z="8"/>

		</StartFrom>

		<CompoundTranslate name="compound_translate" randomize_order="false" allow_overlap="false">

			<Translate chain="X" distribution="uniform" angstroms="2.0" cycles="50"/>

		</CompoundTranslate>

		<Rotate name="rotate_X" chain="X" distribution="uniform" degrees="360" cycles="500"/>

		<SlideTogether name="slide_together" chains="X"/>

		<HighResDocker name="high_res_docker" cycles="6" repack_every_Nth="3" scorefxn="ligand_soft_rep" movemap_builder="docking"/>

		<FinalMinimizer name="final" scorefxn="hard_rep" movemap_builder="final"/>

		<InterfaceScoreCalculator name="add_scores" chains="X" scorefxn="hard_rep"/>

	compound movers

		<ParsedProtocol name="low_res_dock">

			<Add mover_name="start_from_X"/>

			<Add mover_name="compound_translate"/>

			<Add mover_name="rotate_X"/>

			<Add mover_name="slide_together"/>

		</ParsedProtocol>

		<ParsedProtocol name="high_res_dock">

			<Add mover_name="high_res_docker"/>

			<Add mover_name="final"/>

		</ParsedProtocol>

	</MOVERS>

	<PROTOCOLS>

		<Add mover_name="low_res_dock"/>

		<Add mover_name="high_res_dock"/>

		<Add mover_name="add_scores"/>

	</PROTOCOLS>

</ROSETTASCRIPTS>

对应的flag options文件：

-in:file:s inputs_GAD/GALS_TPP_GAD.pdb

-in:file:extra_res_fa inputs_GAD/TPP.params inputs_GAD/GAD.params

-packing

	-ex1

	-ex2aro

	-ex2

	-no_optH false

	-flip_HNQ true

	-ignore_ligand_chi true

-parser

	-protocol inputs_GAD/ligand_dock.xml

-out

	-path:all outputs_GAD

	-nstruct 1000

	-overwrite

high_resolution_dock_only

This protocol will simply do high-resolution docking.

It will also report the binding energy (ddg) and buried-surface area (sasa) in the score file.

<ROSETTASCRIPTS>

	<SCOREFXNS>

		<ligand_soft_rep weights="ligand_soft_rep">

		</ligand_soft_rep>

		<hard_rep weights="ligand">

		</hard_rep>

	</SCOREFXNS>

	<LIGAND_AREAS>

		<docking_sidechain_X chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true" minimize_ligand="10"/>

		<final_sidechain_X chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true"/>

		<final_backbone_X chain="X" cutoff="7.0" add_nbr_radius="false" all_atom_mode="true" Calpha_restraints="0.3"/>

	</LIGAND_AREAS>

	<INTERFACE_BUILDERS>

		<side_chain_for_docking ligand_areas="docking_sidechain_X"/>

		<side_chain_for_final ligand_areas="final_sidechain_X"/>

		<backbone ligand_areas="final_backbone_X" extension_window="3"/>

	</INTERFACE_BUILDERS>

	<MOVEMAP_BUILDERS>

		<docking sc_interface="side_chain_for_docking" minimize_water="true"/>

		<final sc_interface="side_chain_for_final" bb_interface="backbone" minimize_water="true"/>

	</MOVEMAP_BUILDERS>

	<SCORINGGRIDS ligand_chain="X" width="15">

		<classic grid_type="ClassicGrid" weight="1.0"/>

	</SCORINGGRIDS>

	<MOVERS>

		<Transform name="transform" chain="X"box_size="7.0" move_distance="0.2" angle="20" cycles="500" repeats="1" temperature="5"/>

		<HighResDocker name="high_res_docker" cycles="6" repack_every_Nth="3" scorefxn="ligand_soft_rep" movemap_builder="docking"/>

		<FinalMinimizer name="final" scorefxn="hard_rep" movemap_builder="final"/>

		<InterfaceScoreCalculator name="add_scores" chains="X" scorefxn="hard_rep"/>

	</MOVERS>

	<PROTOCOLS>

		<Add mover_name="transform"/>

		<Add mover_name="high_res_docker"/>

		<Add mover_name="final"/>

		<Add mover_name="add_scores"/>

	</PROTOCOLS>

</ROSETTASCRIPTS>

对应的flag_options文件：

-in:file:s inputs_IMA/GALS_IMA.pdb

-in:file:extra_res_fa inputs_IMA/IMA.params

-packing

	-ex1

	-ex2aro

	-ex2

	-no_optH false

	-flip_HNQ true

	-ignore_ligand_chi true

-parser

	-protocol inputs_IMA/ligand_dock.xml

-out

	-path:all outputs_IMA

	-nstruct 1000

	-overwrite

多底物对接（两个小分子的链分别对应“X”和“F”）

1. 手动指定两个小分子（“X”，“F”）的位置（使用StartFrom mover指定，对应的Coordinates）
This protocol will simply do low-resolution followed by high-resolution docking.

It will also report the binding energy (ddg) and buried-surface area (sasa) in the score file.

<ROSETTASCRIPTS>

	<SCOREFXNS>

		<ligand_soft_rep weights="ligand_soft_rep">

		</ligand_soft_rep>

		<hard_rep weights="ligand">

		</hard_rep>

	</SCOREFXNS>

	<LIGAND_AREAS>

		<docking_sidechain_X chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true" minimize_ligand="10"/>

		<final_sidechain_X chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true"/>

		<final_backbone_X chain="X" cutoff="7.0" add_nbr_radius="false" all_atom_mode="true" Calpha_restraints="0.3"/>

		<docking_sidechain_F chain="F" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true" minimize_ligand="10"/>

		<final_sidechain_F chain="F" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true"/>

		<final_backbone_F chain="F" cutoff="7.0" add_nbr_radius="false" all_atom_mode="true" Calpha_restraints="0.3"/>

	</LIGAND_AREAS>

	<INTERFACE_BUILDERS>

		<side_chain_for_docking ligand_areas="docking_sidechain_X,docking_sidechain_F"/>

		<side_chain_for_final ligand_areas="final_sidechain_X,final_sidechain_F"/>

		<backbone ligand_areas="final_backbone_X,final_backbone_F" extension_window="3"/>

	</INTERFACE_BUILDERS>

	<MOVEMAP_BUILDERS>

		<docking sc_interface="side_chain_for_docking" minimize_water="true"/>

		<final sc_interface="side_chain_for_final" bb_interface="backbone" minimize_water="true"/>

	</MOVEMAP_BUILDERS>

	<MOVERS>

	single movers_X

		<StartFrom name="start_from_X" chain="X">

			<Coordinates x="-31" y="-35.5" z="8"/>

		</StartFrom>

		<StartFrom name="start_from_F" chain="F">

			<Coordinates x="-27.26" y="-34.56" z="4.87"/>

		</StartFrom>

		<CompoundTranslate name="compound_translate" randomize_order="false" allow_overlap="false">

			<Translate chain="X" distribution="uniform" angstroms="2.0" cycles="50"/>

			<Translate chain="F" distribution="uniform" angstroms="2.0" cycles="50"/>

		</CompoundTranslate>

		<Rotate name="rotate_X" chain="X" distribution="uniform" degrees="360" cycles="500"/>

		<Rotate name="rotate_F" chain="F" distribution="uniform" degrees="360" cycles="500"/>

		<SlideTogether name="slide_together" chains="X,F"/>

		<HighResDocker name="high_res_docker" cycles="6" repack_every_Nth="3" scorefxn="ligand_soft_rep" movemap_builder="docking"/>

		<FinalMinimizer name="final" scorefxn="hard_rep" movemap_builder="final"/>

		<InterfaceScoreCalculator name="add_scores" chains="X,F" scorefxn="hard_rep"/>

	compound movers

		<ParsedProtocol name="low_res_dock">

			<Add mover_name="start_from_X"/>

			<Add mover_name="start_from_F"/>

			<Add mover_name="compound_translate"/>

			<Add mover_name="rotate_X"/>

			<Add mover_name="rotate_F"/>

			<Add mover_name="slide_together"/>

		</ParsedProtocol>

		<ParsedProtocol name="high_res_dock">

			<Add mover_name="high_res_docker"/>

			<Add mover_name="final"/>

		</ParsedProtocol>

	</MOVERS>

	<PROTOCOLS>

		<Add mover_name="low_res_dock"/>

		<Add mover_name="high_res_dock"/>

		<Add mover_name="add_scores"/>

	</PROTOCOLS>

</ROSETTASCRIPTS> 
2. 其中一个小分子（“F”）位置已知，另一个小分子（“X”）的位置需要指定（首先需使用SCORINGGRIDS为已知位置的小分子指定盒子，然后使用Transform指定其动作模式）
This protocol will simply do low-resolution followed by high-resolution docking.

It will also report the binding energy (ddg) and buried-surface area (sasa) in the score file.

<ROSETTASCRIPTS>

	<SCOREFXNS>

		<ligand_soft_rep weights="ligand_soft_rep">

		</ligand_soft_rep>

		<hard_rep weights="ligand">

		</hard_rep>

	</SCOREFXNS>

	<LIGAND_AREAS>

		<docking_sidechain_X chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true" minimize_ligand="10"/>

		<final_sidechain_X chain="X" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true"/>

		<final_backbone_X chain="X" cutoff="7.0" add_nbr_radius="false" all_atom_mode="true" Calpha_restraints="0.3"/>

		<docking_sidechain_F chain="F" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true" minimize_ligand="10"/>

		<final_sidechain_F chain="F" cutoff="6.0" add_nbr_radius="true" all_atom_mode="true"/>

		<final_backbone_F chain="F" cutoff="7.0" add_nbr_radius="false" all_atom_mode="true" Calpha_restraints="0.3"/>

	</LIGAND_AREAS>

	<INTERFACE_BUILDERS>

		<side_chain_for_docking ligand_areas="docking_sidechain_X,docking_sidechain_F"/>

		<side_chain_for_final ligand_areas="final_sidechain_X,final_sidechain_F"/>

		<backbone ligand_areas="final_backbone_X,final_backbone_F" extension_window="3"/>

	</INTERFACE_BUILDERS>

	<MOVEMAP_BUILDERS>

		<docking sc_interface="side_chain_for_docking" minimize_water="true"/>

		<final sc_interface="side_chain_for_final" bb_interface="backbone" minimize_water="true"/>

	</MOVEMAP_BUILDERS>

	<SCORINGGRIDS ligand_chain="F" width="15">

		<classic grid_type="ClassicGrid" weight="1.0"/>

	</SCORINGGRIDS>

	<MOVERS>

	single movers_X

		<StartFrom name="start_from_X" chain="X">

			<Coordinates x="-31" y="-35.5" z="8"/>

		</StartFrom>

		<Transform name="transform_F" chain="F" box_size="7.0" move_distance="0.2" angle="20" cycles="500" repeats="1" temperature="5"/>

		<CompoundTranslate name="compound_translate" randomize_order="false" allow_overlap="false">

			<Translate chain="X" distribution="uniform" angstroms="2.0" cycles="50"/>

			<Translate chain="F" distribution="uniform" angstroms="2.0" cycles="50"/>

		</CompoundTranslate>

		<Rotate name="rotate_X" chain="X" distribution="uniform" degrees="360" cycles="500"/>

		<Rotate name="rotate_F" chain="F" distribution="uniform" degrees="360" cycles="500"/>

		<SlideTogether name="slide_together" chains="X,F"/>

		<HighResDocker name="high_res_docker" cycles="6" repack_every_Nth="3" scorefxn="ligand_soft_rep" movemap_builder="docking"/>

		<FinalMinimizer name="final" scorefxn="hard_rep" movemap_builder="final"/>

		<InterfaceScoreCalculator name="add_scores" chains="X,F" scorefxn="hard_rep"/>

	compound movers

		<ParsedProtocol name="low_res_dock">

			<Add mover_name="start_from_X"/>

			<Add mover_name="transform_F"/>

			<Add mover_name="compound_translate"/>

			<Add mover_name="rotate_X"/>

			<Add mover_name="rotate_F"/>

			<Add mover_name="slide_together"/>

		</ParsedProtocol>

		<ParsedProtocol name="high_res_dock">

			<Add mover_name="high_res_docker"/>

			<Add mover_name="final"/>

		</ParsedProtocol>

	</MOVERS>

	<PROTOCOLS>

		<Add mover_name="low_res_dock"/>

		<Add mover_name="high_res_dock"/>

		<Add mover_name="add_scores"/>

	</PROTOCOLS>

</ROSETTASCRIPTS>

　　

ROSETTA使用技巧随笔--RosettaLigand Docking的更多相关文章

  1. ROSETTA使用技巧随笔--蛋白蛋白对接

    先写简略版,以后再详细写. 1. 对输入结构进行预处理(refine) $> relax.default.linuxgccrelease -in:file:s input_files/from_ ...

  2. ROSETTA使用技巧随笔--relax使用

    Purpose: 主要说目的,relax的作用就是对一个给定的蛋白进行构象搜索,寻找与WT相似并能量低于WT的结构,既包含packer又包含minimizer.主要的应用在对一个结构构象进行取样,获得 ...

  3. ROSETTA使用技巧随笔--PyMOL实时观测ROSETTA模拟过程中的结构变化

    没有梦想的人,就是一只咸鱼,像我,就有一个梦想,就是让蛋白模拟过程变成动画,动起来! 虽然MD中有很多方法可以方模拟过程像动画一样播放出来,但是我一直想在ROSETTA中也找一个这样的功能,这不,我发 ...

  4. ROSETTA使用技巧随笔--控制Log输出等级

    一般运行ROSETTA,屏幕上的Log很多,而且很复杂,让我们看着眼晕,现在我们可以通过控制Log等级来控制屏幕上输出的东西.  Integer Level 0 Fatal 100 Error 200 ...

  5. ROSETTA使用技巧随笔--Full Atom Representation和Centroid Representation

    Full Atom Representation vs Centroid Representation Full Atom Representation即全原子标识,氨基酸残基的所有相关原子,均原封不 ...

  6. ROSETTA使用技巧随笔--score.sc处理

    对score.sc的便利操作: $ sort -n -k2 example_score_file.sc 此命令会以score文件的第二列进行排序(sort -n -k2,-n表示对数值型数据排序,可用 ...

  7. RF常用技巧随笔

    合并xml输出文件: pybot --outputdir d:\test --log log1.html --report report1.html --output output1.xml -T d ...

  8. angular小技巧随笔

    1. 重新刷新页面 同页面切换状态: $state.go('tab.index', {inviteId:inviteId}); self.location.reload();

  9. JavaScript小技巧随笔

    1. 在用||做条件判断时,如果情况较多,我们可以考虑是否能够用Array.includes()方法代替 var conditionArray = [ test1, test2, test3 ]; i ...

随机推荐

  1. java整形数据和浮点型数据

    在定义数据的时候,如果有不是整形的数据(单精度,双精度之类的),那么运算的时候要在其所在式子里的数据后面加上.0,否则会出现(1/2=0.1.0/2=0.5)的现象,比如如下代码 public cla ...

  2. POJ 1102 - LC-Display

    Description A friend of you has just bought a new computer. Until now, the most powerful computer he ...

  3. 文本分类-TensorRT优化结果对比图

    做的文本二分类,使用tensorRT进行图优化和加速,输出预测概率结果对比如下: 从结果对比来看,概率值有微小的变化,但不影响最终的分类

  4. hash_map

    点开一道第是自己oj的第440大关,想a了,一直想却无果,学长一句点醒,开始写hash. 关于这道题呢很无语了,两天卡在这上面,而且有些dalao不到20min就a了.我太菜了. 所以要深入讨论这道题 ...

  5. DbGridEh 一个单元格的值改变时另一单元格的值随之改变

    你可以为每个字段设置OnSetText事件,这样在输入完后回车会移动时就会触发,或者在adoquery的beforepost中或afterpost中都可以grid也提供了一些事件,也可以在某些条件下做 ...

  6. 转：Java 异常结构体系

    原文地址:Java 异常结构体系 保存一份资料 前几天在参加网易和360公司的在线考试的时候,都出了一道关于java中异常类的多项选择题.这几天翻看了相关书籍和网上一些资料,结合自己的理解与思考,将自 ...

  7. Linux svn服务器搭建

    一.安装 yum -y install subversion 二.建立版本库目录 mkdir /var/svn/svnrepos 三.创建svn版本库 svnadmin create /var/svn ...

  8. 浏览器数据库 IndexedDB 入门教程

    一.概述 随着浏览器的功能不断增强,越来越多的网站开始考虑,将大量数据储存在客户端,这样可以减少从服务器获取数据,直接从本地获取数据. 现有的浏览器数据储存方案,都不适合储存大量数据:Cookie 的 ...

  9. java 线程 (一) Thread

    package cn.sasa.demo1; public class Test { public static void main(String[] args) throws Interrupted ...

  10. cd 命令

    [root@localhost ~]# cd # 进入当前用户的家目录 [root@localhost ~]# cd ~ # 进入当前用户的家目录 [root@localhost ~]# cd /da ...