Conformationally controlled pK-switching in membrane proteins: one more mechanism specific to the enzyme catalysis?

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

Modeling of interaction between cytochrome c and the WD domains of Apaf-1: bifurcated salt bridges underlying apoptosome assembly

BACKGROUND

Binding of cytochrome c, released from the damaged mitochondria, to the apoptotic protease activating factor 1 (Apaf-1) is a key event in the apoptotic signaling cascade. The binding triggers a major domain rearrangement in Apaf-1, which leads to oligomerization of Apaf-1/cytochrome c complexes into an apoptosome. Despite the availability of crystal structures of cytochrome c and Apaf-1 and cryo-electron microscopy models of the entire apoptosome, the binding mode of cytochrome c to Apaf-1, as well as the nature of the amino acid residues of Apaf-1 involved remain obscure.

RESULTS

We investigated the interaction between cytochrome c and Apaf-1 by combining several modeling approaches. We have applied protein-protein docking and energy minimization, evaluated the resulting models of the Apaf-1/cytochrome c complex, and carried out a further analysis by means of molecular dynamics simulations. We ended up with a single model structure where all the lysine residues of cytochrome c that are known as functionally-relevant were involved in forming salt bridges with acidic residues of Apaf-1. This model has revealed three distinctive bifurcated salt bridges, each involving a single lysine residue of cytochrome c and two neighboring acidic resides of Apaf-1. Salt bridge-forming amino acids of Apaf-1 showed a clear evolutionary pattern within Metazoa, with pairs of acidic residues of Apaf-1, involved in bifurcated salt bridges, reaching their highest numbers in the sequences of vertebrates, in which the cytochrome c-mediated mechanism of apoptosome formation seems to be typical.

CONCLUSIONS

The reported model of an Apaf-1/cytochrome c complex provides insights in the nature of protein-protein interactions which are hard to observe in crystallographic or electron microscopy studies. Bifurcated salt bridges can be expected to be stronger than simple salt bridges, and their formation might promote the conformational change of Apaf-1, leading to the formation of an apoptosome. Combination of structural and sequence analyses provides hints on the evolution of the cytochrome c-mediated apoptosis.

Modeling of interaction between cytochrome c and the WD domains of Apaf-1: bifurcated salt bridges underlying apoptosome assembly

Background

Binding of cytochrome c, released from the damaged mitochondria, to the apoptotic protease activating factor 1 (Apaf-1) is a key event in the apoptotic signaling cascade. The binding triggers a major domain rearrangement in Apaf-1, which leads to oligomerization of Apaf-1/cytochrome c complexes into an apoptosome. Despite the availability of crystal structures of cytochrome c and Apaf-1 and cryo-electron microscopy models of the entire apoptosome, the binding mode of cytochrome c to Apaf-1, as well as the nature of the amino acid residues of Apaf-1 involved remain obscure.

Results

We investigated the interaction between cytochrome c and Apaf-1 by combining several modeling approaches. We have applied protein-protein docking and energy minimization, evaluated the resulting models of the Apaf-1/cytochrome c complex, and carried out a further analysis by means of molecular dynamics simulations. We ended up with a single model structure where all the lysine residues of cytochrome c that are known as functionally-relevant were involved in forming salt bridges with acidic residues of Apaf-1. This model has revealed three distinctive bifurcated salt bridges, each involving a single lysine residue of cytochrome c and two neighboring acidic resides of Apaf-1. Salt bridge-forming amino acids of Apaf-1 showed a clear evolutionary pattern within Metazoa, with pairs of acidic residues of Apaf-1, involved in bifurcated salt bridges, reaching their highest numbers in the sequences of vertebrates, in which the cytochrome c-mediated mechanism of apoptosome formation seems to be typical.

Conclusions

The reported model of an Apaf-1/cytochrome c complex provides insights in the nature of protein-protein interactions which are hard to observe in crystallographic or electron microscopy studies. Bifurcated salt bridges can be expected to be stronger than simple salt bridges, and their formation might promote the conformational change of Apaf-1, leading to the formation of an apoptosome. Combination of structural and sequence analyses provides hints on the evolution of the cytochrome c-mediated apoptosis.

Reviewers

This article was reviewed by Andrei L. Osterman, Narayanaswamy Srinivasan, Igor N. Berezovsky, and Gerrit Vriend (nominated by Martijn Huynen).

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0059-4) contains supplementary material, which is available to authorized users.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

Divergence and convergence in enzyme evolution

Comparative analysis of the sequences of enzymes encoded in a variety of prokaryotic and eukaryotic genomes reveals convergence and divergence at several levels. Functional convergence can be inferred when structurally distinct and hence non-homologous enzymes show the ability to catalyze the same biochemical reaction. In contrast, as a result of functional diversification, many structurally similar enzyme molecules act on substantially distinct substrates and catalyze diverse biochemical reactions. Here, we present updates on the ATP-grasp, alkaline phosphatase, cupin, HD hydrolase, and N-terminal nucleophile (Ntn) hydrolase enzyme superfamilies and discuss the patterns of sequence and structural conservation and diversity within these superfamilies. Typically, enzymes within a superfamily possess common sequence motifs and key active site residues, as well as (predicted) reaction mechanisms. These observations suggest that the strained conformation (the entatic state) of the active site, which is responsible for the substrate binding and formation of the transition complex, tends to be conserved within enzyme superfamilies. The subsequent fate of the transition complex is not necessarily conserved and depends on the details of the structures of the enzyme and the substrate. This variability of reaction outcomes limits the ability of sequence analysis to predict the exact enzymatic activities of newly sequenced gene products. Nevertheless, sequence-based (super)family assignments and generic functional predictions, even if imprecise, provide valuable leads for experimental studies and remain the best approach to the functional annotation of uncharacterized proteins from new genomes.

