rosetta_ddg 使用-rosetta 2020版

参考https://www.rosettacommons.org/docs/latest/application_documentation/analysis/ddg-monomer
https://zhuanlan.zhihu.com/p/26071679

参考作者的进行了下操作，有些部分需要修改

先处理PDB

python ~/Downloads/rosetta_src_2020.50.61505_bundle/main/tools/protein_tools/scripts/clean_pdb.py protein.pdb A > cleanpdb.log

protein.pdb为自己的WT蛋白结构文件，A表示只用处理A链，处理好的蛋白结构文件为protein_A.pdb，

能量最小化

~/Downloads/rosetta_src_2020.50.61505_bundle/main/source/bin/minimize_with_cst.mpi.linuxgccrelease -in:file:s protein_A.pdb  -in:file:fullatom -ignore_unrecognized_res -fa_max_dis 9.0 -database ~/Downloads/rosetta_src_2020.50.61505_bundle/main/database/ -ddg::harmonic_ca_tether 0.5 -ddg::constraint_weight 1.0 -ddg::out_pdb_prefix min_cst_0.5 -ddg::sc_min_only false > mincst.log

获得距离限制文件

./convert_to_cst_file.sh mincst.log > input.cst

这里会出现input.cst文档为空的情况，是因为rosetta中的convert_to_cst_file.sh文档不太符合自己的，因此需要修改convert_to_cst_file.sh 文档，该文档目的就是提取mincst.log 匹配的字符重定向生产input.cst

这里为修改后converst_to_cst_file.sh的文档
cat mincst.log | grep ‘^apps.public.ddg.minimize_with_cst: (0) c-alpha’ | awk ‘{print "AtomPair CA “$8” CA “$10” HARMONIC “$12” "$15}’ > input.cst

首先就是匹配到apps.public.ddg.minimize_with_cst: (0) c-alpha字符，然后提取字符的第8、10、12、15列，重定向到位置限制input.cst文档中。

编写突变文件flags，这里一定要与已重新编号的pdb 文件相符,比如这里我需要突变两个点，分别为N13P/V193F，那么6d3u.mutfile 文件内容如下：
total 2
1
N 13 P
1
V 193 F

执行ddg_monomer程序

~/Downloads/rosetta_src_2020.50.61505_bundle/main/source/bin/ddg_monomer.mpi.linuxgccrelease -in:file:s min_cst_0.5.protein_A_0001.pdb @flags

这里flags 文件编写参考了rosetta tests 文件
-ddg:weight_file standard_plus_score12
-ddg::iterations 50
-ddg:dump_pdbs true
-ddg::mut_file ./6d3u.mutfile
-database /home/tust/Downloads/rosetta_src_2020.50.61505_bundle/main/database
-ignore_unrecognized_res
-ddg:local_opt_only true
-ddg:min_cst true
-ddg:suppress_checkpointing true
-in:file:fullatom
-ddg:mean false
-ddg:min true
-ddg:sc_min_only false
-ddg:ramp_repulsive true
-unmute core.optimization.LineMinimizer
-constraints:cst_file ./input.cst
-ddg:output_silent true
-ddg:opt_radius 12.0
-score:fa_max_dis 9.0

-in:file:s 为能量最小化后的文件
第一个打分文件使用了standard_plus_score12 ，所以在本目录下有该文档，文献中也有使用soft_rep_design
flags要修改的参数为mut_file文件，也就是编写的需要突变文件。

rosetta_ddg 使用-rosetta 2020版相关推荐

怎么将两个datatable 合并_视频音频合并大师免费版-视频音频合并大师2020版下载(附使用方法) v2020.02.08...
视频音频合并大师免费版是一款好用的视频合并软件,包括视频合并和音频合并,可以一次性将多个视频格式文件或音频文件合并成一个单独的视频文件,进行转换合并成你想要的视频,并且有多种自适应屏幕尺寸可以选择, ...
会议论文影响因子多少_中国学术期刊影响因子年报(2020版)发布｜中国心血管杂志影响力指数3年连升...
中国学术期刊影响因子年报(2020版)发布 <中国心血管杂志>影响力指数3年连续增长近日,中国知网与中国科学文献计量评价研究中心联合发布了2020年<中国学术期刊影响因子年报> ...
英特尔530和535哪个好_详细介绍2020版联想小新Air 14酷睿版和锐龙版怎么选，哪款好...
最近在后台收到很多朋友的咨询,问联想小新Air 14 2020酷睿版和联想小新Air 14 2020锐龙版哪个比较好,怎么选呢.下面我们一起从这两个版本的价格对比,配置参数对比方面了解下.供需要的朋友 ...
北大核刊最新版2020_爱尔兰地图（2020版）
友情提示:公众号首页输入关键字"分国地图"或查阅菜单"分地查 - 分国地图"均可以查阅其他国家. 自然地理:爱尔兰是岛国,位于欧洲爱尔兰岛中南部.西濒大西洋,东 ...
2021年东港二中高考成绩查询,辽宁省东港二中2021届英语高考模拟试卷1（新课改原创2020版，供2021年课改省份考生使用）...
ID:12674887 资源大小:288KB 资料简介: 说明:新高考变化:2021届考生取消听力,4个阅读每个2.5分共37.5分,七选五阅读每个2.5分共12.5分,完形填空由20个变成15个每个 ...
建筑专业规范大全 2020版_房屋建筑工程现行规范标准目录汇编(2020版)—建筑电气...
房屋建筑工程现行规范标准目录汇编(2020版) 建筑电气规范编号规范名称 GB 50034-2013 建筑照明设计标准 GB 50052-2009 供配电系统设计规范 GB 50053-2013 ...
javaweb不同用户需要几张表_程序员必备2020版：JavaWeb快速进阶全套教程
Java Web应用由一组Servlet.HTML页.类.以及其它可以被绑定的资源构成.它可以在各种供应商提供的实现Servlet规范的Servlet容器中运行. JavaWeb项目简单来说就是一个应 ...
dumpbin发现没有入口函数_JavaScript基础之入口函数-2020版
JavaScript基础之入口函数-2020版 1.入口函数 window.onload = function(){ 内部放js} 这个函数的意思就是说,当我们页面加载完毕之后(就是说等页面的结构样 ...
小马虎想用计算机计算396乘19,2020版苏教版数学四年级下册第四单元《用计算器计算》单元测试卷C卷...
2020版苏教版数学四年级下册第四单元<用计算器计算>单元测试卷C卷姓名:________ 班级:________ 成绩:________ 小朋友,带上你一段时间的学习成果,一起来做个自 ...
【报告分享】中国物联网平台产业市场研究报告（2020版）.pdf（附下载链接）...
今天给大家分享中国物联网产业应用联盟.深圳市物联网产业协会和深圳市物联传媒有限公司联合发布的报告<中国物联网平台产业市场研究报告(2020版).pdf>,报告共包含如下七大章节: 1.云平 ...

rosetta_ddg 使用-rosetta 2020版

rosetta_ddg 使用-rosetta 2020版相关推荐

最新文章

热门文章