Variation of Ser-L223 Hydrogen Bonding with the QB Redox State in Reaction Centers from Rhodobacter sphaeroides

Mechanism of Asparagine-Mediated Proton Transfer in Photosynthetic Reaction Centers

In photosynthetic reaction centers from purple bacteria (PbRCs), light-induced charge separation leads to the reduction of the terminal electron acceptor quinone, QB. The reduction of QB to QB•- is followed by protonation via Asp-L213 and Ser-L223 in PbRC from Rhodobacter sphaeroides. However, Asp-L213 is replaced with nontitratable Asn-L222 and Asn-L213 in PbRCs from Thermochromatium tepidum and Blastochloris viridis, respectively. Here, we investigated the energetics of proton transfer along the asparagine-involved H-bond network using a quantum mechanical/molecular mechanical approach. The potential energy profile for the H-bond between H3O+ and the carbonyl O site of Asn-L222 shows that the proton is predominantly localized at the Asn-L222 moiety in the T. tepidum PbRC protein environment, easily forming the enol species. The release of the proton from the amide -NH2 site toward Ser-L232 via tautomerization suffers from the energy barrier. Upon reorientation of Asn-L222, the enol -OH site forms a short low-barrier H-bond with Ser-L232, facilitating protonation of QB•- in a Grotthuss-like mechanism. This is a basis of how asparagine or glutamine side chains function as acceptors/donors in proton transfer pathways.

Tools for analyzing protonation states and for tracing proton transfer pathways with examples from the Rb. sphaeroides photosynthetic reaction centers

Protons participate in many reactions. In proteins, protons need paths to move in and out of buried active sites. The vectorial movement of protons coupled to electron transfer reactions establishes the transmembrane electrochemical gradient used for many reactions, including ATP synthesis. Protons move through hydrogen bonded chains of waters and hydroxy side chains via the Grotthuss mechanism and by proton binding and release from acidic and basic residues. MCCE analysis shows that proteins exist in a large number of protonation states. Knowledge of the equilibrium ensemble can provide a rational basis for setting protonation states in simulations that fix them, such as molecular dynamics (MD). The proton path into the Q_B site in the bacterial reaction centers (RCs) of Rb. sphaeroides is analyzed by MD to provide an example of the benefits of using protonation states found by the MCCE program. A tangled web of side chains and waters link the cytoplasm to Q_B. MCCE analysis of snapshots from multiple trajectories shows that changing the input protonation state of a residue in MD biases the trajectory shifting the proton affinity of that residue. However, the proton affinity of some residues is more sensitive to the input structure. The proton transfer networks derived from different trajectories are quite robust. There are some changes in connectivity that are largely restricted to the specific residues whose protonation state is changed. Trajectories with Q_B^•- are compared with earlier results obtained with Q_B [Wei et. al Photosynthesis Research volume 152, pages153-165 (2022)] showing only modest changes. While introducing new methods the study highlights the difficulty of establishing the connections between protein conformation.

Electron Transfer Route between Quinones in Type-II Reaction Centers

In photosynthetic reaction centers from purple bacteria (PbRCs) and photosystem II (PSII), the photoinduced charge separation is terminated by an electron transfer between the primary (QA) and secondary (QB) quinones. Here, we investigate the electron transfer route, calculating the superexchange coupling (HQA-QB) for electron transfer from QA to QB in the protein environment. HQA-QB is significantly larger in PbRC than in PSII. In superexchange electron tunneling, the electron transfer via unoccupied molecular orbitals of the nonheme Fe complex (QA → Fe → QB) is pronounced in PbRC, whereas the electron transfer via occupied molecular orbitals (Fe → QB followed by QA → Fe) is pronounced in PSII. The significantly large HQA-QB is caused by a water molecule that donates the H-bond to the ligand Glu-M234 in PbRC. The corresponding water molecule is absent in PSII due to the existence of D1-Tyr246. HQA-QB increases in response to the Ser-L223···QB H-bond formation caused by an extension of the H-bond network, which facilitates charge delocalization over the QB site. This explains the observed discrepancy in the QA-to-QB electron transfer between PbRC and PSII, despite their structural similarity.

Conformational Changes and H-Bond Rearrangements during Quinone Release in Photosystem II

In photosystem II (PSII) and photosynthetic reaction centers from purple bacteria (PbRC), the electron released from the electronically excited chlorophyll is transferred to the terminal electron acceptor quinone, Q_B. Q_B accepts two electrons and two protons before leaving the protein. We investigated the molecular mechanism of quinone exchange in PSII, conducting molecular dynamics (MD) simulations and quantum mechanical/molecular mechanical (QM/MM) calculations. MD simulations suggest that the release of Q_B leads to the transformation of the short helix (D1-Phe260 to D1-Ser264), which is adjacent to the stromal helix de (D1-Asn247 to D1-Ile259), into a loop and to the formation of a water-intake channel. Water molecules enter the Q_B binding pocket via the channel and form an H-bond network. QM/MM calculations indicate that the H-bond network serves as a proton-transfer pathway for the reprotonation of D1-His215, the proton donor during Q_BH^-/Q_BH₂ conversion. Together with the absence of the corresponding short helix but the presence of Glu-L212 in PbRC, it seems likely that the two type-II reaction centers undergo quinone exchange via different mechanisms.

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

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

Mechanism of the formation of proton transfer pathways in photosynthetic reaction centers

The crystal structures of photosynthetic reaction centers from purple bacteria (PbRCs) and photosystem II show large structural similarity. However, the proposed mechanisms of proton transfer toward the terminal electron acceptor quinone (Q_B) are not consistent. In particular, not His-L190, which is an H-bond partner of Q_B, but rather Glu-L212, which is ∼6 Å away from Q_B, was assumed to be the direct proton donor for Q_B. We demonstrate that the H-bond between His-L190 and Q_B is a low-barrier H-bond, which facilitates proton transfer from singly protonated His-L190 to Q_B. Furthermore, Glu-L212 is not a direct H-bond donor for Q_B. However, it facilitates proton transfer toward deprotonated His-L190 via water molecules after Q_BH₂ forms and leaves the PbRC.

In photosynthetic reaction centers from purple bacteria (PbRCs) from Rhodobacter sphaeroides, the secondary quinone Q_B accepts two electrons and two protons via electron-coupled proton transfer (PT). Here, we identify PT pathways that proceed toward the Q_B binding site, using a quantum mechanical/molecular mechanical approach. As the first electron is transferred to Q_B, the formation of the Grotthuss-like pre-PT H-bond network is observed along Asp-L213, Ser-L223, and the distal Q_B carbonyl O site. As the second electron is transferred, the formation of a low-barrier H-bond is observed between His-L190 at Fe and the proximal Q_B carbonyl O site, which facilitates the second PT. As Q_BH₂ leaves PbRC, a chain of water molecules connects protonated Glu-L212 and deprotonated His-L190 forms, which serves as a pathway for the His-L190 reprotonation. The findings of the second pathway, which does not involve Glu-L212, and the third pathway, which proceeds from Glu-L212 to His-L190, provide a mechanism for PT commonly used among PbRCs.

