Cytochrome bd is a distinctive family of terminal oxidases present in the respiratory chains of many prokaryotes. Despite their biological importance, the redox chemistry of these proteins remains poorly understood, largely due to the presence of two b-type hemes and one d-type heme. Here, we report the first computational study of interheme electron transfer in the cytochrome bd family. We performed 10 μs of molecular dynamics simulations of E. coli cytochrome bd-I embedded in realistic membranes, combined with quantum chemical calculations to estimate the thermodynamic parameters of electron transfer from heme b595 to heme d within the framework of Marcus theory. We further identify the respective contributions of the hemes, protein scaffold, lipid bilayer, water, and counterions to the driving force and reorganization energy. The interheme electronic coupling was calculated using the Projected Orbital Diabatization (POD) method in a hybrid Quantum Mechanics/Molecular Mechanics scheme and rationalized through electron transfer pathway analysis. This study provides fundamental insights into how electron transfer steps are orchestrated in the catalytic cycle of E. coli cytochrome bd-I.
Siddeeque et al. (Thu,) studied this question.