Proton Release Pathways Research Articles

Theoretical elucidation of the structure, bonding, and reactivity of the CaMn4Ox clusters in the whole Kok cycle for water oxidation embedded in the oxygen evolving center of photosystem II. New molecular and quantum insights into the mechanism of the O-O bond formation.

This paper reviews our historical developments of broken-symmetry (BS) and beyond BS methods that are applicable for theoretical investigations of metalloenzymes such as OEC in PSII. The BS hybrid DFT (HDFT) calculations starting from high-resolution (HR) XRD structure in the most stable S1 state have been performed to elucidate structure and bonding of whole possible intermediates of the CaMn4Ox cluster (1) in the Si (i = 0 ~ 4) states of the Kok cycle. The large-scale HDFT/MM computations starting from HR XRD have been performed to elucidate biomolecular system structures which are crucial for examination of possible water inlet and proton release pathways for water oxidation in OEC of PSII. DLPNO CCSD(T0) computations have been performed for elucidation of scope and reliability of relative energies among the intermediates by HDFT. These computations combined with EXAFS, XRD, XFEL, and EPR experimental results have elucidated the structure, bonding, and reactivity of the key intermediates, which are indispensable for understanding and explanation of the mechanism of water oxidation in OEC of PSII. Interplay between theory and experiments have elucidated important roles of four degrees of freedom, spin, charge, orbital, and nuclear motion for understanding and explanation of the chemical reactivity of 1 embedded in protein matrix, indicating the participations of the Ca(H2O)n ion and tyrosine(Yz)-O radical as a one-electron acceptor for the O-O bond formation. The Ca-assisted Yz-coupled O-O bond formation mechanisms for water oxidation are consistent with recent XES and very recent time-resolved SFX XFEL and FTIR results.

Proton Release Reactions in the Inward H+ Pump NsXeR.

Directional ion transport across biological membranes plays a central role in many cellular processes. Elucidating the molecular determinants for vectorial ion transport is key to understanding the functional mechanism of membrane-bound ion pumps. The extensive investigation of the light-driven proton pump bacteriorhodopsin from Halobacterium salinarum(HsBR) enabled a detailed description of outward proton transport. Although the structure of inward-directed proton pumping rhodopsins is very similar to HsBR, little is known about their protonation pathway, and hence, the molecular reasons for the vectoriality of proton translocation remain unclear. Here, we employ a combined experimental and theoretical approach to tracking protonation steps in the light-driven inward proton pump xenorhodopsin from Nanosalina sp. (NsXeR). Time-resolved infrared spectroscopy reveals the transient deprotonation of D220 concomitantly with deprotonation of the retinal Schiff base. Our molecular dynamics simulations support a proton release pathway from the retinal Schiff base via a hydrogen-bonded water wire leading to D220 that could provide a putative gating point for the proton release and with allosteric interactions to the retinal Schiff base. Our findings support the key role of D220 in mediating proton release to the cytoplasmic side and provide evidence that this residue is not the primary proton acceptor of the proton transiently released by the retinal Schiff base.

Identification of amino acid residues in a proton release pathway near the bacteriochlorophyll dimer in reaction centers from Rhodobacter sphaeroides.

Insight into control of proton transfer, a crucial attribute of cellular functions, can be gained from investigations of bacterial reaction centers. While the uptake of protons associated with the reduction of the quinone is well characterized, the release of protons associated with the oxidized bacteriochlorophyll dimer has been poorly understood. Optical spectroscopy and proton release/uptake measurements were used to examine the proton release characteristics of twelve mutant reaction centers, each containing a change in an amino acid residue near the bacteriochlorophyll dimer. The mutant reaction centers had optical spectra similar to wild-type and were capable of transferring electrons to the quinones after light excitation of the bacteriochlorophyll dimer. They exhibited a large range in the extent of proton release and in the slow recovery of the optical signal for the oxidized dimer upon continuous illumination. Key roles were indicated for six amino acid residues, Thr L130, Asp L155, Ser L244, Arg M164, Ser M190, and His M193. Analysis of the results points to a hydrogen-bond network that contains these residues, with several additional residues and bound water molecules, forming a proton transfer pathway. In addition to proton transfer, the properties of the pathway are proposed to be responsible for the very slow charge recombination kinetics observed after continuous illumination. The characteristics of this pathway are compared to proton transfer pathways near the secondary quinone as well as those found in photosystem II and cytochrome c oxidase.