Proton transfer pathway from the oxygen-evolving complex in photosystem II substantiated by extensive mutagenesis

We report a structure-based biological approach to identify the proton-transfer pathway in photosystem II. First, molecular dynamics (MD) simulations were conducted to analyze the H-bond network that may serve as a Grotthuss-like proton conduit. MD simulations show that D1-Asp61, the H-bond acceptor of H₂O at the Mn₄CaO₅ cluster (W1), forms an H-bond via one water molecule with D1-Glu65 but not with D2-Glu312. Then, D1-Asp61, D1-Glu65, D2-Glu312, and the adjacent residues, D1-Arg334, D2-Glu302, and D2-Glu323, were thoroughly mutated to the other 19 residues, i.e., 114 Chlamydomonas chloroplast mutant cells were generated. Mutation of D1-Asp61 was most crucial. Only the D61E and D61C cells grew photoautotrophically and exhibit O₂-evolving activity. Mutations of D2-Glu312 were less crucial to photosynthetic growth than mutations of D1-Glu65. Quantum mechanical/molecular mechanical calculations indicated that in the PSII crystal structure, the proton is predominantly localized at D1-Glu65 along the H-bond with D2-Glu312, i.e., pK_a(D1-Glu65) > pK_a(D2-Glu312). The potential-energy profile shows that the release of the proton from D1-Glu65 leads to the formation of the two short H-bonds between D1-Asp61 and D1-Glu65, which facilitates downhill proton transfer along the Grotthuss-like proton conduit in the S₂ to S₃ transition. It seems possible that D1-Glu65 is involved in the dominant pathway that proceeds from W1 via D1-Asp61 toward the thylakoid lumen, whereas D2-Glu312 and D1-Arg334 may be involved in alternative pathways in some mutants.

Chemical rescue of H+ delivery in proton transfer mutants of reaction center of photosynthetic bacteria

In the native and most mutant reaction centers of bacterial photosynthesis, the electron transfer is coupled to proton transfer and is rate limiting for the second reduction of Q_B^- → Q_BH₂. In the presence of divalent metal ions (e.g. Cd²⁺) or in some ("proton transfer") mutants (L210DN/M17DN or L213DN), the proton delivery to Q_B^- is made rate limiting and the properties of the proton pathway can be directly examined. We found that small weak acids and buffers in large concentrations (up to 1 M) were able to rescue the severely impaired proton transfer capability differently depending on the location of the defects: lesions at the protein surface (proton gate H126H/H128H + Cd²⁺), beneath the surface (M17DN + Cd²⁺, L210DN/M17DN) or deep inside the protein (L213DN) could be completely, partially or to very small extent recovered, respectively. Small zwitterionic acids (azide/hydrazoic acid) and buffers (tricine) proved to be highly effective rescuers consistent with their enhanced binding affinity and access to any of the proton acceptors (including Q_B^- itself) in the pathway. As a consequence, back titration of the protons at L212Glu could be observed as a pH-dependence of the rate constant of the charge recombination in the presence of azide or formate. Model calculations support the collective influence of the acid cluster on the change of the protonation states upon extension of the cluster with the bound small acid. In proton transfer mutants, the rescuing agents decreased the free energy of activation together with their enthalpic and entropic components. This is in agreement with the hypothesis that they function as protein-penetrating protonophores delivering protons into the chain and select dominating paths out of many alternate routes. We estimate that the proton delivery will be accelerated in one pathway out of 100-200 alternate pathways. The implications for design of the chemical recovery of impaired intra-protein proton transfer pathways in proton transfer mutants are discussed.

Redox potentials of ubiquinone, menaquinone, phylloquinone, and plastoquinone in aqueous solution

Quinones serve as redox active cofactors in bacterial photosynthetic reaction centers: photosystem I, photosystem II, cytochrome bc ₁, and cytochrome b ₆ f. In particular, ubiquinone is ubiquitous in animals and most bacteria and plays a key role in several cellular processes, e.g., mitochondrial electron transport. Their experimentally measured redox potential values for one-electron reduction E _m(Q/Q^·-) were already reported in dimethylformamide (DMF) versus saturated calomel electrode but not in water versus normal hydrogen electrode (NHE). We calculated E _m(Q/Q^·-) of 1,4-quinones using a quantum chemical approach. The calculated energy differences of reduction of Q to Q^·- in DMF and water for 1,4-quinone derivatives correlated highly with the experimentally measured E _m(Q/Q^·-) in DMF and water, respectively. E _m(Q/Q^·-) were calculated to be -163 mV for ubiquinone, -260 mV for menaquinone and phylloquinone, and -154 mV for plastoquinone in water versus NHE.

Terminal Electron-Proton Transfer Dynamics in the Quinone Reduction of Respiratory Complex I

Complex I functions as a redox-driven proton pump in aerobic respiratory chains. By reducing quinone (Q), complex I employs the free energy released in the process to thermodynamically drive proton pumping across its membrane domain. The initial Q reduction step plays a central role in activating the proton pumping machinery. In order to probe the energetics, dynamics, and molecular mechanism for the proton-coupled electron transfer process linked to the Q reduction, we employ here multiscale quantum and classical molecular simulations. We identify that both ubiquinone (UQ) and menaquinone (MQ) can form stacking and hydrogen-bonded interactions with the conserved Q-binding-site residue His-38 and that conformational changes between these binding modes modulate the Q redox potentials and the rate of electron transfer (eT) from the terminal N2 iron-sulfur center. We further observe that, while the transient formation of semiquinone is not proton-coupled, the second eT process couples to a semiconcerted proton uptake from conserved tyrosine (Tyr-87) and histidine (His-38) residues within the active site. Our calculations indicate that both UQ and MQ have low redox potentials around -260 and -230 mV, respectively, in the Q-binding site, respectively, suggesting that release of the Q toward the membrane is coupled to an energy transduction step that could thermodynamically drive proton pumping in complex I.

Redox potential tuning through differential quinone binding in the photosynthetic reaction center of Rhodobacter sphaeroides

