Unraveling the mechanism of proton translocation in the extracellular half‐channel of bacteriorhodopsin

Bacteriorhodopsin, a light activated protein that creates a proton gradient in halobacteria, has long served as a simple model of proton pumps. Within bacteriorhodopsin, several key sites undergo protonation changes during the photocycle, moving protons from the higher pH cytoplasm to the lower pH extracellular side. The mechanism underlying the long‐range proton translocation between the central (the retinal Schiff base SB216, D85, and D212) and exit clusters (E194 and E204) remains elusive. To obtain a dynamic view of the key factors controlling proton translocation, a systematic study using molecular dynamics simulation was performed for eight bacteriorhodopsin models varying in retinal isomer and protonation states of the SB216, D85, D212, and E204. The side‐chain orientation of R82 is determined primarily by the protonation states of the residues in the EC. The side‐chain reorientation of R82 modulates the hydrogen‐bond network and consequently possible pathways of proton transfer. Quantum mechanical intrinsic reaction coordinate calculations of proton‐transfer in the methyl guanidinium‐hydronium‐hydroxide model system show that proton transfer via a guanidinium group requires an initial geometry permitting proton donation and acceptance by the same amine. In all the bacteriorhodopsin models, R82 can form proton wires with both the CC and the EC connected by the same amine. Alternatively, rare proton wires for proton transfer from the CC to the EC without involving R82 were found in an O′ state where the proton on D85 is transferred to D212. Proteins 2016; 84:639–654. © 2016 Wiley Periodicals, Inc.     

The photosynthetic water oxidase: its proton pumping activity is short-circuited within the protein by DCCD

The photosynthetic water oxidase is composed of ˜15 polypeptides which are grouped around two functional parts: photosystem II and the catalytic manganese centre. Photochemically driven vectorial electron transfer between the manganese centre and bound plastoquinone causes deprotonation–protonation reactions at opposite sides of the thylakoid membrane. Thereby the water oxidase acts as a proton pump. Incubation of stacked thylakoids with N,N'-dicyclohexylcarbodiimide (DCCD) short-circuited its proton pumping activity. Under flashing light, the extent of both proton release into the lumen by water oxidation and of proton uptake from the medium by reduced quinone was diminished. Instead there was a rapid electrogenic backreaction with a strong H/D-isotope effect. Apparently protons which were produced by water oxidation were channelled across the transmembrane protein to the bound quinone. A more rapid protonation of the reduced quinone was evident from a shortening of the time lag for the reduction of photosystem I. These effects were paralleled by the preferential labelling with [¹⁴C]DCCD in stacked thylakoids of two polypeptides with 20 and 24 kd apparent molecular mass. These may be capping the oxidizing and the reducing terminus of the water oxidase to control proton extrusion and proton uptake respectively.     

A model for conformational coupling of membrane potential and proton translocation to ATP synthesis and to active transport.

Acceptance of a membrane potential and/or a proton gradient as a possible means of transmitting energy from oxidations to ATP synthesis rests in part on a satisfactory hypothesis for how the potential or proton gradient could drive ATP synthesis. Recognition that energy input may drive ATP synthesis by change in binding of reactants at the catalytic site has led to the suggestions presented in this paper. These are that in oxidative phosphorylation and photophosphorylation, the requisite conformational changes may be coupled to exposure of charged groups to different sides of the membrane. The cycle of charged group exposure or movement may be driven by the membrane potential or, through protonation and deprotonation, may be coupled to proton translocation across the membrane. Effects of proton gradient and membrane potential may be additive. Similar conformational coupling suggestions may explain proton translocation coupled to ATP cleavage and active transport of metabolites coupled to membrane potential, proton gradients of ATP cleavage.     

Charge requirements for proton gradient-driven translocation of anthrax toxin

