Storage of an excess proton in the hydrogen-bonded network of the d-pathway of cytochrome C oxidase: identification of a protonated water cluster

Water, Protons, and the Gating of Voltage-Gated Potassium Channels

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

Proton transfer in the D-channel of cytochrome c oxidase modeled by a transition network approach

BACKGROUND

Determination of proton uptake pathways in Cytochrome c Oxidase is difficult due to the complexity of the system. The transition networks approach allows sampling of proton transfer pathways without predefined reaction coordinate.

METHODS

Computation of the proton transfer pathways in a model of the D-channel of cytochrome c oxidase has been performed by a transition network approach that combines discrete, optimisation based and molecular dynamics based sampling.

RESULTS

The optimal pathway involves an opening of the so-called asparagine gate, hydration of the asparagine region, the formation of a hydrogen-bonded chain, and finally concerted proton hole transport along this chain. The optimal pathway finds the protonation of residue H26 close to the channel entrance favourable for lowering the transition energies of subsequent steps, in particular, opening of the Asn gate and formation of a hydrogen-bonded chain. Residue Y33 plays an important role in shuttling the transferred proton hole.

CONCLUSIONS

The optimal pathway found by the transition network approach shows the same important characteristics as pathways determined earlier by other methods. The computed barrier and reaction energies are also in good agreement with previous studies. The transition network approach provides an alternative to explore pathways in complex systems.

GENERAL SIGNIFICANCE

The correct function of the enzyme as oxidase and proton pump depends on the interplay of several redox and proton transport steps. Understanding the proton transport mechanism is therefore key to understanding the protein's function. The complex nature of long- distances proton transfer through a protein requires a non-trivial simulation strategy.

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⁺.

Spin Interconversion of Heme-Peroxo-Copper Complexes Facilitated by Intramolecular Hydrogen-Bonding Interactions

Synthetic peroxo-bridged high-spin (HS) heme-(μ-η2:η1-O22-)-Cu(L) complexes incorporating (as part of the copper ligand) intramolecular hydrogen-bond (H-bond) capabilities and/or steric effects are herein demonstrated to affect the complex's electronic and geometric structure, notably impacting the spin state. An H-bonding interaction with the peroxo core favors a low-spin (LS) heme-(μ-η1:η1-O22-)-Cu(L) structure, resulting in a reversible temperature-dependent interconversion of spin state (5 coordinate HS to 6 coordinate LS). The LS state dominates at low temperatures, even in the absence of a strong trans-axial heme ligand. Lewis base addition inhibits the H-bond facilitated spin interconversion by competition for the H-bond donor, illustrating the precise H-bonding interaction required to induce spin-crossover (SCO). Resonance Raman spectroscopy (rR) shows that the H-bonding pendant interacts with the bridging peroxide ligand to stabilize the LS but not the HS state. The H-bond (to the Cu-bound O atom) acts to weaken the O-O bond and strengthen the Fe-O bond, exhibiting ν(M-O) and ν(O-O) values comparable to analogous known LS complexes with a strong donating trans-axial ligand, 1,5-dicyclohexylimidazole, (DCHIm)heme-(μ-η1:η1-O22-)-Cu(L). Variable-temperature (-90 to -130 °C) UV-vis and 2H NMR spectroscopies confirm the SCO process and implicate the involvement of solvent binding. Examining a case of solvent binding without SCO, thermodynamic parameters were obtained from a van't Hoff analysis, accounting for its contribution in SCO. Taken together, these data provide evidence for the H-bond group facilitating a core geometry change and allowing solvent to bind, stabilizing a LS state. The rR data, complemented by DFT analysis, reveal a stronger H-bonding interaction with the peroxo core in the LS compared to the HS complexes, which enthalpically favors the LS state. These insights enhance our fundamental understanding of secondary coordination sphere influences in metalloenzymes.

Hydrogen-Bonded Network and Water Dynamics in the D-channel of Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase (CcO) from the cellular inside to the binuclear redox centre as well as proton pumping through the membrane takes place through proton entrance via two distinct pathways, the D- and K-channel. Both channels show a dependence of their hydration level on the protonation states of their key residues, K362 for the K-channel, and E286 or D132 for the D-channel. In the oxidative half of CcO's catalytic cycle the D-channel is the proton-conducting path. For this channel, an interplay of protonation state of the D-channel residues with the water and hydrogen-bond dynamics has been observed in molecular dynamics simulations of the CcO protein, embedded in a lipid bi-layer, modelled in different protonation states. Protonation of residue E286 at the end of the D-channel results in a hydrogen-bonded network pointing from E286 to N139, that is against proton transport, and favouring N139 conformations which correspond to a closed asparagine gate (formed by residues N121 and N139). Consequently, the hydration level is lower than with unprotonated E286. In those models, the Asn gate is predominantly open, allowing water molecules to pass and thus increase the hydration level. The hydrogen-bonded network in these states exhibits longer life times of the Asn residues with water than other models and shows the D-channel to be traversable from the entrance, D132, to exit, E286. The D-channel can thus be regarded as auto-regulated with respect to proton transport, allowing proton passage only when required, that is the proton is located at the lower part of the D-channel (D132 to Asn gate) and not at the exit (E286).

Orientation of non-spherical protonated water clusters revealed by infrared absorption dichroism

Infrared continuum bands that extend over a broad frequency range are a key spectral signature of protonated water clusters. They are observed for many membrane proteins that contain internal water molecules, but their microscopic mechanism has remained unclear. Here we compute infrared spectra for protonated and unprotonated water chains, discs, and droplets from ab initio molecular dynamics simulations. The continuum bands of the protonated clusters exhibit significant anisotropy for chains and discs, with increased absorption along the direction of maximal cluster extension. We show that the continuum band arises from the nuclei motion near the excess charge, with a long-ranged amplification due to the electronic polarizability. Our experimental, polarization-resolved light–dark difference spectrum of the light-driven proton pump bacteriorhodopsin exhibits a pronounced dichroic continuum band. Our results suggest that the protonated water cluster responsible for the continuum band of bacteriorhodopsin is oriented perpendicularly to the membrane normal.