Ubiquinone forms an integral part of the electron transport chain in cellular respiration and photosynthesis across a vast number of organisms. Prior experimental results have shown that the photosynthetic reaction center (RC) from Rhodobacter sphaeroides is only fully functional with a limited set of methoxy-bearing quinones, suggesting that specific interactions with this substituent are required to drive electron transport and the formation of quinol. The nature of these interactions has yet to be determined. Through parameterization of a CHARMM-compatible quinone force field and subsequent molecular dynamics simulations of the quinone-bound RC, we have investigated and characterized the interactions of the protein with the quinones in the Q(A) and Q(B) sites using both equilibrium simulation and thermodynamic integration. In particular, we identify a specific interaction between the 2-methoxy group of ubiquinone in the Q(B) site and the amide nitrogen of GlyL225 that we implicate in locking the orientation of the 2-methoxy group, thereby tuning the redox potential difference between the quinones occupying the Q(A) and Q(B) sites. Disruption of this interaction leads to weaker binding in a ubiquinone analogue that lacks a 2-methoxy group, a finding supported by reverse electron transfer electron paramagnetic resonance experiments of the Q(A)⁻Q(B)⁻ biradical and competitive binding assays.

Proton-Coupled Electron Transfer During the S-State Transitions of the Oxygen-Evolving Complex of Photosystem II

The oxygen-evolving complex (OEC) of photosystem II (PSII) is a unique Mn4O5Ca cluster that catalyzes water oxidation via four photoactivated electron transfer steps. As the protein influence on the redox and protonation chemistry of the OEC remains an open question, we present a classical valence model of the OEC that allows the redox state of each Mn and the protonation state of bridging μ-oxos and terminal waters to remain in equilibrium with the PSII protein throughout the redox cycle. We find that the last bridging oxygen loses its proton during the transition from S0 to S1. Two possible S2 states are found depending on the OEC geometry: S2 has Mn4(IV) with a proton lost from a terminal water (W1) trapped by the nearby D1-D61 if O5 is closer to Mn4, or Mn1(IV), with partial deprotonation of D1-H337 and D1-E329 if O5 is closer to Mn1. In S3, the OEC is Mn4(IV) with W2 deprotonated. The estimated OEC Em's range from +0.7 to +1.3 V, enabling oxidation by P680(+), the primary electron donor in PSII. In chloride-depleted PSII, the proton release increases during the S1 to S2 transition, leaving the OEC unable to properly advance through the water-splitting cycle.

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.

Mechanism of proton-coupled quinone reduction in Photosystem II

Photosystem II uses light to drive water oxidation and plastoquinone (PQ) reduction. PQ reduction involves two PQ cofactors, Q(A) and Q(B), working in series. Q(A) is a one-electron carrier, whereas Q(B) undergoes sequential reduction and protonation to form Q(B)H(2). Q(B)H(2) exchanges with PQ from the pool in the membrane. Based on the atomic coordinates of the Photosystem II crystal structure, we analyzed the proton transfer (PT) energetics adopting a quantum mechanical/molecular mechanical approach. The potential-energy profile suggests that the initial PT to Q(B)(•-) occurs from the protonated, D1-His252 to Q(B)(•)(-) via D1-Ser264. The second PT is likely to occur from D1-His215 to Q(B)H(-) via an H-bond with an energy profile with a single well, resulting in the formation of Q(B)H(2) and the D1-His215 anion. The pathway for reprotonation of D1-His215(-) may involve bicarbonate, D1-Tyr246 and water in the Q(B) site. Formate ligation to Fe(2+) did not significantly affect the protonation of reduced Q(B), suggesting that formate inhibits Q(B)H(2) release rather than its formation. The presence of carbonate rather than bicarbonate seems unlikely because the calculations showed that this greatly perturbed the potential of the nonheme iron, stabilizing the Fe(3+) state in the presence of Q(B)(•-), a situation not encountered experimentally. H-bonding from D1-Tyr246 and D2-Tyr244 to the bicarbonate ligand of the nonheme iron contributes to the stability of the semiquinones. A detailed mechanistic model for Q(B) reduction is presented.

Slow dissociation of a charged ligand: analysis of the primary quinone Q(A) site of photosynthetic bacterial reaction centers

Reaction centers (RCs) are integral membrane proteins that undergo a series of electron transfer reactions during the process of photosynthesis. In the Q(A) site of RCs from Rhodobacter sphaeroides, ubiquinone-10 is reduced, by a single electron transfer, to its semiquinone. The neutral quinone and anionic semiquinone have similar affinities, which is required for correct in situ reaction thermodynamics. A previous study showed that despite similar affinities, anionic quinones associate and dissociate from the Q(A) site at rates ≈10(4) times slower than neutral quinones indicating that anionic quinones encounter larger binding barriers (Madeo, J.; Gunner, M. R. Modeling binding kinetics at the Q(A) site in bacterial reaction centers. Biochemistry 2005, 44, 10994-11004). The present study investigates these barriers computationally, using steered molecular dynamics (SMD) to model the unbinding of neutral ground state ubiquinone (UQ) and its reduced anionic semiquinone (SQ(-)) from the Q(A) site. In agreement with experiment, the SMD unbinding barrier for SQ(-) is larger than for UQ. Multi Conformational Continuum Electrostatics (MCCE), used here to calculate the binding energy, shows that SQ(-) and UQ have comparable affinities. In the Q(A) site, there are stronger binding interactions for SQ(-) compared to UQ, especially electrostatic attraction to a bound non-heme Fe(2+). These interactions compensate for the higher SQ(-) desolvation penalty, allowing both redox states to have similar affinities. These additional interactions also increase the dissociation barrier for SQ(-) relative to UQ. Thus, the slower SQ(-) dissociation rate is a direct physical consequence of the additional binding interactions required to achieve a Q(A) site affinity similar to that of UQ. By a similar mechanism, the slower association rate is caused by stronger interactions between SQ(-) and the polar solvent. Thus, stronger interactions for both the unbound and bound states of charged and highly polar ligands can slow their binding kinetics without a conformational gate. Implications of this for other systems are discussed.

The measured and calculated affinity of methyl- and methoxy-substituted benzoquinones for the Q(A) site of bacterial reaction centers