Structural dynamics in the water and proton channels of photosystem II during the S2 to S3 transition

Light-driven oxidation of water to molecular oxygen is catalyzed by the oxygen-evolving complex (OEC) in Photosystem II (PS II). This multi-electron, multi-proton catalysis requires the transport of two water molecules to and four protons from the OEC. A high-resolution 1.89 Å structure obtained by averaging all the S states and refining the data of various time points during the S2 to S3 transition has provided better visualization of the potential pathways for substrate water insertion and proton release. Our results indicate that the O1 channel is the likely water intake pathway, and the Cl1 channel is the likely proton release pathway based on the structural rearrangements of water molecules and amino acid side chains along these channels. In particular in the Cl1 channel, we suggest that residue D1-E65 serves as a gate for proton transport by minimizing the back reaction. The results show that the water oxidation reaction at the OEC is well coordinated with the amino acid side chains and the H-bonding network over the entire length of the channels, which is essential in shuttling substrate waters and protons.

Structural Basis for the Function of the C-Terminal Proton Release Pathway in the Calcium Pump.

The calcium pump (sarco/endoplasmic reticulum Ca2+-ATPase, SERCA) plays a major role in calcium homeostasis in muscle cells by clearing cytosolic Ca2+ during muscle relaxation. Active Ca2+ transport by SERCA involves the structural transition from a low-Ca2+ affinity E2 state toward a high-Ca2+ affinity E1 state of the pump. This structural transition is accompanied by the countertransport of protons to stabilize the negative charge and maintain the structural integrity of the transport sites and partially compensate for the positive charges of the two Ca2+ ions passing through the membrane. X-ray crystallography studies have suggested that a hydrated pore located at the C-terminal domain of SERCA serves as a conduit for proton countertransport, but the existence and function of this pathway have not yet been fully characterized. We used atomistic simulations to demonstrate that in the protonated E2 state and the absence of initially bound water molecules, the C-terminal pore becomes hydrated in the nanosecond timescale. Hydration of the C-terminal pore is accompanied by the formation of water wires that connect the transport sites with the cytosol. Water wires are known as ubiquitous proton-transport devices in biological systems, thus supporting the notion that the C-terminal domain serves as a conduit for proton release. Additional simulations showed that the release of a single proton from the transport sites induces bending of transmembrane helix M5 and the interaction between residues Arg762 and Ser915. These structural changes create a physical barrier against full hydration of the pore and prevent the formation of hydrogen-bonded water wires once proton transport has occurred through this pore. Together, these findings support the notion that the C-terminal proton release pathway is a functional element of SERCA and also provide a mechanistic model for its operation in the catalytic cycle of the pump.

Transmembrane helix 6b links proton and metal release pathways and drives conformational change in an Nramp-family transition metal transporter

Transmembrane helix 6b links proton and metal release pathways and drives conformational change in an Nramp-family transition metal transporter

Transmembrane helix 6b links proton and metal release pathways and drives conformational change in an Nramp-family transition metal transporter.

