[gmx-users] best parameters for energy determination with GBSA
schreib at compugen.co.il
Thu Oct 7 12:29:49 CEST 2010
Dear GROMACS users,
I try to compute energies of a small peptide-domain complex and of its
constituents using implicit solvent, so that I can get the binding Delta
E. I need the accuracy to be high enough so that I can reliably measure
the Delta Delta E between the wildtype peptide and a mutation of it.
I would like to get advice as to the best setting and parameters to use;
are the ones given below the correct and best ones (most in an mdp file,
some run control parameters of GROMACS commands)? Are other defaults in
need to be changed?
As the GBSA method necessitates using the inexact cut-off interactions,
and my system is rather small, I consider using no cut-off at all, i.e.
an all-vs.-all approach. For that I understand the parameters should be
-pd ; in mdrun
Then, there are the implicit solvent parameters themselves:
Finally, regarding the simulation box, I see no need for periodic
boundary conditions, and I pick some big margin around the system (as
there are no water molecules, no need to spare):
-d 5.0 ; in editconf
What do you say?
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gromacs.org_gmx-users