输出体系的整个带电量:
measure sumweights $all weight charge给PDB文件设置周期边界条件:
pbc set {54 54 24 } -all将此晶胞内原子脱除周期性
pbc unwrap显示晶胞边界
pbc box体系总质量:
measure sumweights $all weight mass体系总电荷:
measure sumweights $all weight charge 添加周期性水环境package require solvatesolvate ubq.psf ubq.pdb -t 5 -o ubq wb-t (override with any of the following)-o (data will be written to output.psf/output.pdb)添加离子中和多余电荷autoionize -psf file.psf -pdb file.pdb [options]可以直接使用VMD中extension>modeling>add ions的autoionize完成TK console中选择全部原子,保存psf与pdb文件参考:Autoionize Plugin命令模式运行ubq.pgn> vmd -dispdev text -e ubq.pgn測量周期最大最小边界set everyone [atomselect top all]measure minmax $everyone測量周期中心点set everyone [atomselect top all]measure center $everyone測量质量中心measure center $sel weight mass位置移动atomselect0 moveby {1 1 6} # 把所选原子向1,1,6向量方向和距离上移动atomselect0 moveto { 3 6 5} # 把所选内容移动到3,6,5位置原文:VMD的TK Console中的内置命令设置resname/chain/residmol load pdb fileA.pdbset sel [atomselect top "serial <=9"]$sel set resname CD$sel set chain X$sel set resid 8set all [atomselect top all]$all writepdb fileB.pdbrefer the VMD Graphical Representations>Selections>Keyword for more valuestop文件IC A B C D [bond(AB)] [angle(ABC)] [dihedral(ABCD)] [angle(BCD)] [bond(CD)]IC A B *C D [bond(AC)] [angle(BCA)] [improper(ABCD)] [angle(BCD)] [bond(CD)].The * next to the C atom indicates that it is at the center of an improper angle definition. No * indicates that the 4 atoms do not have an improper topology.Specifying IC’s is not necessary if you already have all the atoms in your PDB file. The entries are simply thereto have a way to create missing atoms from the positions of present ones.参考:NAMD topology-tutorialVMD键长单位number in the “Value” field corresponds to the length of the bond in ?ngstroms.reference:Using VMD - An Introductory Tutorialpsfgen命令pdbalias residue HOH TIP3 #aliasing residue HOH to TIP3pdbalias atom TIP3 O OH2 #aliasing residue TIP3 atom O to OH2vmd中tk命令打开文件mol new filename.psf #打开一个文件mol addfile filename.pdb #打开另一个文件,并叠加到第一个文件结构上最适分子数Presently, the maximum number of atoms one can expect to realistically compute with high level quantum chemistry is approximately 120 atoms平衡态模拟常见的模拟思路是,先在NVT下约束住你的溶质(剂)做限制性模拟,这是一个升温的过程,当温度达到你的设定后, 接着做NPT模拟,此过程将调整体系的压强进而使体系密度收敛。经过一段时间的平衡模拟,在确定系统弛豫已经完全消除之后,就可以开始取数据了。如何判断体系达到平衡,简单的讲可以通过以 下几种方式,一,看能量(势能,动能和总能)是否收敛;二,看系统的压强,密度等等是否收敛;三看系统的RMSD是否达到你能接受的范围,等等。升温/降温控制# IF HeatingreassignFreq 2000reassignTemp 100reassignIncr 1reassignHold 300#从100K升温到300K,每2000步升高1K# IF Cooling
reassignFreq 2000reassignTemp 300reassignIncr -1reassignHold 100 固定分子all atoms with a value of 1 (or a number different of 0) in a predetermined column will be fixed; atoms with a value of 0 in the same column will not be affected.load pdb file into the psf file in vmdin TKconset allatoms [atomselect top all]$allatoms set beta 0set fixedatom [atomselect top "resid 1 and name CA"]$fixedatom set beta 1$allatoms writepdb fileName.ref拉动原子to set which atom is to be pulled (SMD atom).uses the occupancy column of the pdb file to distinguish itset allatoms [atomselect top all]$allatoms set occupancy 0set smdatom [atomselect top "resid 76 and name CA"]$smdatom set occupancy 1$allatoms writepdb fileName.ref牵拉方向set smdpos [lindex [$smdatom get {x y z}] 0]set fixedpos [lindex [$fixedatom get {x y z}] 0]vecnorm [vecsub $smdpos $fixedpos] VMD计算SASA(solvent-accessible surface area )set all [atomselect top "all"]
set some [atomselect top "resid 1 to 5"]measure sasa 1.4 $all -restrict $some -points sasapointsforeach pt $sasapoints { draw point $pt}以上可以在窗口中显示出1-5号氨基酸的溶剂可及表面,以point绘出。当然如果不需要浏览仅仅计算的话,只需要:set all [atomselect top "all"]set some [atomselect top "resid 1 to 5"]measure sasa 1.4 $all -restrict $some如果想看一下准不准可以算一下除了1-5号氨基酸的是多少set someothers [atomselect top "not(resid 1 to 5)"]measure sasa 1.4 $all -restrict $someothers将两个值相加,看是否是measure sasa 1.4 $all得到的结果。原文摘录自:关于溶剂可及表面积(Solvent Accessible Surface)- Biocheming
try VMD timeline Tool to calculate the change in SASA over time.Load the molecule and trajectory and then go to:Extensions -> Analysis -> Timeline.SASA is one of the many parameters you can calculate in the Timeline window.2014-9-10原文摘录自:How to calculate solvent accessible surface (SASA)?-Uttam Pal[笔记]VMD/NAMD命令/规则此博文包含图片 (2013-10-10 16:18:42)转载▼标签:vmdnamd笔记分类: 分子模拟添加周期性水环境package require solvatesolvate ubq.psf ubq.pdb -t 5 -o ubq wb-t (override with any of the following)-o (data will be written to output.psf/output.pdb)添加离子中和多余电荷autoionize -psf file.psf -pdb file.pdb [options]可以直接使用VMD中extension>modeling>add ions的autoionize完成TK console中选择全部原子,保存psf与pdb文件参考:Autoionize Plugin命令模式运行ubq.pgn> vmd -dispdev text -e ubq.pgn測量周期最大最小边界set everyone [atomselect top all]measure minmax $everyone測量周期中心点set everyone [atomselect top all]measure center $everyone測量质量中心measure center $sel weight mass位置移动atomselect0 moveby {1 1 6} # 把所选原子向1,1,6向量方向和距离上移动atomselect0 moveto { 3 6 5} # 把所选内容移动到3,6,5位置原文:VMD的TK Console中的内置命令设置resname/chain/residmol load pdb fileA.pdbset sel [atomselect top "serial <=9"]$sel set resname CD$sel set chain X$sel set resid 8set all [atomselect top all]$all writepdb fileB.pdbrefer the VMD Graphical Representations>Selections>Keyword for more valuestop文件IC A B C D [bond(AB)] [angle(ABC)] [dihedral(ABCD)] [angle(BCD)] [bond(CD)]IC A B *C D [bond(AC)] [angle(BCA)] [improper(ABCD)] [angle(BCD)] [bond(CD)].The * next to the C atom indicates that it is at the center of an improper angle definition. No * indicates that the 4 atoms do not have an improper topology.Specifying IC’s is not necessary if you already have all the atoms in your PDB file. The entries are simply thereto have a way to create missing atoms from the positions of present ones.参考:NAMD topology-tutorialVMD键长单位number in the “Value” field corresponds to the length of the bond in ?ngstroms.reference:Using VMD - An Introductory Tutorialpsfgen命令pdbalias residue HOH TIP3 #aliasing residue HOH to TIP3pdbalias atom TIP3 O OH2 #aliasing residue TIP3 atom O to OH2vmd中tk命令打开文件mol new filename.psf #打开一个文件mol addfile filename.pdb #打开另一个文件,并叠加到第一个文件结构上最适分子数Presently, the maximum number of atoms one can expect to realistically compute with high level quantum chemistry is approximately 120 atoms平衡态模拟常见的模拟思路是,先在NVT下约束住你的溶质(剂)做限制性模拟,这是一个升温的过程,当温度达到你的设定后, 接着做NPT模拟,此过程将调整体系的压强进而使体系密度收敛。经过一段时间的平衡模拟,在确定系统弛豫已经完全消除之后,就可以开始取数据了。如何判断体系达到平衡,简单的讲可以通过以 下几种方式,一,看能量(势能,动能和总能)是否收敛;二,看系统的压强,密度等等是否收敛;三看系统的RMSD是否达到你能接受的范围,等等。升温/降温控制# IF HeatingreassignFreq 2000reassignTemp 100reassignIncr 1reassignHold 300#从100K升温到300K,每2000步升高1K# IF Cooling
reassignFreq 2000reassignTemp 300reassignIncr -1reassignHold 100 固定分子all atoms with a value of 1 (or a number different of 0) in a predetermined column will be fixed; atoms with a value of 0 in the same column will not be affected.load pdb file into the psf file in vmdin TKconset allatoms [atomselect top all]$allatoms set beta 0set fixedatom [atomselect top "resid 1 and name CA"]$fixedatom set beta 1$allatoms writepdb fileName.ref拉动原子to set which atom is to be pulled (SMD atom).uses the occupancy column of the pdb file to distinguish itset allatoms [atomselect top all]$allatoms set occupancy 0set smdatom [atomselect top "resid 76 and name CA"]$smdatom set occupancy 1$allatoms writepdb fileName.ref牵拉方向set smdpos [lindex [$smdatom get {x y z}] 0]set fixedpos [lindex [$fixedatom get {x y z}] 0]vecnorm [vecsub $smdpos $fixedpos]run a NAMD simulationnamd2 +p[procs] configfile > outfile力场中二面角倒序正序区别软件里的二面角有正负是得考虑4号原子在123号原子平面的哪一侧,举MS里二面角的来说,按顺序点1234共4个原子,二面角指的是123平面和234 平面的二面角,沿着23的轴看过去,1号原子在4号原子的逆时针方向(左手边),则取正值,反之取负值;绝对值仍在[0,180]。原文摘录自:二面角和扭转角区别Reversing the atoms in a proper dihedral does not have any mathematical or practical consequences. Conversely, the ordering of the atoms is very important for improper dihedrals.原文摘录自:Parameters for a Protein-RNA covalent bond力场参数图示
[笔记]VMD/NAMD命令/规则2014-06-18图片来源:Brève introduction à la mécanique (MM) et à la dynamique moléculaire (DM)improper
dihedralBy convention, the first atom of an improper
dihedral (type A-X-X-B or A-B-C-D) is usually the central atom. This hadbeen a general rule in the past.2014-06-19
原文摘录自:CHARMM c32b2 parmfile.doc
拓扑文件AUTOGENERATE ANGLES DIHEDRALAUTOgenerate default options to be used when building a sturcture.AUTO ANGLes specifies that all possible angles and DIHEdral specifies that all possible dihedral angles be generated when building a structure. If these options are not included the angles and/or dihedrals must be listed explicitly in the topology file2014-06-19原文转自:CHARMM c38b1 rtop.docnamd开始NPT计算过程报错Periodic cell has become too small for original patch grid解决方法
I have a different approach to fixing this problem. I believe that the underlying cause of this error is that the system is resizing too quickly (due to the system pressure being far from the specified pressure). To slow the resizing down I increase langevinPistonPeriod to 1000 and langevinPistonDecay to 500 for a short time when I start the NPT simulation, and then switch to the normal values of 100 and 50.For example, in your case I would do one input file for 100 ps with a constraintscaling of 10, langevinPistonPeriod of 1000, and a langevinPistonDecay of 500. Then I would load the restart files into a new simulation and proceed as you were above. Also, by restarting the simulation you are forcing NAMD to re-assign the patches, which will help prevent this error.
2014-08-26原文摘录自:Re: Periodic cell has become too small for original patch grid!FATAL ERROR: Periodic cell has become too small for original patch grid!问题出现在nvt计算完成后第一次进行npt的时候.通过在namd的邮件列表寻找,最后发现是周期性盒子设置的时候四舍五入数值,造成盒子在npt计算 过程中溶剂环境密度过小(盒子体积过大,而分子总数一定,最终水溶液密度过低)。解决方法,周期盒子大小不能直接四舍五入,要考虑到盒子体积与分子总数得 到溶液的密度是否恰当好处。改用该方法后,顺利计算,为提示错误。2014-9-10该内容原创计算namd体系密度
首先计算体系的质量,然后根据namd输出的体积,两者相除即可。过程如下:例如:C 11H 1596O 790N 2上面表示的是每种原子的个数则mass=(11*12+1596*1+790*16+2*14)/NA 单位为克。NA为阿伏伽德罗常数:6.02214129(27)×1023那么密度为mass/V2014-08-26原文摘录自:heating过程中水盒子形状发生变化 VMD计算SASA(solvent-accessible surface area )控制台下用命令即可得出:set all [atomselect top "all"]set some [atomselect top "resid 1 to 5"]measure sasa 1.4 $all -restrict $some -points sasapointsforeach pt $sasapoints { draw point $pt}以上可以在窗口中显示出1-5号氨基酸的溶剂可及表面,以point绘出。当然如果不需要浏览仅仅计算的话,只需要:set all [atomselect top "all"]set some [atomselect top "resid 1 to 5"]measure sasa 1.4 $all -restrict $some如果想看一下准不准可以算一下除了1-5号氨基酸的是多少set someothers [atomselect top "not(resid 1 to 5)"]measure sasa 1.4 $all -restrict $someothers将两个值相加,看是否是measure sasa 1.4 $all得到的结果。2014-9-10原文摘录自:关于溶剂可及表面积(Solvent Accessible Surface)- Biochemingtry VMD timeline Tool to calculate the change in SASA over time.Load the molecule and trajectory and then go to:Extensions -> Analysis -> Timeline.SASA is one of the many parameters you can calculate in the Timeline window.2014-9-10