I am attempting to work through what I believe is a fairly well-trodden
path but have come a little unstuck with RESP charge-fitting so I hope
somebody can assist me...
I have a large protein which contains a non-standard residue (Asparagine
with a GlcNAc sugar attached). In order to prepare the protein for MD
simulations in Amber, I am attempting to calculate the charge
distribution of the non-standard residue then fit partial RESP charges.
To do this I first capped the Asn+GlcNac fragment with an ACE and an
NME in xleap then ran a Gaussian03 energy calculation using the input
#hf/6-31g(d) geom=connectivity Pop=MK iop(6/33=2) iop(6/42=6)
I now have a .log output file that I need to feed into RESP to get my
partial charges. I have used antechamber to do this, but I am unsure
how to restrain the partial charges on the capping residues to the
standard values as given in the Cornell 1995 forcefield paper. So
essentially my question is: How do you calculate a RESP fit whilst
imposing fixed partial charge values for specific atoms? Can
antechamber do this?
Any help greatfully appreciated.
On a related note: Does anybody know if the 1995 RESP paper by Cieplak,
Cornell and Bayly [ J. Comp. Chem. 16 (11) 1357-1377 ] exists
electronically anywhere?
Many thanks,
Thomas E. Turner
Imperial College London.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Tue Feb 17 2004 - 16:53:00 PST