Protein-bound water clusters play a key role for proton transport and storage in molecular biology. Here, the authors show by simulations and experiments that the orientation of non-spherical protonated water clusters in bacteriorhodopsin is unveiled by polarization-resolved infrared spectroscopy.

Protonation-State-Dependent Communication in Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase from the cellular inside to the binuclear redox center (BNC) can occur through two distinct pathways, the D- and K-channels. For the protein to function as both a redox enzyme and a proton pump, proton transfer into the protein toward the BNC or toward a proton loading site (and ultimately through the membrane) must be highly regulated. The P_R → F transition is the first step in a catalytic cycle that requires proton transfer from the bulk at the N-side to the BNC. Molecular dynamics simulations of the P_R → F intermediate of this transition, with 16 different combinations of protonation states of key residues in the D- and K-channel, show the impact of the K-channel on the D-channel to be protonation-state dependent. Strength as well as means of communication, correlations in positions, or communication along the hydrogen-bonded network depends on the protonation state of the K-channel residue K362. The conformational and hydrogen-bond dynamics of the D-channel residue N139 is regulated by an interplay of protonation in the D-channel and K362. N139 thus assumes a gating function by which proton passage through the D-channel toward E286 is likely facilitated for states with protonated K362 and unprotonated E286. In contrast, proton passage through the D-channel is hindered by N139's preference for a closed conformation in situations with protonated E286.

Distal Proton Shuttle Mechanism of Ribosome Catalysed Peptide Bond Formation-A Theoretical Study

In this work, we have investigated a novel distal proton shuttle mechanism of ribosome catalyzed peptide bond formation reaction. The reaction was found to follow a two-step mechanism. A distal water molecule located about 6.1 Å away from the attacking amine plays as a proton acceptor and results in a charge-separated intermediate that is stabilized by the N terminus of L27 and the A-site A76 5′-phosphate. The ribose A2451 bridges the proton shuttle pathway, thus plays critical role in the reaction. The calculated 27.64 kcal·mol⁻¹ free energy barrier of the distal proton shuttle mechanism is lower than that of eight-membered ring transition state. The distal proton shuttle mechanism studied in this work can provide new insights into the important biological peptide synthesis process.

Monte Carlo Simulations of Glu-242 in Cytochrome c Oxidase

Monte Carlo (MC) simulations of conformational changes and protonation of Glu-242, a key residue that shuttles protons in cytochrome c oxidase (CcO), are reported. Previous studies suggest that this residue may play a role of the valve of the enzyme proton pump. Here we examine how sensitive the results of simulations are to the computational method used. We applied both molecular mechanic (MM) and hybrid quantum mechanic:molecular mechanic (QM:MM) methods and find that the results are qualitatively different. The results indicate that the mechanism for proton gating in CcO is still an open issue.

Molecular Basis for Differential Anion Binding and Proton Coupling in the Cl(-)/H(+) Exchanger ClC-ec1

Cl–/H+ transporters of the CLC superfamily form a ubiquitous class of membrane proteins that catalyze stoichiometrically coupled exchange of Cl– and H+ across biological membranes. CLC transporters exchange H+ for halides and certain polyatomic anions, but exclude cations, F–, and larger physiological anions, such as PO43– and SO42–. Despite comparable transport rates of different anions, the H+ coupling in CLC transporters varies significantly depending on the chemical nature of the transported anion. Although the molecular mechanism of exchange remains unknown, studies on bacterial ClC-ec1 transporter revealed that Cl– binding to the central anion-binding site (Scen) is crucial for the anion-coupled H+ transport. Here, we show that Cl–, F–, NO3–, and SCN– display distinct binding coordinations at the Scen site and are hydrated in different manners. Consistent with the observation of differential bindings, ClC-ec1 exhibits markedly variable ability to support the formation of the transient water wires, which are necessary to support the connection of the two H+ transfer sites (Gluin and Gluex), in the presence of different anions. While continuous water wires are frequently observed in the presence of physiologically transported Cl–, binding of F– or NO3– leads to the formation of pseudo-water-wires that are substantially different from the wires formed with Cl–. Binding of SCN–, however, eliminates the water wires altogether. These findings provide structural details of anion binding in ClC-ec1 and reveal a putative atomic-level mechanism for the decoupling of H+ transport to the transport of anions other than Cl–.

Microscopic basis for kinetic gating in Cytochrome c oxidase: insights from QM/MM analysis

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps.

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps. For the proton pump cytochrome c oxidase, while the regulation of thermodynamic driving force for key proton transfers has been discussed in great detail, the microscopic basis for the control of proton transfer kinetics has been poorly understood. Here we carry out extensive QM/MM free energy simulations to probe the kinetics of relevant proton transfer steps and analyze the effects of local structure and hydration level. We show that protonation of the proton loading site (PLS, taken to be a propionate of heme a₃) requires a concerted process in which a key glutamic acid (Glu286H) delivers the proton to the PLS while being reprotonated by an excess proton coming from the D-channel. The concerted nature of the mechanism is a crucial feature that enables the loading of the PLS before the cavity containing Glu286 is better hydrated to lower its pK_a to experimentally measured range; the charged rather than dipolar nature of the process also ensures a tight coupling with heme a reduction, as emphasized by Siegbahn and Blomberg. In addition, we find that rotational flexibility of the PLS allows its protonation before that of the binuclear center (the site where oxygen gets reduced to water). Together with our recent study (P. Goyal, et al., Proc. Natl. Acad. Sci. U. S. A., 2013, 110, 18886–18891) that focused on the modulation of Glu286 pK_a, the current work suggests a mechanism that builds in a natural sequence for the protonation of the PLS prior to that of the binuclear center. This provides microscopic support to the kinetic constraints revealed by kinetic network analysis as essential elements that ensure an efficient vectorial proton transport in cytochrome c oxidase.

Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site

Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O₂) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O₂ to travel from the solvent/membrane to the heme a₃-Cu_B binuclear center (BNC). Moreover, the identification of all possible O₂ channels as well as the atomic details of O₂ diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O₂ distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O₂ channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O₂ molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O₂ inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the BNC. Both alternative O₂ channels start in the membrane spanning region and terminate close to Y288_I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O₂ to reach the BNC are the alternative channels, rather than the X-ray inferred pathway.

