Monte Carlo simulations of proton pumps: On the working principles of the biological valve that controls proton pumping in cytochrome c oxidase

Gaining a detailed understanding of the proton-pumping process in cytochrome c oxidase (COX) is one of the challenges of modern biophysics. Recent mutation experiments have highlighted this challenge by showing that a single mutation (the N139D mutation) blocks the overall pumping while continuing to channel protons to the binuclear center without inhibiting the oxidase activity. Rationalizing this result has been a major problem because the mutation is quite far from E286, which is believed to serve as the branching point for the proton transport in the pumping process. In the absence of a reasonable explanation for this important observation, we have developed a Monte Carlo simulation method that can convert mutation and structural information to pathways for proton translocation and simulate the pumping process in COX on a millisecond and even subsecond time scale. This tool allows us to reproduce and propose a possible explanation to the effect of the N139D mutation and to offer a consistent model for the origin of the valve effect in COX, which is crucial for maintaining uphill proton pumping. Furthermore, obtaining the first structure-based simulation of proton pumping in COX, or in any other protein, indicates that our approach should provide a powerful tool for verification of mechanistic hypotheses about the action of proton transport proteins.