Search

Cytochrome c oxidase is a membrane protein of the respiratory chain that consumes protons and molecular oxygen to produce water and uses the resulting energy to pump protons across the membrane. Our molecular dynamics simulations with an excess proton located at different positions in one of the proton-conducting channels, the K-channel, show a clear dependence of the number of water molecules inside the channel on the proton position. A higher hydration level facilitates the formation of hydrogen-bonded chains along which proton transfer can occur. However, a sufficiently high hydration level for such proton transport is observed only when the excess proton is located above S365, i.e., the lower third of the channel. From the channel entrance up to this point, proton transport is via water molecules as proton carriers. These hydronium ions move with their surrounding water molecules, up to K362, filling and widening the channel. The conformation of K362 depends on its own protonation state and on the hydration level, suggesting its role to be proton transport from a hydronium ion at the height of K362 to the upper part of the channel via a conformational change. The protonation-dependent conformational dynamics of E101 at the bottom of the channel renders proton transfer via E101 unlikely. Instead, its role is rather that of an amplifier of H96’s proton affinity, suggesting H96 as the initial proton acceptor.

Proton transfer reactions are one of the most fundamental processes in biochemistry. We present a simplistic approach for estimating proton transfer probabilities in a membrane protein, cytochrome c oxidase. We combine short molecular dynamics simulations at discrete protonation states with a Monte Carlo approach to exchange between those states. Requesting for a proton transfer the existence of a hydrogen-bonded connection between the two source and target residues of the exchange, restricts the acceptance of transfers to only those in which a proton-relay is possible. Together with an analysis of the hydrogen-bonded connectivity in one of the proton-conducting channels of cytochrome c oxidase, this approach gives insight into the protonation dynamics of the hydrogen-bonded networks. The connectivity and directionality of the networks are coupled to the conformation of an important protein residue in the channel, K362, rendering proton transfer in the entire channel feasible in only one of the two major conformations. Proton transport in the channel can thus be regulated by K362 not only through its possible role as a proton carrier itself, but also by allowing or preventing proton transport via water residues.