[gmx-users] Choosing terminals with AMBER ff
jalemkul at vt.edu
Wed Nov 8 12:57:39 CET 2017
On 11/8/17 6:25 AM, Simon Kit Sang Chu wrote:
> Dear Mark,
> Thank you very much. I am worried with the ends attracting each other for a
> long peptide. May I ask for the general strategy simulating segments of a
> protein? I believe simulating part of the whole protein might not be
> something new.
You have to do like Mark suggested and manually build an ACE cap on the
N-terminus if you want it to be uncharged. There are also NME/NHE
termini for the C-terminus to be modeled as an amide. Indeed, this isn't
a new concept. But AMBER defaults to assuming specific N- and C-terminal
versions of their amino acids. The user has to go the extra step of
building proper termini and naming them as such. All force fields
support capping groups, and they are routinely used, but some force
fields are easier to deal with via normal .n.tdb and .c.tdb entries.
Justin A. Lemkul, Ph.D.
Virginia Tech Department of Biochemistry
303 Engel Hall
340 West Campus Dr.
Blacksburg, VA 24061
jalemkul at vt.edu | (540) 231-3129
More information about the gromacs.org_gmx-users