I use DFT and also HF,too.There exits the repulsive reaction between them,
while force field can not show this repulsive reaction.
I just calculate the binding energy between protein and drug.
I'm going to do a systemic study on this to show the E = E(R) profile.
I just wanna to ask whether it may be a problem of Amber force field
according to other users' experience, and whether it is a good project
otherwise it will waste my time.
Thank for anyone's suggestion.
He Xiao
I would need to know the van der Waals overlap in your system.
Neither DFT nor HF should be used for systems with dispersion
such as stacking, intercalation, etc.. They fail completely.
cf. for example.
Structure, energetics, and dynamics of the nucleic acid base pairs: Nonempirical ab initio calculations
Hobza P, Sponer J
Chemical Reviews
99, 3247-3276, 1999
Base stacking in cytosine dimer. A comparison of correlated ab initio
calculations with three empirical potential models and density functional theory calculations
Sponer J, Leszczynski J, Hobza P
J. Comput. Chem.
17, 841-850 1996
If this is a pure H-bonding than the force field result would
be strange, but I have never seen a failure of this kind for AMBER.
It would be worth to analyze. What is the Lennard Jones component
of your AMBER interaction energy?
Best wishes, Jiri Sponer
-------------------------------------------------------
Jiri Sponer
Institute of Biophysics
Academy of Sciences of the Czech Republic
Kralovopolska 135
CZ-61265 Brno
Czech Republic
e-mail: sponer.ncbr.chemi.muni.cz
fax: 420 5412 12179
phone: 420 5415 17133
http://www.ibp.cz/labs/LSDNA/
http://ncbr.chemi.muni.cz
Senior Wellcome Trust International Research Fellow for Biomedical Science
--------------------------------------------------------
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sat May 22 2004 - 15:53:01 PDT
