[gmx-users] Free energy perturbation protocol

[Anton Feenstra]

bmbbl at leeds.ac.uk wrote:

> Hi there,
> 
> Since nobody has replied my post, I am sending it around once more.
> 
> I want to confirm my protocol of free energy difference calculation using FEP
> with GROMACS.
> 
> A-B and A-B' complex
> 
> 1. Full MD on A-B complex. input 1.tpr; output: 1.xtc, 1.gro
> 2. perturbation topology indicated as: B --> B' names as newtp.itp
>     grompp -f FEP.mdp -o 2.tpr -c 1.gro -p newtp.itp
> 3. mdrun -c 1.tpr -rerun 1.xtc -o AB.edr
>    mdrun -c 2.tpr -rerun 1.xtc -o AB'.edr
> 4. g_energy -f AB.edr -f2 AB'.edr -fee -o FE.xvg
> 
> My question is on step 2 and 3. Do I have to run full MD using 2.tpr and 1.gro
> to generate a new 2.xtc file or I can simply apply 1.xtc to the -rerun?

No. The Free Energy Perturbation (FEP) assumes that the perturbation 
does not (significantly) affect the enemble, only the energies. So, you 
must in both cases (2 and 3) use the same trajectory (1.xtc).

However, as an alternate check, you can apply the perturbation the other 
way around as well, and use mdrun -c 2.tpr to generate '2.xtc', and then 
rerun 'mdrun -rerun 2.xtc -c 1.tpr -o 2AB.edr' and 'mdrun -rerun 2.xtc 
-c 2.tpr -o 2AB'.edr'. In principle, this must give the same answer, but 
in practice you will see a (small?) difference, which is a good 
indication of the actual error you make in your method.

By the way, I'm not sure if you have all the file options correct in 
your examples above. Check the manuals of mdrun, grompp and g_energy.

Also, as a side remark, it may not be a good idea to actually use a 
prime (') in a filename (as in AB'.edr in your example).


-- 
Groetjes,

Anton

* NOTE: New Phone & Fax numbers (below) *

  _____________ _______________________________________________________
|             |                                                       |
|  _   _  ___,| K. Anton Feenstra                                     |
| / \ / \'| | | Dept. of Pharmacochem. - Vrije Universiteit Amsterdam |
|(   |   )| | | De Boelelaan 1083 - 1081 HV Amsterdam - Netherlands   |
| \_/ \_/ | | | Tel: +31 20 59 87608 - Fax: +31 20 59 87610           |
|             | Feenstra at few.vu.nl - www.few.vu.nl/~feenstra/         |
|             | "If You See Me Getting High, Knock Me Down" (RHCP)    |
|_____________|_______________________________________________________|