Anthrax lethal toxin is used as a model system to study protein translocation. The toxin is composed of a translocase channel, called protective antigen (PA), and an enzyme, called lethal factor (LF). A proton gradient (ΔpH) can drive LF unfolding and translocation through PA channels; however, the mechanism of ΔpH-mediated force generation, substrate unfolding, and establishment of directionality are poorly understood. One recent hypothesis suggests that the ΔpH may act through changes in the protonation state of residues in the substrate. Here we report the charge requirements of LF's amino-terminal binding domain (LF(N)) using planar lipid bilayer electrophysiology. We found that acidic residues are required in LF(N) to utilize a proton gradient for translocation. Constructs lacking negative charges in the unstructured presequence of LF(N) translocate independently of the ΔpH driving force. Acidic residues markedly increase the rate of ΔpH-driven translocation, and the presequence is optimized in its natural acidic residue content for efficient ΔpH-driven unfolding and translocation. We discuss a ΔpH-driven charge state Brownian ratchet mechanism for translocation, where glutamic and aspartic acid residues in the substrate are the "molecular teeth" of the ratchet. Our Brownian ratchet model includes a mechanism for unfolding and a novel role for positive charges, which we propose chaperone negative charges through the PA channel during ΔpH translocation.     

Two conformational states of Glu242 and pKas in bovine cytochrome c oxidase.

Cytochrome c oxidase (CcO) is the terminal enzyme in the respiratory electron transport chain of aerobic organisms. It catalyses the reduction of atmospheric oxygen to water, and couples this reaction to proton pumping across the membrane; this process generates the electrochemical gradient that subsequently drives the synthesis of ATP. The molecular details of the mechanism by which electron transfer is coupled to proton pumping in CcO is poorly understood. Recent calculations from our group indicate that His291, a ligand of the Cu(B) center of the enzyme, may play the role of the pumping element. In this paper we describe calculations in which a DFT/continuum electrostatic method is used to explore the coupling of the conformational changes of Glu242 residue, the main proton donor of both chemical and pump protons, to its pKa, and the pKa of His291, a putative proton loading site of our pumping model. The computations are done for several redox states of metal centers, different protonation states of Glu242 and His291, and two well-defined conformations of the Glu242 side chain. Thus, in addition to equilibrium redox/protonation states of the catalytic cycle, we also examine the transient and intermediate states. Different dielectric models are employed to investigate the robustness of the results, and their viability in the light of the proposed proton pumping mechanism of CcO. The main results are in agreement with the experimental measurements and support the proposed pumping mechanism. Additionally, the present calculations indicate a possibility of gating through conformational changes of Glu242; namely, in the pumping step, we find that Glu242 needs to be reprotonated before His291 can eject a proton to the P-site of membrane. As a result, the reprotonation of Glu can control proton release from the proton loading site.     

Proton coupling and the multiscale kinetic mechanism of a peptide transporter

Proton-coupled peptide transporters (POTs) are crucial for the uptake of di- and tripeptides as well as drug and prodrug molecules in prokaryotes and eukaryotic cells. We illustrate from multiscale modeling how transmembrane proton flux couples within a POT protein to drive essential steps of the full functional cycle: 1) protonation of a glutamate on transmembrane helix 7 (TM7) opens the extracellular gate, allowing ligand entry; 2) inward proton flow induces the cytosolic release of ligand by varying the protonation state of a second conserved glutamate on TM10; 3) proton movement between TM7 and TM10 is thermodynamically driven and kinetically permissible via water proton shuttling without the participation of ligand. Our results, for the first time, give direct computational confirmation for the alternating access model of POTs, and point to a quantitative multiscale kinetic picture of the functioning protein mechanism.     

The protonation state of the Glu-71/Asp-80 residues in the KcsA potassium channel: a first-principles QM/MM molecular dynamics study

Although a few x-ray structures of the KcsA K(+) channel have been crystallized several issues concerning the mechanisms of the ionic permeation and the protonation state of the selectivity filter ionizable side chains are still open. Using a first-principles quantum mechanical/molecular mechanical simulation approach, we have investigated the protonation state of Glu-71 and Asp-80, two important residues located in the vicinity of the selectivity filter. Results from the dynamics show that a proton is shared between the two residues, with a slight preference for Glu-71. The proton is found to exchange on the picosecond timescale, an interesting phenomenon that cannot be observed in classical molecular dynamics. Simulations of different ionic loading states of the filter show that the probability for the proton transfer is correlated with the filter occupancy. In addition, the Glu-71/Asp-80 pair is able to modulate the potential energy profile experienced by a K(+) ion as it translates along the pore axis. These theoretical predictions, along with recent experimental results, suggest that changes of the filter structure could be associated with a shift in the Glu-Asp protonation state, which in turn would influence the ion translocation.     

