[gmx-users] oplsaa force field
jalemkul at vt.edu
Mon Dec 17 13:07:14 CET 2012
On 12/16/12 11:55 PM, Shine A wrote:
> I studying the dynamics of membrane proteins.First I did a 20 ns
> simulation using GROMOS96 53a6.This force field cause some problems in the
> helical part of the protein.Now I am trying to do the same simulation with
> opls-aa force field.Is the changes I should have to do in ffnonbonded.itp
> and ffbonded,itp is same as in the justin manual?
The concept is the same, but the actual information will be very different. As
I recall, there are OPLS-AA lipid parameters somewhere out there, but I've never
used them personally. One can also implement a hybrid approach wherein the
Berger lipids can be made compatible with OPLS-AA. There's a link to that
method in my tutorial and it has been discussed on the mailing list at length.
For what it's worth, the helical instability is almost certainly due to the use
of 53A6, which is known to over-stabilize extended configurations and
Justin A. Lemkul, Ph.D.
Department of Biochemistry
jalemkul[at]vt.edu | (540) 231-9080
More information about the gromacs.org_gmx-users