The distal residue-CO interaction in carbonmonoxy myoglobins: a molecular dynamics study of two distal histidine tautomers.

Four independent 90 ps molecular dynamics simulations of sperm-whale wild-type carbonmonoxy myoglobin (MbCO) have been calculated using a new AMBER force field for the haem prosthetic group. Two trajectories have the distal 64N delta nitrogen protonated, and two have the 64N epsilon nitrogen protonated; all water molecules within 16 A of the carbonyl O are included. In three trajectories, the distal residue remains part of the haem pocket, with the protonated distal nitrogen pointing into the active site. This is in contrast with the neutron diffraction crystal structure, but is consistent with the solution phase CO stretching frequencies (upsilon CO) of MbCO and various of its mutants. There are significant differences in the "closed" pocket structures found for each tautomer: the 64N epsilon H trajectories both show stable distal-CO interactions, whereas the 64N delta H tautomer) has a weaker interaction resulting in a more mobile distal side chain. One trajectory (a 64N delta H tautomer) has the distal histidine moving out into the "solvent", leaving the pocket in an "open" structure, with a large unhindered entrance to the active site. These trajectories suggest that the three upsilon CO frequencies observed for wild-type MbCO in solution, rather than representing significantly different Fe-C-O geometries as such, arise from three different haem pocket structures, each with different electric fields at the ligand. Each pocket structure corresponds to a different distal histidine conformer: the A3 band to the 64N epsilon H tautomer, the A1,2 band to the 64N delta H tautomer, and the A0 band to the absence of any significant interaction with the distal side chain.