Abstract
The protein contribution to the relative binding affinity of the ligands CO and O2 toward myoglobin (Mb) has been simulated using free energy perturbation calculations. The tautomers of the His E7 residue are different for the oxymyoglobin (MbO2) and carboxymyoglobin (MbCO) systems. This was modeled by performing two-step calculations that mutate the ligand and mutate the His E7 tautomers in separate steps. Differences in hydrogen bonding to the O2 and CO ligands were incorporated into the model. The O2 complex was calculated to be 2-3 kcal/mol more stable than the corresponding CO complex when compared to the same difference in an isolated heme control. This value agrees well with the experimental value of 2.0 kcal/mol. In qualitative agreement with experiments, the Fe-C-O bond is found to be bent (theta = 159.8 degrees) with a small tilt (theta = 6.2 degrees). The contributions made by each of the 29 residues--within the 9.0-A radius of the iron atom--to the free energy difference are separated into van der Waals and electrostatic contributions; the latter contributions are dominant. Aside from the proximal histidine and the heme group, the residues having the largest difference in free energy in mutating MbO2-->MbCO are His E7, Phe CD1, Phe CD4, Val E11, and Thr E10.