A modeling study of the interaction and electron transfer between cytochrome b5 and some oxidized haemoglobins

Using Brownian motion simulations we have studied the formation of docked complexes of reduced cytochrome b5 and oxidized haemoglobin. Our results indicate that the presence of molecular electrostatic fields has a significant role to play in the formation of these complexes. In contrast to previous modeling studies on this system, we clearly identify electron transfer within an ensemble of similarly docked complexes rather than the formation of a single complex. Docking involves a number of acidic residues surrounding the exposed haem edge of cytochrome b5 and a set of basic residues surrounding the exposed haem edge of the globins. Although amino acids from the partner globin proteins are involved to a small extent in the binding of some of the complexes, the reactivity of any particular globin is essentially independent of the nature of its partner globin chain within the haemoglobin molecule. Comparison of results from adult and embryonic haemoglobins indicates a significant difference in complex formation. Application of electron tunneling analysis to the complexes allows us to predict the rates of electron transfer within each ensemble of complexes. These data provide a theoretical insight into the important process of re-reduction of oxidized haemoglobins as well as explaining the general inability to produce crystalline forms of many docked electron transfer complexes.