AMBER分子动力学模拟之分子模拟-- HIV蛋白酶-抑制剂复合物(2)

news2024/12/22 19:04:12

AMBER分子动力学模拟之分子模拟-- HIV蛋白酶-抑制剂复合物(2)

体系平衡

体系的优化分为两步,第一步固定蛋白和小分子,对添加的水盒优化。第二步则是对整个体系的优化。

对水盒子优化

vim min1.in

&cntrl
imin=1,
maxcyc=10000, 
ncyc=5000,
ntb=1,
ntr=1,
ntpr=500,
cut=10,
ioutfm=1,
ntxo=1,
restraint_wt=500.0,
restraintmask=':1-199',
/

对复合物优化

vim min2.in

&cntrl
  imin = 1,
  maxcyc = 50000, 
  ncyc = 25000, 
  ntb = 1, 
  ntr = 0, 
  ntpr = 500, 
  cut = 10, 
  ioutfm = 1, 
  ntxo = 1, 
/

升温

由于动力学模拟的是真实的生物体环境,因此必须使研究对象升温升压到临界值,体系达到平衡,才能做实际的动力学模拟

vim md0.in

&cntrl
imin=0, 
irest =0, 
ntx=1, 
ntb=1, 
cut =10, 
ntr=1, 
ntc =2, 
ntf=2, 
tempi=0.0, 
temp0=300.0, 
ntt =3, 
gamma_ln=1.0, 
nstlim=150000, 
dt=0.002, 
ntpr=500, 
ntwx=1000, 
ntwr=500, 
restraint_wt=10.0, 
restraintmask=':1-199',
/

动力学模拟

一般来说第一步的MD模拟时间要长一些,确保体系构象已经稳定,之后进行第二步的MD用于能量、相互作用分析等。
vim md1.in

&cntrl
imin=0, 
irest=1, 
ntx=5, 
ntb=2, 
pres0=1.0, 
ntp=1, 
taup=2.0, 
cut=10, 
ntr=0, 
ntc=2, 
ntf=2, 
tempi=300.0, 
temp0=300.0, 
ntt=3, 
gamma_ln=1.0, 
nstlim=20000000, 
dt=0.002, 
ntpr=500, 
ntwx=500, 
ntwr=500,
/

vim md2.in

&cntrl
imin=0, 
irest=1, 
ntx=5, 
ntb=2, 
pres0=1.0, 
ntp=1, 
taup=2.0, 
cut=10, 
ntr=0,
ntc=2, 
ntf=2, 
tempi=300.0, 
temp0=300.0, 
ntt=3, 
gamma_ln=1.0, 
nstlim=5000000, 
dt=0.002, 
ntpr=500, 
ntwx=500, 
ntwr=500, 
ntwprt=3198,
/

mdrun 执行文件

在mdrun.exe文件中,调用amber的pmemd.cuda模块进行动力学模拟
vi mdrun.exe写入

### CPU 
sander -O -i min1.in -p pep.top -c pep.crd -o min1.out -r min1.rst -ref pep.crd < /dev/null
sander -O -i min2.in -p pep.top -c min1.rst -o min2.out -r min2.rst < /dev/null
sander -O -i md0.in -p pep.top -c min2.rst -o md0.out -x md0.crd -r md0.rst -ref min2.rst < /dev/null
sander -O -i md1.in -p pep.top -c md0.rst -o md1.out -x md1.crd -r md1.rst < /dev/null
sander -O -i md2.in -p pep.top -c md1.rst -o md1.out -x md2.crd -r md2.rst < /dev/null


### cpu_mpi
mpirun -np 2 sander -O -i min1.in -p pep.top -c pep.crd -o min1.out -r min1.rst -ref pep.crd < /dev/null
mpirun -np 2 sander -O -i min2.in -p pep.top -c min1.rst -o min2.out -r min2.rst < /dev/null
mpirun -np 8 sander -O -i md0.in -p pep.top -c min2.rst -o md0.out -x md0.crd -r md0.rst -ref min2.rst < /dev/null
mpirun -np 8 sander -O -i md1.in -p pep.top -c md0.rst -o md1.out -x md1.crd -r md1.rst < /dev/null
mpirun -np 8 sander -O -i md2.in -p pep.top -c md1.rst -o md1.out -x md2.crd -r md2.rst < /dev/null