Structural origins of high apparent dielectric constants experienced by ionizable groups in the hydrophobic core of a protein

The side chains of Lys66, Asp66, and Glu66 in staphylococcal nuclease are fully buried and surrounded mainly by hydrophobic matter, except for internal water molecules associated with carboxylic oxygen atoms. These ionizable side chains titrate with pK(a) values of 5.7, 8.8, and 8.9, respectively. To reproduce these pK(a) values with continuum electrostatics calculations, we treated the protein with high dielectric constants. We have examined the structural origins of these high apparent dielectric constants by using NMR spectroscopy to characterize the structural response to the ionization of these internal side chains. Substitution of Val66 with Lys66 and Asp66 led to increased conformational fluctuations of the microenvironments surrounding these groups, even under pH conditions where Lys66 and Asp66 are neutral. When Lys66, Asp66, and Glu66 are charged, the proteins remain almost fully folded, but resonances for a few backbone amides adjacent to the internal ionizable residues are broadened. This suggests that the ionization of the internal groups promotes a local increase in dynamics on the intermediate timescale, consistent with either partial unfolding or increased backbone fluctuations of helix 1 near residue 66, or, less likely, with increased fluctuations of the charged side chains at position 66. These experiments confirm that the high apparent dielectric constants reported by internal Lys66, Asp66, and Glu66 reflect localized changes in conformational fluctuations without incurring detectable global structural reorganization. To improve structure-based pK(a) calculations in proteins, we will need to learn how to treat this coupling between ionization of internal groups and local changes in conformational fluctuations explicitly.

On the origin of life in the zinc world: 1. Photosynthesizing, porous edifices built of hydrothermally precipitated zinc sulfide as cradles of life on Earth

BACKGROUND

The complexity of the problem of the origin of life has spawned a large number of possible evolutionary scenarios. Their number, however, can be dramatically reduced by the simultaneous consideration of various bioenergetic, physical, and geological constraints.

RESULTS

This work puts forward an evolutionary scenario that satisfies the known constraints by proposing that life on Earth emerged, powered by UV-rich solar radiation, at photosynthetically active porous edifices made of precipitated zinc sulfide (ZnS) similar to those found around modern deep-sea hydrothermal vents. Under the high pressure of the primeval, carbon dioxide-dominated atmosphere ZnS could precipitate at the surface of the first continents, within reach of solar light. It is suggested that the ZnS surfaces (1) used the solar radiation to drive carbon dioxide reduction, yielding the building blocks for the first biopolymers, (2) served as templates for the synthesis of longer biopolymers from simpler building blocks, and (3) prevented the first biopolymers from photo-dissociation, by absorbing from them the excess radiation. In addition, the UV light may have favoured the selective enrichment of photostable, RNA-like polymers. Falsification tests of this hypothesis are described in the accompanying article (A.Y. Mulkidjanian, M.Y. Galperin, Biology Direct 2009, 4:27).

CONCLUSION

The suggested "Zn world" scenario identifies the geological conditions under which photosynthesizing ZnS edifices of hydrothermal origin could emerge and persist on primordial Earth, includes a mechanism of the transient storage and utilization of solar light for the production of diverse organic compounds, and identifies the driving forces and selective factors that could have promoted the transition from the first simple, photostable polymers to more complex living organisms.

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.

An extended dead-end elimination algorithm to determine gap-free lists of low energy states

Proteins are flexible systems and commonly populate several functionally important states. To understand protein function, these states and their energies have to be identified. We introduce an algorithm that allows the determination of a gap-free list of the low energy states. This algorithm is based on the dead-end elimination (DEE) theorem and is termed X-DEE (extended DEE). X-DEE is applicable to discrete systems whose state energy can be formulated as pairwise interaction between sites and their intrinsic energies. In this article, the computational performance of X-DEE is analyzed and discussed. X-DEE is implemented to determine the lowest energy protonation states of proteins, a problem to which DEE has not been applied so far. We use X-DEE to calculate a list of low energy protonation states for two bacteriorhodopsin structures that represent the first proton transfer step of the bacteriorhodopsin photocycle.

Chance and design—Proton transfer in water, channels and bioenergetic proteins

Proton transfer and transport in water, gramicidin and some selected channels and bioenergetic proteins are reviewed. An attempt is made to draw some conclusions about how Nature designs long distance, proton transport functionality. The prevalence of water rather than amino acid hydrogen bonded chains is noted, and the possible benefits of waters as the major component are discussed qualitatively.

Proton in the well and through the desolvation barrier

The concept of the membrane proton well was suggested by Peter Mitchell to account for the energetic equivalence of the chemical (DeltapH) and electrical (Deltapsi) components of the proton-motive force. The proton well was defined as a proton-conducting crevice passing down into the membrane dielectric and able to accumulate protons in response to the generation either of Deltapsi or of DeltapH. In this review, the concept of proton well is contrasted to the desolvation penalty of > 500 meV for transferring protons into the membrane core. The magnitude of the desolvation penalty argues against deep proton wells in the energy-transducing enzymes. The shallow DeltapH- and Deltapsi-sensitive proton traps, mechanistically linked to the functional groups in the membrane interior, seem more realistic. In such constructs, the draw of a trapped proton into the membrane core can happen at the expense of some exergonic reaction, e.g., release of another proton from the membrane into the aqueous phase. It is argued that the proton transfer in the ATP synthase and the cytochrome bc complex could proceed in this way.

