Hello,

I want to derive the atomic partial charges for heme in chloroperoxidase,
since these parameters cannot be found in the GROMOS53a6 force field.
Chloroperoxidase is a cysteine-ligated heme protein (high spin) much like
P450. In my work, I want to study the interactions between the protein
residues of chloroperoxidase and its substrates. Thus, I used the QM
software ORCA to optimize the heme-thiolate structure (TPSS/def2-SPVD), and
then I used CHELPG method to derive the heme charges. Could anyone tell me
if the method I used is suitable for my purpose?

In addition, the charges I obtained for Fe (0.931877) and nitrogen
(-0.44572, -0.34427, -0.39467, and -0.39467) are largely consistent with
two reference. The reference are one in 1995 (Fe: 1.0, N: -0.4; Helms et
al. Thermodynamics of water mediating protein-ligand interactions in
cytochrome P450cam: a molecular dynamics study, Biophys. J. 69: 810-24) and
one in recent (Fe: 0.847, N: -0.37, -0.423, -0.504, and -0.528; Favia et
al. (2006) Three-dimensional model of the human aromatase enzyme and
density functional parameterization of the iron-containing protoporphyrin
IX for a molecular dynamics study of heme-cysteinato cytochromes, Proteins
62: 1073-84). I will greatly appreciate any comment and suggestions.

Best Regards,

Rui
-- 
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to