## GPU
pmemd.cuda  -O -i min1.in -p pep.top -c pep.crd -o min1.out -r min1.rst -ref pep.crd < /dev/null
pmemd.cuda  -O -i min2.in -p pep.top -c min1.rst -o min2.out -r min2.rst < /dev/null
pmemd.cuda  -O -i md0.in -p pep.top -c min2.rst -o md0.out -x md0.crd -r md0.rst -ref min2.rst < /dev/null
pmemd.cuda  -O -i md1.in -p pep.top -c md0.rst -o md1.out -x md1.crd -r md1.rst < /dev/null
pmemd.cuda  -O -i md2.in -p pep.top -c md1.rst -o md2.out -x md2.crd -r md2.rst < /dev/null

目前准备已经准备好进行md的参数文件(带有水盒子信息的拓扑和坐标文件)、md的输入文件,min1.in min2.in md0.in md1.in md2.in,以及MD执行文件mdrun.exe

将mdrun.exe的权限设置为可执行文件
chmod +xmdrun.exe
nohup./mdrun.exe&

输出文件

在md的输出文件中,所有的*out文件均含有每一步相对应的Amber版本,工作日期,工作路径、输入文件的详细参数,包含未曾修改的默认参数、引用文献、我们重点关注和分析的是每一步输出的能量项等等。
vi min1.out

vi min2.out

vi md0.out

vi md1.out

vi md2.out

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/531207.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

杂记——22.强弱类型语言的区别

这篇文章我们来讲一下强类型语言和弱类型语言&#xff0c;并分析一下二者的区别 目录 1.强类型语言 2.弱类型语言 3.区别 3.1编译型和解释型 3.2动态语言和静态语言 3.3 强类型定义语言和弱类型定义语言 4.静态类型语言的优势 1.强类型语言 强类型语言也称为强类型定义…

error C2143: 语法错误 : 缺少“)”(在“常量”的前面)

【背景】最近做一个项目&#xff0c;需要使用到凌华的运动控制卡和IO卡&#xff0c;在测试工程中&#xff0c;使用的类和引用的库文件都一切正常。 测试完成后&#xff0c;本着可重复移植的考虑&#xff0c;将整个文件夹添加到现有工程中&#xff0c;却出现了非常奇怪的编译错误…

JVM系列-第9章-StringTable(字符串常量池)

StringTable&#xff08;字符串常量池&#xff09; String的基本特性 String&#xff1a;字符串&#xff0c;使用一对 “” 引起来表示 String s1 "atguigu" ; // 字面量的定义方式String s2 new String("hello"); // new 对象的方式String被…

基于html+css的图展示73

准备项目 项目开发工具 Visual Studio Code 1.44.2 版本: 1.44.2 提交: ff915844119ce9485abfe8aa9076ec76b5300ddd 日期: 2020-04-16T16:36:23.138Z Electron: 7.1.11 Chrome: 78.0.3904.130 Node.js: 12.8.1 V8: 7.8.279.23-electron.0 OS: Windows_NT x64 10.0.19044 项目…

【fly-iot飞凡物联】():开源项目Apache IoTDB,开源项目学习,原来还有这样的项目,关于IOT的几个开源项目汇总下

目录 前言1&#xff0c;关于&#xff1a;开源项目Apache IoTDB2&#xff0c;还有个admin后台3&#xff0c;thinglinks项目4&#xff0c;thingsboard-ui-vue项目5&#xff0c;apache pulsar项目6&#xff0c;ActorCloud项目 前言 本文的原文连接是: https://blog.csdn.net/freew…

ASP.NET MVC-WebAPI请求

一、无参数Get请求 一般的get请求我们可以使用jquery提供的$.get() 或者$.ajax({type:"get"}) 来实现&#xff1a; 请求的后台Action方法仍为上篇文章中的GetUser() : 也可以用$.ajax({type:"get"}) 方式&#xff0c;正确的获得了返回数据&#xff1a; 二、…

Grafana功能菜单介绍(04)

Grafana的功能菜单设计为侧边栏(sidebar)形式,可以折叠隐藏,便于我们更加专注数据的可视化。现将菜单栏各项功能进行编号讲解,如下图所示。 ① Grafana的logo,即当前页为Grafana的Home page,在任何页面点击Grafana的logo,都会跳到Home Page。 ② 新建按钮,用于创建Da…

Camtasia Studio2023新功能及下载安装使用教程

Camtasia Studio2023新版本包含了屏幕录像、视频剪辑和编辑、视频录音配音、视频菜单制作、视频剧场和视频播放等功能。它可以在任何颜色模式下记录屏幕动作&#xff0c;包括影像、 音效、鼠标移动轨迹&#xff0c;解说声音等等&#xff1b;它支持视频播放和编辑压缩&#xff0…

Mars3d实现单击按钮高亮矢量对象

调用高亮的关键代码&#xff1a; let graphictruck new mars3d.graphic.ModelEntity({ position: [116.327881, 31.018378, 5000], style: { url: //data.mars3d.cn/gltf/mars/qiche.gltf, heading: 90, scale: 0.9, …

