Investigating the Many Roles of Internal Water in Cytochrome c Oxidase

Searching for proton transfer channels in respiratory complex I

We have explored a strategy to identify potential proton transfer channels using computational analysis of a protein structure based on Voronoi partitioning and applied it for the analysis of proton transfer pathways in redox-driven proton-pumping respiratory complex I. The analysis results in a network of connected voids/channels, which represent the dual structure of the protein; we then hydrated the identified channels using our water placement program Dowser++. Many theoretical water molecules found in the channels perfectly match the observed experimental water molecules in the structure; some other predicted water molecules have not been resolved in the experiments. The channels are of varying cross sections. Some channels are big enough to accommodate water molecules that are suitable to conduct protons; others are too narrow to hold water but require only minor conformational changes to accommodate proton transfer. We provide a preliminary analysis of the proton conductivity of the network channels, classifying the proton transfer channels as open, closed, and partially open, and discuss possible conformational changes that can modulate, i.e., open and close, the channels.

Interplay of Hydration and Protonation Dynamics in the K-Channel of Cytochrome c Oxidase

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.

Comparison of proton transfer paths to the QA and QB sites of the Rb. sphaeroides photosynthetic reaction centers

The photosynthetic bacterial reaction centers from purple non-sulfur bacteria use light energy to drive the transfer of electrons from cytochrome c to ubiquinone. Ubiquinone bound in the Q_A site cycles between quinone, Q_A, and anionic semiquinone, Q_A^·-, being reduced once and never binding protons. In the Q_B site, ubiquinone is reduced twice by Q_A^·-, binds two protons and is released into the membrane as the quinol, QH₂. The network of hydrogen bonds formed in a molecular dynamics trajectory was drawn to investigate proton transfer pathways from the cytoplasm to each quinone binding site. Q_A is isolated with no path for protons to enter from the surface. In contrast, there is a complex and tangled network requiring residues and waters that can bring protons to Q_B. There are three entries from clusters of surface residues centered around HisH126, GluH224, and HisH68. The network is in good agreement with earlier studies, Mutation of key nodes in the network, such as SerL223, were previously shown to slow proton delivery. Mutational studies had also shown that double mutations of residues such as AspM17 and AspL210 along multiple paths in the network presented here slow the reaction, while single mutations do not. Likewise, mutation of both HisH126 and HisH128, which are at the entry to two paths reduce the rate of proton uptake.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima

Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λ_max) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λ_max of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λ_max from the experimental λ_max. Although the major effect of water molecules on λ_max originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λ_max mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.

Conformational Motions and Water Networks at the α/β Interface in E. coli Ribonucleotide Reductase

Ribonucleotide reductases (RNRs) catalyze the conversion of all four ribonucleotides to deoxyribonucleotides and are essential for DNA synthesis in all organisms. The active form of E. coli Ia RNR is composed of two homodimers that form the active α₂β₂ complex. Catalysis is initiated by long-range radical translocation over a ∼32 Å proton-coupled electron transfer (PCET) pathway involving Y356β and Y731α at the interface. Resolving the PCET pathway at the α/β interface has been a long-standing challenge due to the lack of structural data. Herein, molecular dynamics simulations based on a recently solved cryogenic-electron microscopy structure of an active α₂β₂ complex are performed to examine the structure and fluctuations of interfacial water, as well as the hydrogen-bonding interactions and conformational motions of interfacial residues along the PCET pathway. Our free energy simulations reveal that Y731 is able to sample both a flipped-out conformation, where it points toward the interface to facilitate interfacial PCET with Y356, and a stacked conformation with Y730 to enable collinear PCET with this residue. Y356 and Y731 exhibit hydrogen-bonding interactions with interfacial water molecules and, in some conformations, share a bridging water molecule, suggesting that the primary proton acceptor for PCET from Y356 and from Y731 is interfacial water. The conformational flexibility of Y731 and the hydrogen-bonding interactions of both Y731 and Y356 with interfacial water and hydrogen-bonded water chains appear critical for effective radical translocation along the PCET pathway. These simulations are consistent with biochemical and spectroscopic data and provide previously unattainable atomic-level insights into the fundamental mechanism of RNR.

The structure of the oxidized state of cytochrome c oxidase - experiments and theory compared

Cytochrome c oxidase (CcO), the terminal enzyme in the respiratory chain, reduces molecular oxygen to water. Experimental data on the midpoint potentials of the heme iron/copper active site cofactors do not match the overall reaction energetics, and are also in conflict with the observed efficiency of energy conservation in CcO. Therefore it has been postulated that the ferric/cupric intermediate (the oxidized state) exists in two forms. One form, labelled O_H, is presumably involved during catalytic turnover, and should have a high Cu_B midpoint potential due to a metastable high energy structure. When no more electrons are supplied, the O_H state supposedly relaxes to the resting form, labelled O, with a lower energy and a lower midpoint potential. It has been suggested that there is a pure geometrical difference between the O_H and O states, obtained by moving a water molecule inside the active site. It is shown here that the difference between the two forms of the oxidized state must be of a more chemical nature. The reason is that all types of geometrically relaxed structures of the oxidized intermediate have similar energies, all with a high proton coupled reduction potential in accordance with the postulated O_H state. One hypothesized chemical modification of the O_H state is the transfer of an extra proton, possibly internal, into the active site. Such a protonated state has several properties that agree with experimental data on the relaxed oxidized state, including a decreased midpoint potential.

Allosteric Cooperativity in Proton Energy Conversion in A1-Type Cytochrome c Oxidase

Cytochrome c oxidase (CcO), the Cu_A, heme a, heme a₃, Cu_B enzyme of respiratory chain, converts the free energy released by aerobic cytochrome c oxidation into a membrane electrochemical proton gradient (ΔμH⁺). ΔμH⁺ derives from the membrane anisotropic arrangement of dioxygen reduction to two water molecules and transmembrane proton pumping from a negative (N) space to a positive (P) space separated by the membrane. Spectroscopic, potentiometric, and X-ray crystallographic analyses characterize allosteric cooperativity of dioxygen binding and reduction with protonmotive conformational states of CcO. These studies show that allosteric cooperativity stabilizes the favorable conformational state for conversion of redox energy into a transmembrane ΔμH⁺.