The natural resistance-associated macrophage protein (Nramp) family encompasses transition metal and proton cotransporters that are present in many organisms from bacteria to humans. Recent structures of Deinococcus radiodurans Nramp (DraNramp) in multiple conformations revealed the intramolecular rearrangements required for alternating access of the metal-binding site to the external or cytosolic environment. Here, using recombinant proteins and metal transport and cysteine accessibility assays, we demonstrate that two parallel cytoplasm-accessible networks of conserved hydrophilic residues in DraNramp, one lining the wide intracellular vestibule for metal release and the other forming a narrow proton transport pathway, are essential for metal transport. We further show that mutagenic or posttranslational modifications of transmembrane helix (TM) 6b, which structurally links these two pathways, impede normal conformational cycling and metal transport. TM6b contains two highly conserved histidines, His232 and His237 We found that different mutagenic perturbations of His232, just below the metal-binding site along the proton exit route, differentially affect DraNramp's conformational state, suggesting that His232 serves as a pivot point for conformational changes. In contrast, any replacement of His237, lining the metal exit route, locked the transporter in a transport-inactive outward-closed state. We conclude that these two histidines, and TM6b more broadly, help trigger the bulk rearrangement of DraNramp to the inward-open state upon metal binding and facilitate return of the empty transporter to an outward-open state upon metal release.

Settling the Long-Standing Debate on the Proton Storage Site of the Prototype Light-Driven Proton Pump Bacteriorhodopsin.

Despite decades of research, the location and molecular identity of the proton release group together with the subsequent proton release pathway remain controversial even for the simplest light-driven proton pump, bacteriorhodopsin, according to the most recent experiments and simulations. Yet despite this nagging lack of knowledge for the generic case, even more complex pumps are currently under investigation. The proton release group disclosed by our large-scale simulations satisfies available experimental results, especially the broad Zundel continuum absorption subject to a striking anisotropy observed only recently. Moreover, our simulations delineate the seamless pathway by which the excess proton (being stored in an ultrastrong centered H-bond involving two glutamates) is finally translocated into the extracellular medium.

The cytochrome b Zn binding amino acid residue histidine 291 is essential for ubihydroquinone oxidation at the Qo site of bacterial cytochrome bc1

The ubiquinol:cytochrome (cyt) c oxidoreductase (or cyt bc1) is an important membrane protein complex in photosynthetic and respiratory energy transduction. In bacteria such as Rhodobacter capsulatus it is constituted of three subunits: the iron-sulfur protein, cyt b and cyt c1, which form two catalytic domains, the Qo (hydroquinone (QH2) oxidation) and Qi (quinone (Q) reduction) sites. At the Qo site, the pathways of bifurcated electron transfers emanating from QH2 oxidation are known, but the associated proton release routes are not well defined. In energy transducing complexes, Zn2+ binding amino acid residues often correlate with proton uptake or release pathways. Earlier, using combined EXAFS and structural studies, we identified Zn coordinating residues of mitochondrial and bacterial cyt bc1. In this work, using the genetically tractable bacterial cyt bc1, we substituted each of the proposed Zn binding residues with non-protonatable side chains. Among these mutants, only the His291Leu substitution destroyed almost completely the Qo site catalysis without perturbing significantly the redox properties of the cofactors or the assembly of the complex. In this mutant, which is unable to support photosynthetic growth, the bifurcated electron transfer reactions that result from QH2 oxidation at the Qo site, as well as the associated proton(s) release, were dramatically impaired. Based on these findings, on the putative role of His291 in liganding Zn, and on its solvent exposed and highly conserved position, we propose that His291 of cyt b is critical for proton release associated to QH2 oxidation at the Qo site of cyt bc1.

Crystal structure of Cruxrhodopsin-3 from Haloarcula vallismortis.

