modeller建模
off999 2024-11-24 20:12 29 浏览 0 评论
蛋白质建模方法 Modeller 和Swiss model
第一种方法: 最简便的方法,swiss_model建模
材料准备:一段已知的氨基酸序列:(本篇文章采用从modeller的基础教程里的TvLDH)
打开swiss model官网 (https://swissmodel.expasy.org/)
开始:一键建模 Builde Model
图1
下一步:将氨基酸序列粘贴到如下文本框
图2
第三步:点击 Build Model开始建模,swiss model 网站会将你的刚刚复制到文本框的上传到后台,后台建模器会自动选择一个最为相似的蛋白质序列,作为模板进行建模。这个过程大概需要5到10分钟
图3
第四步:得到结果(づ ●─● )づ,嘻嘻
不过大家看红框的地方序列相似度100%,说明这个模型已经被建立过来,而且这个被建立过的模型已经上传的世界蛋白质数据库了,如果你没有特殊的需要,这个模型(如图中4uum)的可靠性非常高了,你就不用再建了,当然这里我们仍然使用这段序列,因为它是modeller官网给的序列,为了方便大家按照modeller官网的教程步调一致。
图4
关于Modeller建模
前言: swiss model建模简单方便,可靠性高,但是它不能进行多模板建模,而且也不能修改其中你想优化的序列片段,那modeller就提供了另一种方法,灵活性高,可进行多模板建模,和序列片段优化,虽然操作相对有些繁琐,不过这对于理解建立蛋白质模型的思路和历程都很有帮助。
本篇教程和其他教程略有不同的地方:
1:采用在Pycharm直接运行脚本文件,而不同于如图5,在Modeller终端运行脚本。
2:在搜索模板的步骤,我们采用了Swiss model搜索模板的方法。
这两种操作都是官网提供的思路,也是modeller建模,脚本运用的灵活之处,也为了大家更方便的学习和操作。
图5
第二种:运用modeller建模
第一步:材料准备
材料准备:一段已知的氨基酸序列:(从modeller的基础教程里找 , 以官网给的TvDHL为例 )注意要修改为 .ali 的格式哦,如 TvLDH.ali ,直接重命名修改后缀名即可。
MSEAAHVLITGAAGQIGYILSHWIASGELYGDRQVYLHLLDIPPAMNRLTALTMELEDCAFPHLAGFVATTDPKA
AFKDIDCAFLVASMPLKPGQVRADLISSNSVIFKNTGEYLSKWAKPSVKVLVIGNPDNTNCEIAMLHAKNLKPEN
FSSLSMLDQNRAYYEVASKLGVDVKDVHDIIVWGNHGESMVADLTQATFTKEGKTQKVVDVLDHDYVFDTFFKKI
GHRAWDILEHRGFTSAASPTKAAIQHMKAWLFGTAPGEVLSMGIPVPEGNPYGIKPGVVFSFPCNVDKEGKIHVV
EGFKVNDWLREKLDFTEKDLFHEKEIALNHLAQGG
第二步:搜索模板
运行:搜寻模板:在Modeller通过用build_profile.py和compare.py两个python脚本,来筛选与目标序列相似的模板。耗时较长,而且对比的数据库pdb_95.pir有限,在这里我们采用swiss model网站的搜寻模板功能
图6
结果:找到模板,为了和官网教程保持一致,这里还是按照官网给的1bdmA模板
第三步:模板与目标序列对齐
运行:用align2d.py脚本使目标序列这里用的是TvDHL和模板对齐,模板用.pdb格式
结果:生成TvLDH-1bdmA.pap文件
第四步:建立模型
运行:用model-single.py进行建模
结果:TvLDH.B99990001.pdb
TvLDH.B99990002.pdb
TvLDH.B99990003.pdb
TvLDH.B99990004.pdb
TvLDH.B99990005.pdb
注意:可以根据自己的需要选择要建模的个数,同时建模个数越多,耗时也越长,但是建出好模型的概率也会越大。
筛选最优建模结果:从model-single.log筛选最优的模板,一般看molpdf分数最低的。
第五步:评测画图
运行:用evaluate_model.py脚本对建立的蛋白质模型进行评测
结果:得到TvLDH.profile文件
用gnuplot或者 plot_profiles.py画图
那在这里我们直接用 plot_profiles.py画图了, plot_profiles.py在官网给的basic example文件里有。
当然如果你习惯用gunplot软件画图,也可以达到同样的效果。如果你没有这个软件的话,这里提供一个下载链接https://jaist.dl.sourceforge.net/project/gnuplot/gnuplot/5.4.3/gp543-win64-mingw.exe。使用命令行语句,'plot "TvLDH.profile" using 1:42 with lines'
我们这里采用plot_profiles.py
图7
这个pylab模块在Matplotlib库下,所以要提前安装Matplotlib 在终端Terminal 执行命令 pip install matplotlib即可
图8
在这里要注意一个事项,evaluate.py第一次运行一次,那么我只能得到一个画图的数据文件TvLDH.profile,但是我们还需要模板Template所用的数据,也就是图中的绿线(Template)所用的数据,所以我们需要修改一下evaluate.py脚本,把下图中圈红线的修改一下,运行一下evaluate.py得到画模板线的数据1bdmA.profile文件。
如果大家有Python基础的话,理解起来还是很简单的,当然对这些脚本的简单理解和修改,也能很好的帮助大家理解modeller建模过程,如果大家感兴趣的话,Ctrl+左键,可以直接查看一些方法的源代码,进一步了解建模所运用的数学算法,如果不了解Python的话也没关系啦,按照我介绍的步骤也能得到结果,同时我也会上传一些步骤视频。
图9
修改为下图:
图10
最后一步:运行plot_profiles.py,生成dope_profile.png如下图,学会制作这个图很重要哦,接下来的多模板建模还要用到画图。
加油,各位!
写的不好的地方,多多见谅(づ ●─● )づ,有什么疑问的地方,欢迎评论留言,别忘了点赞哦
图11
相关推荐
- windows server 2008激活
-
CAD不注册激活的话,一般只能试用30天,而且每次打开都会弹出激活窗口。为了能够永久使用,我们需要一个注册机来帮我们生成一个激活码激活它。下载CAD2008注册机(注册机必须与所要激活的CAD版本相对...
- 邮件服务器(hmailserver搭建邮件服务器)
-
电子邮件服务器名称:比如添加的是网易邮箱帐号在“接收邮件(pop、IMAP或HTTP)服务器:”字段中输入pop.163.com。在“发送邮件服务器(SMTP):”字段中输入smtp.163.com...
- win7蓝屏0x0000000a(win7蓝屏代码0x000000f4)
-
电脑蓝屏代码0x0000000a的原因以及解决办法如下:1、在BIOS界面内,进入“Intogratedperipherals”选项里把“SATAConfiguration”项的值改成IDE。&...
- 手机看nwd格式的3d图软件(手机看3d图片)
-
NWD只有具备NavisWorksPublisher许可才能保存NWD文件。这种文件格式存储NWF文件格式存储的所有NavisWorks特定数据,外加模型的几何图形。NWD文件一般比原始的CAD文件...
- windows系统正版(win10正版系统)
-
WINDOWS激活了不一定是正版,可以使用以下方法验证:"开始"菜单——"运行"中输入:1、slmgr.vbs-dli显示:操作系统版本、部分产品密钥、许可证状态...
- 创新声卡驱动安装教程(创新声卡7.1驱动安装)
-
1、准备工作:准备好声卡驱动安装文件,根据自己计算机操作系统的不同而选择不同的版本2、双击安装文件,根据弹出的提示进行操作,进行安装,可以选择安装至默认位置,安装过程中会出现若干设置,最好按照制造厂商...
- win10进入安全模式黑屏(win10进安全模式黑屏什么都不显示)
-
正常情况下进入“安全模式”屏幕是黑色的,没有壁纸。并且,与显卡,显示器没有任何关系。正常模式比安全模式多加载了很多启动,视频模式有所改变。如果出现异常,则有可能是以下问题引起:1、正常模式下...
- 问7升级win10(win 7升win 10)
-
Windows7系统如何升级为win10。嗯,这个是需要重新安装操作系统的。请在安装前先预备好备份好windows7系统下面你自己的个人数据备份到自己的移动硬盘上面。然后通过USB导的方式进行安装,...
- windows中文叫什么(win的中文是什么)
-
windows是窗口的意思,翻译到中文可以是窗户的意思。现在windows的意思大多用在PC领域,泛指微软出品的电脑、手机操作系统。windows操作系统业可以叫做视窗操作系统、可视化图形界面操作系统...
- win7系统激活工具下载免费(win7激活工具免费版)
-
KMSpico是一个非常好用的Win7系统激活工具。它简单易用,只需一步操作即可激活系统,无需复杂的设置和操作。同时,它支持离线激活和在线激活两种方式,让用户可以根据自己的需求选择使用。此外,KMSp...
- tenda路由器怎么设置网速快(tendawifi怎么设置网速快)
-
tenda设置网速最快方法如下:1.登入无线路由器;在浏览器地址输入192.168.0.1;无线路由器进入无线路由器控制界面几乎都是这个地址;原始密码admin。2.接着点击“高级设置”或者页面右上角...
- 功能最全的pe系统(pe系统功能介绍)
-
1、Windows预安装环境,是带有有限服务的最小Win32子系统,基于以保护模式运行的WindowsXPProfessional及以上内核。它包括运行Windows安装程序及脚本、连接网络共享、...
- 华为鸿蒙操作系统下载(华为鸿蒙电脑操作系统下载)
-
鸿蒙系统可以下载拼多多,但需要注意的是,目前拼多多官方并未推出面向鸿蒙系统的专门版本,所以需要通过安装第三方应用市场或者通过APK安装包的方式才能在鸿蒙系统上安装使用。此外,由于鸿蒙系统和拼多多应用的...
- 系统iso怎么安装(系统iso怎么安装应用)
-
ISO系统安装详细教程步骤如下:1.准备所需材料:一台符合系统要求的计算机、ISO系统文件、可启动的USB驱动器或光盘、系统激活密钥(如果需要的话)。2.制作启动盘:如果你使用USB驱动器作为启动...
欢迎 你 发表评论:
- 一周热门
-
-
抖音上好看的小姐姐,Python给你都下载了
-
全网最简单易懂!495页Python漫画教程,高清PDF版免费下载
-
飞牛NAS部署TVGate Docker项目,实现内网一键转发、代理、jx
-
Python 3.14 的 UUIDv6/v7/v8 上新,别再用 uuid4 () 啦!
-
python入门到脱坑 输入与输出—str()函数
-
Python三目运算基础与进阶_python三目运算符判断三个变量
-
(新版)Python 分布式爬虫与 JS 逆向进阶实战吾爱分享
-
失业程序员复习python笔记——条件与循环
-
系统u盘安装(win11系统u盘安装)
-
Python 批量卸载关联包 pip-autoremove
-
- 最近发表
- 标签列表
-
- python计时 (73)
- python安装路径 (56)
- python类型转换 (93)
- python进度条 (67)
- python吧 (67)
- python的for循环 (65)
- python格式化字符串 (61)
- python静态方法 (57)
- python列表切片 (59)
- python面向对象编程 (60)
- python 代码加密 (65)
- python串口编程 (77)
- python封装 (57)
- python写入txt (66)
- python读取文件夹下所有文件 (59)
- python操作mysql数据库 (66)
- python获取列表的长度 (64)
- python接口 (63)
- python调用函数 (57)
- python多态 (60)
- python匿名函数 (59)
- python打印九九乘法表 (65)
- python赋值 (62)
- python异常 (69)
- python元祖 (57)