Cytochrome c oxidases (Ccoxs), the terminal enzymes of the respiratory electron transport chain in eukaryotes and many prokaryotes, are key enzymes in aerobic respiration. These proteins couple the reduction of molecular dioxygen to water with the creation of a transmembrane electrochemical proton gradient. Over the last decades, most of the Ccoxs research focused on the mechanisms and energetics of reduction and/or proton pumping, and little emphasis has been given to the pathways used by dioxygen to reach the binuclear center, where dioxygen reduction takes place. In particular, the existence and the characteristics of the channel(s) used by O₂ to travel from the solvent/membrane to the binuclear site are still unclear. In this work, we combine all-atom molecular dynamics simulations and implicit ligand sampling calculations in order to identify and characterize the O₂ delivery channels in the Ccox from Rhodobacter sphaeroides. Altogether, our results suggest that, in this Ccox, O₂ can diffuse via three well-defined channels that start in membrane region (where O₂ solubility is higher than in the water). One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the binuclear center.

Mutations in the D-channel of cytochrome c oxidase causes leakage of the proton pump

It has experimentally been found that certain mutations close to the entry point of the proton transfer channel in cytochrome c oxidase stop proton translocation but not the oxygen reduction chemistry. This effect is termed uncoupling. Since the mutations are 20Å away from the catalytic center, this is very surprising. A new explanation for this phenomenon is suggested here, involving a local effect at the entry point of the proton channel, rather than the long range effects suggested earlier.

Detection of an intermediary, protonated water cluster in photosynthetic oxygen evolution

In photosynthesis, photosystem II evolves oxygen from water by the accumulation of photooxidizing equivalents at the oxygen-evolving complex (OEC). The OEC is a Mn4CaO5 cluster, and its sequentially oxidized states are termed the Sn states. The dark-stable state is S1, and oxygen is released during the transition from S3 to S0. In this study, a laser flash induces the S1 to S2 transition, which corresponds to the oxidation of Mn(III) to Mn(IV). A broad infrared band, at 2,880 cm(-1), is produced during this transition. Experiments using ammonia and (2)H2O assign this band to a cationic cluster of internal water molecules, termed "W5(+)." Observation of the W5(+) band is dependent on the presence of calcium, and flash dependence is observed. These data provide evidence that manganese oxidation during the S1 to S2 transition results in a coupled proton transfer to a substrate-containing, internal water cluster in the OEC hydrogen-bonded network.

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N₂O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666–70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb₃ oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.

Denitrification is an anaerobic process performed by several bacteria as an alternative to aerobic respiration. A key intermediate step is catalyzed by the nitric oxide reductase (NOR) enzyme, which is situated in the cytoplasmic membrane. Proton delivery to the catalytic site inside NOR is an important part of its functioning. In this work we use molecular dynamics simulations to describe water distribution and to identify proton transfer pathways in cNOR. Our results reveal two channels from the periplasmic side of the membrane and none from the cytoplasmic side, indicating that cNOR is not a proton pump. It is our hope that these results will provide a basis for further experimental and computational studies aimed to understand details of the NOR mechanism. Furthermore, this work sheds light on the molecular evolution of respiratory enzymes.

Expanding the view of proton pumping in cytochrome c oxidase through computer simulation

In cytochrome c oxidase (CcO), a redox-driven proton pump, protons are transported by the Grotthuss shuttling via hydrogen-bonded water molecules and protonatable residues. Proton transport through the D-pathway is a complicated process that is highly sensitive to alterations in the amino acids or the solvation structure in the channel, both of which can inhibit proton pumping and enzymatic activity. Simulations of proton transport in the hydrophobic cavity showed a clear redox state dependence. To study the mechanism of proton pumping in CcO, multi-state empirical valence bond (MS-EVB) simulations have been conducted, focusing on the proton transport through the D-pathway and the hydrophobic cavity next to the binuclear center. The hydration structures, transport pathways, effects of residues, and free energy surfaces of proton transport were revealed in these MS-EVB simulations. The mechanistic insight gained from them is herein reviewed and placed in context for future studies.

Gating and regulation of the cytochrome c oxidase proton pump

As a consumer of 95% of the oxygen we breathe, cytochrome c oxidase plays a major role in the energy balance of the cell. Regulation of its oxygen reduction and proton pumping activity is therefore critical to physiological function in health and disease. The location and structure of pathways for protons that are required to support cytochrome c oxidase activity are still under debate, with respect to their requirements for key residues and fixed waters, and how they are gated to prevent (or allow) proton backflow. Recent high resolution structures of bacterial and mammalian forms reveal conserved lipid and steroid binding sites as well as redox-linked conformational changes that provide new insights into potential regulatory ligands and gating modes. Mechanistic interpretation of these findings and their significance for understanding energy regulation is discussed.

Glu-286 rotation and water wire reorientation are unlikely the gating elements for proton pumping in cytochrome C oxidase