Ubiquinone reduction in the photosynthetic reaction centre of Rhodobacter sphaeroides: interplay between electron transfer, proton binding and flips of the quinone ring

This review is focused on reactions that gate (control) the electron transfer between the primary quinone QA and secondary quinone QB in the photosynthetic reaction centre of Rhodobacter sphaeroides. The results on electron and proton transfer are discussed in relation to structural information and to the steered molecular dynamics simulations of the QB ring flip in its binding pocket. Depending on the initial position of QB in the pocket and on certain conditions, the rate of electron transfer is suggested to be limited either by the quinone ring flip or by the charge-compensating proton equilibration between the surface and the buried QB site.

Ubiquinol oxidation in the cytochrome bc1 complex: Reaction mechanism and prevention of short-circuiting

This review is focused on the mechanism of ubiquinol oxidation by the cytochrome bc1 complex (bc1). This integral membrane complex serves as a "hub" in the vast majority of electron transfer chains. The bc1 oxidizes a ubiquinol molecule to ubiquinone by a unique "bifurcated" reaction where the two released electrons go to different acceptors: one is accepted by the mobile redox active domain of the [2Fe-2S] iron-sulfur Rieske protein (FeS protein) and the other goes to cytochrome b. The nature of intermediates in this reaction remains unclear. It is also debatable how the enzyme prevents short-circuiting that could happen if both electrons escape to the FeS protein. Here, I consider a reaction mechanism that (i) agrees with the available experimental data, (ii) entails three traits preventing the short-circuiting in bc1, and (iii) exploits the evident structural similarity of the ubiquinone binding sites in the bc1 and the bacterial photosynthetic reaction center (RC). Based on the latter congruence, it is suggested that the reaction route of ubiquinol oxidation by bc1 is a reversal of that leading to the ubiquinol formation in the RC. The rate-limiting step of ubiquinol oxidation is then the re-location of a ubiquinol molecule from its stand-by site within cytochrome b into a catalytic site, which is formed only transiently, after docking of the mobile redox domain of the FeS protein to cytochrome b. In the catalytic site, the quinone ring is stabilized by Glu-272 of cytochrome b and His-161 of the FeS protein. The short circuiting is prevented as long as: (i) the formed semiquinone anion remains bound to the reduced FeS domain and impedes its undocking, so that the second electron is forced to go to cytochrome b; (ii) even after ubiquinol is fully oxidized, the reduced FeS domain remains docked to cytochrome b until electron(s) pass through cytochrome b; (iii) if cytochrome b becomes (over)reduced, the binding and oxidation of further ubiquinol molecules is hampered; the reason is that the Glu-272 residue is turned towards the reduced hemes of cytochrome b and is protonated to stabilize the surplus negative charge; in this state, this residue cannot participate in the binding/stabilization of a ubiquinol molecule.

Specific Radiation Damage Illustrates Light-Induced Structural Changes in the Photosynthetic Reaction Center

The photosynthetic reaction center of the purple non-sulfur bacterium Blastochloris viridis was frozen in the presence and absence of illumination. Differences in the resulting datasets are monitored using the difference Fourier method. Radiation damage is localized to those parts of the protein that are significant for electron transfer, and show changes that are sensitive to oxidation and protonation state.

Elemental sulfur: toxicity in vivo and in vitro to bacterial luciferase, in vitro yeast alcohol dehydrogenase, and bovine liver catalase

The aim of this research was to analyze the effects and the modes of action of elemental sulfur (S(0)) in bioluminescence and respiration of Vibrio fischeri cells and the enzymes crude luciferase, pure catalase, and alcohol dehydrogenase (ADH). Metallic copper removed sulfur and reduced the toxicity of acetone extracts of sediment samples analyzed in the bioluminescence test. The sulfur inhibition of cell bioluminescence was noncompetitive with decanal, the luciferase substrate; reversible, with maximum toxicity after 15 min (EC(50) = 11.8 microg/L); and almost totally recovered after 2 h. In vitro preincubation of crude luciferase extract with sulfur (0.28 ppm) weakly inhibited bioluminescence at 5 min, but at 30 min the inhibition reached 60%. Increasing the concentration of sulfur in the parts per million concentration range in vitro decreased bioluminescence, which was not constant, but depended on exposure time, and no dead-end/total inhibition was observed. The redox state of enzymes in the in vitro system significantly affected inhibition. Hydrogen peroxide restored fully and the reducing agent dithiothreitol, itself toxic, restored only partially luciferase activity in the presence of sulfur. Sulfur (5.5 ppm) slightly inhibited ADH and catalase, and dithiothreitol enhanced sulfur inhibition. High sulfur concentrations (2.2 ppm) inhibited the bioluminescence and enhanced the respiration rate of V. fischeri cells. Elemental sulfur data were interpreted to show that sulfur acted on at least a few V. fischeri cell sites: reversibly modifying luciferase at sites sensitive to/protected by oxidative and reducing agents and by affecting electron transport processes, resulting in enhanced oxygen consumption. Sulfur together with an enzyme reducing agent inhibited the oxidoreductive enzymes ADH and catalase, which have --SH groups, metal ion cofactors, or heme, respectively, in their active centers.