Quinones play important roles in mitochondrial and photosynthetic energy conversion acting as intramembrane, mobile electron, and proton carriers between catalytic sites in various electron transfer proteins. They display different affinity, selectivity, functionality, and exchange dynamics in different binding sites. The computational analysis of quinone binding sheds light on the requirements for quinone affinity and specificity. The affinities of 10 oxidized, neutral benzoquinones were measured for the high affinity Q(A) site in the detergent-solubilized Rhodobacter sphaeroides bacterial photosynthetic reaction center. Multiconformation Continuum Electrostatics was then used to calculate their relative binding free energies by grand canonical Monte Carlo sampling with a rigid protein backbone, flexible ligand, and side chain positions and protonation states. Van der Waals and torsion energies, Poisson-Boltzmann continuum electrostatics, and accessible surface area-dependent ligand-solvent interactions are considered. An initial, single cycle of GROMACS backbone optimization improves the match with experiment as do coupled-ligand and side-chain motions. The calculations match experiment with an root mean square deviation (RMSD) of 2.29 and a slope of 1.28. The affinities are dominated by favorable protein-ligand van der Waals rather than electrostatic interactions. Each quinone appears in a closely clustered set of positions. Methyl and methoxy groups move into the same positions as found for the native quinone. Difficulties putting methyls into methoxy sites are observed. Calculations using a solvent-accessible surface area-dependent implicit van der Waals interaction smoothed out small clashes, providing a better match to experiment with a RMSD of 0.77 and a slope of 0.97.

Hydration of Watson-Crick base pairs and dehydration of Hoogsteen base pairs inducing structural polymorphism under molecular crowding conditions

It has been revealed recently that molecular crowding, which is one of the largest differences between in vivo and in vitro conditions, is a critical factor determining the structure, stability, and function of nucleic acids. However, the effects of molecular crowding on Watson-Crick and Hoogsteen base pairs remain unclear. In order to investigate directly and quantitatively the molecular crowding effects on base pair types in nucleic acids, we designed intramolecular parallel- and antiparallel-stranded DNA duplexes consisting of Hoogsteen and Watson-Crick base pairs, respectively, as well as an intramolecular parallel-stranded triplex containing both types of base pairs. Thermodynamic analyses demonstrated that the values of free energy change at 25 degrees C for Hoogsteen base-pair formations decreased from +1.45 +/- 0.15 to +1.09 +/- 0.13 kcal mol(-1), and from -1.89 +/- 0.13 to -2.71 +/- 0.11 kcal mol(-1) in the intramolecular duplex and triplex, respectively, when the concentration of PEG 200 (polyethylene glycol with average molecular weight 200) increased from 0 to 20 wt %. However, corresponding values for Watson-Crick formation in the duplex and triplex increased from -10.2 +/- 0.2 to -8.7 +/- 0.1 kcal mol(-1), and from -10.8 +/- 0.2 to -9.2 +/- 0.2 kcal mol(-1), respectively. Furthermore, it was revealed that the opposing effects of molecular crowding on the Hoogsteen and Watson-Crick base pairs were due to different behaviors of water molecules binding to the DNA strands.

Modification of quinone electrochemistry by the proteins in the biological electron transfer chains: examples from photosynthetic reaction centers

Quinones such as ubiquinone are the lipid soluble electron and proton carriers in the membranes of mitochondria, chloroplasts and oxygenic bacteria. Quinones undergo controlled redox reactions bound to specific sites in integral membrane proteins such as the cytochrome bc(1) oxidoreductase. The quinone reactions in bacterial photosynthesis are amongst the best characterized, presenting a model to understand how proteins modulate cofactor chemistry. The free energy of ubiquinone redox reactions in aqueous solution and in the Q(A) and Q(B) sites of the bacterial photosynthetic reaction centers (RCs) are compared. In the primary Q(A) site ubiquinone is reduced only to the anionic semiquinone (Q(*-)) while in the secondary Q(B) site the product is the doubly reduced, doubly protonated quinol (QH(2)). The ways in which the protein modifies the relative energy of each reduced and protonated intermediate are described. For example, the protein stabilizes Q(*-) while destabilizing Q(=) relative to aqueous solution through electrostatic interactions. In addition, kinetic and thermodynamic mechanisms for stabilizing the intermediate semiquinones are compared. Evidence for the protein sequestering anionic compounds by slowing both on and off rates as well as by binding the anion more tightly is reviewed.

Investigating the mechanisms of photosynthetic proteins using continuum electrostatics

Computational methods based on continuum electrostatics are widely used in theoretical biochemistry to analyze the function of proteins. Continuum electrostatic methods in combination with quantum chemical and molecular mechanical methods can help to analyze even very complex biochemical systems. In this article, applications of these methods to proteins involved in photosynthesis are reviewed. After giving a short introduction to the basic concepts of the continuum electrostatic model based on the Poisson-Boltzmann equation, we describe the application of this approach to the docking of electron transfer proteins, to the comparison of isofunctional proteins, to the tuning of absorption spectra, to the analysis of the coupling of electron and proton transfer, to the analysis of the effect of membrane potentials on the energetics of membrane proteins, and to the kinetics of charge transfer reactions. Simulations as those reviewed in this article help to analyze molecular mechanisms on the basis of the structure of the protein, guide new experiments, and provide a better and deeper understanding of protein functions.

ENDOR spectroscopy reveals light induced movement of the H-bond from Ser-L223 upon forming the semiquinone (Q(B)(-)(*)) in reaction centers from Rhodobacter sphaeroides

Proton ENDOR spectroscopy was used to monitor local conformational changes in bacterial reaction centers (RC) associated with the electron-transfer reaction DQB --> D+*QB-* using mutant RCs capable of photoreducing QB at cryogenic temperatures. The charge separated state D+*QB-* was studied in mutant RCs formed by either (i) illuminating at low temperature (77 K) a sample frozen in the dark (ground state protein conformation) or (ii) illuminating at room temperature prior to and during freezing (charge separated state protein conformation). The charge recombination rates from the two states differed greatly (>10(6) fold) as shown previously, indicating a structural change (Paddock et al. (2006) Biochemistry 45, 14032-14042). ENDOR spectra of QB-* from both samples (35 GHz, 77 K) showed several H-bond hyperfine couplings that were similar to those for QB-* in native RCs indicating that in all RCs, QB-* was located at the proximal position near the metal site. In contrast, one set of hyperfine couplings were not observed in the dark frozen samples but were observed only in samples frozen under illumination in which the protein can relax prior to freezing. This flexible H-bond was assigned to an interaction between the Ser-L223 hydroxyl and QB-* on the basis of its absence in Ser L223 --> Ala mutant RCs. Thus, part of the protein relaxation, in response to light induced charge separation, involves the formation of an H-bond between the OH group of Ser-L223 and the anionic semiquinone QB-*. These results show the flexibility of the Ser-L223 H-bond, which is essential for its function in proton transfer to reduced QB.

pH modulates the quinone position in the photosynthetic reaction center from Rhodobacter sphaeroides in the neutral and charge separated states

