Autodock的基本使用步骤

2023-11-10 05:21
文章标签 使用 步骤 基本 autodock

本文主要是介绍Autodock的基本使用步骤,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

 

 1.文件准备

        在进行对接前,我们必须保证在工作目录中有这五个文件。其中hsg表示受体,ind表示配体。

2.导入作为受体分子的pdb文件、去游离O原子(水)、加H原子,保存为pdbqt文件

        1>选择File-->Read Molecule        将工作目录下的pdb文件导入

        2>选择Edit-->Delete Water        将水分子去除

        3>选择Edit-->Hydrogens-->Add-->OK        对残基加上H原子

        4>选择Grid-->Macromolecule-->Choose        将处理好的构象选择为受体,保存为作为受体的pdbqt文件,在此步骤中auto dock对受体分子进行电荷计算

3.处理小分子配体pdb文件,保存为pdbqt文件

        1>Edit-->Hydrogens-->Add-->OK        进行加H原子操作

        2>ligand-->input-->choose        保存配体的pdbqt文件(这一步完成自动添加gasteiger电荷(修改原子类型)、合并非极性的H原子和检查结构的可扭转键。(在预览窗里结构的H数量减少)

        3>Ligand-->Torsion Tree--> Detect Root 以及Ligand-->Torsion Tree--> Choose Torsion       检查配体分子的所有扭转键,方便之后查找配体的不同构象

前者找到这个结构的中心root,后者使用后,不可扭转的键为红色,可扭转的键为绿色,表示酰胺键为粉色

        4>按住shift键,对粉色键进行点击,可将其转变为绿色的可旋转键,对绿色键可以做这个操作,红色键不行        调整可旋转键的位置和个数

        5>Ligand-->Output-->Save as PDBQT        进行pdbqt文件保存

4.导入受体和配体的pdbqt文件

        1>选择Grid-->Macromolecule-->Open        打开准备好的受体的pdbqqt文件,对弹出选择Yes

        2>选择Grid-->Set Map Types-->Open Ligand        打开准备的配体的pdbqt文件

 5.设置工作盒子

        1>在预览窗里出现一个三色盒子,在弹窗来修改盒子的3轴,让盒子包围受体和配体,也可以只是包围配体和活性位点将小分子拖到盒子外

        2>选择DejaVu GUI 打开新弹窗

        2>在上面的红框内,取消mouse transforms apply to "root" object only 选择,点击下方红框内的ligand将指定鼠标的操作对象

        3>鼠标左键进行旋转,右键进行移动,滑轮进行缩放。右键预览窗里配体移出盒子,原先取消mouse transforms apply to "root" object only 重新勾选,关闭弹窗

6.保存调整好的盒子为gpf文件

        1>选择File-->Close saving current

        2>选择Gird-->Output-->Save GPF        进行文件命名和保存

        3>选择run--Run AutoGrid        打开新的弹窗,注意工作目录,在Parameter Filename选择wd里的保存的gpf文件,会在下面的Log Filename建立同名glg文件

7.进行分子对接计算,结果查看和保存

        1>进行半柔性对接,所以作为受体的大分子是刚性的(Rigid)

        2>选择作为对接配体的pdbqt文件,弹窗中选择accept

        3>选择合适算法

        Docking-->Search Parameters -->Genetic Algorithm -->打开的弹窗中选择Accept

        Docking-->Docking Parameters-->打开的弹窗里选择Accept

        Docking-->Output-->Lamarckian GA 4.2,,-->定义dpf文件名,作为对接结果的保存文件

Number of GA Runs    表示我们对接多少次,默认为0次,官方建议对接50次以上。
Maximum Number of evals 和Maximum Number of generations 值可以设置大一些,从而获得较优构象,但值越大,后面运行AutoDock的时间会越长

        4>Run-->RunAutoDock-->弹窗里确保工作目录正确,brower定义的dpf文件,选择launch

同样等待弹窗自动消失        运行Dock计算

        5.查看计算完成得到的对接构象

         Analyze-->Dockings-->Open...,打开计算完成的dlg文件,读取获取的前10个对接构象

        Analyze-->Macromolecule-->Open..自动打开对应的受体大分子

        Analyze-->Conformations-->Played,ranked  by energy...,打开一个播放10个构象的弹窗      

                ①点击播放10个构象,构象的顺序是按结合能从小到大排列

                

                ②点击打开查看H键等信息的弹窗

                         ①选择show info        查看该构象结合能的大小和各个能量的贡献和形成H键的个数和对应的残基原子

                         ②选择build H-bonds        勾选后在结构显示窗里可视化H键

         6>对对接结果进行保存

                ①选择结合能最小的对接结果

                ②点击Write Complex按钮,对第1个对接构象保存为pdbqt文件

8.对接结果是否可靠的判断

         用文本文档软件打开我们的dlg文件

1>文件头部是使用的软件版本等信息,然后是程序运行的一些参数记录
2>FINAL GENETIC ALGORITHM DOCKED STATE往下是我们的对接结果,如果不是用GENETIC ALGORITHM DOCKED,该出名称不一样,结果中,Run = 1是第一次对接的结果
3>在CLUSTERING HISTOGRAM处会看见一个表。10个运行结果统计在这里表示

       上面一直说的结合能的单位是kcal/mol,kj/mol是标准单位,1KCal等于4.184KJ。kj/mol是焓的单位。因为焓是一个状态函数,所以某物质所含的能量是随摩尔质量的改变而改变的,物质的量的单位是mol,能量的单位是J,所以焓的单位就是J/mol或者kJ/mol。如果结合能小于-1.2kcal/mol或者小于-5kj/mol,那么我们认为对接结果是可行的

这篇关于Autodock的基本使用步骤的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


原文地址:https://blog.csdn.net/weixin_72518947/article/details/131881848
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.chinasem.cn/article/380599

相关文章

使用jenv工具管理多个JDK版本的方法步骤

《使用jenv工具管理多个JDK版本的方法步骤》jenv是一个开源的Java环境管理工具,旨在帮助开发者在同一台机器上轻松管理和切换多个Java版本,:本文主要介绍使用jenv工具管理多个JD... 目录一、jenv到底是干啥的?二、jenv的核心功能(一)管理多个Java版本(二)支持插件扩展(三)环境隔

SQL中JOIN操作的条件使用总结与实践

《SQL中JOIN操作的条件使用总结与实践》在SQL查询中,JOIN操作是多表关联的核心工具,本文将从原理,场景和最佳实践三个方面总结JOIN条件的使用规则,希望可以帮助开发者精准控制查询逻辑... 目录一、ON与WHERE的本质区别二、场景化条件使用规则三、最佳实践建议1.优先使用ON条件2.WHERE用

Java中Map.Entry()含义及方法使用代码

《Java中Map.Entry()含义及方法使用代码》:本文主要介绍Java中Map.Entry()含义及方法使用的相关资料,Map.Entry是Java中Map的静态内部接口,用于表示键值对,其... 目录前言 Map.Entry作用核心方法常见使用场景1. 遍历 Map 的所有键值对2. 直接修改 Ma

MySQL 衍生表(Derived Tables)的使用

《MySQL衍生表(DerivedTables)的使用》本文主要介绍了MySQL衍生表(DerivedTables)的使用,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学... 目录一、衍生表简介1.1 衍生表基本用法1.2 自定义列名1.3 衍生表的局限在SQL的查询语句select

MyBatis ResultMap 的基本用法示例详解

《MyBatisResultMap的基本用法示例详解》在MyBatis中,resultMap用于定义数据库查询结果到Java对象属性的映射关系,本文给大家介绍MyBatisResultMap的基本... 目录MyBATis 中的 resultMap1. resultMap 的基本语法2. 简单的 resul

Mybatis Plus Join使用方法示例详解

《MybatisPlusJoin使用方法示例详解》:本文主要介绍MybatisPlusJoin使用方法示例详解,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,... 目录1、pom文件2、yaml配置文件3、分页插件4、示例代码:5、测试代码6、和PageHelper结合6

MySQL分区表的具体使用

《MySQL分区表的具体使用》MySQL分区表通过规则将数据分至不同物理存储,提升管理与查询效率,本文主要介绍了MySQL分区表的具体使用,具有一定的参考价值,感兴趣的可以了解一下... 目录一、分区的类型1. Range partition(范围分区)2. List partition(列表分区)3. H

使用SpringBoot整合Sharding Sphere实现数据脱敏的示例

《使用SpringBoot整合ShardingSphere实现数据脱敏的示例》ApacheShardingSphere数据脱敏模块,通过SQL拦截与改写实现敏感信息加密存储,解决手动处理繁琐及系统改... 目录痛点一:痛点二:脱敏配置Quick Start——Spring 显示配置:1.引入依赖2.创建脱敏

Python使用smtplib库开发一个邮件自动发送工具

《Python使用smtplib库开发一个邮件自动发送工具》在现代软件开发中,自动化邮件发送是一个非常实用的功能,无论是系统通知、营销邮件、还是日常工作报告,Python的smtplib库都能帮助我们... 目录代码实现与知识点解析1. 导入必要的库2. 配置邮件服务器参数3. 创建邮件发送类4. 实现邮件

Go语言中Recover机制的使用

《Go语言中Recover机制的使用》Go语言的recover机制通过defer函数捕获panic,实现异常恢复与程序稳定性,具有一定的参考价值,感兴趣的可以了解一下... 目录引言Recover 的基本概念基本代码示例简单的 Recover 示例嵌套函数中的 Recover项目场景中的应用Web 服务器中