Proton release upon oxidation of tyrosine in reaction centers from Rhodobacter sphaeroides

Markedly different light-induced protonational changes were measured in two reaction center mutants of Rhodobacter sphaeroides. A quadruple mutant containing alterations, at residues L131, M160, M197, and M210, that elevate the midpoint potential of the bacteriochlorophyll dimer was compared to the Y(M) mutant, which contains these alterations plus a tyrosine at M164 serving as a secondary electron donor [Kálmán et al., Nature 402 (1999) 696]. In the quadruple mutant, a proton uptake of 0.1-0.3 H(+)/reaction center between pH 6 and 10 resulted from formation of the oxidized bacteriochlorophyll donor and reduced primary quinone. In the Y(M) mutant, a maximal proton release of -0.5 H(+)/reaction center at pH 8 was attributed to formation of the tyrosyl radical and modeled using electrostatic and direct proton-releasing mechanisms.

The position of QB in the photosynthetic reaction center depends on pH: a theoretical analysis of the proton uptake upon QB reduction

Electrostatics-based calculations have been performed to examine the proton uptake upon reduction of the terminal electron acceptor Q(B) in the photosynthetic reaction center of Rhodobacter sphaeroides as a function of pH and the associated conformational equilibrium. Two crystal structures of the reaction center were considered: one structure was determined in the dark and the other under illumination. In the two structures, the Q(B) was found in two different positions, proximal or distal to the nonheme iron. Because Q(B) was found mainly in the distal position in the dark and only in the proximal position under illumination, the two positions have been attributed mostly to the oxidized and the reduced forms of Q(B), respectively. We calculated the proton uptake upon Q(B) reduction by four different models. In the first model, Q(B) is allowed to equilibrate between the two positions with either oxidation state. This equilibrium was allowed to vary with pH. In the other three models the distribution of Q(B) between the proximal position and the distal position was pH-independent, with Q(B) occupying only the distal position or only the proximal position or populating the two positions with a fixed ratio. Only the first model, which includes the pH-dependent conformational equilibrium, reproduces both the experimentally measured pH dependence of the proton uptake and the crystallographically observed conformational equilibrium at pH 8. From this model, we find that Q(B) occupies only the distal position below pH 6.5 and only the proximal position above pH 9.0 in both oxidation states. Between these pH values both positions are partially occupied. The reduced Q(B) has a higher occupancy in the proximal position than the oxidized Q(B). In summary, the present results indicate that the conformational equilibrium of Q(B) depends not only on the redox state of Q(B), but also on the pH value of the solution.

Is the glutamate residue Glu-373 the proton acceptor of the excitatory amino acid carrier 1?

Glutamate transport by the neuronal excitatory amino acid carrier (EAAC1) is accompanied by the coupled movement of one proton across the membrane. We have demonstrated previously that the cotransported proton binds to the carrier in the absence of glutamate and, thus, modulates the EAAC1 affinity for glutamate. Here, we used site-directed mutagenesis together with a rapid kinetic technique that allows one to generate sub-millisecond glutamate concentration jumps to locate possible binding sites of the glutamate transporter for the cotransported proton. One candidate for this binding site, the highly conserved glutamic acid residue Glu-373 of EAAC1, was mutated to glutamine. Our results demonstrate that the mutant transporter does not catalyze net transport of glutamate, whereas Na(+)/glutamate homoexchange is unimpaired. Furthermore, the voltage dependence of the rates of Na(+) binding and glutamate translocation are unchanged compared with the wild-type. In contrast to the wild-type, however, homoexchange of the E373Q transporter is completely pH-independent. In line with these findings the transport kinetics of the mutant EAAC1 show no deuterium isotope effect. Thus, we suggest a new transport mechanism, in which Glu-373 forms part of the binding site of EAAC1 for the cotransported proton. In this model, protonation of Glu-373 is required for Na(+)/glutamate translocation, whereas the relocation of the carrier is only possible when Glu-373 is negatively charged. Interestingly, the Glu-373-homologous amino acid residue is glutamine in the related neutral amino acid transporter alanine-serine-cysteine transporter. The function of alanine-serine-cysteine transporter is neither potassium- nor proton-dependent. Consequently, our results emphasize the general importance of glutamate and aspartate residues for proton transport across membranes.

Electrostatics and proton transfer in photosynthetic water oxidation

Photosystem II (PSII) oxidizes two water molecules to yield dioxygen plus four protons. Dioxygen is released during the last out of four sequential oxidation steps of the catalytic centre (S(0) --> S(1), S(1) --> S(2), S(2) --> S(3), S(3) --> S(4) --> S(0)). The release of the chemically produced protons is blurred by transient, highly variable and electrostatically triggered proton transfer at the periphery (Bohr effect). The extent of the latter transiently amounts to more than one H(+)/e(-) under certain conditions and this is understood in terms of electrostatics. By kinetic analyses of electron-proton transfer and electrochromism, we discriminated between Bohr-effect and chemically produced protons and arrived at a distribution of the latter over the oxidation steps of 1 : 0 : 1 : 2. During the oxidation of tyr-161 on subunit D1 (Y(Z)), its phenolic proton is not normally released into the bulk. Instead, it is shared with and confined in a hydrogen-bonded cluster. This notion is difficult to reconcile with proposed mechanisms where Y(Z) acts as a hydrogen acceptor for bound water. Only in manganese (Mn) depleted PSII is the proton released into the bulk and this changes the rate of electron transfer between Y(Z) and the primary donor of PSII P(+)(680) from electron to proton controlled. D1-His190, the proposed centre of the hydrogen-bonded cluster around Y(Z), is probably further remote from Y(Z) than previously thought, because substitution of D1-Glu189, its direct neighbour, by Gln, Arg or Lys is without effect on the electron transfer from Y(Z) to P(+)(680) (in nanoseconds) and from the Mn cluster to Y(ox)(Z).