Connectivity of the retinal Schiff base to Asp85 and Asp96 during the bacteriorhodopsin photocycle: the local-access model.

In the recently proposed local-access model for proton transfers in the bacteriorhodopsin transport cycle (Brown et al. 1998. Biochemistry. 37:3982-3993), connection between the retinal Schiff base and Asp85 (in the extracellular direction) and Asp96 (in the cytoplasmic direction)is maintained as long as the retinal is in its photoisomerized state. The directionality of the proton translocation is determined by influences in the protein that make Asp85 a proton acceptor and, subsequently, Asp96 a proton donor. The idea of concurrent local access of the Schiff base in the two directions is now put to a test in the photocycle of the D115N/D96N mutant. The kinetics had suggested that there is a single sequence of intermediates, L<-->M1<-->M2<-->N, and the M2-->M1 reaction depends on whether a proton is released to the extracellular surface. This is now confirmed. We find that at pH 5, where proton release does not occur, but not at higher pH, the photostationary state created by illumination with yellow light contains not only the M1 and M2 states, but also the L and the N intermediates. Because the L and M1 states decay rapidly, they can be present only if they are in equilibrium with later intermediates of the photocycle. Perturbation of this mixture with a blue flash caused depletion of the M intermediate, followed by its partial recovery at the expense of the L state. The change in the amplitude of the C=O stretch band at 1759 cm-1 demonstrated protonation of Asp85 in this process. Thus, during the reequilibration the Schiff base lost its proton to Asp85. Because the N state, also present in the mixture, arises by protonation of the Schiff base from the cytoplasmic surface, these results fulfill the expectation that under the conditions tested the extracellular access of the Schiff base would not be lost at the time when there is access in the cytoplasmic direction. Instead, the connectivity of the Schiff base flickers rapidly (with the time constant of the M1<-->M2 equilibration) between the two directions during the entire L-to-N segment of the photocycle.     

A possible role for iron-sulfur cluster N2 in proton translocation by the NADH: ubiquinone oxidoreductase (complex I)

The proton-pumping NADH:ubiquinone oxidoreductase, the respiratory complex I, couples the transfer of electrons from NADH to ubiquinone with the translocation of protons across the membrane. The enzyme mechanism is still unknown due to the lack of a high-resolution structure and its complicated composition. The complex from Escherichia coli is made up of 13 subunits called NuoA through NuoN and contains one FMN and nine iron-sulfur (Fe/S) clusters as redox groups. The pH dependence of the midpoint redox potential of the Fe/S cluster named N2 and its spin-spin interaction with ubiquinone radicals made it an ideal candidate for a key component in redox-driven proton translocation. During the past years we have assigned the subunit localization of cluster N2 to subunit NuoB by site-directed mutagenesis and predicted its ligation by molecular simulation. Redox-induced FT-IR spectroscopy has shown that its redox reaction is accompanied by the protonation and deprotonation of individual amino acid residues. These residues have been identified by site-directed mutagenesis. The enzyme catalytic activity depends on the presence of cluster N2 and is coupled with major conformational changes. From these data a model for redox-induced conformation-driven proton translocation has been derived.     

A review of the binding-change mechanism for proton-translocating transhydrogenase