Cruxrhodopsin-3 (cR3), a retinylidene protein found in the claret membrane of Haloarcula vallismortis, functions as a light-driven proton pump. In this study, the membrane fusion method was applied to crystallize cR3 into a crystal belonging to space group P321. Diffraction data at 2.1 Å resolution show that cR3 forms a trimeric assembly with bacterioruberin bound to the crevice between neighboring subunits. Although the structure of the proton-release pathway is conserved among proton-pumping archaeal rhodopsins, cR3 possesses the following peculiar structural features: 1) The DE loop is long enough to interact with a neighboring subunit, strengthening the trimeric assembly; 2) Three positive charges are distributed at the cytoplasmic end of helix F, affecting the higher order structure of cR3; 3) The cytoplasmic vicinity of retinal is more rigid in cR3 than in bacteriorhodopsin, affecting the early reaction step in the proton-pumping cycle; 4) the cytoplasmic part of helix E is greatly bent, influencing the proton uptake process. Meanwhile, it was observed that the photobleaching of retinal, which scarcely occurred in the membrane state, became significant when the trimeric assembly of cR3 was dissociated into monomers in the presence of an excess amount of detergent. On the basis of these observations, we discuss structural factors affecting the photostabilities of ion-pumping rhodopsins.

Theoretical modelling of biomolecular systems I. Large-scale QM/MM calculations of hydrogen-bonding networks of the oxygen evolving complex of photosystem II

Quantum mechanical (QM)/molecular mechanics (MM) calculations by the use of a large-scale QM model (QM Model V) have been performed to elucidate hydrogen-bonding networks and proton wires for proton release pathways (PRP) of water oxidation reaction in the oxygen evolving complex (OEC) of photosystem II (PSII). Full geometry optimisations of PRP by the QM/MM model have been carried out starting from the geometry of heavy atoms determined by the recent high-resolution X-ray diffraction (XRD) experiment of PSII refined to 1.9 Å resolution. Computational results by the QM/MM calculations have elucidated the hydrogen-bonding O···O(N) and O···H distances and O(N)–H···O angles in PRP, together with the Cl–O(N) and Cl···H distances and O(N)–H···Cl angles for chloride anions. The optimised hydrogen-bonding networks are well consistent with the XRD results and available experiments such as extended X-ray absorption fine structure, showing the reliability of channel structures of OEC of PSII revealed by the XRD experiment. The QM/MM computations have elucidated possible roles of chloride anions in the OEC of PSII. The QM/MM computational results have provided useful information for understanding and explanation of accumulated mutation experiments of key amino acid residues in the OEC of PSII. Implications of the present results are discussed in relation to three steps for theoretical modelling of water oxidation in the OEC of PSII and bio-inspired working hypotheses for developments of artificial water oxidation systems by use of 3d transition-metal complexes.

Structure of archaerhodopsin-2 at 1.8 Å resolution.

Archaerhodopsin-2 (aR2), the sole protein found in the claret membrane of Halorubrum sp. Aus-2, functions as a light-driven proton pump. In this study, structural analysis of aR2 was performed using a novel three-dimensional crystal prepared by the successive fusion of claret membranes. The crystal is made up of stacked membranes, in each of which aR2 trimers are arranged on a hexagonal lattice. This lattice structure resembles that found in the purple membrane of H. salinarum, except that lipid molecules trapped within the trimeric structure are not distributed with perfect threefold symmetry. Nonetheless, diffraction data at 1.8 Å resolution provide accurate structural information about functionally important residues. It is shown that two glutamates in the proton-release channel form a paired structure that is maintained by a low-barrier hydrogen bond. Although the structure of the proton-release pathway is highly conserved among proton-pumping archaeal rhodopsins, aR2 possesses the following peculiar structural features: (i) the motional freedom of the tryptophan residue that makes contact with the C13 methyl group of retinal is restricted, affecting the formation/decay kinetics of the L state, and (ii) the N-terminal polypeptide folds into an Ω-loop, which may play a role in organizing the higher-order structure.

Crystal structure of deltarhodopsin‐3 from Haloterrigena thermotolerans

