Computational studies on imidazole heme conformations

Density functional theory computations of heme with ionized propionic acid groups, axially coordinated with two imidazoles, were performed for different mutual orientations of the imidazole planes. Environmental influences from water or protein were considered with a continuum dielectric medium by solving the Poisson equation. In vacuum, optimized geometries yielded imidazole-heme conformations where the NH groups of imidazoles are oriented toward the heme propionic groups in agreement with data from crystal structures of heme proteins. Conformational free-energy dependencies of the mutual orientation of axially ligated imidazoles calculated in protein (epsilon=10) and water (epsilon=80) environments confirmed the vacuum results, albeit the energy difference between the preferred and the 180 degrees opposite orientations of the imidazole ligand decreased from 3.84 kcal/mol in vacuum to 2.35 and 2.40 kcal/mol in protein and water, respectively. Two main factors determine the imidazole orientation: (1) the direct intramolecular electrostatic interactions of propionic groups with the polar NH groups of imidazole and (2) the electrostatic interaction of the total dipole moment of the imidazole-heme complex with the reaction field. In vacuum, only the first type of interaction is present, while in a dielectric medium the latter effect becomes competitive at high dielectric constant, resulting in a decrease of the orientational preference. Interestingly, the orientational preference of the imidazole axially ligated to heme becomes even more pronounced, if the negatively charged propionates are neutralized by counter charges that mimic salt bridges or protonation of the propionates.