Proton-translocating transhydrogenase is found in the inner membranes of animal mitochondria, and in the cytoplasmic membranes of many bacteria. It catalyses hydride transfer from NADH to NADP(+) coupled to inward proton translocation. Evidence is reviewed suggesting the enzyme operates by a "binding-change" mechanism. Experiments with Escherichia coli transhydrogenase indicate the enzyme is driven between "open" and "occluded" states by protonation and deprotonation reactions associated with proton translocation. In the open states NADP(+)/NADPH can rapidly associate with, or dissociate from, the enzyme, and hydride transfer is prevented. In the occluded states bound NADP(+)/NADPH cannot dissociate, and hydride transfer is allowed. Crystal structures of a complex of the nucleotide-binding components of Rhodospirillum rubrum transhydrogenase show how hydride transfer is enabled and disabled at appropriate steps in catalysis, and how release of NADP(+)/NADPH is restricted in the occluded state. Thermodynamic and kinetic studies indicate that the equilibrium constant for hydride transfer on the enzyme is elevated as a consequence of the tight binding of NADPH relative to NADP(+). The protonation site in the translocation pathway must face the outside if NADP(+) is bound, the inside if NADPH is bound. Chemical shift changes detected by NMR may show where alterations in protein conformation resulting from NADP(+) reduction are initiated. This article is part of a Special Issue entitled: 17th European Bioenergetics Conference (EBEC 2012).     

Protonation changes upon ligand binding to trypsin and thrombin: structural interpretation based on pK(a) calculations and ITC experiments

The protonation states of a protein and a ligand can be altered upon complex formation. Such changes can be detected experimentally by isothermal titration calorimetry (ITC). For a series of ligands binding to the serine proteases trypsin and thrombin, we previously performed an extensive ITC and crystallographic study and were able to identify protonation changes for four complexes. However, since ITC measures only the overall proton exchange, it does not provide structural insights into the functional groups involved in the proton transfer. Using Poisson-Boltzmann calculations based on our recently developed PEOE_PB charges, we compute pK(a) values for all complexes of our former study in order to reveal the residues with altered protonation states. The results indicate that His57, a member of the catalytic triad, is responsible for the most relevant pK(a) shifts leading to the experimentally detected protonation changes. This finding is in contrast to our previous assumption that the observed protonation changes occur at the carboxylic group of the ligands. The newly detected proton acceptor is used for a revised factorization of the ITC data, which is necessary whenever the protonation inventory changes upon complexation. The pK(a) values of complexes showing no protonation change in the ITC experiment are reliably predicted in most cases, whereas predictions of strongly coupled systems remain problematic.     

Identifying the proton loading site cluster in the ba3 cytochrome c oxidase that loads and traps protons

Cytochrome c Oxidase (CcO) is the terminal electron acceptor in aerobic respiratory chain, reducing O₂ to water. The released free energy is stored by pumping protons through the protein, maintaining the transmembrane electrochemical gradient. Protons are held transiently in a proton loading site (PLS) that binds and releases protons driven by the electron transfer reaction cycle. Multi-Conformation Continuum Electrostatics (MCCE) was applied to crystal structures and Molecular Dynamics snapshots of the B-type Thermus thermophilus CcO. Six residues are identified as the PLS, binding and releasing protons as the charges on heme b and the binuclear center are changed: the heme a₃ propionic acids, Asp287, Asp372, His376 and Glu126B. The unloaded state has one proton and the loaded state two protons on these six residues. Different input structures, modifying the PLS conformation, show different proton distributions and result in different proton pumping behaviors. One loaded and one unloaded protonation states have the loaded/unloaded states close in energy so the PLS binds and releases a proton through the reaction cycle. The alternative proton distributions have state energies too far apart to be shifted by the electron transfers so are locked in loaded or unloaded states. Here the protein can use active states to load and unload protons, but has nearby trapped states, which stabilize PLS protonation state, providing new ideas about the CcO proton pumping mechanism. The distance between the PLS residues Asp287 and His376 correlates with the energy difference between loaded and unloaded states.     

Interplay between local protein interactions and water bridging of a proton antenna carboxylate cluster

Proteins that bind protons at cell membrane interfaces often expose to the bulk clusters of carboxylate and histidine sidechains that capture protons transiently and, in proton transporters, deliver protons to an internal site. The protonation-coupled dynamics of bulk-exposed carboxylate clusters, also known as proton antennas, is poorly described. An essential open question is how water-mediated bridges between sidechains of the cluster respond to protonation change and facilitate transient proton storage. To address this question, here I studied the protonation-coupled dynamics at the proton-binding antenna of PsbO, a small extrinsinc subunit of the photosystem II complex, with atomistic molecular dynamics simulations and systematic graph-based analyses of dynamic protein and protein-water hydrogen-bond networks. The protonation of specific carboxylate groups is found to impact the dynamics of their local protein-water hydrogen-bond clusters. Regardless of the protonation state considered for PsbO, carboxylate pairs that can sample direct hydrogen bonding, or bridge via short hydrogen-bonded water chains, anchor to nearby basic or polar protein sidechains. As a result, carboxylic sidechains of the hypothesized antenna cluster are part of dynamic hydrogen bond networks that may rearrange rapidly when the protonation changes.     

