Effects of mutations on the molecular dynamics of oxygen escape from the dimeric hemoglobin of Scapharca inaequivalvis

Like many hemoglobins, the structure of the dimeric hemoglobin from the clam Scapharca inaequivalvis is a “closed bottle” since there is no direct tunnel from the oxygen binding site on the heme to the solvent. The proximal histidine faces the dimer interface, which consists of the E and F helicies. This is significantly different from tetrameric vertebrate hemoglobins and brings the heme groups near the subunit interface. The subunit interface is also characterized by an immobile, hydrogen-bonded network of water molecules. Although there is data which is consistent with the histidine gate pathway for ligand escape, these aspects of the structure would seem to make that pathway less likely. Locally enhanced sampling molecular dynamics are used here to suggest alternative pathways in the wild-type and six mutant proteins. In most cases the point mutations change the selection of exit routes observed in the simulations. Exit via the histidine gate is rarely seem although oxygen molecules do occasionally cross over the interface from one subunit to the other. The results suggest that changes in flexibility and, in some cases, creation of new cavities can explain the effects of the mutations on ligand exit paths.