Calculation of atom-centered partial charges for heme

Atom-centered partial charges are calculated for the Fe-heme in cytochrome P450cam for use in molecular dynamics simulations of polar substrates bound in the active site of the enzyme. Charges are fit to the electrostatic potential produced by ab initio UHF wavefunctions for an Fe-porphine model. Basis set dependence of these charges is observed using the LANL1DZ, LANL2DZ and augmented 6–31G levels of theory. Upon geometry optimization of the enzyme, these charge sets cause varying degrees of distortion of the porphyrin from its crystallographically observed conformation. Scaling the charges calculated from the augmented 6–31G basis by 75% reduces the heme distortion while preserving reasonable interactions with a polar substrate. A comparison of the calculated charges with other published values is presented.