[gmx-users] parameter for repeating units with cgenff
jalemkul at vt.edu
Fri Dec 11 19:33:27 CET 2015
On 12/11/15 1:28 PM, gromacs query wrote:
> Hi Justin
> I am bit lost I think. For e.g. in the amino acid residues library we have
> -CO-CH(R)-NH- which we can combine these amino acids in any way and in the
> rtp file it defines only these atoms and connecting atoms are mentioned
> Now lets say I have some new amino acid -CO-CH(X)-NH- which I want to
> combine with some existing amino acids. If I want to derive charges from
> Cgenff what kind of residue I should prepare?
> Is it NH2-CO-CH(X)-NH-COOH or may be CH3-NH-CO-CH(X)-NH-COCH3. But
> ultimately I will be requiring only -CO-CH(X)-NH- part. So not considering
> charges for the caps in topology will give non-integral charge or in other
> words this part -CO-CH(X)-NH- will always be non-integral. Sorry am
The backbone group always has a charge of zero, so you only deal with the side
chain when deriving charges.
Initial charge fitting would be done on a side chain analog that terminates in a
-CH3 that is analogous to CB. Once that work is done, the model compound is
merged with the backbone and CB (which is -CH2, so generally the third H charge
is just lumped into the C charge) and torsional fitting is done for
chi1/chi2/etc using a dipeptide model of the amino acid. But at that point the
charges are done.
I also wouldn't use CGenFF for a modified amino acid; derive parameters by
analogy from the parent CHARMM force field. CGenFF is generalized, at the
expense of some accuracy. Mixing CGenFF into a polypeptide chain is not the
most robust approach. CHARMM has excellent coverage for most chemical moieties.
Simply assigning charges and types by analogy should be quite straightforward.
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
More information about the gromacs.org_gmx-users