The structure of the photosynthetic reaction-center from Rhodobacter sphaeroides has been determined at four different pH values (6.5, 8.0, 9.0, 10.0) in the neutral and in charge separated states. At pH 8.0, in the neutral state, we obtain a resolution of 1.87 A, which is the best ever reported for the bacterial reaction center protein. Our crystallographic data confirm the existence of two different binding positions of the secondary quinone (QB). We observe a new orientation of QB in its distal position, which shows no ring-flip compared to the orientation in the proximal position. Datasets collected for the different pH values show a pH-dependence of the population of the proximal position. The new orientation of QB in the distal position and the pH-dependence could be confirmed by continuum electrostatics calculations. Our calculations are in agreement with the experimentally observed proton uptake upon charge separation. The high resolution of our crystallographic data allows us to identify new water molecules and external residues being involved in two previously described hydrogen bond proton channels. These extended proton-transfer pathways, ending at either of the two oxo-groups of QB in its proximal position, provide additional evidence that ring-flipping is not required for complete protonation of QB upon reduction.

Protonation States of Ammonia/Ammonium in the Hydrophobic Pore of Ammonia Transporter Protein AmtB

The crystal structure of the ammonia transport (Amt) protein AmtB at 1.4 Angstrom resolution revealed four ammonia/ammonium (NH(3)/NH(4)(+)) binding sites along the approximately 20 Angstrom narrow pore. It is an open question whether the bound NH(3)/NH(4)(+) are neutral (NH(3)) or cationic (NH(4)(+)). On the basis of the AmtB crystal structure, we calculated the pK(a) of these four NH(3)/NH(4)(+) by solving the Poisson-Boltzmann equation. Except for one NH(3)/NH(4)(+) binding site (Am1) at the entry point of the Amt pore, binding sites are occupied by NH(3) due to lack of energy contributions from solvation, eliminating an existence of charged form NH(4)(+) and, inevitably, its potential cation-pi interaction. The only two titratable residues in the pore, His168 and His318, are in the neutral charge state. The NH(4)(+) charge state at the Am1 site is stabilized by Ser219 functioning as an H-bond acceptor. However, when involving explicit crystal water nearby, the NH(3) charge state is stabilized by the reorientation of Ser219-OH group. This H-bond donor Ser219 significantly decreases the pK(a) of NH(3)/ NH(4)(+) at the Am1 site to approximately 1. The flip/flop H-bond of Ser219 may play a dual role first in binding and subsequently in deprotonating NH(4)(+), which is a prerequisite to conduct NH(3) through the Amt pore across the membrane.

Monitoring the pH Dependence of IR Carboxylic Acid Signals upon QB- Formation in the Glu-L212 → Asp/Asp-L213 → Glu Swap Mutant Reaction Center from Rhodobacter sphaeroides

In the photosynthetic reaction center (RC) from the purple bacterium Rhodobacter sphaeroides, proton-coupled electron-transfer reactions occur at the secondary quinone (QB) site. Involved in the proton uptake steps are carboxylic acids, which have characteristic infrared vibrations in the 1770-1700 cm-1 spectral range that are sensitive to 1H/2H isotopic exchange. With respect to the native RC, a novel protonation pattern for carboxylic acids upon QB photoreduction has been identified in the Glu-L212 --> Asp/Asp-L213 --> Glu mutant RC using light-induced FTIR difference spectroscopy (Nabedryk, E., Breton, J., Okamura, M. Y., and Paddock, M. L. (2004) Biochemistry 43, 7236-7243). These carboxylic acids are structurally close and have been implicated in proton transfer to reduced QB. In this work, we extend previous studies by measuring the pH dependence of the QB-/QB FTIR difference spectra of the mutant in 1H2O and 2H2O. Large pH dependent changes were observed in the 1770-1700 cm-1 spectral range between pH 8 and pH 4. The IR fingerprints of the protonating carboxylic acids upon QB- formation were obtained from the calculated double-difference spectra 1H2O minus 2H2O. These IR fingerprints are specific for each pH, indicative of the contribution of different titrating groups. In particular, the 1752 cm-1 signal indicates that Glu-L213 protonates upon QB- formation at pH >or= 5, whereas the 1746 cm-1 signal indicates protonation of Asp-L212 even at pH 4. An unidentified carboxylic acid absorbing at approximately 1765 cm-1 could be the proton donor between pH 8 and 5. The observation that in the swap mutant there are several uniquely behaving carboxylic acids shows that electrostatic interactions occurring between them are sufficiently modified from the native RC to reveal their IR signatures.

Trapped conformational states of semiquinone (D+*QB-*) formed by B-branch electron transfer at low temperature in Rhodobacter sphaeroides reaction centers

The reaction center (RC) from Rhodobacter sphaeroides captures light energy by electron transfer between quinones QA and QB, involving a conformational gating step. In this work, conformational states of D+*QB-* were trapped (80 K) and studied using EPR spectroscopy in native and mutant RCs that lack QA in which QB was reduced by the bacteriopheophytin along the B-branch. In mutant RCs frozen in the dark, a light induced EPR signal due to D+*QB-* formed in 30% of the sample with low quantum yield (0.2%-20%) and decayed in 6 s. A small signal with similar characteristics was also observed in native RCs. In contrast, the EPR signal due to D+*QB-* in mutant RCs illuminated while freezing formed in approximately 95% of the sample did not decay (tau >107 s) at 80 K (also observed in the native RC). In all samples, the observed g-values were the same (g = 2.0026), indicating that all active QB-*'s were located in a proximal conformation coupled with the nonheme Fe2+. We propose that before electron transfer at 80 K, the majority (approximately 70%) of QB, structurally located in the distal site, was not stably reducible, whereas the minority (approximately 30%) of active configurations was in the proximal site. The large difference in the lifetimes of the unrelaxed and relaxed D+*QB-* states is attributed to the relaxation of protein residues and internal water molecules that stabilize D+*QB-*. These results demonstrate energetically significant conformational changes involved in stabilizing the D+*QB-* state. The unrelaxed and relaxed states can be considered to be the initial and final states along the reaction coordinate for conformationally gated electron transfer.

How proteins trigger excitation energy transfer in the FMO complex of green sulfur bacteria

