| 查看: 3169 | 回复: 6 | |||
[交流]
关于用amber跑REMD计算的流程及输入文件的设置
|
|||
|
最近看了一篇论文，http://www.ncbi.nlm.nih.gov/pubmed/22339436 作者使用amber做了REMD（replica exchange molecular dynamic）的计算。 自己正好也有类似需要，想用REMD模拟。就想按论文相同的条件实现一下，不过论文里的过程写的不是很详细。下面是原文的method部分： In this work, the temperature setting used for sheep PrP 125-230 such that MD simulations were exchanged at 320.0, 322.0,324.0, 326.0, 328.1, 330.2, 332.3, 334.4, 336.5, 338.6, 340.8, 343.0, 345.2, 347.4, 349.6, 351.8, 354.0, 356.2, 358.5, 360.8, 363.1, 365.4, 367.7, and 370.0 K was used. All of the replicas were equilibrated for 20 ns without exchanging temperatures and then extended for 65 ns of REMD simulation. The generalized Born model used in this study modified the calculation of Born radii and improved the accuracy in the solvent polarization for macromolecules. The combinational use of the all-atom point-charge force-field (also known as ff03) and the generalized Born model led to successful folding of several proteins. The AMBER 11 simulation package 26 was used in both REMD simulation and data analysis. The melted huPrP 121-230 was computed starting from an extended huPrP. To generate the initial extended structure, a heating method was applied to a known NMR structure (PDBcode 1hjn,15Figure 1A), enabling it to unfold at 600 K for 40 ns of MD simulation to result in an extended conformation (Figure 1B) as described previously. During this simulation, the disulfide covalent bond between residues 179 and 214 was preserved. In total, 24 replicas with duration of 65 ns and with an integration time step of 2 fs were computed based on the extended huPrP with different random number seeds to generate the initial conditions. A 16 Å force-shifted non-bonded cutoff and generalized Born solvent models with salt concentration of 0.2 M were applied. 根据文章的描述我自己总结了一下，并按文章描述写了各个环节的输入文件。想和大家交流一下，看看我写的有没有不合适的问题。 Simulation Procedure: 1. system building 2. system minimization 3. heating system 4. generate the extend conformation 5. local minimization after heating system 6. equilibrate the every replica 7. REMD simulation 1.leap.inp for system building source pdb: 1ag2 use the the ff03 (Duan et al.) force field leap.inp source leaprc.ff03.r1 # load pdb file 1ag2 = loadPdb input.pdb # solvation # solvatebox 1ag2 TIP3PBOX 18.0 # save 1ag2 to pdb file savePdb 1ag2 1ag2.pdb # add countinos # addions 1ag2 Cl- 0 # addions 1ag2 Na+ # s-s bond bond 1ag2.179.SG 1ag2.214.SG # save 1ag2 to prmtop and inpcrd files saveAmberParm 1ag2 1ag2.prmtop 1ag2.inpcrd # finish quit 2. system minimization minimisation for heated system &cntrl imin=1, maxcyc=1000, ncyc=500, cut=999., rgbmax=999.,igb=1, ntb=0, ntpr=100 / ~ ~ 3. heat the system heating system from 0 K to 600K. &cntrl nstlim = 50000, dt = 0.002, ntt = 1, tautp = 1.0, tempi = 0, temp0 = 600, ntc =2, ntf = 2, ntpr =100, ntwx = 100, ntb = 0, igb = 1, #这里设置和后面remd的设置不同（igb = 5）结果可能受影响。 cut = 999.0,rgbmax = 999.0, / 4. generate the full unfolded conformation To generate the initial extend structure, a heating method was used to a known NMR structure (PDB code:1ag2), enabling it to unfold at 600 K for 40 ns of LD simulation to result in an extended conformation. 40nsld.inp enabling the heated NMR structure to unfold at 600 K for 40ns of LD simulation &cntrl irest = 1, ntx = 5, nstlim = 20000000, dt = 0.002, ntt = 3, gamma_ln = 1.0, tempi = 600,temp0 = 600, ntb = 0, igb = 2, ntpr = 500, ntwx = 1000, ntwr = 2000000, ntc = 2, ntf = 2, cut = 999.0,rgbmax = 999.0, / | RE_POSITION Moving by -2.485812 -1.096141 0.618430 NSTEP = 20000000 TIME(PS) = 40060.000 TEMP(K) = 577.57 PRESS = 0.0 Etot = 1351.4737 EKtot = 2407.9819 EPtot = -1056.5082 BOND = 626.4085 ANGLE = 1429.4823 DIHED = 1419.4285 1-4 NB = 376.4803 1-4 EEL = 3772.0639 VDWAALS = -368.8173 EELEC = -5949.1848 EGB = -2362.3695 RESTRAINT = 0.0000 ------------------------------------------------------------------------- A V E R A G E S O V E R ******* S T E P S NSTEP = 20000000 TIME(PS) = 40060.000 TEMP(K) = 600.20 PRESS = 0.0 Etot = 1514.3153 EKtot = 2502.2899 EPtot = -987.9746 BOND = 644.7390 ANGLE = 1506.2954 DIHED = 1362.4635 1-4 NB = 397.8849 1-4 EEL = 3727.1078 VDWAALS = -361.6498 EELEC = -5967.0181 EGB = -2297.7973 RESTRAINT = 0.0000 ------------------------------------------------------------------------- 结果Eptot 的plot图 600kmd.inp enabling the heated NMR structure to unfold at 600 K for 40ns of LD simulation &cntrl irest = 0, ntx = 1, nstlim = 20000000, dt = 0.002, ntt = 3, gamma_ln = 1.0, tempi = 600,temp0 = 600, ntb = 0, igb = 2, ntpr = 500, ntwx = 1000, ntwr = 2000000, ntc = 2, ntf = 2, cut = 999.0,rgbmax = 999.0, / 600kmd.out ... | RE_POSITION Moving by -0.033675 -0.790281 -1.085325 NSTEP = 2000000 TIME(PS) = 4060.000 TEMP(K) = 601.00 PRESS = 0.0 Etot = 1554.3820 EKtot = 2505.6360 EPtot = -951.2540 BOND = 654.9106 ANGLE = 1485.1469 DIHED = 1383.1021 1-4 NB = 422.1188 1-4 EEL = 3774.2620 VDWAALS = -379.3759 EELEC = -6038.4682 EGB = -2252.9505 RESTRAINT = 0.0000 ------------------------------------------------------------------------- A V E R A G E S O V E R 2000000 S T E P S NSTEP = 2000000 TIME(PS) = 4060.000 TEMP(K) = 599.90 PRESS = 0.0 Etot = 1513.2011 EKtot = 2501.0782 EPtot = -987.8771 BOND = 644.5795 ANGLE = 1506.1933 DIHED = 1361.8286 1-4 NB = 398.1118 1-4 EEL = 3724.2906 VDWAALS = -363.4705 EELEC = -5968.1047 EGB = -2291.3058 RESTRAINT = 0.0000 irest = 1, ntx = 5,和irest = 0, ntx = 1 模拟结过差别还是有的，从EPtot的结果来看，后者结构变化比较激烈，前者较稳定。但是不知道在这个模拟里用哪个合适？ 5. local minimization after heating system 使用第三步minimization的输入文件。 6. equilibrate the every replica equilibrate.mdin equilibration 20 ns, every 10ps save output. equilibration &cntrl irest=0, ntx=1, nstlim=10000000, dt=0.002, irest=0, ntt=3, gamma_ln=1.0, temp0=XXXXX, ig=RANDOM_NUMBER, ntc=2, ntf=2, nscm=1000, ntb=0, igb=5, cut=999.0, rgbmax=999.0, ntpr=5000, ntwx=5000, ntwr=10000000, nmropt=1, / &wt TYPE='END' / DISANG=system_chir.dat 7. REMD simulation remd 40ns exchange every 2ps <- 这里一直有困惑，到底间隔多久交换一下比较合适？ remd.mdin remd 40ns exchange every 2ps &cntrl irest=0, ntx=1, nstlim=1000, dt=0.002, irest=0, ntt=3, gamma_ln=1.0, temp0=XXXXX, ig=RANDOM_NUMBER, ntc=2, ntf=2, nscm=1000, ntb=0, igb=5, cut=999.0, rgbmax=999.0, ntpr=100, ntwx=1000, ntwr=100000, nmropt=1, numexchg=20000, / &wt TYPE='END' / DISANG=system_chir.dat |
回复此楼» 本帖已获得的红花（最新10朵）
超人与小木虫» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
amber跑蛋白复合物动力学前，小分子用高斯如何处理
已经有7人回复- AMBER MMPBSA 计算时建立prmtop文件问题
已经有6人回复
- 新手想知道如何用Amber做蛋白质分子固定某两个或多个氨基酸距离的MD计算
已经有4人回复
- 【转帖】用Amber+Gaussian做小分子力场
已经有16人回复
» 抢金币啦！回帖就可以得到:
-
坐标深圳，诚征女友
+1/158
-
坐标北京 真诚脱单
+1/89
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收各类博士生
+1/85
-
华东师范大学 程义云 课题组招2026年博士研究生 - 有机化学、材料化学、高分子合成等
+1/82
-
东北大学数字钢铁全国重点实验室刘振宇教授课题组拟招收2026级入学博士研究生1~2名
+2/82
-
湘潭大学化学学院陈华杰教授课题组招收有机/高分子方向的博士研究生
+1/81
-
华中科技大学2026级申请考核制博士生1名-新型共轭MOF/COF合成及其超电应用方向
+1/75
-
因为雪而勾起的一些往事
+1/65
-
北京理工大学博士招生（Micro-LED、钙钛矿喷墨打印等）
+1/40
-
青岛大学 丁欣 课题组 招收2026秋化学博士1名
+1/30
-
浙江工业大学国家优青朱艺涵团队在固态电池解构与设计方向招收2026年博士生2名
+1/27
-
2026年度智能交通课题组诚招理工科背景博士
+1/24
-
长江大学武汉校区诚招新能源博士-2025
+1/9
-
【陕西师范大学】催化化学课题组2026年招收博士后/讲师/副高
+1/9
-
伦敦大学学院（University College London）机械工程系博士招生/CSC联培招生
+1/7
-
想替换掉环状DNA中心通道中的金属离子 如何替换才是正确操作
+1/4
-
法赫德国王石油与矿业大学（King Fahd University of Petroleum and Minerals-KFUPM）
+1/3
-
有没有一款可以听文献的APP
+1/2
-
上海交通大学AIMS-Lab招收AI for Science方向2026级博士生
+1/1
-
长春理工大学和西安工业大学主动光电探测成像技术重点实验室王春阳教授招收博士生
+1/1
2楼2013-07-03 09:54:58
3楼2013-07-04 09:18:13
4楼2013-07-09 14:25:43
6楼2013-07-09 14:27:19
7楼2013-10-13 15:59:59
简单回复
jackyma5楼
2013-07-09 14:25
回复
