[gmx-users] Higher that usual salt bridge occupancies in Amber99sb*-ildn forcefield despite adding nonbonded corrections
priyeshmohanty at yahoo.in
Mon Apr 17 16:06:37 CEST 2017
Well, yes. I have been preparing .tpr files in the exact same way for the last 6 months or so when I was getting the correct occupancies with these corrections. There were no errors whatsoever. In fact, when I delete the non-bonded salt bridge corrections from ffnonbonded.itp, I get considerably stronger salt bridges as one would expect from this particular forcefield. Hence, I really don't understand how grompp would suddenly start messing up when it comes to reading the topology file. The salt bridges occupancies that I get as of now are essentially stronger than they used to be but nowhere near as strong as they are without the corrections. This subtle increase however appears to considerably alter the dynamics for some multidomain proteins for which I previously performed simulations with these corrections.
On Monday, 17 April 2017 5:13 PM, Justin Lemkul <jalemkul at vt.edu> wrote:
On 4/17/17 6:09 AM, Priyesh Mohanty wrote:
> I've also tried increasing the sigma to higher values beyond 0.325 which should in principle considerably weaken salt bridges but there appears to be no decrease in salt bridge occupancy despite the increase in sigma. It is as though mdrun is not reading the nonbonded_param addition correctly.
mdrun doesn't read these files at all; grompp does when building your .tpr. So
if you're making adjustments to ffnonbonded.itp after writing a .tpr file, there
will be no effect on the simulation because you haven't actually changed the
physical model. Are you correctly re-creating the .tpr file each time?
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
jalemkul at outerbanks.umaryland.edu | (410) 706-7441
Gromacs Users mailing list
* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.
More information about the gromacs.org_gmx-users