基于Gromacs模拟软件分析小分子配体与蛋白结合之后的稳定性

背景：随着分子病理学的发展，人们对疾病发生和药物发挥药效的分子机制有了更深入的认识，伴随着蛋白质晶体学的发展，越来越多的蛋白质晶体被解析出来，更多的药物与其靶标相互作用的三维结构被揭示，对于理解药物发挥作用的机制起到关键作用。基于药物靶点结构的合理药物设计和筛选愈来愈受到重视，药物靶点和配体的亲和力以及相互作用的结构是十分重要的研究数据。但是 Protein Data Bank（PDB）数据库并不直接提供蛋白质和配体的亲和力数据，大量的亲和力数据需要用户手动在文献中寻找；与此同时，在蛋白质纯化和结晶的过程中需要加入大量的添加剂（如去污剂、盐离子、缓冲剂、金属离子、溶剂分子等），其将会成为污染物出现在解析后的晶体结构中，这些污染物并不能与蛋白质特异性结合，并将在许多自动化的处理中造成干扰。因此收集靶点蛋白质和受体亲和力数据，评估 PDB 结构中配体的生物相关性对药物虚拟筛选、分子对接、蛋白质功能注释及相关软件设计具有重要意义。

具体步骤：

1. 蛋白和小分子的准备

如下列途径：

2. 对接

对接软件目前市面上有许多软件都能做docking

如 DS、Chimera、Autodock系列

下面用autodock vina演示对接过程

a.file-read molecule

删除水分子和其他配体，常规操作不用解释

b. file-read molecule

同样的处理方式，加极性氢，计算电荷和添加原子类型

c.处理完蛋白质和小分子后，开始设置对接盒子grid-box，最后使用vina输出参数，最后run autodock vina即可，具体的参照《autodock vina的安装及使用》。

对接完成后，可以分析小分子的极性接触。

如：

3. Gromacs 模拟步骤

将配体pdb文件加载至文件夹中

3.1 蛋白质拓扑结构生成

gmx pdb2gmx -f protein.pdb -o protein_processed.gro -water spc -ignh

3.2 配体拓扑结构处理

根据你所需要的力场修改小分子拓扑。

3.3 复合物构建

复制protein_processed.gro文件并命名为complex.gro，将每个配体gro文件的坐标部分复制到complex.gro的蛋白质原子坐标末尾，最后更新complex.gro开头的原子总数。

3.3 周期性边界条件

4. # 定义盒子尺寸

5. gmx editconf -f complex.gro -o newbox.gro -bt cubic -box 8.0 8.0 8.0

6. # 添加溶剂

7. gmx solvate -cp newbox.gro -cs spc216.gro -p topol.top -o solv.gro

8. # 加入离子

9. gmx grompp -f em.mdp -c solv.gro -r solv.gro -p topol.top -o ions.tpr -maxwarn 2

10. gmx genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL -np 5

3.5 能量最小化，预平衡NVT，NPT

gmx grompp -f em_real.mdp -c solv_ions.gro -r solv_ions.gro -p topol.top -n index.ndx -o em.tpr -maxwarn 2

得到em.tpr文件之后通过以下命令运行能量最小化：

gmx mdrun -v -deffnm em

gmx grompp -f nvt.mdp -c em.gro -r em.gro -p topol.top -n index.ndx -o nvt.tpr -maxwarn 2

gmx mdrun -deffnm nvt

gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -t nvt.cpt -p topol.top -n index.ndx -o npt.tpr -maxwarn 2

gmx mdrun -deffnm npt

3.6 生产模拟

首先运行如下命令生成 tpr 文件：

gmx grompp -f md.mdp -c npt.gro -r npt.gro -t npt.cpt -p topol.top -n index.ndx -o md_results_1.tpr -maxwarn 2

运行主体模拟：

gmx mdrun -deffnm md_results_1

模拟完成后，可以计算蛋白稳定性，通过RMSD、RMSF

最后通过PBSA计算结合能

。

总结：通过上述对接，模拟，分析，我们可以了解蛋白与小分子结合的稳定性，来判断结合的好坏。也可以利用模拟来提前预测小分子能否在实验中得到好的验证。

最后,如有MD相关需求,可以联系我们.

微信公众号:320科技工作室

VX: CAE320

基于Gromacs模拟软件分析小分子配体与蛋白结合之后的稳定性相关推荐

javplayer 使用教程_药物设计软件Sybyl教程（一）：基于Tripos力场对小分子配体进行结构优化...
文 / 利刃君微信ID/ ziyuanliren666全文共1104字,推荐阅读时间6分钟. 教程内容: 以SYBYL-X 2.0软件为例,对小分子配体进行基于Tripos力场的能量最小化计算,优化分 ...
【Amber】带小分子配体分子动力学模拟
目录适用于Amber18\20版本一.文件准备 1.文件检查与拆分 2.小分子文件处理 3.蛋白质文件处理 4.复合体文件处理二.分子动力学模拟 1.能量最小化 2.体系加热 3.恒压平衡 4. ...
分子模拟软件amber_使用Amber创建小分子与蛋白质复合蛋白的坐标和拓扑文件
复合蛋白amber坐标和拓扑文件的创建作者:朱宁来源:大科研小分享前言分子动力学(Molecular Dynamics, MD)是一门结合物理,数学和化学的综合技术.目前主流分子动力学软 ...
火车售票系统（基于UML的软件分析与设计模型建模实验）
火车售票系统(基于UML的软件分析与设计模型建模试验) 文章目录火车售票系统(基于UML的软件分析与设计模型建模试验) 一.需求二.UML用例图三.UML类图四.UML时序图一.需求实现一 ...
小分子化合物的蛋白靶点（蛋白质、酶、受体） + IC50、EC50、Ki 是什么？
受体拮抗剂(英语:receptor antagonist),也叫阻断剂(英语:blocker),是药理学术语,指能与受体结合,并能阻止激动剂产生效应的一类配体物质. 抑制剂 & 拮抗剂等小分子 ...
基于NLP的软件分析和程序漏洞挖掘（二）
找到NLP在软件分析和漏洞挖掘的契合点之后,我们需要去调研这方面现有的研究情况.发现已经有人使用NLP来对二进制代码进行分析,主要有3种方向上使用NLP分析二进制代码,但大多数都是对二进制代码进行静态 ...
软件设计实验（一）C#火车售票系统-基于UML的软件分析与设计模型建模实验（用例图、类图、时序图）
主要从用户管理模块分析: (1)用例图: (2)类图: (3)时序图:
Amber小分子-蛋白复合体分子动力学模拟
Amber小分子-蛋白复合体分子动力学模拟以前经常用GROMACS进行分子动力学模拟,后来试了一下Amber后发现,在我当前配置的GPU资源上,果然还是Amber更快一些,GROMACS太吃CPU资 ...
全球与中国铸造模拟软件市场深度研究分析报告（2022）
[报告篇幅]:85 [报告图表数]:123 [报告出版时间]:2022年1月报告摘要本文研究全球及中国市场铸造模拟软件现状及未来发展趋势,侧重分析全球及中国市场的主要企业,同时对比北美.欧洲.中国 ...

基于Gromacs模拟软件分析小分子配体与蛋白结合之后的稳定性

基于Gromacs模拟软件分析小分子配体与蛋白结合之后的稳定性相关推荐

最新文章

热门文章