A simple electrostatic method for the calculation of optical transition energies of pigments in protein environments is presented and applied to the Fenna-Matthews-Olson (FMO) complex of Prosthecochloris aestuarii and Chlorobium tepidum. The method, for the first time, allows us to reach agreement between experimental optical spectra and calculations based on transition energies of pigments that are calculated in large part independently, rather than fitted to the spectra. In this way it becomes possible to understand the molecular mechanism allowing the protein to trigger excitation energy transfer reactions. The relative shift in excitation energies of the seven bacteriochlorophyll-a pigments of the FMO complex of P. aestuarii and C. tepidum are obtained from calculations of electrochromic shifts due to charged amino acids, assuming a standard protonation pattern of the protein, and by taking into account the three different ligand types of the pigments. The calculations provide an explanation of some of the earlier results for the transition energies obtained from fits of optical spectra. In addition, those earlier fits are verified here by using a more advanced theory of optical spectra, a genetic algorithm, and excitonic couplings obtained from electrostatic calculations that take into account the influence of the dielectric protein environment. The two independent calculations of site energies strongly favor one of the two possible orientations of the FMO trimer relative to the photosynthetic membrane, which were identified by electron microscopic studies and linear dichroism experiments. Efficient transfer of excitation energy to the reaction center requires bacteriochlorophylls 3 and 4 to be the linker pigments. The temporal and spatial transfer of excitation energy through the FMO complex is calculated to proceed along two branches, with transfer times that differ by an order of magnitude.

Electrostatic role of the non-heme iron complex in bacterial photosynthetic reaction center

To elucidate the role of the non-heme iron complex (Fe-complex) in the electron transfer (ET) events of bacterial photosynthetic reaction centers (bRC), we calculated redox potentials of primary/secondary quinones Q(A/B) (E(m)(Q(A/B))) in the Fe-depleted bRC. Removing the Fe-complex, the calculated E(m)(Q(A/B)) are downshifted by approximately 220 mV/ approximately 80 mV explaining both the 15-fold decrease in ET rate from bacteriopheophytin (H(A)(-)) to Q(A) and triplet state occurrence in Fe-depleted bRC. The larger downshift in E(m)(Q(A)) relative to E(m)(Q(B)) increases the driving-energy for ET from Q(A) to Q(B) by 140 meV, in agreement with approximately 100 meV increase derived from kinetic studies.

Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations

A protein structure should provide the information needed to understand its observed properties. Significant progress has been made in developing accurate calculations of acid/base and oxidation/reduction reactions in proteins. Current methods and their strengths and weaknesses are discussed. The distribution and calculated ionization states in a survey of proteins is described, showing that a significant minority of acidic and basic residues are buried in the protein and that most of these remain ionized. The electrochemistry of heme and quinones are considered. Proton transfers in bacteriorhodopsin and coupled electron and proton transfers in photosynthetic reaction centers, 5-coordinate heme binding proteins and cytochrome c oxidase are highlighted as systems where calculations have provided insight into the reaction mechanism.

Ubiquinone (coenzyme Q10) binding sites: Low dielectric constant of the gate prevents the escape of the semiquinone

The photosynthetic reaction center (RC) from purple bacteria is frequently used as a model for the interaction of ubiquinones (coenzyme Q) with membrane proteins. Single-turnover flash activation of RC leads to formation of the semiquinone (SQ) of the secondary acceptor quinone after odd flashes and quinol after even flashes. The ubiquinol escapes the binding site in 1 ms, while the SQ does not leave the binding site for at least 5 min. Observed difference between these times suggests a large energetic barrier for the SQ. However, high apparent dielectric constant in the vicinity of the quinone ring (>or=25) results in a relatively small electrostatic energy of SQ stabilization. To resolve this apparent contradiction I suggest that a significant part of the kinetic stabilization of the SQ is achieved by the special topology of the binding site in which quinone can exit the binding site only by moving its headgroup toward the center of the membrane. The large energetic penalty of transferring the charged headgroup to the membrane dielectric can explain the observed kinetic stability of the SQ.

Structural and Computational Analysis of the Quinone-binding Site of Complex II (Succinate-Ubiquinone Oxidoreductase)

The transfer of electrons and protons between membrane-bound respiratory complexes is facilitated by lipid-soluble redox-active quinone molecules (Q). This work presents a structural analysis of the quinone-binding site (Q-site) identified in succinate:ubiquinone oxidoreductase (SQR) from Escherichia coli. SQR, often referred to as Complex II or succinate dehydrogenase, is a functional member of the Krebs cycle and the aerobic respiratory chain and couples the oxidation of succinate to fumarate with the reduction of quinone to quinol (QH(2)). The interaction between ubiquinone and the Q-site of the protein appears to be mediated solely by hydrogen bonding between the O1 carbonyl group of the quinone and the side chain of a conserved tyrosine residue. In this work, SQR was co-crystallized with the ubiquinone binding-site inhibitor Atpenin A5 (AA5) to confirm the binding position of the inhibitor and reveal additional structural details of the Q-site. The electron density for AA5 was located within the same hydrophobic pocket as ubiquinone at, however, a different position within the pocket. AA5 was bound deeper into the site prompting further assessment using protein-ligand docking experiments in silico. The initial interpretation of the Q-site was re-evaluated in the light of the new SQR-AA5 structure and protein-ligand docking data. Two binding positions, the Q(1)-site and Q(2)-site, are proposed for the E. coli SQR quinone-binding site to explain these data. At the Q(2)-site, the side chains of a serine and histidine residue are suitably positioned to provide hydrogen bonding partners to the O4 carbonyl and methoxy groups of ubiquinone, respectively. This allows us to propose a mechanism for the reduction of ubiquinone during the catalytic turnover of the enzyme.

Modulation of Fragmental Charge Transfer via Hydrogen Bonds. Direct Measurement of Electronic Contributions†

Hydrogen bonds play an important role in an overwhelming variety of fields from biology to surface and supramolecular chemistry. The term "hydrogen bond" refers to a wide range of interactions with various covalent and polar contributions. In particular, hydrogen bonds have an important role in the folding and packing of peptides and nucleic acids. Recent studies also point to the importance of hydrogen bonding in the context of second-shell interactions, in metal binding and selectivity in metalloproteins, and in controlling the dynamics of membrane proteins. In this study, we demonstrate and quantify the modulation of fragmental charge transfer from hydrogen-bonded ligands to a metal center, by employing our recently introduced molecular potentiometer. The molecular details that affect this type of fragmental charge transfer are presented and a path for transferring chemical information is demonstrated. We found that H-bond interactions in the extended positions of axial ligands provide an effective means of modulating the amount of fragmental charge transfer to a metal center, thereby dramatically influencing the electronic properties of the ligand, the binding affinity, and the binding of additional ligands. The magnitude of fragmental charge-transfer modulation induced by a single ligand-solvent H-bond interaction is comparable to those induced by covalent substitution, although H-bond enthalpy is only on the order of several kilojoules per mole. Importantly, we find a significant change in the ligand electronic properties, even for weak C-H...O=C H-bond formation, where the bond enthalpy is substantially lower than for conventional H-bond interactions. The excess fragmental charge transferred to the metal center, deduced from the spectroscopic measurements, correlates well with the computationally determined values. Our findings underscore the importance of second-shell interactions in the active sites of enzymes, beyond the structural and electrostatic importance that is widely recognized today.

