NAMD靶向分子动力学模拟Target molecular dynamic stimulation（附conf文件）

下面附上靶向分子动力学的conf文件，具体问题具体分析，参数根据自己体系进行修改。TMD比常规MD的conf不同就在于EXTRA PARAMETERS 这部分多了一些参数。具体设置可参考上述TMD介绍。

conf:

#############################################################
## JOB DESCRIPTION ##
#############################################################

# Min. and Eq. of KcsA
# embedded in POPC membrane, ions and water.
# Protein constrained. PME, Constant Pressure.

#############################################################
## ADJUSTABLE PARAMETERS ##
#############################################################

structure p-gpoutward_popc_water_ion.psf
coordinates p-gpoutward_popc_water_ion.pdb
outputName p-gpoutward_popc_water_ion_tmd

set temperature 310

# Continuing a job from the restart files
if {1} {
set inputname p-gpoutward_popc_water_ion_eqm2
binCoordinates $inputname.restart.coor
binVelocities $inputname.restart.vel ;# remove the "temperature" entry if you use this!
extendedSystem $inputname.restart.xsc
}

firsttimestep 1775000

#############################################################
## SIMULATION PARAMETERS ##
#############################################################

# Input
paraTypeCharmm on
parameters par_all27_prot_lipidNBFIX.prm

# NOTE: Do not set the initial velocity temperature if you
# have also specified a .vel restart file!
# temperature $temperature

# Periodic Boundary Conditions
# NOTE: Do not set the periodic cell basis if you have also
# specified an .xsc restart file!
if {0} {
cellBasisVector1 164. 0. 0.
cellBasisVector2 0. 167. 0.
cellBasisVector3 0. 0. 180.
cellOrigin -0.40401795506477356 -1.9578770399093628 7.172184467315674
}
wrapWater on
wrapAll on

# Force-Field Parameters
exclude scaled1-4
1-4scaling 1.0
cutoff 12.
switching on
switchdist 10.
pairlistdist 13.5

# Integrator Parameters
timestep 1.0 ;# 2fs/step
# rigidBonds all ;# needed for 2fs steps
nonbondedFreq 1
fullElectFrequency 2
stepspercycle 20

#PME (for full-system periodic electrostatics)
if {1} {
PME yes
PMEGridSizeX 145
PMEGridSizeY 145
PMEGridSizeZ 183
}

# Constant Temperature Control
langevin on ;# do langevin dynamics
langevinDamping 1 ;# damping coefficient (gamma) of 5/ps
langevinTemp $temperature

# Constant Pressure Control (variable volume)
if {1} {
useGroupPressure yes ;# needed for 2fs steps
useFlexibleCell yes ;# no for water box, yes for membrane
useConstantArea yes ;# no for water box, yes for membrane

langevinPiston on
langevinPistonTarget 1.01325 ;# in bar -> 1 atm
langevinPistonPeriod 400.
langevinPistonDecay 200.
langevinPistonTemp $temperature
}

restartfreq 500 ;# 1000steps = every 2ps
dcdfreq 500
xstFreq 500
outputEnergies 500
outputPressure 500

# Fixed Atoms Constraint (set PDB beta-column to 1)
if {0} {
fixedAtoms on
fixedAtomsFile nottails.fix.pdb
fixedAtomsCol B
fixedAtomsForces on
}

#############################################################
## EXTRA PARAMETERS ##
#############################################################

# Put here any custom parameters that are specific to
# this job (e.g., SMD, TclForces, etc...)

# constraints on
# consexp 2
# consref p-gpoutward_popc_water_ion.pdb
# conskfile p-gpoutward_popc_water_ion_ca.pdb
# conskcol B
# margin 3

tclforces on
set waterCheckFreq 100
set lipidCheckFreq 100
set allatompdb p-gpoutward_popc_water_ion.pdb
tclForcesScript keep_water_out.tcl

# eFieldOn yes
# eField 0 0 -0.155

TMD on
TMDFile p-gpinward_tmd_nbd.pdb
TMDk 515
TMDOutputFreq 500
TMDFinalRMSD 0
TMDFirstStep 1775000
TMDLastStep 6775000
#############################################################
## EXECUTION SCRIPT ##
#############################################################

# Minimization
if {0} {
minimize 20000
reinitvels $temperature
}

run 5000000 ;# 5 ns

NAMD靶向分子动力学模拟Target molecular dynamic stimulation（附conf文件）相关推荐

NAMD跑分子动力学模拟出现的一些问题（更新中）
1.能量最小化平衡的时候,如果没有控制好NVT,会出现盒子变小,蛋白跑出水盒子.但是如果仅仅是控制体积不变,又会造成盒子内溶剂密度太小,产生'Periodic cell has become too ...
干货分享 | 分子对接与分子动力学模拟在药物研发中的应用
前言分子对接(Molecular docking)与分子动力学模拟(Molecular dynamics simulation)是计算生物学中重要的一部分,在生物学研究中不断发挥着重要的作用.分子对 ...
[分子动力学模拟资料]几种常用力场参数网站及LJ参数计算
原创一直陪着你的 LAMMPS交流站我们在使用分子动力学模拟时,常常需要力场参数或势函数文件,本文就推荐几个常用的EAM势＼MEAM势及LJ势的力场参数网站,欢迎大家补充. 力场参数常用网站 ...
vasp 模拟退火_[编] 使用vasp进行分子动力学模拟
vasp做分子动力学的好处,由于vasp是近些年开发的比较成熟的软件,在做电子scf速度方面有较好的优势.缺点:可选系综太少.尽管如此,对于大多数有关分子动力学的任务还是可以胜任的.主要使用的系综是 ...
vasp 模拟退火_科学网—vasp的分子动力学模拟 - 王达的博文
vasp做分子动力学的好处,由于vasp是近些年开发的比较成熟的软件,在做电子scf速度方面有较好的优势. 缺点:可选系综太少. 尽管如此,对于大多数有关分子动力学的任务还是可以胜任的. 主要使用的系 ...
分子动力学模拟软件_机器学习模拟1亿原子：中美团队获2020「超算诺贝尔奖」戈登贝尔奖...
在前世界第一超算 Summit 上,研究人员在保持「从头算」精度的前提下成功模拟了 1 亿原子的运动轨迹,将超大系统的分子动力学模拟带进了一个新时代. 机器之心报道,作者:泽南.张倩.小舟. 有超算界 ...
AMBER分子动力学模拟之结果分析(构象分析)-- HIV蛋白酶-抑制剂复合物(6)
AMBER分子动力学模拟之结果分析(构象分析)-- HIV蛋白酶-抑制剂复合物(6) RMSD RMSF b-facto计算 RMSD RMSD measures the deviation of a ...
计算机辅助药物设计 fda,计算机辅助药物设计中的分子动力学模拟.pdf
计算机辅助药物设计中的分子动力学模拟计算机辅助药物设计中的分子动力学模拟肖旭东一.分子动力学概述随着生命科学理论和计算分析方法的快速发展,新药研究已经进入一个全新的时期,计算机辅助药物设计正 ...
amber软件的分子动力学模拟步骤
进行分子动力学模拟的步骤: 1.准备模拟需要的文件: autodock对接完成得分高的复合物构象. 文本编辑器内对构象中分子分离,得到: 1): 带有受体和配体的pdb文件 2): 不含有配体的受体p ...

NAMD靶向分子动力学模拟Target molecular dynamic stimulation（附conf文件）

NAMD靶向分子动力学模拟Target molecular dynamic stimulation（附conf文件）相关推荐

最新文章

热门文章