rmsf和PCA的分析怎么做？

胡姐姐

我做的是关于蛋白质的模拟，现在轨迹跑出来了。老师让我做一下rmsf和PCA分析，请问，大家都是怎么做出来的？纠结好久了。看网上有关于vmdICE可以进行rmsf的计算，但是window  size 死活搞不懂  
求大神指点啊

霜晨月

我一般用Gromacs，rmsf 可以用 gmx rmsf 命令，PCA的话，先用gmx covar生成本征向量和本征值，然后用gmx anaeig 进行分析。

胡姐姐

霜晨月 发表于 2018-1-5 16:52
我一般用Gromacs，rmsf 可以用 gmx rmsf 命令，PCA的话，先用gmx covar生成本征向量和本征值，然后用gmx an ...

谢谢你啦

胡姐姐

霜晨月 发表于 2018-1-5 16:52
我一般用Gromacs，rmsf 可以用 gmx rmsf 命令，PCA的话，先用gmx covar生成本征向量和本征值，然后用gmx an ...

请问，文件格式不对，怎么办？用什么转化呢？还是用gromacs重做一遍吗

霜晨月

用catdcd可以将dcd转换成gromacs所用的trr格式

霜晨月

我的方法可能太笨了，论坛上应该有高人知道如何直接分析NAMD的轨迹

胡姐姐

霜晨月 发表于 2018-1-5 22:31
我的方法可能太笨了，论坛上应该有高人知道如何直接分析NAMD的轨迹

我觉得这个方法已经很好啦，还有更好的方法

胡姐姐

霜晨月 发表于 2018-1-5 22:29
用catdcd可以将dcd转换成gromacs所用的trr格式

我现在已经将轨迹文件转换为trr形式，可是，如何将xst转化为gromacs的xtc？

fhh2626

分子模拟软件不是黑箱，做科研也不是学会写个NAMD的输入文件就可以了。
你算一个量之前，至少要明白这个量是什么意思：
When a dynamical system fluctuates about some well-defined average position, the RMSD from the average over time can be referred to as the RMSF or root mean square fluctuation.
你如果知道RMSF就是一种RMSD，就知道用VMD自带的RMSD Visualizer点两下鼠标就可以算

PCA也一样，用NMWiz就可以算

霜晨月

胡姐姐 发表于 2018-1-10 10:31
我现在已经将轨迹文件转换为trr形式，可是，如何将xst转化为gromacs的xtc？

Gromacs是用不到xst文件的，有trr和pdb就够了。
但是你将dcd转成trr之前，要先修正了PBC。否则，没有tpr文件，gromacs是无法修正PBC的。

胡姐姐

霜晨月 发表于 2018-1-10 11:58
Gromacs是用不到xst文件的，有trr和pdb就够了。
但是你将dcd转成trr之前，要先修正了PBC。否则，没有tpr ...

请问，你有什么资料或文献推荐吗？大神，我还是没有听懂您说的

霜晨月

胡姐姐 发表于 2018-1-10 12:54
请问，你有什么资料或文献推荐吗？大神，我还是没有听懂您说的

虽然我不是什么大神，但已经说得很明白了啊
用 gmx trjconv 可以把trr转成xtc。用不到xst文件的。详情可以看看gmx trjconv -h 和gmx covar、gmx anaeig的帮助内容。

不知道什么是PBC的话，可能就得详细看看NAMD或GMX的manual了。

胡姐姐

霜晨月 发表于 2018-1-10 13:23
虽然我不是什么大神，但已经说得很明白了啊
用 gmx trjconv 可以把trr转成xtc。用不到xst文件的 ...

在我这里，您就是大神的！我知道大神说的很明白了，由于我平时知识积累的太少，所以没办法理解大神说的  ，希望大神可以多多帮忙解惑，拜托啦

胡姐姐

fhh2626 发表于 2018-1-10 10:43
分子模拟软件不是黑箱，做科研也不是学会写个NAMD的输入文件就可以了。
你算一个量之前，至少要明白这个量 ...

那我是在算rmsf之前，先算一下整个轨迹的平均结构是吗？

胡姐姐

胡姐姐 发表于 2018-1-10 13:32
在我这里，您就是大神的！我知道大神说的很明白了，由于我平时知识积累的太少，所以没办法理解大神说的{: ...

gmx grompp -f npt-nopr-md.mdp -c npt-pr.gro -p fws.top -o npt-nopr.tpr。  这是教程中指出的生成tpr文件。
1.我知道PBC是周期性边界条件，改的是边长和中心坐标吗？
2.npt.mdp是把续跑的这个文件
structure              all_wb_ion.psf
coordinates            all_wb_ion.pdb

set init_temp         100
set desired_temp      300
#temperature          $init_temp
set outputname        all-2nd
# Continuing a job from the restart files
binCoordinates     all-1st.restart.coor
binVelocities      all-1st.restart.vel  ;# remove the "temperature" entry if you use this!
extendedSystem          all-1st.xsc
firsttimestep      21520000
paraTypeCharmm         on
parameters             par_all36_prot.prm
parameters             toppar_water_ions.str
parameters             sub.final.str

exclude                scaled1-4
1-4scaling             1.0
cutoff                 12.0
switching              on
switchdist             10.0
pairlistdist           14.0

timestep               2.0
rigidBonds             all
nonbondedFreq          1
fullElectFrequency     1
stepspercycle          20

wrapAll                on
wrapWater              on
wrapNearest            on

PME                    yes
PMEGridSizeX           72
PMEGridSizeY           72
PMEGridSizeZ           72

outputName             $outputname

restartfreq            2000     ;# 2000steps = every 1ps
dcdfreq                2000
xstFreq                2000
outputEnergies         500
outputPressure         500

langevin               on   ;# do langevin dynamics
langevinDamping        5    ;# damping coefficient (gamma) of 5/ps
langevinTemp           $desired_temp
langevinHydrogen       on    ;# don't couple langevin bath to hydrogens

useGroupPressure       yes ;# needed for rigidBonds
useFlexibleCell        no
useConstantArea        no
langevinPiston         on
langevinPistonTarget   1.01325 ;#  in bar -> 1 atm
langevinPistonPeriod   200
langevinPistonDecay    500
langevinPistonTemp     $desired_temp
run                    20000000

保存为mdp吗？
3.coord坐标文件改为gro？
4.拓扑文件是将 all_wb_ion.psf  改为top  对嘛？