One of the key unresolved issues regarding proton pumping in cytochrome c oxidase (CcO) is the identity of the gating element that prevents the backflow of protons. In this study, we analyze two popular proposals for this element: isomerization of the key branching residue (Glu-286) and (re)orientation of water molecules in the hydrophobic cavity. Using a multifaceted set of computational analyses that involve CcO embedded in either an implicit or explicit treatment of lipid membrane, we show that neither Glu-286 nor active-site water likely constitutes the gating element. Detailed energetic and structural analyses of the simulation results indicate that the gating-relevant properties of these structural motifs observed in previous work are likely a result of the simplified computational models employed in those studies.

The mechanism for proton pumping in cytochrome c oxidase from an electrostatic and quantum chemical perspective

The mechanism for proton pumping in cytochrome c oxidase in the respiratory chain, has for decades been one of the main unsolved problems in biochemistry. However, even though several different suggested mechanisms exist, many of the steps in these mechanisms are quite similar and constitute a general consensus framework for discussing proton pumping. When these steps are analyzed, at least three critical gating situations are found, and these points are where the suggested mechanisms in general differ. The requirements for gating are reviewed and analyzed in detail, and a mechanism is suggested, where solutions for all the gating situations are formulated. This mechanism is based on an electrostatic analysis of a kinetic experiment fior the O to E transition. The key component of the mechanism is a positively charged transition state. An electron on heme a opens the gate for proton transfer from the N-side to a pump loading site (PLS). When the negative charge of the electron is compensated by a chemical proton, the positive transition state prevents backflow from the PLS to the N-side at the most critical stage of the pumping process. The mechanism has now been tested by large model DFT calculations, and these calculations give strong support for the suggested mechanism.

Energetics and dynamics of proton transfer reactions along short water wires

Proton transfer (pT) reactions in biochemical processes are often mediated by chains of hydrogen-bonded water molecules. We use hybrid density functional calculations to study pT along quasi one-dimensional water arrays that connect an imidazolium-imidazole proton donor-acceptor pair. We characterize the structures of intermediates and transition states, the energetics, and the dynamics of the pT reactions, including vibrational contributions to kinetic isotope effects. In molecular dynamics simulations of pT transition paths, we find that for short water chains with four water molecules, the pT reactions are semi-concerted. The formation of a high-energy hydronium intermediate next to the proton-donating group is avoided by a simultaneous transfer of a proton from the donor to the first water molecule, and from the first water molecule into the water chain. Lowering the dielectric constant of the environment and increasing the water chain length both reduce the barrier for pT. We study the effect of the driving force on the energetics of the pT reaction by changing the proton affinity of the donor and acceptor groups through halogen and methyl substitutions. We find that the barrier of the pT reaction depends linearly on the proton affinity of the donor but is nearly independent of the proton affinity of the acceptor, corresponding to Brønsted slopes of one and zero, respectively.

Exploration of the cytochrome c oxidase pathway puzzle and examination of the origin of elusive mutational effects

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is a problem of great current interest. Despite promising mechanistic proposals, so far, a physically consistent model that would reproduce all the relevant barriers needed to create a working pump has not been presented. In addition, there are major problems in elucidating the origin of key mutational effects and in understanding the nature of the apparent pK(a) values associated with the pH dependencies of specific proton transfer (PT) reactions in CcO. This work takes a key step in resolving the above problems, by considering mutations, such as the Asn139Asp replacement, that blocks proton pumping without affecting PT to the catalytic site. We first introduce a formulation that makes it possible to relate the apparent pK(a) of Glu286 to different conformational states of this residue. We then use the new formulation along with the calculated pK(a) values of Glu286 at these different conformations to reproduce the experimentally observed apparent pK(a) of the residue. Next, we take the X-ray structures of the native and Asn139Asp mutant of the Paracoccus denitrificans CcO (N131D in this system) and reproduce for the first time the change in the primary PT pathways (and other key features) based on simulations that start with the observed structural changes. We also consider the competition between proton transport to the catalytic site and the pump site, as a function of the bulk pH, as well as the H/D isotope effect, and use this information to explore the relative height of the two barriers. The paper emphasizes the crucial role of energy-based considerations that include the PT process, and the delicate control of PT in CcO.

Crystallographic and online spectral evidence for role of conformational change and conserved water in cytochrome oxidase proton pump

Crystal structures in both oxidized and reduced forms are reported for two bacterial cytochrome c oxidase mutants that define the D and K proton paths, showing conformational change in response to reduction and the loss of strategic waters that can account for inhibition of proton transfer. In the oxidized state both mutants of the Rhodobacter sphaeroides enzyme, D132A and K362M, show overall structures similar to wild type, indicating no long-range effects of mutation. In the reduced state, the mutants show an altered conformation similar to that seen in reduced wild type, confirming this reproducible, reversible response to reduction. In the strongly inhibited D132A mutant, positions of residues and waters in the D pathway are unaffected except in the entry region close to the mutation, where a chloride ion replaces the missing carboxyl and a 2-Å shift in N207 results in loss of its associated water. In K362M, the methionine occupies the same position as the original lysine, but K362- and T359-associated waters in the wild-type structure are missing, likely accounting for the severe inhibition. Spectra of oxidized frozen crystals taken during X-ray radiation show metal center reduction, but indicate development of a strained configuration that only relaxes to a native form upon annealing. Resistance of the frozen crystal to structural change clarifies why the oxidized conformation is observable and supports the conclusion that the reduced conformation has functional significance. A mechanism is described that explains the conformational change and the incomplete response of the D-path mutant.

Intricate role of water in proton transport through cytochrome c oxidase