Key role of proline L209 in connecting the distant quinone pockets in the reaction center of Rhodobacter sphaeroides

Photosynthetic bacterial reaction centers convert light excitation into chemical free energy. The initial electron transfer leads to the consecutive semireductions of the primary (Q(A)) and secondary (Q(B)) quinone acceptors. The Q(A)(-) and Q(B)(-) formations induce proton uptake from the bulk. Their magnitudes (H(+)/Q(A)(-) and H(+)/Q(B)(-), respectively) probe the electrostatic interactions within the complex. The pH dependence of H(+)/Q(A)(-) and H(+)/Q(B)(-) were studied in five single mutants modified at the L209 site (L209P-->F,Y,W,E,T). This residue is situated at the border of a continuous chain of water molecules connecting Q(B) to the bulk. In the wild type (WT), a proton uptake band is present at high pH in the H(+)/Q(A)(-) and H(+)/Q(B)(-) curves and is commonly attributed to a cluster of acidic groups situated nearby Q(B). In the H(+)/Q(A)(-) curves of the L209 variants, this band is systematically absent but remains in the H(+)/Q(B)(-) curves. Moreover, notable increase of H(+)/Q(B)(-) is observed in the L209 mutants at neutral pH as compared with the WT. The large effects observed in all L209 mutants are not associated with significant structural changes (Kuglstatter, A., Ermler, U., Michel, H., Baciou, L. & Fritzsch, G. Biochemistry (2001) 40, 4253-4260). Our data suggest that, in the L209 mutants, the Q(B) cluster does not respond to the Q(A)(-) formation as observed in the WT. We propose that, in the mutants, removal of the rigid proline L209 breaks a necessary hydrogen bonding connection between the quinone sites. These findings suggest an important role for structural rigidity in ensuring a functional interaction between quinone binding sites.

Identification of the proton pathway in bacterial reaction centers: decrease of proton transfer rate by mutation of surface histidines at H126 and H128 and chemical rescue by imidazole identifies the initial proton donors

The pathway for proton transfer to Q(B) was studied in the reaction center (RC) from Rhodobacter sphaeroides. The binding of Zn(2+) or Cd(2+) to the RC surface at His-H126, His-H128, and Asp-H124 inhibits the rate of proton transfer to Q(B), suggesting that the His may be important for proton transfer [Paddock, M. L., Graige, M. S., Feher, G. and Okamura, M. Y. (1999) Proc. Natl. Acad. Sci. U.S.A. 96, 6183-6188]. To assess directly the role of the histidines, mutant RCs were constructed in which either one or both His were replaced with Ala. In the single His mutant RCs, no significant effects were observed. In contrast, in the double mutant RC at pH 8.5, the observed rates of proton uptake associated with both the first and the second proton-coupled electron-transfer reactions k(AB)(()(1)()) [Q(A)(-)(*)Q(B)-Glu(-) + H(+) --> Q(A)(-)(*)Q(B)-GluH --> Q(A)Q(B)(-)(*)-GluH] and k(AB)(()(2)()) [Q(A)(-)(*)Q(B)(-)(*) + H(+) --> Q(A)(-)(*)(Q(B)H)(*) --> Q(A)(Q(B)H)(-)], were found to be slowed by factors of approximately 10 and approximately 4, respectively. Evidence that the observed changes in the double mutant RC are due to a reduction in the proton-transfer rate constants are provided by the observations: (i) k(AB)(1) at pH approximately pK(a) of GluH became biphasic, indicating that proton transfer is slower than electron transfer and (ii) k(AB)(2) became independent of the driving force for electron transfer, indicating that proton transfer is the rate-limiting step. These changes were overcome by the addition of exogenous imidazole which acts as a proton donor in place of the imidazole groups of His that were removed in the double mutant RC. Thus, we conclude that His-H126 and His-H128 facilitate proton transfer into the RC, acting as RC-bound proton donors at the entrance of the proton-transfer pathways.

Photosystem II of peas: effects of added divalent cations of Mn, Fe, Mg, and Ca on two kinetic components of P(+)(680) reduction in Mn-depleted core particles

