[gmx-users] Umbrella sampling restraining potential
agnivogromacs14 at gmail.com
Sun Feb 15 01:56:52 CET 2015
The PMF generated by umbrella sampling of my system appears to be converged
and the calculated free energy of binding is close to that in published
I believe a necessary check for the convergence is the overlap and
smoothness of the histograms generated from sampling of the reaction
coordinate. Earlier I discussed in the forum regarding my histogram spacing
Dr. Lemkul had then pointed out the rocky nature of many histograms and the
previous post :
For my system there is no experimental study and I am resorting to trial
and error. My pull force constant is 1000 kJ/mol/nm^2 whereas the only
literature doing a MD study of my system ( with slightly different pulling
coordinates) uses a very high force constant of 12600 kJ/mol/nm^2 and they
have provided no explanation for this.
I have increased simulation time of my windows from 1 ns to 10 ns and I am
thinking of using a higher force constant in one half of my windows where
the pulled group is closer to the reference group and use a lesser
potential in the rest of the windows.
However I am declined to use the 12600 units of force constant as it is
very high and there is no validation. I found some literature using 1000 -
10,000 kJ/mol/nm^2 force constants for proteins.
Any suggestions please ?
Thanks & Regards
Grad Student, Iowa State University.
More information about the gromacs.org_gmx-users