Induced conformational changes upon Cd2+ binding at photosynthetic reaction centers

Cd(2+) binding at the bacterial photosynthetic reaction center (bRC) from Rhodobacter sphaeroides is known to inhibit proton transfer (PT) from bulk solvent to the secondary quinone Q(B). To elucidate this mechanism, we calculated the pK(a) for residues along the water channels connecting Q(B) with the stromal side based on the crystal structures of WT-bRC and Cd(2+)-bound bRC. Upon Cd(2+) binding, we observed the release of two protons from His-H126/128 at the Cd(2+) binding site and significant pK(a) shifts for residues along the PT pathways. Remarkably, Asp-L213 near Q(B), which is proposed to play a significant role in PT, resulted in a decrease in pK(a) upon Cd(2+) binding. The direct electrostatic influence of the Cd(2+)-positive charge on these pK(a) shifts was small. Instead, conformational changes of amino acid side chains induced electrostatically by Cd(2+) binding were the main mechanism for these pK(a) shifts. The long-range electrostatic influence over approximately 12 A between Cd(2+) and Asp-L213 is likely to originate from a set of Cd(2+)-induced successive reorientations of side chains (Asp-H124, His-H126, His-H128, Asp-H170, Glu-H173, Asp-M17, and Asp-L210), which propagate along the PT pathways as a "domino" effect.

Oxidation of the Non-Heme Iron Complex in Photosystem II

In photosystem II (PSII), the redox properties of the non-heme iron complex (Fe complex) are sensitive to the redox state of quinones (Q(A/)(B)), which may relate to the electron/proton transfer. We calculated the redox potentials for one-electron oxidation of the Fe complex in PSII [E(m)(Fe)] based on the reference value E(m)(Fe) = +400 mV at pH 7 in the Q(A)(0)Q(B)(0) state, considering the protein environment in atomic detail and the associated changes in protonation pattern. Our model yields the pH dependence of E(m)(Fe) with -60 mV/pH as observed in experimental redox titration. We observed significant deprotonation at D1-Glu244 in the hydrophilic loop region upon Fe complex oxidation. The calculated pK(a) value for D1-Glu244 depends on the Fe complex redox state, yielding a pK(a) of 7.5 and 5.5 for Fe(2+) and Fe(3+), respectively. To account for the pH dependence of E(m)(Fe), a model involving not only D1-Glu244 but also the other titratable residues (five Glu in the D-de loops and six basic residues near the Fe complex) seems to be needed, implying the existence of a network of residues serving as an internal proton reservoir. Reduction of Q(A/B) yields +302 mV and +268 mV for E(m)(Fe) in the Q(A)(-)Q(B)(0) and Q(A)(0)Q(B)(-) states, respectively. Upon formation of the Q(A)(0)Q(B)(-) state, D1-His252 becomes protonated. Forming Fe(3+)Q(B)H(2) by a proton-coupled electron transfer process from the initial state Fe(2+)Q(B)(-) results in deprotonation of D1-His252. The two EPR signals observed at g = 1.82 and g = 1.9 in the Fe(2+)Q(A)(-) state of PSII may be attributed to D1-His252 with variable and fixed protonation, respectively.

An Isotope-Edited FTIR Investigation of the Role of Ser-L223 in Binding Quinone (QB) and Semiquinone (QB-) in the Reaction Center from Rhodobacter sphaeroides

In the photosynthetic reaction center (RC) from the purple bacterium Rhodobacter sphaeroides, proton-coupled electron-transfer reactions occur at the secondary quinone (Q(B)) site. Several nearby residues are important for both binding and redox chemistry involved in the light-induced conversion from Q(B) to quinol Q(B)H(2). Ser-L223 is one of the functionally important residues located near Q(B). To obtain information on the interaction between Ser-L223 and Q(B) and Q(B)(-), isotope-edited Q(B)(-)/Q(B) FTIR difference spectra were measured in a mutant RC in which Ser-L223 is replaced with Ala and compared to the native RC. The isotope-edited IR fingerprint spectra for the C=O [see text] and C=C [see text] modes of Q(B) (Q(B)(-)) in the mutant are essentially the same as those of the native RC. These findings indicate that highly equivalent interactions of Q(B) and Q(B)(-) with the protein occur in both native and mutant RCs. The simplest explanation of these results is that Ser-L223 is not hydrogen bonded to Q(B) or Q(B)(-) but presumably forms a hydrogen bond to a nearby acid group, preferentially Asp-L213. The rotation of the Ser OH proton from Asp-L213 to Q(B)(-) is expected to be an important step in the proton transfer to the reduced quinone. In addition, the reduced quinone remains firmly bound, indicating that other distinct hydrogen bonds are more important for stabilizing Q(B)(-). Implications on the design features of the Q(B) binding site are discussed.

Control of Quinone Redox Potentials in Photosystem II: Electron Transfer and Photoprotection

In O(2)-evolving complex Photosystem II (PSII), an unimpeded transfer of electrons from the primary quinone (Q(A)) to the secondary quinone (Q(B)) is essential for the efficiency of photosynthesis. Recent PSII crystal structures revealed the protein environment of the Q(A/B) binding sites. We calculated the plastoquinone (Q(A/B)) redox potentials (E(m)) for one-electron reduction with a full account of the PSII protein environment. We found two different H-bond patterns involving Q(A) and D2-Thr217, resulting in an upshift of E(m)(Q(A)) by 100 mV if the H bond between Q(A) and Thr is present. The formation of this H bond to Q(A) may be the origin of a photoprotection mechanism, which is under debate. At the Q(B) side, the formation of a H bond between D2-Ser264 and Q(B) depends on the protonation state of D1-His252. Q(B) adopts the high-potential form if the H bond to Ser is present. Conservation of this residue and H-bond pattern for Q(B) sites among bacterial photosynthetic reaction centers (bRC) and PSII strongly indicates their essential requirement for electron transfer function.

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.

Redox potentials of chlorophylls in the photosystem II reaction center

Water oxidation generating atmospheric oxygen occurs in photosystem II (PSII), a large protein-pigment complex located in the thylakoid membrane. The recent crystal structures at 3.2 and 3.5 A resolutions provide novel details on amino acid side chains, especially in the D1/D2 subunits. We calculated the redox potentials for one-electron oxidation of the chlorophyll a (Chla) molecules in PSII, considering the protein environment in atomic detail. The calculated redox potentials for the dimer Chla (P(D1/D2)) and accessory Chla (Chl(D1/D2)) were 1.11-1.30 V relative to the normal hydrogen electrode at pH 7, which is high enough for water oxidation. The D1/D2 proteins and their cofactors contribute approximately 390 mV to the enormous upshift of 470 mV compared to the redox potential of monomeric Chla in dimethylformamide. The other subunits are responsible for the remaining 80 mV. The high redox potentials of the two accessory Chla Chl(D1/D2) suggests that they also participate in the charge separation process.