The catalytic Mn cluster of the photosynthetic oxygen-evolving system is oxidized via a tyrosine, Y(Z), by a photooxidized chlorophyll a moiety, P(+)(680). The rapid reduction of P(+)(680) by Y(Z) in nanoseconds requires the intactness of an acid/base cluster around Y(Z) with an apparent functional pK of <5. The removal of Mn (together with bound Ca) shifts the pK of the acid/base cluster from the acid into the neutral pH range. At alkaline pH the electron transfer (ET) from Y(Z) to P(+)(680) is still rapid (<1 micros), whereas at acid pH the ET is much slower (10-100 micros) and steered by proton release. In the intermediate pH domain one observes a mix of these kinetic components (see R. Ahlbrink, M. Haumann, D. Cherepanov, O. Bögershausen, A. Mulkidjanian, W. Junge, Biochemistry 37 (1998)). The overall kinetics of P(680)(+) reduction by Y(Z) in Mn-depleted photosystem II (PS II) has been previously shown to be slowed down by divalent cations (added at >10 microM), namely: Mn(2+), Co(2+), Ni(2+), Cu(2+), Zn(2+) (C.W. Hoganson, P.A. Casey, O. Hansson, Biochim. Biophys. Acta 1057 (1991)). Using Mn-depleted PS II core particles from pea as starting material, we re-investigated this phenomenon at nanosecond resolution, aiming at the effect of divalent cations on the particular kinetic components of P(+)(680) reduction. To our surprise we found only the slower, proton steered component retarded by some added cations (namely Co(2+)/Zn(2+)>Fe(2+)>Mn(2+)). Neither the fast component nor the apparent pK of the acid/base cluster around Y(Z) was affected. Apparently, the divalent cations acted (electrostatically) on the proton release channel that connects the oxygen-evolving complex with the bulk water, but not on the ET between Y(Z) and P(+)(680), proper. Contrastingly, Ca(2+) and Mg(2+), when added at >5 mM, accelerated the slow component of P(+)(680) reduction by Y(Z) and shifted the apparent pK of Y(Z) from 7.4 to 6.6 and 6.7, respectively. It was evident that the binding site(s) for added Ca(2+) and Mg(2+) were close to Y(Z) proper. The data obtained are discussed in relation to the nature of the metal-binding sites in photosystem II.

Proton transfer in Azotobacter vinelandii ferredoxin I: entatic Lys84 operates as elastic counterbalance for the proton-carrying Asp15

In ferredoxin I from Azotobacter vinelandii, the reduction of a [3Fe-4S] iron-sulphur cluster is coupled with the protonation of the mu2S sulphur atom that is approx. 6 A away from the protein boundary. The recent study of the site-specific mutants of ferredoxin I led to the conclusion that a particular surface aspartic residue (Asp15) is solely responsible for the proton transfer to the mu2S atom by 'rapid penetrative excursions' (K. Chen, J. Hirst, R. Camba, C.A. Bonagura, C.D. Stout, B.K. Burgess, F.A. Armstrong, Nature 405 (2000) 814-817). In the same paper it has been reported that the replacement of Asp15 by glutamate led to the blockage of the enzyme, although glutamate, with its longer and more flexible side chain, should apparently do even better as a mobile proton carrier than aspartate. We tackled this puzzling incompetence of Glu15 by molecular dynamics simulations. It was revealed that the conformational alterations of Asp15 are energetically balanced by the straining of the nearby Lys84 side chain in wild-type ferredoxin I but not in the Asp15-->Glu mutant. Lys84 in ferredoxin I of A. vinelandii seems to represent the first case where the strained (entatic) conformation of a particular amino acid side chain could be directly identified in the ground state of an enzyme and assigned to a distinct mechanism of energy balance during the catalytic transition.

Photosynthetic electron transfer controlled by protein relaxation: analysis by Langevin stochastic approach

Relaxation processes in proteins range in time from picoseconds to seconds. Correspondingly, biological electron transfer (ET) could be controlled by slow protein relaxation. We used the Langevin stochastic approach to describe this type of ET dynamics. Two different types of kinetic behavior were revealed, namely: oscillating ET (that could occur at picoseconds) and monotonically relaxing ET. On a longer time scale, the ET dynamics can include two different kinetic components. The faster one reflects the initial, nonadiabatic ET, whereas the slower one is governed by the medium relaxation. We derived a simple relation between the relative extents of these components, the change in the free energy (DeltaG), and the energy of the slow reorganization Lambda. The rate of ET was found to be determined by slow relaxation at -DeltaG < or = Lambda. The application of the developed approach to experimental data on ET in the bacterial photosynthetic reaction centers allowed a quantitative description of the oscillating features in the primary charge separation and yielded values of Lambda for the slower low-exothermic ET reactions. In all cases but one, the obtained estimates of Lambda varied in the range of 70-100 meV. Because the vast majority of the biological ET reactions are only slightly exothermic (DeltaG > or = -100 meV), the relaxationally controlled ET is likely to prevail in proteins.

Equilibrium and kinetic parameters for the binding of inhibitors to the QB pocket in bacterial chromatophores: dependence on the state of QA

The equilibrium and kinetic parameters for the binding of various inhibitors to the Q(B) pocket of the bacterial reaction center were investigated in chromatophores from Rhodobacter capsulatus and Rhodobacter sphaeroides. By monitoring the near-IR absorption changes specific to Q(A)(-) and Q(B)(-), we measured the fraction of inhibited centers in the dark and the kinetics and extent of inhibitor displacement after one flash due to the formation of the Q(A)Q(B)(-) state. The inhibitor release rate was much faster for triazines and o-phenanthroline (t(1/2) in the 50 ms to 1 s range) than for stigmatellin (t(1/2) approximately 20 s). For inhibitors with a rapid release rate, the fast phase of P(+) decay observed in the absence of secondary donor reflects the competition between P(+)Q(A)(-) recombination and inhibitor release: it is thus faster than the P(+)Q(A)(-) recombination, and its relative extent is smaller than the fraction of initially inhibited centers. At appropriate inhibitor concentrations, one can have almost total binding in the dark and almost total inhibitor displacement after one flash. Under such conditions, a pair of closely spaced flashes resets the two-electron gate in a single state (Q(A)Q(B)(-)), irrespective of the initial state. The apparent dissociation constant of terbutryn was significantly increased (by a factor of 4-7) in the presence of Q(A)(-), in agreement with the conclusion of Wraight and co-workers [Stein, R. R., et al. (1984) J. Cell. Biochem. 24, 243-259]. We suggest that this effect is essentially due to a tighter binding of ubiquinone in the Q(A)(-) state.

