[gmx-users] g_covar + g_anaeig entropy calculation
mailmd2011 at gmail.com
Thu Dec 4 11:07:04 CET 2014
I am trying to calculate ligand binding energy for protein/ligand
system. The enthapy part was done by GMXPBSA. Now I would like to
evalutate the entropy part with g_covar + g_anaeig.
Here I got some questions:
--- is it necessary using double precison Gromacs for entropy evaluation?
--- I googled a lot concering on this topic from Gromacs mailist. If I
understand correctly, we should
(1) generate eigenvec.trr as input file for g_anaeig
(2) use g_anaeig calculate the entropy
is it correct?
---I am trying to run command:
g_covar -f input.xtc -s input.tpr -n -o eigenval.xvg -v eigenvec.trr -av
average.pdb -l covar.log -ascii covar.dat -xpm covar.xpm -xpma covara.xpm
The trajecotry I am using for above commend was extracted from MD
simulations and water/ions were removed from the system. When g_covar
pop up a dialoug asking which compoenent I am going to analyze, I
selected "System" which contains protien+ligand
Is it correct to do so?
---I am planning to run g_anaeig with following command for entropy
g_anaeig -v eigenvec.trr -entropy -s input.tpr -comp
Shall I speicify additional options for this purpose?
Thank you very much.
More information about the gromacs.org_gmx-users