应用软件系统架构设计的“七种武器”

对于软件架构这一概念&#xff0c;有太多的版本&#xff0c;目前在业界由大师级人物或组织提出的对这一概念的阐述就超过十种以上&#xff0c;我个人比较赞同RUP(Rational Unified Process)中对软件架构的定义&#xff0c;即软件架构包含了关于以下问题的重要决策&#xff1a; …

剑指offer(C++)-JZ48:最长不含重复字符的子字符串(算法-动态规划)

作者&#xff1a;翟天保Steven 版权声明&#xff1a;著作权归作者所有&#xff0c;商业转载请联系作者获得授权&#xff0c;非商业转载请注明出处 题目描述&#xff1a; 请从字符串中找出一个最长的不包含重复字符的子字符串&#xff0c;计算该最长子字符串的长度。 数据范围…

c++ 11标准模板(STL) std::set(四)

定义于头文件 <set> template< class Key, class Compare std::less<Key>, class Allocator std::allocator<Key> > class set;(1)namespace pmr { template <class Key, class Compare std::less<Key>> using se…

让 ChatGPT 来担任这个召之即来挥之即去的私人健身教练

健身教练 和ChatGPT聊天可以运用在各种正式、非正式&#xff0c;工作、休闲场合&#xff0c;让我们再来看一个场景&#xff0c;也是年轻人群体中非常热门的话题&#xff1a;健身。健身已经是年轻人最流行的活动&#xff0c;既可以是私密的个人健身&#xff0c;也可以是呼朋唤友…

JAVA集合,复杂度,泛型

泛型: 1.存储数据的时候 可以帮我们进行自动的类型检查 2.获取元素的时候&#xff0c;可以帮我们进行类型转换 获取数据时&#xff0c;不需要强转类型转换。 类型形参一般使用一个大写字母表示&#xff0c;常用的名称有&#xff1a; E 表示 Element K 表示 Key V 表示 …

多目标跟踪:视觉联合检测和跟踪

国内头部以自动驾驶全站技术为主线的交流学习社区&#xff08;感知、归控等&#xff09;&#xff0c;包含大量前沿论文解读、工程实践(源代码)、视频课程&#xff0c;热招岗位。欢迎加入&#xff01; 点击上方“迈微AI研习社”&#xff0c;选择“星标★”公众号 重磅干货&#…

Python每日一练(20230516) 打家劫舍 I\II\III\IV HouseRobber

目录 1. 打家劫舍 I House Robber i 2. 打家劫舍 II House Robber ii 3. 打家劫舍 III House Robber iii 4. 打家劫舍 IV House Robber iv &#x1f31f; 每日一练刷题专栏 &#x1f31f; Golang每日一练 专栏 Python每日一练 专栏 C/C每日一练 专栏 Java每日一练 专栏…

除了OA以外,只用点击按钮就能实现协同办公,靠一个工具就能实现

OA不同时期&#xff0c;不同功能 OA是办公自动化系统&#xff0c;与互联网技术的发展有着密切的联系。 很多办公场景都离不开OA协同系统。OA是办公系统&#xff0c;借助与互联网的交互性&#xff0c;开创出的一种新型的办公方式。 所以&#xff0c;说到OA软件&#xff0c;在不…

Linux 网络基础(1)基础知识、IP地址、端口、协议、网络字节序

基础知识 网络发展背景&#xff1a; 网络的划分&#xff1a;局域网&#xff08;覆盖范围在1000m以内&#xff09;、城域网&#xff08;覆盖范围在20km以内&#xff09;、广域网&#xff08;更大范围&#xff09; 组网方式&#xff1a;以太网、令牌环网.... 日常名词&#xff1a…

ChatGPT「最强竞争对手」Claude有了史诗级更新~~~百页数据一分钟读完还能做出总结

号称ChatGPT「最强竞争对手」的Claude&#xff0c;迎来史诗级更新—— 模型记忆力原地起飞&#xff0c;现在1分钟看完一本数万字的小说&#xff0c;已经是一件很理所当然的事情了。 消息一出&#xff0c;评论区直接炸了&#xff0c;网友们蜂拥而至&#xff0c;「wocao连连」&am…

报表解决方案Telerik Reporting - 轻松完成嵌入式报表开发

假设您是一名开发人员&#xff0c;思考一下您的用户将从Web报表设计器中得到什么吗&#xff1f;他们能突出异常数据吗?他们能够从头开始创建一个报告吗?让我们站在商业用户的角度&#xff0c;看看这到底是怎么回事。 Telerik Reporting Web Designer是一个多功能且用户友好的…