[gmx-users] Fwd: free energy calculation.

neelam wafa neelam.wafa at gmail.com
Tue Aug 28 01:38:32 CEST 2018

Hi everyone.
Kindlt need your help urgently.

---------- Forwarded message ---------
From: neelam wafa <neelam.wafa at gmail.com>
Date: Mon, 27 Aug 2018, 3:41 pm
Subject: free energy calculation.
To: <gmx-users at gromacs.org>, <gromacs.org_gmx-users at maillist.sys.kth.se>

Dear gmx users,
I have md simmulations of a protein with four different ligands. now I want
to calculate the free energy of binding. which method does suit after md
simmulation?  either g_mmpbsa or gmx BAR  method.  If I use BAR method is
it necessary to set charges to zero as in the tutorial given at

I s necessary to turn off couloumb interactions for calculating free energy
of ligand protein complex. if i want to use coloumb interactions also then
should i set the

coul_lambdas values in the same way as vdw_lambdas.

Another question is that there is another tutorial available at

In this tutorial the system has been equilliberated before generating
the lambda states.

while in the justin tutorial each lambda state is equiliberated separately.

kindly tell me which one is more suitable for a protein ligand complex.

looking forward for your cooperation.

sorry if the questions are trivial.



Ph.D scholar

More information about the gromacs.org_gmx-users mailing list