Dear All,
In order to calculate the binding free energy of a ligand to a protein
I am using Thermodynamic Integration, i.e. making the ligand "disappear".
Because this ligand is charged (-1), at lambdas different than zero
(when the ligand is already fading) the total system charge increases until
it is +1 when lambda=1.
I am using 1 nm cutoff for the non-bonded interactions and PME. Should
I make something to account for this total system charge? Won't it induce
artifacts in my calculation?
thanks in advance,
Leandro Bortot
--
gmx-users mailing list [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [email protected].
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
