Hello
I am trying to parametrize an unnatural amino acid-phenylglycine
(PHG)- in the OPLS force field in GROMACS for an MD simulation. In
the ffoplsaa.hdb file, the following is the format for the hydrogen
addition that I have put in:
PHG
7
1 1 N -C CA
1 5 CA N C CG
1 1 CD1 CG CE1
1 1 CD2 CG CE2
1 1 CE1 CD1 CZ
1 1 CE2 CD2 CZ
1 1 CZ CE1 CE2
As can be seen, I have added the correct tetrahedral hydrogen of
type-5 in the second line (CA N C CG). Yet during the MD, the carbon
fails to remain tetrahedral, with especially the hydrogen moving
around a lot (it is staying bonded but the C CA H angle is getting
highly distorted). What can I do to remedy this?
Thanks very much
Ashutosh Jogalekar
Emory University Chemistry Department
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
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/mailing_lists/users.php