Cytochrome c oxidase (CytcO), the final electron acceptor in the respiratory chain, catalyzes the reduction of O(2) to H(2)O while simultaneously pumping protons across the inner mitochondrial or bacterial membrane to maintain a transmembrane electrochemical gradient that drives, for example, ATP synthesis. In this work mutations that were predicted to alter proton translocation and enzyme activity in preliminary computational studies are characterized with extensive experimental and computational analysis. The mutations were introduced in the D pathway, one of two proton-uptake pathways, in CytcO from Rhodobacter sphaeroides . Serine residues 200 and 201, which are hydrogen-bonded to crystallographically resolved water molecules halfway up the D pathway, were replaced by more bulky hydrophobic residues (Ser200Ile, Ser200Val/Ser201Val, and Ser200Val/Ser201Tyr) to query the effects of changing the local structure on enzyme activity as well as proton uptake, release, and intermediate transitions. In addition, the effects of these mutations on internal proton transfer were investigated by blocking proton uptake at the pathway entrance (Asp132Asn replacement in addition to the above-mentioned mutations). Even though the overall activities of all mutant CytcO's were lowered, both the Ser200Ile and Ser200Val/Ser201Val variants maintained the ability to pump protons. The lowered activities were shown to be due to slowed oxidation kinetics during the P(R) → F and F → O transitions (P(R) is the "peroxy" intermediate formed at the catalytic site upon reaction of the four-electron-reduced CytcO with O(2), F is the oxoferryl intermediate, and O is the fully oxidized CytcO). Furthermore, the P(R) → F transition is shown to be essentially pH independent up to pH 12 (i.e., the apparent pK(a) of Glu286 is increased from 9.4 by at least 3 pK(a) units) in the Ser200Val/Ser201Val mutant. Explicit simulations of proton transport in the mutated enzymes revealed that the solvation dynamics can cause intriguing energetic consequences and hence provide mechanistic insights that would never be detected in static structures or simulations of the system with fixed protonation states (i.e., lacking explicit proton transport). The results are discussed in terms of the proton-pumping mechanism of CytcO.

Secondary water pore formation for proton transport in a ClC exchanger revealed by an atomistic molecular-dynamics simulation

Several prokaryotic ClC proteins have been demonstrated to function as exchangers that transport both chloride ions and protons simultaneously in opposite directions. However, the path of the proton through the ClC exchanger, and how the protein brings about the coupled movement of both ions are still unknown. In this work, we use an atomistic molecular dynamics (MD) simulation to demonstrate that a previously unknown secondary water pore is formed inside an Escherichia coli ClC exchanger. The secondary water pore is bifurcated from the chloride ion pathway at E148. From the systematic simulations, we determined that the glutamate residue exposed to the intracellular solution, E203, plays an important role as a trigger for the formation of the secondary water pore, and that the highly conserved tyrosine residue Y445 functions as a barrier that separates the proton from the chloride ion pathways. Based on our simulation results, we conclude that protons in the ClC exchanger are conducted via a water network through the secondary water pore, and we propose a new mechanism for the coupled transport of chloride ions and protons. It has been reported that several members of ClC proteins are not just channels that simply transport chloride ions across lipid bilayers; rather, they are exchangers that transport both the chloride ion and proton in opposite directions. However, the ion transit pathways and the mechanism of the coupled movement of these two ions have not yet been unveiled. In this article, we report a new finding (to our knowledge) of a water pore inside a prokaryotic ClC protein as revealed by computer simulation. This water pore is bifurcated from the putative chloride ion, and water molecules inside the new pore connect two glutamate residues that are known to be key residues for proton transport. On the basis of our simulation results, we conclude that the water wire that is formed inside the newly found pore acts as a proton pathway, which enables us to resolve many problems that could not be addressed by previous experimental studies.

Simulating proton transport through a simplified model for trans-membrane proteins

Ab initio MD simulations on a polyglycine helix and water-wire expressed under periodic boundary conditions have created a channel that supports proton transfer up to distances of 10.5 A. The effect of varying the density of water molecules in the channel has been investigated. A range of cationic states are identified with widely varying lifetimes. The mechanism of proton transport in this model shares some features with the simulations reported for bulk water, with, e.g., the hydrogen bond distance shortening in the time period leading up to successful proton transfer. However, there are also some important differences such as the observation of a heightened number of proton rattling events. We also observe that the helix plays an important role in directing the behavior of the water wire: the most active proton transport regions of the water-wire are found in areas where the helix is most tightly coiled. Finally, we report on the effects of different DFT functionals to model a water-wire and on the importance of including dispersion corrections to stabilize the alpha-helical structure.

Partial steps of charge translocation in the nonpumping N139L mutant of Rhodobacter sphaeroides cytochrome c oxidase with a blocked D-channel

The N139L substitution in the D-channel of cytochrome oxidase from Rhodobacter sphaeroides results in an approximately 15-fold decrease in the turnover number and a loss of proton pumping. Time-resolved absorption and electrometric assays of the F --> O transition in the N139L mutant oxidase result in three major findings. (1) Oxidation of the reduced enzyme by O(2) shows approximately 200-fold inhibition of the F --> O step (k approximately 2 s(-1) at pH 8) which is not compatible with enzyme turnover ( approximately 30 s(-1)). Presumably, an abnormal intermediate F(deprotonated) is formed under these conditions, one proton-deficient relative to a normal F state. In contrast, the F --> O transition in N139L oxidase induced by single-electron photoreduction of intermediate F, generated by reaction of the oxidized enzyme with H(2)O(2), decelerates to an extent compatible with enzyme turnover. (2) In the N139L mutant, the protonic phase of Deltapsi generation coupled to the flash-induced F --> O transition greatly decreases in rate and magnitude and can be assigned to the movement of a proton from E286 to the binuclear site, required for reduction of heme a(3) from the Fe(4+) horizontal lineO(2-) state to the Fe(3+)-OH(-) state. Electrogenic reprotonation of E286 from the inner aqueous phase is missing from the F --> O step in the mutant. (3) In the N139L mutant, the KCN-insensitive rapid electrogenic phase may be composed of two components with lifetimes of approximately 10 and approximately 40 mus and a magnitude ratio of approximately 3:2. The 10 mus phase matches vectorial electron transfer from Cu(A) to heme a, whereas the 40 mus component is assigned to intraprotein proton displacement across approximately 20% of the membrane dielectric thickness. This proton displacement might be triggered by rotation of the charged K362 side chain coupled to heme a reduction. The two components of the rapid electrogenic phase have been resolved subsequently with other D-channel mutants as well as with cyanide-inhibited wild-type oxidase. The finding helps to reconcile the unusually high relative contribution of the microsecond electrogenic phase in the bacterial enzyme ( approximately 30%) with the net electrogenicity of the F --> O transition coupled to transmembrane transfer of two charges per electron.

