Definition of the interaction domain for cytochrome c on cytochrome c oxidase. III. Prediction of the docked complex by a complete, systematic search

The electron transfer complex between bovine cytochrome c oxidase and horse cytochrome c has been predicted with the docking program DOT, which performs a complete, systematic search over all six rotational and translational degrees of freedom. Energies for over 36 billion configurations were calculated, providing a free-energy landscape showing guidance of positively charged cytochrome c to the negative region on the cytochrome c oxidase surface formed by subunit II. In a representative configuration, the solvent-exposed cytochrome c heme edge is within 4 A of the indole ring of subunit II residue Trp(104), indicating a likely electron transfer path. These two groups are surrounded by a small, hydrophobic contact region, which is surrounded by electrostatically complementary hydrophilic interactions. Cytochrome c/cytochrome c oxidase interactions of Lys(13) with Asp(119) and Lys(72) with Gln(103) and Asp(158) are the most critical polar interactions due to their proximity to the hydrophobic region and exclusion from bulk solvent. The predicted complex matches previous mutagenesis, binding, and time-resolved kinetics studies that implicate Trp(104) in electron transfer and show the importance of specific charged residues to protein affinity. Electrostatic forces not only enhance long range protein/protein association; they also predominate in short range alignment, creating the transient interaction needed for rapid turnover.