Insights into functions of the H channel of cytochrome c oxidase from atomistic molecular dynamics simulations [Biophysics and Computational Biology]

Proton pumping A-type cytochrome c oxidase (CcO) terminates the respiratory chains of mitochondria and many bacteria. Three possible proton transfer pathways (D, K, and H channels) have been identified based on structural, functional, and mutational data. Whereas the D channel provides the route for all pumped protons in bacterial A-type CcOs, studies of bovine mitochondrial CcO have led to suggestions that its H channel instead provides this route. Here, we have studied H-channel function by performing atomistic molecular dynamics simulations on the entire, as well as core, structure of bovine CcO in a lipid-solvent environment. The majority of residues in the H channel do not undergo large conformational fluctuations. Its upper and middle regions have adequate hydration and H-bonding residues to form potential proton-conducting channels, and Asp51 exhibits conformational fluctuations that have been observed crystallographically. In contrast, throughout the simulations, we do not observe transient water networks that could support proton transfer from the N phase toward heme a via neutral His413, regardless of a labile H bond between Ser382 and the hydroxyethylfarnesyl group of heme a. In fact, the region around His413 only became sufficiently hydrated when His413 was fixed in its protonated imidazolium state, but its calculated pK_a is too low for this to provide the means to create a proton transfer pathway. Our simulations show that the electric dipole moment of residues around heme a changes with the redox state, hence suggesting that the H channel could play a more general role as a dielectric well.