B3LYP study on reduction mechanisms from O2 to H2O at the catalytic sites of fully reduced and mixed-valence bovine cytochrome c oxidases

Reduction mechanisms of oxygen molecule to water molecules in the fully reduced (FR) and mixed-valence (MV) bovine cytochrome c oxidases (CcO) have been systematically examined based on the B3LYP calculations. The catalytic cycle using four electrons and four protons has been also shown consistently. The MV CcO catalyses reduction to produce one water molecule, while the FR CcO catalyses to produce two water molecules. One water molecule is added into vacant space between His240 and His290 in the catalytic site. This water molecule constructs the network of hydrogen bonds of Tyr244, farnesyl ethyl, and Thr316 that is a terminal residue of the K-pathway. It plays crucial roles for the proton transfer to the dioxygen to produce the water molecules in both MV and FR CcOs. Tyr244 functions as a relay of the proton transfer from the K-pathway to the added water molecule, not as donors of a proton and an electron to the dioxygen. The reduction mechanisms of MV and FR CcOs are strictly distinguished. In the FR CcO, the Cu atom at the Cu(B) site maintains the reduced state Cu(I) during the process of formation of first water molecule and plays an electron storage. At the final stage of formation of first water molecule, the Cu(I) atom releases an electron to Fe-O. During the process of formation of second water molecule, the Cu atom maintains the oxidized state Cu(II). In contrast with experimental proposals, the K-pathway functions for formation of first water molecule, while the D-pathway functions for second water molecule. The intermediates, P(M), P(R), F, and O, obtained in this work are compared with those proposed experimentally.

Microscopic pKa analysis of Glu286 in cytochrome c oxidase (Rhodobacter sphaeroides): toward a calibrated molecular model

As stringent tests for the molecular model and computational protocol, microscopic pK(a) calculations are performed for the key residue, Glu286, in cytochrome c oxidase (CcO) using a combined quantum mechanical/molecular mechanical (QM/MM) potential and a thermodynamic integration protocol. The impact of the number of water molecules in the hydrophobic cavity and protonation state of several key residues (e.g., His334, Cu(B)-bound water, and PRD(a3)) on the computed microscopic pK(a) values of Glu286 has been systematically examined. To help evaluate the systematic errors in the QM/MM-based protocol, microscopic pK(a) calculations have also been carried out for sites in a soluble protein (Asp70 in T4 lysozyme) and a better-characterized membrane protein (Asp85 in bacteriorhodopsin). Overall, the results show a significant degree of internal consistency and reproducibility that support the effectiveness of the computational framework. Although the number of water molecules in the hydrophobic cavity does not greatly influence the computed pK(a) of Glu286, the protonation states of several residues, some of which are rather far away, have more significant impacts. Adopting the standard protonation state for all titratable residues leaves a large net charge on the system and a significantly elevated pK(a) for Glu286, highlighting that any attempt to address the energetics of proton transfers in CcO at a microscopic level should carefully select the protonation state of residues, even those not in the immediate neighborhood of the active site. The calculations indirectly argue against the deprotonation of His334 for the proton pumping process, although further studies that explicitly compute its pK(a) are required for a more conclusive statement. Finally, the deprotonated Glu286 is found to be in a stable water-mediated connection with PRD(a3) for at least several nanoseconds when this presumed pumping site is protonated. This does not support the proposed role of Glu286 as a robust gating valve that prevents proton leakage, although a conclusive statement awaits a more elaborate characterization of the Glu286-PRD(a3) connectivity with free energy simulations and a protonated PRD(a3). The large sets of microscopic simulations performed here have provided useful guidance to the establishment of a meaningful molecular model and effective computational protocol for explicitly analyzing the proton transfer kinetics in CcO, which is required for answering key questions regarding the pumping function of this fascinating and complex system.

A pathogenic mutation in cytochrome c oxidase results in impaired proton pumping while retaining O(2)-reduction activity

In this work we have investigated the effect of a pathogenic mitochondrial DNA mutation found in human colon cells, at a functional-molecular level. The mutation results in the amino-acid substitution Tyr19His in subunit I of the human CytcO and it is associated with respiratory deficiency. It was introduced into Rhodobacter sphaeroides, which carries a cytochrome c oxidase (cytochrome aa(3)) that serves as a model of the mitochondrial counterpart. The residue is situated in the middle of a pathway that is used to transfer substrate protons as well as protons that are pumped across the membrane. The Tyr33His (equivalent residue in the bacterial CytcO) structural variant of the enzyme was purified and its function was investigated. The results show that in the structurally altered CytcO the activity decreased due to slowed proton transfer; proton transfer from an internal proton donor, the highly-conserved Glu286, to the catalytic site was slowed by a factor of approximately 5, while reprotonation of the Glu from solution was slowed by a factor of approximately 40. In addition, in the structural variant proton pumping was completely impaired. These results are explained in terms of introduction of a barrier for proton transfer through the D pathway and changes in the coordination of water molecules surrounding the Glu286 residue. The study offers an explanation, at the molecular level, to the link between a specific amino-acid substitution and a pathogenic phenotype identified in human colon cells.