Deltarhodopsin, a new member of the microbial rhodopsin family, functions as a light-driven proton pump. Here, we report the three-dimensional structure of deltarhodopsin (dR3) from Haloterrigena thermotolerans at 2.7 Å resolution. A crystal belonging to space group R32 (a, b = 111.71 Å, c = 198.25 Å) was obtained by the membrane fusion method. In this crystal, dR3 forms a trimeric structure as observed for bacteriorhodopsin (bR). Structural comparison of dR with bR showed that the inner part (the proton release and uptake pathways) is highly conserved. Meanwhile, residues in the protein-protein contact region are largely altered so that the diameter of the trimeric structure at the cytoplasmic side is noticeably larger in dR3. Unlike bR, dR3 possesses a helical segment at the C-terminal region that fills the space between the AB and EF loops. A significant difference is also seen in the FG loop, which is one residue longer in dR3. Another peculiar property of dR3 is a highly crowded distribution of positively charged residues on the cytoplasmic surface, which may be relevant to a specific interaction with some cytoplasmic component.

Theoretical insight in to hydrogen-bonding networks and proton wire for the CaMn4O5 cluster of photosystem II. Elongation of Mn–Mn distances with hydrogen bonds

Quantum mechanical (QM) and QM/molecular mechanics (MM) calculations of three different cluster models have been performed to shed light on hydrogen bonding networks and proton wires for proton release pathways (PRP) of water oxidation reactions in the oxygen evolving complex (OEC) of photosystem II (PSII). Positions of all the hydrogen atoms in an extended QM Model III including the second coordination sphere for the active-site CaMn4O5 complex of OEC of PSII have been optimized assuming the geometry of heavy atoms determined by the recent high-resolution X-ray diffraction (XRD) experiment of PSII refined to 1.9 A resolution. Full geometry optimizations of the first coordination sphere model (QM Model I) embedded in the Model III and QM (QM Model I plus seven water molecules, namely QM Model II)/MM models, together with full QM Model III, have also been conducted to elucidate confinement effects for geometrical parameters of the CaMn4O5 cluster by proteins. Computational results by these methods have elucidated the O⋯O(N), O⋯H distances and O(N)–H⋯O angles for hydrogen bonds in proton release paths (PRP) I and II that construct a proton wire from Asp61 toward His190. The hydrogen-bonding structures revealed have also been examined in relation to the possibilities of protonation of bridge oxygen dianions within the CaMn4O5 cluster. The optimized inter-atomic distances by QM Models I and III, together with QM(Model II)/MM, have elucidated the elongation of the Mn–Mn distances with hydrogen bonds and variations of the Mnd–O(5) length with confinement effects by protein. Implications of the computational results are discussed in relation to the available EXAFS experiments, and internal, semi-internal and external reductions of Mn ions and long Mn–Mn distances of the high-resolution (SP8) XRD, and rational design of artificial catalysts for water oxidation that are current topics in the field of OEC of PSII.

Water oxidation mechanism in photosystem II, including oxidations, proton release pathways, O―O bond formation and O2 release

The present status of DFT studies on water oxidation in photosystem II is described. It is argued that a full understanding of all steps is close. In each S-transition, the manganese that is oxidized and the proton released are strongly implicated, and structures of all intermediates have been determined. For the S2-state, recent important experimental findings support key elements of the structure and the mechanism. In this mechanism, the O―O bond is formed between an oxyl radical in the center of the cluster and an Mn-bridging μ-oxo ligand, which was suggested already in 2006. The DFT structure of the oxygen evolving complex, suggested in 2008, is very similar to the recent high-resolution X-ray structure. Some new aspects of the interaction between P680 and the OEC are suggested. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

Possible mechanisms of water splitting reaction based on proton and electron release pathways revealed for CaMn4O5 cluster of PSII refined to 1.9 Å X‐ray resolution

