Dear Amber users,
I am using MCPB to generate prep and frcmod files for my protein with Zn2+
atom coordinated by HieHieCymCym motif.
I have successfully generated the prep and frcmod files with restrains on
all the backbone atoms (CA, H, HA, N, HN, C, O) to the AMBER parm99.dat
force field values (ChgModC).
When I check the charges generated finally, it is different from the ones
seen in $AMBERHOME/dat/leap/lib/all_amino94.lib.
When I compare the charges for CYM from
$AMBERHOME/dat/leap/lib/all_amino94.lib and
file:///usr/local/amber11/dat/mtkpp/amino94.xml, it is completely different.
amino94.xml (MCPB) all_amino94.lib (leap)
N -0.46300 -0.415700
HN 0.25200 0.271900
CA 0.03500 -0.035100
HA 0.04800 0.050800
CB 0.17552 -0.241300
HB3 0.00265 0.112200
HB2 0.00265 0.112200
SG -0.65632 -0.884400
C 0.61600 0.597300
O -0.50400 -0.567900
Why is this difference observed? Someone suggest me which values should be
used for the CYM backbone atoms (including Hs).
Also I have attached my chg2.prep file obtained from MCPB. Please check the
charges in the file and suggest if they are fine to use in simulations.
--
Thank you,
With regards,
Sangeetha
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
- application/octet-stream attachment: chg2.prep
Received on Mon Oct 01 2012 - 01:00:03 PDT
