Monte Carlo Simulations of Glu-242 in Cytochrome c Oxidase

Electric fields control water-gated proton transfer in cytochrome c oxidase

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O₂ reduction into proton pumping. Here we show that CcO creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a, leads to dissociation of an arginine (Arg438)-heme a₃ D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å^-1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.

Cytochrome c oxidase structures suggest a four-state stochastic pump mechanism

A simple stochastic model for a cytochrome c oxidase proton pump.

Investigating the Many Roles of Internal Water in Cytochrome c Oxidase

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain. As part of its catalytic cycle, C cO transfers protons to its Fe-Cu binuclear center (BNC) to reduce oxygen, and in addition, it pumps protons across the mitochondrial inner, or bacterial, membrane where it is located. It is believed that this proton transport is facilitated by a network of water chains inside the enzyme. Here we present an analysis of the hydration of C cO, including the BNC region, using a semi-empirical hydration program, Dowser++, recently developed in our group. Using high-resolution X-ray data, we show that Dowser++ predictions match very accurately the water molecules seen in the D- and K-channels of C cO, as well as in the vicinity of its BNC. Moreover, Dowser++ predicts many more internal water molecules than is typically seen in the experiment. However, no significant hydration of the catalytic cavity in C cO described recently in the literature is observed. As Dowser++ itself does not account for structural changes of the protein, this result supports the earlier assessment that the proposed wetting transition in the catalytic cavity can only either be due to structural rearrangements of BNC, possibly induced by the charges during the catalytic cycle, or occur transiently, in concert with the proton transfer. Molecular dynamics simulations were performed to investigate the global dynamic nature of Dowser++ waters in C cO, and the results suggest a consistent explanation as to why some predicted water molecules would be missing in the experimental structures. Furthermore, in light of the significant protein hydration predicted by Dowser++, the dielectric constant of the hydrated cavities in C cO was also investigated using the Fröhlich-Kirkwood model; the results indicate that in the cavities where water is packed sufficiently densely the dielectric constant can approach values comparable even to that of bulk water.

Dewetting transitions coupled to K-channel activation in cytochrome c oxidase

Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane.

Cytochrome c oxidase (CcO) drives aerobic respiratory chains in all organisms by transducing the free energy from oxygen reduction into an electrochemical proton gradient across a biological membrane. CcO employs the so-called D- and K-channels for proton uptake, but the molecular mechanism for activation of the K-channel has remained elusive for decades. We show here by combining large-scale atomistic molecular simulations with graph-theoretical water network analysis, and hybrid quantum/classical (QM/MM) free energy calculations, that the K-channel is activated by formation of a reactive oxidized intermediate in the binuclear heme a₃/Cu_B active site. This state induces electrostatic, hydration, and conformational changes that lower the barrier for proton transfer along the K-channel by dewetting pathways that connect the D-channel with the active site. Our combined results reconcile previous experimental findings and indicate that water dynamics plays a decisive role in the proton pumping machinery in CcO.

Cavity hydration dynamics in cytochrome c oxidase and functional implications

Cytochrome c oxidase (CcO) is a transmembrane protein that uses the free energy of O₂ reduction to generate the proton concentration gradient across the membrane. The regulation of competitive proton transfer pathways has been established to be essential to the vectorial transport efficiency of CcO, yet the underlying mechanism at the molecular level remains lacking. Recent studies have highlighted the potential importance of hydration-level change in an internal cavity that connects the proton entrance channel, the site of O₂ reduction, and the putative proton exit route. In this work, we use atomistic molecular dynamics simulations to investigate the energetics and timescales associated with the volume fluctuation and hydration-level change in this central cavity. Extensive unrestrained molecular dynamics simulations (accumulatively [Formula: see text]4 [Formula: see text]s) and free energy computations for different chemical states of CcO support a model in which the volume and hydration level of the cavity are regulated by the protonation state of a propionate group of heme a₃ and, to a lesser degree, the redox state of heme a and protonation state of Glu286. Markov-state model analysis of [Formula: see text]2-[Formula: see text]s trajectories suggests that hydration-level change occurs on the timescale of 100-200 ns before the proton-loading site is protonated. The computed energetic and kinetic features for the cavity wetting transition suggest that reversible hydration-level change of the cavity can indeed be a key factor that regulates the branching of proton transfer events and therefore contributes to the vectorial efficiency of proton transport.

Building a More Predictive Protein Force Field: A Systematic and Reproducible Route to AMBER-FB15

The increasing availability of high-quality experimental data and first-principles calculations creates opportunities for developing more accurate empirical force fields for simulation of proteins. We developed the AMBER-FB15 protein force field by building a high-quality quantum chemical data set consisting of comprehensive potential energy scans and employing the ForceBalance software package for parameter optimization. The optimized potential surface allows for more significant thermodynamic fluctuations away from local minima. In validation studies where simulation results are compared to experimental measurements, AMBER-FB15 in combination with the updated TIP3P-FB water model predicts equilibrium properties with equivalent accuracy, and temperature dependent properties with significantly improved accuracy, in comparison with published models. We also discuss the effect of changing the protein force field and water model on the simulation results.

Insights into proton translocation in cbb3 oxidase from MD simulations

Heme-copper oxidases are membrane protein complexes that catalyse the final step of the aerobic respiration, namely the reduction of oxygen to water. The energy released during catalysis is coupled to the active translocation of protons across the membrane, which contributes to the establishment of an electrochemical gradient that is used for ATP synthesis. The distinctive C-type (or cbb₃) cytochrome c oxidases, which are mostly present in proteobacteria, exhibit a number of unique structural and functional features, including high catalytic activity at low oxygen concentrations. At the moment, the functioning mechanism of C-type oxidases, in particular the proton transfer/pumping mechanism presumably via a single proton channel, is still poorly understood. In this work we used all-atom molecular dynamics simulations and continuum electrostatics calculations to obtain atomic-level insights into the hydration and dynamics of a cbb₃ oxidase. We provide the details of the water dynamics and proton transfer pathways for both the "chemical" and "pumped" protons, and show that formation of protonic connections is strongly affected by the protonation state of key residues, namely H243, E323 and H337.