Tuning electron transfer by ester-group of chlorophylls in bacterial photosynthetic reaction center

Accessory chlorophylls (B(A/B)) in bacterial photosynthetic reaction center play a key role in charge-separation. Although light-exposed and dark-adapted bRC crystal structures are virtually identical, the calculated B(A) redox potentials for one-electron reduction differ. This can be traced back to different orientations of the B(A) ester-group. This tuning ability of chlorophyll redox potentials modulates the electron transfer from SP* to B(A).

Hydrogen Bonds Not Only Provide a Structural Scaffold to Assemble Donor and Acceptor Moieties of Zinc Porphyrin−Quinone Dyads but Also Control the Photoinduced Electron Transfer to Afford the Long-Lived Charge-Separated States

A series of zinc porphyrin-quinone linked dyads [ZnP-CONH-Q, ZnP-NHCO-Q, and ZnP-n-Q (n = 3, 6, 10)] were designed and synthesized to investigate the effects of hydrogen bonds which can not only provide a structural scaffold to assemble donor and acceptor moieties but also control the photoinduced electron-transfer process. In the case of ZnP-CONH-Q and ZnP-NHCO-Q, the hydrogen bond between the N-H proton and the carbonyl oxygen of Q results in the change in the reduction potential of Q. The strong hydrogen bond between the N-H proton and the carbonyl oxygen of Q*- in ZnP-CONH-Q*-,ZnP-NHCO-Q*-, and ZnP-n-Q*- (n = 3, 6, 10) generated by the chemical reduction has been confirmed by the ESR spectra, which exhibit hyperfine coupling constants in agreement those predicted by the density functional calculations. In the case of ZnP-n-Q (n = 3, 6, 10), on the other hand, the hydrogen bond between two amide groups provides a structural scaffold to assemble the donor (ZnP) and the acceptor (Q) moiety together with the hydrogen bond between the N-H proton and the carbonyl oxygen of Q, leading to attainment of the charge-separated state with a long lifetime up to a microsecond.

Energetics of Proton Transfer Pathways in Reaction Centers from Rhodobacter sphaeroides

Electron transfer between the primary and secondary quinones (Q(A), Q(B)) in the bacterial photosynthetic reaction center (bRC) is coupled with proton uptake at Q(B). The protons are conducted from the cytoplasmic side, probably with the participation of two water channels. Mutations of titratable residues like Asp-L213 to Asn (inhibited mutant) or the double mutant Glu-L212 to Ala/Asp-L213 to Ala inhibit these electron transfer-coupled proton uptake events. The inhibition of the proton transfer (PT) process in the single mutant can be restored by a second mutation of Arg-M233 to Cys or Arg-H177 to His (revertant mutant). These revertant mutants shed light on the location of the main proton transfer pathway of wild type bRC. In contrast to the wild type and inhibited mutant bRC, the revertant mutant bRC showed notable proton uptake at Glu-H173 upon formation of the Q(B)- state. In all of these mutants, the pK(a) of Asp-M17 decreased by 1.4-2.4 units with respect to the wild type bRC, whereas a significant pK(a) upshift of up to 5.8 units was observed at Glu-H122, Asp-H170, Glu-H173, and Glu-H230 in the revertant mutants. These residues belonging to the main PT pathway are arranged along water channel P1 localized mainly in subunit H. bRC possesses subunit H, which has no counterpart in photosystem II. Thus, bRC may possess alternative PT pathways involving water channels in subunit H, which becomes active in case the main PT pathway is blocked.

Redox Potentials of Chlorophylls and β-Carotene in the Antenna Complexes of Photosystem II

Electron transfer (ET) processes in reaction centers (RC) of photosystem II (PSII) are prerequisites of oxygen generation. They are promoted by energy transfer from antenna to RC. Here, we calculated the redox potentials of chlorophylla/beta-carotene (Chla/Car) in PSII CP43/CP47 antenna complexes, solving the linearized Poisson-Boltzmann (LPB) equation based on the PSII crystal structure. The majority of antenna Chla redox potentials for reduction/oxidation were lower than those of RC Chla. Hence, ET events with excess electrons remain localized in the RC. Simultaneously antenna Chla can serve as an efficient cation sink to rereduce RC Chla if normal PSII function is inhibited. Especially three antenna Chla (Chl-47, Chl-18, and Chl-12) and two Car bridging the space between Chl(Z(D1)) and cytochrome (cyt) b559 have the same level of oxidation redox potential. Together with Chl(Z(D2)) they form an electron hole transfer pathway and temporary storage device guiding from the oxidized P680(+.) Chla to the cyt b559. This path may play a photoprotective role as efficient electron hole quencher.

Energetics of Quinone-Dependent Electron and Proton Transfers in Rhodobacter sphaeroides Photosynthetic Reaction Centers

Proteins bind redox cofactors, modifying their electrochemistry and affinity by specific interactions of the binding site with each cofactor redox state. Photosynthetic reaction centers from Rhodobacter sphaeroides have three ubiquinone-binding sites, Q(A), and proximal and distal Q(B) sites. Ubiquinones, which can be doubly reduced and bind 2 protons, have 9 redox states. However, only Q and Q(-) are seen in the Q(A) site and Q, Q(-), and QH(2) in the proximal Q(B) site. The distal Q(B) function is uncertain. Multiple conformation continuum electrostatics (MCCE) was used to compare the ubiquinone electrochemical midpoints (E(m)) and pK(a) values at these three sites. At pH 7, the Q(A)/Q(A)(-) E(m) is -40 mV and proximal Q(B)/Q(B)(-) -10 mV in agreement with the experimental values (assuming a solution ubiquinone E(m) of -145 mV). Q(B) reduction requires changes in nearby residue protonation and SerL223 reorientation. The distal Q(B)/Q(B)(-) E(m) is a much more unfavorable -260 mV. Q(A) and proximal Q(B) sites generally stabilize species with a -1 charge, while the distal Q(B) site prefers binding neutral species. In each site, the dianion is destabilized because favorable interactions with the residues and backbone increase with charge (q), while the unfavorable loss of solvation energy increases with q(2). Therefore, proton binding before a second reduction, forming QH and then QH(-), is always preferred to forming the dianion (Q(-)(2)). The final product QH(2) is higher in energy at the proximal Q(B) site than in solution; therefore, it binds poorly, favoring release. In contrast, QH(2) binds more tightly than Q at the distal Q(B) site.