AMBER Archive (2009)

Subject: RE: [AMBER] glycopeptide Parametrization

From: FyD (
Date: Sat Jun 13 2009 - 02:40:00 CDT

Dear Ross & Carlo,

>> Please let me know if this is what you meant. Which charges I will be
>> using with this parametrization ? Charges based on the Connolly surface
>> for the amino-acid residues and charges based on the CHELPG algorithm
>> for the sugar ?
> I would use RESP charges for everything. See: W.D. Cornell, P. Cieplak, C.I.
> Bayly & P.A. Kollman J. Am. Chem. Soc. 1993, 115, 9620-9631.

Just few additional comments:

- For the "peptide" part use RESP charge derivation approach as in
AMBER force fields: This means use the Connolly surface + 2 RESP
stages with qwt=.0005 & qwt=.001.
See also P. Cieplak, W.D. Cornell, C.I. Bayly & P.A. Kollman J.
Comput. Chem. 1995, 16, 1357-1377.

- For the "sugar" part use RESP charge derivation as in the GLYCAM
force field: This means use the CHELPG algo. + 1 RESP stage with qwt =
See R.J. Woods & R. Chappelle. J. Mol. Struct. (THEOCHEM) 2000, 527, 149-156

For additional references, see:

All this is automatically set up in R.E.D. and/or R.E.D. Server. You
just pick up the "RESP-A1A" or "RESP-C2" charge models:

-1- RESP-A1A: HF/6-31G* Connolly surface algo., 2 stage RESP fit
-5- RESP-C2: HF/6-31G* CHELPG algo., 1 stage RESP fit qwt=.01

This means two different charge derivation approaches for the two
different parts of your molecule (NO am1-bcc and NO gaff):

regards, Francois

AMBER mailing list