AbstractRecently, Umena et al. have revealed the X‐ray diffraction structure of the CaMn4O5 cluster in the oxygen evolving complex (OEC) of photosystem II (PSII) refined to 1.9 Å resolution. Their X‐ray structure has first elucidated hydrogen‐bonding networks and proton release pathways at OEC of PSII. Here, several working hypotheses (heuristic principles) for water splitting reaction are derived from their X‐ray structure for theoretical modeling. These hypotheses suggest how water can be oxidized at OEC of PSII: namely possible reaction mechanisms for the reaction. To confirm them, we have also performed broken‐symmetry (BS) UB3LYP calculations for active site models based on their XRD structure. The bond lengths of formal Mn(V)O with labile dπ‐pπ bonds are optimized to clarify possible roles of the species that are often introduced as a key intermediate in the catalytic (Kok) cycle for water splitting reaction at OEC of PSII. Location of the transition structure for the oxygen‐oxygen (OO) bond formation is also performed by the energy optimization technique. The natural orbital (NO) analysis of the UB3LYP solutions has been performed to obtain the natural molecular orbitals and their occupation numbers that have been useful for classification of localized d‐electrons, labile chemical bonds and closed‐shell (valence) orbitals. The localized d‐electrons characterized by the NO analysis are the origins for the magnetism revealed by ENDOR and other magnetic experiments. On the other hand, the nature of labile (soft) dπ‐pπ bonds responsible for the OO bond formation has been investigated on the basis of chemical indices such as effective bond order (b), diradical character (y), and spin density (Q) indices that are calculated using the orbital overlap between broken‐symmetry orbitals. These chemical indices have been calculated for the transition structure of the OO bond formation at OEC of PSII. Implications of present computational results are discussed in relation to the derived hypotheses and available accumulated experimental results. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

Role of Extracellular Glutamic Acids in the Stability and Energy Landscape of Bacteriorhodopsin

Bacteriorhodopsin (BR), a specialized nanomachine, converts light energy into a proton gradient to power Halobacterium salinarum. In this work, we analyze the mechanical stability of a BR triple mutant in which three key extracellular residues, Glu 9, Glu 194, and Glu 204, were mutated simultaneously to Gln. These three Glu residues are involved in a network of hydrogen bonds, in cation binding, and form part of the proton release pathway of BR. Changes in these features and the robust photocycle dynamics of wild-type (WT) BR are apparent when the three extracellular Glu residues are mutated to Gln. It is speculated that such functional changes of proteins go hand in hand with changes in their mechanical properties. Here, we apply single-molecule dynamic force spectroscopy to investigate how the Glu to Gln mutations change interactions, reaction pathways, and the energy barriers of the structural regions of WT BR. The altered heights and positions of individual energy barriers unravel the changes in the mechanical and the unfolding kinetic properties of the secondary structures of WT BR. These changes in the mechanical unfolding energy landscape cause the proton pump to choose unfolding pathways differently. We suggest that, in a similar manner, the changed mechanical properties of mutated BR alter the functional energy landscape favoring different reaction pathways in the light-induced proton pumping mechanism.

Dynamic aspects of extracellular loop region as a proton release pathway of bacteriorhodopsin studied by relaxation time measurements by solid state NMR

Local dynamics of interhelical loops in bacteriorhodopsin (bR), the extracellular BC, DE and FG, and cytoplasmic AB and CD loops, and helix B were determined on the basis of a variety of relaxation parameters for the resolved 13C and 15N signals of [1- 13C]Tyr-, [ 15N]Pro- and [1- 13C]Val-, [ 15N]Pro-labeled bR. Rotational echo double resonance (REDOR) filter experiments were used to assign [1- 13C]Val-, [ 15N]Pro signals to the specific residues in bR. The previous assignments of [1- 13C]Val-labeled peaks, 172.9 or 171.1 ppm, to Val69 were revised: the assignment of peak, 172.1 ppm, to Val69 was made in view of the additional information of conformation-dependent 15N chemical shifts of Pro bonded to Val in the presence of 13C– 15N correlation, although no assignment of peak is feasible for 13C nuclei not bonded to Pro. 13C or 15N spin–lattice relaxation times ( T 1), spin–spin relaxation times under the condition of CP-MAS ( T 2), and cross relaxation times ( T CH and T NH) for 13C and 15N nuclei and carbon or nitrogen-resolved, 1H spin–lattice relaxation times in the rotating flame ( 1H T 1ρ) for the assigned signals were measured in [1- 13C]Val-, [ 15N]Pro-bR. It turned out that V69–P70 in the BC loop in the extracellular side has a rigid β-sheet in spite of longer loop and possesses large amplitude motions as revealed from 13C and 15N conformation-dependent chemical shifts and T 1, T 2, 1H T 1ρ and cross relaxation times. In addition, breakage of the β-sheet structure in the BC loop was seen in bacterio-opsin (bO) in the absence of retinal.