The cytochrome ba3 oxygen reductase from Thermus thermophilus uses a single input channel for proton delivery to the active site and for proton pumping

The heme-copper oxygen reductases are redox-driven proton pumps that generate a proton motive force in both prokaryotes and mitochondria. These enzymes have been divided into 3 evolutionarily related groups: the A-, B- and C-families. Most experimental work on proton-pumping mechanisms has been performed with members of the A-family. These enzymes require 2 proton input pathways (D- and K-channels) to transfer protons used for oxygen reduction chemistry and for proton pumping, with the D-channel transporting all pumped protons. In this work we use site-directed mutagenesis to demonstrate that the ba(3) oxygen reductase from Thermus thermophilus, a representative of the B-family, does not contain a D-channel. Rather, it utilizes only 1 proton input channel, analogous to that of the A-family K-channel, and it delivers protons to the active site for both O2 chemistry and proton pumping. Comparison of available subunit I sequences reveals that the only structural elements conserved within the oxygen reductase families that could perform these functions are active-site components, namely the covalently linked histidine-tyrosine, the Cu(B) and its ligands, and the active-site heme and its ligands. Therefore, our data suggest that all oxygen reductases perform the same chemical reactions for oxygen reduction and comprise the essential elements of the proton-pumping mechanism (e.g., the proton-loading and kinetic-gating sites). These sites, however, cannot be located within the D-channel. These results along with structural considerations point to the A-propionate region of the active-site heme and surrounding water molecules as the proton-loading site.

The chemistry of the CuB site in cytochrome c oxidase and the importance of its unique His-Tyr bond

The CuB metal center is at the core of the active site of the heme-copper oxidases, comprising a copper atom ligating three histidine residues one of which is covalently bonded to a tyrosine residue. Using quantum chemical methodology, we have studied the CuB site in several redox and ligand states proposed to be intermediates of the catalytic cycle. The importance of the His-Tyr crosslink was investigated by comparing energetics, charge, and spin distributions between systems with and without the crosslink. The His-Tyr bond was shown to decrease the proton affinity and increase the electron affinity of both Tyr-244 and the copper. A previously unnoticed internal electronic equilibrium between the copper atom and the tyrosine was observed, which seems to be coupled to the unique structure of the system. In certain states the copper and Tyr-244 compete for the unpaired electron, the localization of which is determined by the oxygenous ligand of the copper. This electronic equilibrium was found to be sensitive to the presence of a positive charge 10 A away from the center, simulating the effect of Lys-319 in the K-pathway of proton transfer. The combined results provide an explanation for why the heme-copper oxidases need two pathways of proton uptake, and why the K-pathway is active only in the second half of the reaction cycle.

Molecular dynamics simulation of water in cytochrome c oxidase reveals two water exit pathways and the mechanism of transport

We have examined the network of connected internal cavities in cytochrome c oxidase along which water produced at the catalytic center is removed from the enzyme. Using combination of structural analysis, molecular dynamics simulations, and free energy calculations we have identified two exit pathways that connect the Mg2+ ion cavity to the outside of the enzyme. Each pathway has a well-defined bottleneck, which determines the overall rate of water traffic along the exit pathway, and a specific cooperative mechanism of passing it. One of the pathways is going via Arg438/439 (in bovine numbering) toward the CuA center, approaching closely its His204B ligand and Lys171B residue; and the other is going toward Asp364 and Thr294. Comparison of the pathways among different aa3-type enzymes shows that they are well conserved. Possible connections of the finding to redox-coupled proton pumping mechanism are discussed. We propose specific mutations near the bottlenecks of the exit pathways that can test some of our hypotheses.

High resolution crystal structure of Paracoccus denitrificans cytochrome c oxidase: new insights into the active site and the proton transfer pathways

The structure of the two-subunit cytochrome c oxidase from Paracoccus denitrificans has been refined using X-ray cryodata to 2.25 A resolution in order to gain further insights into its mechanism of action. The refined structural model shows a number of new features including many additional solvent and detergent molecules. The electron density bridging the heme a(3) iron and Cu(B) of the active site is fitted best by a peroxo-group or a chloride ion. Two waters or OH(-) groups do not fit, one water (or OH(-)) does not provide sufficient electron density. The analysis of crystals of cytochrome c oxidase isolated in the presence of bromide instead of chloride appears to exclude chloride as the bridging ligand. In the D-pathway a hydrogen bonded chain of six water molecules connects Asn131 and Glu278, but the access for protons to this water chain is blocked by Asn113, Asn131 and Asn199. The K-pathway contains two firmly bound water molecules, an additional water chain seems to form its entrance. Above the hemes a cluster of 13 water molecules is observed which potentially form multiple exit pathways for pumped protons. The hydrogen bond pattern excludes that the Cu(B) ligand His326 is present in the imidazolate form.

Functional hydration and conformational gating of proton uptake in cytochrome c oxidase

