Computational Study of Heme <i>b</i> ₅₉₅ to Heme <i>d</i> Electron Transfer in <i>E. coli</i> Cytochrome <i>bd</i> -I Oxidase

Cytochrome <i>bd</i> 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 <i>b</i>-type hemes and one <i>d</i>-type heme. Here, we report the first computational study of interheme electron transfer in the cytochrome <i>bd</i> family. We performed 10 μs of molecular dynamics simulations of <i>E. coli</i> cytochrome <i>bd</i>-I embedded in realistic membranes, combined with quantum chemical calculations to estimate the thermodynamic parameters of electron transfer from heme <i>b</i>₅₉₅ to heme <i>d</i> 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 <i>E. coli</i> cytochrome <i>bd</i>-I.