Cytochrome c oxidase (CcO) reduces oxygen to water and uses the released free energy to pump protons across the membrane. We have used multiscale reactive molecular dynamics simulations to explicitly characterize (with free-energy profiles and calculated rates) the internal proton transport events that enable proton pumping during first steps of oxidation of the fully reduced enzyme. Our results show that proton transport from amino acid residue E286 to both the pump loading site (PLS) and to the binuclear center (BNC) are thermodynamically driven by electron transfer from heme a to the BNC, but that the former (i.e., pumping) is kinetically favored whereas the latter (i.e., transfer of the chemical proton) is rate-limiting. The calculated rates agree with experimental measurements. The backflow of the pumped proton from the PLS to E286 and from E286 to the inside of the membrane is prevented by large free-energy barriers for the backflow reactions. Proton transport from E286 to the PLS through the hydrophobic cavity and from D132 to E286 through the D-channel are found to be strongly coupled to dynamical hydration changes in the corresponding pathways and, importantly, vice versa.cytochrome c oxidase | proton-pumping mechanism | reactive molecular dynamics | bioenergetics | proton transport C ytochrome c oxidase (CcO) (Fig. 1) catalyzes the reduction of O 2 to H 2 O and couples the free energy of this exergonic reaction to the pumping of protons across the membrane, creating a transmembrane proton electrochemical gradient that drives, for example, ATP synthesis (1-3). During each reaction cycle eight protons are taken up from the negatively charged inside (N-side) of the membrane and either react with oxygen (referred to as "chemical" protons below) or are pumped to the positively charged outside (P-side) of the membrane (referred to as "pumped" protons below). In aa 3 -type CcO, the D-channel is responsible for uptake of all four pumped protons and at least one out of four chemical protons (4,5). Protons on the N-side are taken into D-channel via the amino acid residue D132 at the channel entrance, and then transferred to residue E286 in the middle of the membrane (4, 5). By Grotthuss shuttling through the water molecules in the hydrophobic cavity (HC) above E286, each proton is either transferred to react with oxygen in the binuclear center (BNC), consisting of heme a 3 and the Cu B complex, or transferred to the pump loading site (PLS) and then further released to the P-side of the membrane (Fig. 1).Despite decades of study, the CcO proton-pumping mechanism is still incompletely understood at the atomistic level. For example, controversy remains regarding how the internal proton transport (PT) events are coupled with electron transfer (ET) (3,(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16), what features enable pumping while preventing proton backflow from the P-side to the N-side (3,6,10,11,15,17,18), and how hydration influences PT through the . Recent computational work proposed a stepwise pumping mechanism in w...