X-Ray Absorption Studies of Zn 2+ Binding Sites in Bacterial, Avian, and Bovine Cytochrome bc1 Complexes

Binding of Zn 2+ has been shown previously to inhibit the ubiquinol cytochrome c oxidoreductase (cyt bc 1 complex). X-ray diffraction data in Zn-treated crystals of the avian cyt bc 1 complex identified two binding sites located close to the catalytic Q o site of the enzyme. One of them (Zn01) might interfere with the egress of protons from the Q o site to the aqueous phase. Using Zn K-edge x-ray absorption fine-structure spectroscopy, we report here on the local structure of Zn 2+ bound stoichiometrically to noncrystallized cyt bc 1 complexes. We performed a comparative x-ray absorption fine-structure spectroscopy study by examining avian, bovine, and bacterial enzymes. A large number of putative clusters, built by combining information from first-shell analysis and metalloprotein databases, were fitted to the experimental spectra by using ab initio simulations. This procedure led us to identify the binding clusters with high levels of confidence. In both the avian and bovine enzyme, a tetrahedral ligand cluster formed by two His, one Lys, and one carboxylic residue was found, and this ligand attribution fit the crystallographic Zn01 location of the avian enzyme. In the chicken enzyme, the ligands were the His 121, His 268, Lys 270, and Asp 253 residues, and in the homologous bovine enzyme they were the His 121, His 267, Lys 269, and Asp 254 residues. Zn 2+ bound to the bacterial cyt bc 1 complex exhibited quite different spectral features, consistent with a coordination number of 6. The best-fit octahedral cluster was formed by one His, two carboxylic acids, one Gln or Asn residue, and two water molecules. It was interesting that by aligning the crystallographic structures of the bacterial and avian enzymes, this group of residues was found located in the region homologous to that of the Zn01 site. This cluster included the His 276, Asp 278, Glu 295, and Asn 279 residues of the cyt b subunit. The conserved location of the Zn 2+ binding sites at the entrance of the putative proton release pathways, and the presence of His residues point to a common mechanism of inhibition. As previously shown for the photosynthetic bacterial reaction center, zinc would compete with protons for binding to the His residues, thus impairing their function as proton donors/acceptors.

Water dynamics simulation as a tool for probing proton transfer pathways in a heptahelical membrane protein

The proton transfer pathway in a heptahelical membrane protein, the light-driven proton pump bacteriorhodopsin (BR), is probed by a combined approach of structural analysis of recent X-ray models and molecular dynamics (MD) simulations that provide the diffusion pathways of internal and external water molecules. Analyzing the hydrogen-bond contact frequencies of the water molecules with protein groups, the complete proton pathway through the protein is probed. Beside the well-known proton binding sites in the protein interior-the protonated Schiff base, Asp85 and Asp96, and the H(5)O(2) (+) complex stabilized by Glu204 and Glu194-the proton release and uptake pathways to the protein surfaces are described in great detail. Further residues were identified, by mutation of which the proposed pathways can be verified. In addition the diffusion pathway of water 502 from Lys216 to Asp96 is shown to cover the positions of the intruding waters 503 and 504 in the N-intermediate. The transiently established water chain in the N-state provides a proton pathway from Asp96 to the Schiff base in the M- to N-transition in a Grotthus-like mechanism, as concluded earlier from time-resolved Fourier transform infrared experiments [le Coutre et al., Proc Nat Acad Sci USA 1995;92:4962-4966].