DFT Fea3-O/O-O Vibrational Frequency Calculations over Catalytic Reaction Cycle States in the Dinuclear Center of Cytochrome c Oxidase.

Density functional vibrational frequency calculations have been performed on eight geometry optimized cytochrome c oxidase (CcO) dinuclear center (DNC) reaction cycle intermediates and on the oxymyoglobin (oxyMb) active site. The calculated Fe-O and O-O stretching modes and their frequency shifts along the reaction cycle have been compared with the available resonance Raman (rR) measurements. The calculations support the proposal that in state A[Fea33+-O2-•···CuB+] of CcO, O2 binds with Fea32+ in a similar bent end-on geometry to that in oxyMb. The calculations show that the observed 20 cm-1 shift of the Fea3-O stretching mode from the PR to F state is caused by the protonation of the OH- ligand on CuB2+ (PR[Fea34+═O2-···HO--CuB2+] → F[Fea34+═O2-···H2O-CuB2+]), and that the H2O ligand is still on the CuB2+ site in the rR identified F[Fea34+═O2-···H2O-CuB2+] state. Further, the observed rR band at 356 cm-1 between states PR and F is likely an O-Fea3-porphyrin bending mode. The observed 450 cm-1 low Fea3-O frequency mode for the OH active oxidized state has been reproduced by our calculations on a nearly symmetrically bridged Fea33+-OH-CuB2+ structure with a relatively long Fea3-O distance near 2 Å. Based on Badger's rule, the calculated Fea3-O distances correlate well with the calculated νFe-O-2/3 (νFe-O is the Fea3-O stretching frequency) with correlation coefficient R = 0.973.