Revisiting the Putative Ferryl Tilting Mode of Oxidized Cytochrome c Oxidase with Density Functional Vibrational Analyses of Model Complexes

Predicting Nuclear Resonance Vibrational Spectra of [Fe(OEP)(NO)]

Nuclear Resonance Vibrational Spectroscopy (NRVS) is a sensitive vibrational probe for biologically important heme complexes. The exquisite sensitivity of the NRVS data to the electronic structure provides detailed insights into the nature of these interesting compounds, but requires highly accurate computational methods for the mode assignments. To determine the best combinations of density functionals and basis sets, a series of benchmark DFT calculations on the previously characterized complex [Fe(OEP)NO] (OEP(2-)=octaethylporphyrinatio dianion) were performed. A test set of 21 methodology combinations including 8 functionals (BP86, mPWPW91, B3LYP, PBE1PBE, M062X, M06L, LC-BP86 and ωB97X-D) and 5 basis set (VTZ, TZVP, Lanl2DZ for iron and 6-31G*, 6-31+G* for other atoms) was carried out to calculate electronic structures and vibrational frequencies. We also implemented the conversion of frequency calculations into orientation-selective mode composition factors (e(2)), which can used to simulate the Vibrational Density Of States (VDOS) using Gaussian normal distribution functions. These use a series of user-friendly scripts for their application to NRVS. The structures as well as the isotropic and anisotropic NRVS of [Fe(OEP)NO] obtained with the M06L functional with a variety of basis sets are found to best reproduce the available experimental data, followed by B3LYP/LanL2DZ calculations. Other density functionals and basis sets do not produce the same level of accuracy. The noticeably worse agreement between theory and experiment for the out-plane NRVS compared with the excellent performance of the M06L functional for the in-plane prediction is attributed to deficiencies of the physical model rather than the computational methodology.

Vibrational resonances and CuB displacement controlled by proton motion in cytochrome c oxidase

Cytochrome c oxidase (CcO), found in the inner mitochondrial membranes or in many bacteria, catalyzes the four-electron reduction of molecular oxygen to water. Four protons are pumped across the inner mitochondrial membrane through CcO. In this study, quantum mechanics/molecular mechanics and molecular dynamics calculations are used to probe the spectroscopic characteristics of the ferryl intermediates in the aa(3) CcO/O(2) reaction. These highly elaborate calculations, supported by several calculations on smaller model systems, demonstrate the sensitivity of vibrational frequencies on the Coulombic field of heme a(3) and their dependence on the distance of the adjacent Cu(B) to the heme a(3)-Fe atom. This distance seems to be associated with the protonation state of the heme a(3) propionate A, and we propose that it plays a crucial role on the mechanism of action of CcO. In detail, we link proton pumping activity in CcO enzyme (a) to a multiple (1:1:2) resonance among the frequencies of Fe(IV)=O bond stretching, the breathing mode of Histidine 411, and a bending mode of the His411-Fe(IV)=O species (aa(3) from Paracoccus denitrificans numbering) and (b) to Cu(B) displacement by electrostatic interactions toward the heme a(3) iron. We find that the vibrations of the His411-Fe(IV)=O unit become highly coupled depending on the protonation state of the heme a(3) ring A propionate/Asp399 pair, and we propose a mechanism for the resonance Raman enhancement of the bending mode delta(His411-Fe(IV)=O). Calculations on model systems demonstrate that the position of Cu(B) in relation to heme a(3) iron-oxo plays a crucial role in regulating that resonance. We also discuss the origin of the coupling between bending, delta(His411-Fe(IV)=O) and nu(Fe=O) stretching modes, and the role played by such vibrational coupling interactions or Cu(B) position in controlling functional properties of the enzyme, including electron/proton coupling as well as experimental spectra.

The intrinsic stability of the second intermediate following the dioxygen-bound form in the O2 reduction by cytochrome c oxidase

Aeration of a two-electron reduced cytochrome c oxidase provides a species with two Raman bands at 804 and 356 cm(-1), identifying it as the second intermediate following the O2-bound species in the enzymatic O2 reduction process. It degrades directly to the fully oxidized form with a half-life time of 70 min at pH 8.0. The stability suggests an effective insulation for the active site in an extremely high oxidation state (Fe4+ with one oxidative equivalent nearby) against spontaneous electron leaks, which would dissipate proton motive force. The formation and degradation of the second intermediate are pH-dependent.

Resonance Raman characterization of the P intermediate in the reaction of bovine cytochrome c oxidase

Reduced cytochrome c oxidase binds molecular oxygen, yielding an oxygenated intermediate first (Oxy) and then converts it to water via the reaction intermediates of P, F, and O in the order of appearance. We have determined the iron-oxygen stretching frequencies for all the intermediates by using time-resolved resonance Raman spectroscopy. The bound dioxygen in Oxy does not form a bridged structure with Cu(B) and the rate of the reaction from Oxy to P (P(R)) is slower at higher pH in the pH range between 6.8 and 8.0. It was established that the P intermediate has an oxo-heme and definitely not the Fe(a(3))-O-O-Cu(B) peroxy bridged structure. The Fe(a(3))=O stretching (nu(Fe=O)) frequency of the P(R) intermediate, 804/764 cm(-1) for (16)O/(18)O, is distinctly higher than that of F intermediate, 785/750 cm(-1). The rate of reaction from P to F in D(2)O solution is evidently slower than that in H(2)O solution, implicating the coupling of the electron transfer with vector proton transfer in this process. The P intermediate (607-nm form) generated in the reaction of oxidized enzyme with H(2)O(2) gave the nu(Fe=O) band at 803/769 cm(-1) for H(2)(16)O(2)/H(2)(18)O(2) and the simultaneously measured absorption spectrum exhibited the difference peak at 607 nm. Reaction of the mixed valence CO adduct with O(2) provided the P intermediate (P(M)) giving rise to an absorption peak at 607 nm and the nu(Fe=O) bands at 804/768 cm(-1). Thus, three kinds of P intermediates are considered to have the same oxo-heme a(3) structure. The nu(4) and nu(2) modes of heme a(3) of the P intermediate were identified at 1377 and 1591 cm(-1), respectively. The Raman excitation profiles of the nu(Fe=O) bands were different between P and F. These observations may mean the formation of a pi cation radical of porphyrin macrocycle in P.

Time dependence of the catalytic intermediates in cytochrome c oxidase

Cytochrome c oxidase, the terminal enzyme in the electron transfer chain, catalyzes the reduction of oxygen to water in a multiple step process by utilizing four electrons from cytochrome c. To study the reaction mechanism, the resonance Raman spectra of the intermediate states were measured during single turnover of the enzyme after catalytic initiation by photolysis of CO from the fully reduced CO-bound enzyme. By measuring the change in intensity of lines associated with heme a, the electron transfer steps were determined and found to be biphasic with apparent rate constants of approximately 40 x 10(3) s(-1) and approximately 1 x 10(3) s(-1). The time dependence for the oxidation of heme a and for the measured formation and decay of the oxy, the ferryl ("F"), and the hydroxy intermediates could be simulated by a simple reaction scheme. In this scheme, the presence of the "peroxy" ("P") intermediate does not build up a sufficient population to be detected because its decay rate is too fast in buffered H(2)O at neutral pH. A comparison of the change in the spin equilibrium with the formation of the hydroxy intermediate demonstrates that this intermediate is high spin. We also confirm the presence of an oxygen isotope-sensitive line at 355 cm(-1), detectable in the spectrum from 130 to 980 micros, coincident with the presence of the F intermediate.