Charge transport in the ClC-type chloride-proton anti-porter from Escherichia coli

The first chloride transporter identified in the superfamily of ClC chloride channels was from Escherichia coli (EClC) (Accardi, A., and Miller, C. (2004) Nature 427, 803-807). Pathways, energetics, and mechanism of proton and chloride translocation and their coupling are up to now unclear. To bridge the hydrophobic gap of proton transport, we modeled four stable buried waters into both subunits of the WT EClC structure. Together they form a "water wire" connecting Glu-203 with the chloride at the central site, which in turn connects to Glu-148, the hypothetical proton exit site. Assuming the transient production of hydrochloride in the central chloride binding site of EClC, the water wire could establish a transmembrane proton transport pathway starting from Glu-203 all the way downstream onto Glu-148. We demonstrated by electrostatic and quantum chemical computations that protonation of the central chloride is energetically feasible. We characterized all chloride occupancies and protonation states possibly relevant for the proton-chloride transport cycle in EClC and constructed a working model. Accordingly, EClC evolves through states involving up to two excess protons and between one and three chlorides, which was required to fulfill the experimentally observed 2:1 stoichiometry. We show that the Y445F and E203H mutants of EClC can operate similarly, thus explaining why they exhibit almost WT activity levels. The proposed mechanism of coupled chloride-proton transport in EClC is consistent with available experimental data and allows predictions on the importance of specific amino acids, which may be probed by mutation experiments.     

Protonation reactions and their coupling in bacteriorhodopsin

Light-induced changes of the proton affinities of amino acid side groups are the driving force for proton translocation in bacteriorhodopsin. Recent progress in obtaining structures of bacteriorhodopsin and its intermediates with an increasingly higher resolution, together with functional studies utilizing mutant pigments and spectroscopic methods, have provided important information on the molecular architecture of the proton transfer pathways and the key groups involved in proton transport. In the present paper I consider mechanisms of light-induced proton release and uptake and intramolecular proton transport and mechanisms of modulation of proton affinities of key groups in the framework of these data. Special attention is given to some important aspects that have surfaced recently. These are the coupling of protonation states of groups involved in proton transport, the complex titration of the counterion to the Schiff base and its origin, the role of the transient protonation of buried groups in catalysis of the chromophore's thermal isomerization, and the relationship between proton affinities of the groups and the pH dependencies of the rate constants of the photocycle and proton transfer reactions.     

The protonation state of a heme propionate controls electron transfer in cytochrome c oxidase

In cytochrome c oxidase (CcO), exergonic electron transfer reactions from cytochrome c to oxygen drive proton pumping across the membrane. Elucidation of the proton pumping mechanism requires identification of the molecular components involved in the proton transfer reactions and investigation of the coupling between internal electron and proton transfer reactions in CcO. While the proton-input trajectory in CcO is relatively well characterized, the components of the output pathway have not been identified in detail. In this study, we have investigated the pH dependence of electron transfer reactions that are linked to proton translocation in a structural variant of CcO in which Arg481, which interacts with the heme D-ring propionates in a proposed proton output pathway, was replaced with Lys (RK481 CcO). The results show that in RK481 CcO the midpoint potentials of hemes a and a(3) were lowered by approximately 40 and approximately 15 mV, respectively, which stabilizes the reduced state of Cu(A) during reaction of the reduced CcO with O(2). In addition, while the pH dependence of the F --> O rate in wild-type CcO is determined by the protonation state of two protonatable groups with pK(a) values of 6.3 and 9.4, only the high-pK(a) group influences this rate in RK481 CcO. The results indicate that the protonation state of the Arg481 heme a(3) D-ring propionate cluster having a pK(a) of approximately 6.3 modulates the rate of internal electron transfer and may act as an acceptor of pumped protons.     