Cytochrome c oxidase couples the reduction of dioxygen to proton pumping against an electrochemical gradient. The D-channel, a 25-A-long cavity, provides the principal pathway for the uptake of chemical and pumped protons. A water chain is thought to mediate the relay of protons via a Grotthuss mechanism through the D-channel, but it is interrupted at N139 in all available crystallographic structures. We use free-energy simulations to examine the proton uptake pathway in the wild type and in single-point mutants N139V and N139A, in which redox and pumping activities are compromised. We present a general approach for the calculation of water occupancy in protein cavities and demonstrate that combining efficient sampling algorithms with long simulation times (hundreds of nanoseconds) is required to achieve statistical convergence of equilibrium properties in the protein interior. The relative population of different conformational and hydration states of the D-channel is characterized. Results shed light on the role of N139 in the mechanism of proton uptake and clarify the physical basis for inactive phenotypes. The conformational isomerization of the N139 side chain is shown to act as a gate controlling the formation of a functional water chain or "proton wire." In the closed state of N139, the spatial distribution of water in the D-channel is consistent with available crystallographic models. However, a metastable state of N139 opens up a narrow bottleneck in which 50% occupancy by a water molecule establishes a proton pathway throughout the D-channel. Results for N139V suggest that blockage of proton uptake resulting from persistent interruption of the water pathway is the cause of this mutant's marginal oxidase activity. In contrast, results for N139A indicate that the D-channel is a continuously hydrated cavity, implying that the decoupling of oxidase activity from proton pumping measured in this mutant is not due to interruption of the proton relay chain.

Amino acids with an intermolecular proton bond as proton storage site in bacteriorhodopsin

The positions of protons are not available in most high-resolution structural data of biomolecules, thus the identity of proton storage sites in biomolecules that transport proton is generally difficult to determine unambiguously. Using combined quantum mechanical/molecular mechanical computations, we demonstrate that a pair of conserved glutamate residues (Glu 194/204) bonded by a delocalized proton is the proton release group that has been long sought in the proton pump, bacteriorhodopsin. This model is consistent with all available experimental structural and infrared data for both the wild-type bacteriorhodopsin and several mutants. In particular, the continuum infrared band in the 1,800- to 2,000-cm(-1) region is shown to arise due to the partially delocalized nature of the proton between the glutamates in the wild-type bacteriorhodopsin; alternations in the flexibility of the glutamates and electrostatic nature of nearby residues in various mutants modulate the degree of proton delocalization and therefore intensity of the continuum band. The strong hydrogen bond between Glu 194/204 also significantly shifts the carboxylate stretches of these residues well <1,700 cm(-1), which explains why carboxylate spectral shift was not observed experimentally in the typical >1,700-cm(-1) region upon proton release. By contrast, simulations with the proton restrained on the nearby water cluster, as proposed by several recent studies [see, for example, Garezarek K, Gerwert K (2006) Functional waters in intraprotein proton transfer monitored by FTIR difference spectroscopy. Nature 439:109], led to significant structural deviations from available X-ray structures. This study establishes a biological function for strong, low-barrier hydrogen bonds.

Cytochrome c oxidase: exciting progress and remaining mysteries

Cytochrome c oxidase generates a proton motive force by two separate mechanisms. The first mechanism is similar to that postulated by Peter Mitchell, and is based on electrons and protons used to generate water coming from opposite sides of the membrane. The second mechanism was not initially anticipated, but is now firmly established as a proton pump. A brief review of the current state of our understanding of the proton pump of cytochrome oxidase is presented. We have come a long way since the initial observation of the pump by Mårten Wikström in 1977, but a number of essential questions remain to be answered.

Toward a chemical mechanism of proton pumping by the B-type cytochrome c oxidases: application of density functional theory to cytochrome ba3 of Thermus thermophilus

A mechanism for proton pumping by the B-type cytochrome c oxidases is presented in which one proton is pumped in conjunction with the weakly exergonic, two-electron reduction of Fe-bound O 2 to the Fe-Cu bridging peroxodianion and three protons are pumped in conjunction with the highly exergonic, two-electron reduction of Fe(III)- (-)O-O (-)-Cu(II) to form water and the active oxidized enzyme, Fe(III)- (-)OH,Cu(II). The scheme is based on the active-site structure of cytochrome ba 3 from Thermus thermophilus, which is considered to be both necessary and sufficient for coupled O 2 reduction and proton pumping when appropriate gates are in place (not included in the model). Fourteen detailed structures obtained from density functional theory (DFT) geometry optimization are presented that are reasonably thought to occur during the four-electron reduction of O 2. Each proton-pumping step takes place when a proton resides on the imidazole ring of I-His376 and the large active-site cluster has a net charge of +1 due to an uncompensated, positive charge formally associated with Cu B. Four types of DFT were applied to determine the energy of each intermediate, and standard thermochemical approaches were used to obtain the reaction free energies for each step in the catalytic cycle. This application of DFT generally conforms with previously suggested criteria for a valid model (Siegbahn, P. E. M.; Blomberg, M. A. R. Chem. Rev. 2000, 100, 421-437) and shows how the chemistry of O 2 reduction in the heme a 3 -Cu B dinuclear center can be harnessed to generate an electrochemical proton gradient across the lipid bilayer.

A chemically explicit model for the mechanism of proton pumping in heme-copper oxidases

A mechanism for proton pumping is described that is based on chemiosmotic principles and the detailed molecular structures now available for cytochrome oxidases. The importance of conserved water positions and a step-wise gated process of proton translocation is emphasized, where discrete electron transfer events are coupled to proton uptake and expulsion. The trajectory of each pumped proton is the same for all four substrate electrons. An essential role for the His-Tyr cross-linked species is discussed, in gating of the D- and K-channels and as an acceptor/donor of electrons and protons at the binuclear center.

Uncovering channels in photosystem II by computer modelling: current progress, future prospects, and lessons from analogous systems

Even prior to the publication of the crystal structures for photosystem II (PSII), it had already been suggested that water, O(2) and H(+) channels exist in PSII to achieve directed transport of these molecules, and to avoid undesirable side reactions. Computational efforts to uncover these channels and investigate their properties are still at early stages, and have so far only been based on the static PSII structure. The rationale behind the proposals for such channels and the computer modelling studies thus far are reviewed here. The need to take the dynamic protein into account is then highlighted with reference to the specific issues and techniques applicable to the simulation of each of the three channels. In particular, lessons are drawn from simulation studies on other protein systems containing similar channels.