Retardation of proton transfer caused by binding of the transition metal ion to the bacterial reaction center is due to pKa shifts of key protonatable residues

Transition metal ions bind to the reaction center (RC) protein of the photosynthetic bacterium Rhodobacter sphaeroides and slow the light-induced electron and proton transfer to the secondary quinone, Q(B). We studied the properties of the metal ion-RC complex by measuring the pH dependence of the dissociation constant and the stoichiometry of proton release upon ligand formation. We investigated the mechanism of inhibition by measuring the stoichiometry and kinetics of flash-induced proton binding, the transfer of (first and second) electrons to Q(B), and the rate of steady-state turnover of the RC in the absence and presence of Cd(2+) and Ni(2+) on a wide pH range. The following results were obtained. (1) The complexation of transition metal ions Cd(2+) and Ni(2+) with the bacterial RC showed strong pH dependence. This observation was explained by different (pH-dependent) states of the metal-ligand cluster: the complex formation was strong when the ligand (Asp and His residues) was deprotonated and was much weaker if the ligand was partly (or fully) protonated. A direct consequence of the model was the pH-dependent proton release upon complexation. (2) The retardation of transfer of electrons and protons to Q(B) was also strongly pH-dependent. The effect was large in the neutral pH range and decreased toward the acidic and alkaline pH values. (3) Steady-state turnover measurements indicated that the rate of the second proton transfer was much less inhibited than that of the first one, which became the rate-limiting step in continuous turnover of the RC. (4) Sodium azide partly recovered the proton transfer rate. The effect is not due to removal of the bound metal ion by azide but probably by formation of a proton-transporting azide network similarly as water molecules may build up proton pathways. (5) We argue that the inhibition comes mainly from pK(a) shifts of key protonatable residues that control the proton transfer along the H-bond network to Q(B). The electrostatic interaction between the metal ion and these residues may result in acidic pK(a) shifts between 1.5 and 2.0 that account for the observed retardation of the electron and proton transfer.

Effect of metal binding on electrogenic proton transfer associated with reduction of the secondary electron acceptor (QB) in Rhodobacter sphaeroides chromatophores

The influence of metal ion (Cd(2+), Zn(2+), Ni(2+)) binding on the electrogenic phases of proton transfer connected with reduction of quinone Q(B) in chromatophores from Rhodobacter sphaeroides was studied by time-resolved electric potential changes. In the presence of metals, the electrogenic transients associated with proton transfer on first and second flash at pH 8 were found to be slower by factors of 3-6. This is essentially the same effect of metal binding that was observed on optical transients in isolated reaction centers (RC), where the metal ion was shown to inhibit proton transfer [Paddock, M. L., Graige, M. S., Feher, G., and Okamura, M. Y. (1999) Proc. Natl. Acad. Sci. U.S.A. 96, 6183-6188]. The effect of metal binding on the kinetics in chromatophores is, therefore, similarly attributed to inhibition of proton uptake, which becomes rate-limiting. A striking observation was an increase in the amplitude of the electrogenic proton-uptake phase after the first flash with bound metal ion. We attribute this to a loss of internal proton rearrangement, requiring that the protons that stabilize Q(B)(-) come from solution. In mutant RCs, in which His-H126 and His-H128 are replaced with Ala, the apparent binding of Cd(2+) and Ni(2+) was decreased, showing that the binding site of these metal ions is the same as found in RC crystals [Axelrod, H. L., Abresch, E. C., Paddock, M. L., Okamura, M. Y., and Feher, G. (2000) Proc. Natl. Acad. Sci. U.S.A. 97, 1542-1547]. Therefore, the unique proton entry point near His-H126, His-H128, and Asp-M17 that was identified in isolated RCs is also the entry point in chromatophores.

Photosynthetic water oxidation: towards a mechanism

This mini-review outlines the current theories on the mechanism of electron transfer from water to P680, the location and structure of the water oxidising complex and the role of the manganese cluster. We discuss how our data fit in with current theories and put forward our ideas on the location and mechanism of water oxidation.

Fourier transform infrared evidence of proton uptake by glutamate L212 upon reduction of the secondary quinone QB in the photosynthetic reaction center from Rhodobacter capsulatus