The Role of Electrostatic Interactions for Cytochrome c Oxidase Function

In recent years, the enormous increase in high-resolution three-dimensional structures of proteins together with the development of powerful theoretical techniques have provided the basis for a more detailed examination of the role of electrostatics in determining the midpoint potentials of redox-active metal centers and in influencing the protonation behavior of titratable groups in proteins. Based on the coordinates of the Paracoccus denitrificans cytochrome c oxidase, we have determined the electrostatic potential in and around the protein, calculated the titration curves for all ionizable residues in the protein, and analyzed the response of the protein environment to redox changes at the metal centers. The results of this study provide insight into how charged groups can be stabilized within a low-dielectric environment and how the range of their electrostatic effects can be modulated by the protein. A cluster of 18 titratable groups around the heme a ₃–Cu_B binuclear center, including a hydroxide ion bound to the copper, was identified that accounts for most of the proton uptake associated with redox changes at the binuclear site. Predicted changes in net protonation were in reasonable agreement with experimentally determined values. The relevance of these findings in the light of possible mechanisms of redox-coupled proton movement is discussed.     

Cu XAS shows a change in the ligation of CuB upon reduction of cytochrome bo3 from Escherichia coli

Copper X-ray absorption spectroscopy (XAS) has been used to examine the structures of the Cu(II) and Cu(I) forms of the cytochrome bo3 quinol oxidase from Escherichia coli. Cytochrome bo3 is a member of the superfamily of heme-copper respiratory oxidases. Of particular interest is the fact that these enzymes function as redox-linked proton pumps, resulting in the net translocation of one H+ per electron across the membrane. The molecular mechanism of how this pump operates and the manner by which it is linked to the oxygen chemistry at the active site of the enzyme are unknown. Several proposals have featured changes in the coordination of CuB during enzyme turnover that would result in sequential protonation or deprotonation events that are key to the functioning proton pump. This would imply lability of the ligands to CuB. In this work, the structure of the protein in the immediate vicinity of CuB, in both the fully oxidized and fully reduced forms of the enzyme, has been examined by Cu XAS, a technique that is particularly sensitive to changes in metal coordination. The results show that in the oxidized enzyme, CuB(II) is four-coordinate, consistent with three imidazoles and one hydroxyl (or water). Upon reduction of the enzyme, the coordination of CuB(I) is significantly altered, consistent with the loss of one of the histidine imidazole ligands in at least a substantial fraction of the population. These data add to the credibility that changes in the ligation of CuB might occur during catalytic turnover of the enzyme and, therefore, could, in principle, be part of the mechanism of proton pumping.     

Sugar binding and protein conformational changes in lactose permease

Lactose permease is an integral membrane protein that uses the cell membrane's proton gradient for import of lactose. Based on extensive biochemical data and a substrate-bound crystal structure, intermediates involved in lactose/H(+) co-transport have been suggested. Yet, the transport mechanism, especially the coupling of protonation states of essential residues and protein conformational changes involved in the transport, is not understood. Here we report molecular-dynamics simulations of membrane-embedded lactose permease in different protonation states, both in the presence and in the absence of lactose. The results analyzed in terms of pore diameter, salt-bridge formation, and substrate motion, strongly implicate Glu(269) as one of the main proton translocation sites, whose protonation state controls several key steps of the transport process. A critical ion pair (Glu(269) and Arg(144)) was found to keep the cytoplasmic entrance open, but via a different mechanism than the currently accepted model. After protonation of Glu(269), the salt bridge between Glu(269) and Arg(144) was found to break, and Arg(144) to move away from Glu(269), establishing a new salt bridge with Glu(126); furthermore, neutralization of Glu(269) and the displacement of Arg(144) and consequently of water molecules from the interdomain region was seen to initiate the closing of the cytoplasmic half channel (2.6-4.0 A reduction in diameter in the cytoplasmic constriction region in 10 ns) by allowing hydrophobic surfaces of the N- and C-domains to fuse. Charged Glu(269) was found to strongly bind the lactose permeant, indicating that proton transfer from water or another residue to Glu(269) is a prerequisite for unbinding of lactose from the binding pocket.