The photoreduction of the secondary quinone Q(B) in native reaction centers (RCs) of Rhodobacter capsulatus and in RCs from the GluL212 --> Gln and GluL212 --> Ala mutants has been investigated at pH 7 in (1)H(2)O and (2)H(2)O by light-induced Fourier transform infrared (FTIR) difference spectroscopy. The Q(B)(-)/Q(B) FTIR difference spectra reflect changes of quinone-protein interactions and of protonation state of carboxylic acid groups as well as reorganization of the protein upon electron transfer. Comparison of Q(B)(-)/Q(B) spectra of native and mutant RCs indicates that the interactions between Q(B) or Q(B)(-) and the protein are similar in all RCs. A differential signal at approximately 1650/1640 cm(-1), which is common to all the spectra, is associated with a movement of a peptide carbonyl or a side chain following Q(B) reduction. On the other hand, Q(B)(-)/Q(B) spectra of native and mutant RCs display several differences, notably between 1700 and 1650 cm(-1) (amide I and side chains), between 1570 and 1530 cm(-1) (amide II), and at 1728-1730 cm(-1) (protonated carboxylic acid groups). In particular, the latter region in native RCs is characterized by a main positive band at 1728 cm(-1) and a negative signal at 1739 cm(-1). In the L212 mutants, the amplitude of the positive band is strongly decreased leading to a differential signal at 1739/1730 cm(-1) that is insensitive to (1)H/(2)H isotopic exchange. In native RCs, only the 1728 cm(-1) band is affected in (2)H(2)O while the 1739 cm(-1) signal is not. The effects of the mutations and of (1)H/(2)H exchange on the Q(B)(-)/Q(B) spectra concur in the attribution of the 1728 cm(-1) band in native RCs to (partial) proton uptake by GluL212 upon the first electron transfer to Q(B), as previously observed in Rhodobacter sphaeroides RCs [Nabedryk, E., Breton, J., Hienerwadel, R., Fogel, C., Mäntele, W., Paddock, M. L., and Okamura, M. Y. (1995) Biochemistry 34, 14722-14732]. More generally, strong homologies of the Q(B) to Q(B)(-) transition in the RCs from Rb. sphaeroides and Rb. capsulatus are detected by differential FTIR spectroscopy. The FTIR data are discussed in relation with the results from global proton uptake measurements and electrogenic events concomitant with the reduction of Q(B) and with a model of the Q(B) turnover in Rb. sphaeroides RCs [Mulkidjanian, A. Y. (1999) FEBS Lett. 463, 199-204].

Proton-coupled electron transfer at the Q(o) site: what type of mechanism can account for the high activation barrier?

In Rhodobacter sphaeroides, transfer of the first electron in quinol oxidation by the bc(1) complex shows kinetic features (a slow rate (approx. 1.5 x 10(3)/s), high activation energy (approx. 65 kJ/mol) and reorganization energy, lambda (2.5 V)) that are unexpected from Marcus theory and the distances shown by the structures. Reduction of the oxidized iron-sulfur protein occurs after formation of the enzyme-substrate complex, and involves a H-transfer in which the electron transfer occurs through the approx. 7 A of a bridging histidine forming a H-bond with quinol and a ligand to 2Fe-2S. The anomalous kinetic features can be explained by a mechanism in which the electron transfer is constrained by coupled transfer of the proton. We discuss this in the context of mutant strains with modified E(m,7) and pK for the iron-sulfur protein, and Marcus theory for proton-coupled electron transfer. We suggest that transfer of the second proton and electron involve movement of semiquinone in the Q(o) site, and rotation of the Glu of the conserved -PEWY- sequence. Mutational studies show a key role for the domain proximal to heme b(L). The effects of mutation at Tyr-302 (Tyr-279 in bovine sequence) point to a possible linkage between conformational changes in the proximal domain, and changes leading to closure of the iron-sulfur protein access channel at the distal domain.

Reduction and protonation of the secondary quinone acceptor of Rhodobacter sphaeroides photosynthetic reaction center: kinetic model based on a comparison of wild-type chromatophores with mutants carrying Arg-->Ile substitution at sites 207 and 217 in the L-subunit

After the light-induced charge separation in the photosynthetic reaction center (RC) of Rhodobacter sphaeroides, the electron reaches, via the tightly bound ubiquinone QA, the loosely bound ubiquinone Q(B) After two subsequent flashes of light, Q(B) is reduced to ubiquinol Q(B)H2, with a semiquinone anion Q-(B) formed as an intermediate after the first flash. We studied Q(B)H2 formation in chromatophores from Rb. sphaeroides mutants that carried Arg-->Ile substitution at sites 207 and 217 in the L-subunit. While Arg-L207 is 17 A away from Q(B), Arg-L217 is closer (9 A) and contacts the Q(B)-binding pocket. From the pH dependence of the charge recombination in the RC after the first flash, we estimated deltaG(AB), the free energy difference between the Q-(A)Q(B) and Q(A)Q-(B) states, and pK212, the apparent pK of Glu-L212, a residue that is only 4 A away from Q(B). As expected, the replacement of positively charged arginines by neutral isoleucines destabilized the Q-(B) state in the L217RI mutant to a larger extent than in the L207RI one. Also as expected, pK212 increased by approximately 0.4 pH units in the L207RI mutant. The value of pK212 in the L217RI mutant decreased by 0.3 pH units, contrary to expectations. The rate of the Q-(A)Q-(B)-->Q(A)Q(B)H2 transition upon the second flash, as monitored by electrometry via the accompanying changes in the membrane potential, was two times faster in the L207RI mutant than in the wild-type, but remained essentially unchanged in the L217RI mutant. To rationalize these findings, we developed and analyzed a kinetic model of the Q-(A)Q-(B)-->Q(A)Q(B)H2 transition. The model properly described the available experimental data and provided a set of quantitative kinetic and thermodynamic parameters of the Q(B) turnover. The non-electrostatic, 'chemical' affinity of the QB site to protons proved to be as important for the attracting protons from the bulk, as the appropriate electrostatic potential. The mutation-caused changes in the chemical proton affinity could be estimated from the difference between the experimentally established pK2J2 shifts and the expected changes in the electrostatic potential at Glu-L212, calculable from the X-ray structure of the RC. Based on functional studies, structural data and kinetic modeling, we suggest a mechanistic scheme of the QB turnover. The detachment of the formed ubiquinol from its proximal position next to Glu-L212 is considered as the rate-limiting step of the reaction cycle.