Factors That Differentiate the H-bond Strengths of Water Near the Schiff Bases in Bacteriorhodopsin and Anabaena Sensory Rhodopsin

Bacteriorhodopsin (BR) functions as a light-driven proton pump, whereas Anabaena sensory rhodopsin (ASR) is believed to function as a photosensor despite the high similarity in their protein sequences. In Fourier transform infrared (FTIR) spectroscopic studies, the lowest O-D stretch for D₂O was observed at ∼2200 cm⁻¹ in BR but was significantly higher in ASR (>2500 cm⁻¹), which was previously attributed to a water molecule near the Schiff base (W402) that is H-bonded to Asp-85 in BR and Asp-75 in ASR. We investigated the factors that differentiate the lowest O-D stretches of W402 in BR and ASR. Quantum mechanical/molecular mechanical calculations reproduced the H-bond geometries of the crystal structures, and the calculated O-D stretching frequencies were corroborated by the FTIR band assignments. The potential energy profiles indicate that the smaller O-D stretching frequency in BR originates from the significantly higher pK_a(Asp-85) in BR relative to the pK_a(Asp-75) in ASR, which were calculated to be 1.5 and −5.1, respectively. The difference is mostly due to the influences of Ala-53, Arg-82, Glu-194–Glu-204, and Asp-212 on pK_a(Asp-85) in BR and the corresponding residues Ser-47, Arg-72, Ser-188-Asp-198, and Pro-206 on pK_a(Asp-75) in ASR. Because these residues participate in proton transfer pathways in BR but not in ASR, the presence of a strongly H-bonded water molecule near the Schiff base ultimately results from the proton-pumping activity in BR.     

Cationic state distribution over the chlorophyll d-containing P(D1)/P(D2) pair in photosystem II

Most of the chlorophyll (Chl) cofactors in photosystem II (PSII) from Acaryochloris marina are Chld, although a few Chla molecules are also present. To evaluate the possibility that Chla may participate in the P(D1)/P(D2) Chl pair in PSII from A. marina, the P(D1)(?+)/P(D2)(?+) charge ratio was investigated using the PSII crystal structure analyzed at 1.9-? resolution, while considering all possibilities for the Chld-containing P(D1)/P(D2) pair, i.e., Chld/Chld, Chla/Chld, and Chld/Chla pairs. Chld/Chld and Chla/Chld pairs resulted in a large P(D1)(?+) population relative to P(D2)(?+), as identified in Chla/Chla homodimer pairs in PSII from other species, e.g., Thermosynechococcus elongatus PSII. However, the Chld/Chla pair possessed a P(D1)(?+)/P(D2)(?+) ratio of approximately 50/50, which is in contrast to previous spectroscopic studies on A. marina PSII. The present results strongly exclude the possibility that the Chld/Chla pair serves as P(D1)/P(D2) in A. marina PSII. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.     

Redox potential difference between Desulfovibrio vulgaris and Clostridium beijerinckii flavodoxins

The redox potential of the flavin mononucleotide (FMN) hydroquinones for one-electron reduction in the Desulfovibrio vulgaris ( D. vulgaris) flavodoxin ( E sq/hq for FMNH (*)/FMNH (-)) was calculated using the crystal structure of the relevant hydroquinone form and compared to the results of the Clostridium beijerinckii ( C. beijerinckii) flavodoxin. In D. vulgaris and C. beijerinckii flavodoxins, the protein side chain causes significant downshifts of 170 and 240 mV in E sq/hq, respectively. In the C. beijerinckii flavodoxin, the E sq/hq downshift because of the protein side chain is essentially compensated by the counter influence of the protein backbone ( E sq/hq upshift of 260 mV). However, in the D. vulgaris flavodoxin, the corresponding protein backbone influence on E sq/hq is significantly small, i.e., less than half of that in the C. beijerinckii flavodoxin. In particular, there is a significant difference in the influence of the protein backbone of the so-called 60s loop region between the two flavodoxins. The E sq/hq difference can be best explained by the lower compensation of the side chain influence by the backbone influence in the D. vulgaris flavodoxin than in the C. beijerinckii flavodoxin.     

Redox potential of the non-heme iron complex in bacterial photosynthetic reaction center

In bacterial photosynthetic reaction centers (bRC), the electron is transferred from the special pair (P) via accessory bacteriochlorophyll (B_A), bacteriopheopytin (H_A), the primary quinone (Q_A) to the secondary quinone (Q_B). Although the non-heme iron complex (Fe complex) is located between Q_A and Q_B, it was generally supposed not to be redox-active. Involvement of the Fe complex in electron transfer (ET) was proposed in recent FTIR studies [A. Remy and K. Gerwert, Coupling of light-induced electron transfer to proton uptake in photosynthesis, Nat. Struct. Biol. 10 (2003) 637-644]. However, other FTIR studies resulted in opposite results [J. Breton, Steady-state FTIR spectra of the photoreduction of Q_A and Q_B in Rhodobacter sphaeroides reaction centers provide evidence against the presence of a proposed transient electron acceptor X between the two quinones, Biochemistry 46 (2007) 4459-4465]. In this study, we calculated redox potentials of Q_A/B (E_m(Q_A/B)) and the Fe complex (E_m(Fe)) based on crystal structure of the wild-type bRC (WT-bRC), and we investigated the energetics of the system where the Fe complex is assumed to be involved in the ET. E_m(Fe) in WT-bRC is much less pH-dependent than that in PSII. In WT-bRC, we observed significant coupling of ET with Glu-L212 protonation upon oxidation of the Fe complex and a dramatic E_m(Fe) downshift by 230 mV upon formation of Q_A⁻ (but not Q_B⁻) due to the absence of proton uptake of Glu-L212. Changes in net charges of the His ligands of the Fe complex appear to be the nature of the redox event if we assume the involvement of the Fe complex in the ET.     

Function of two β-carotenes near the D1 and D2 proteins in photosystem II dimers

The antenna proteins in photosystem II (PSII) not only promote energy transfer to the photosynthetic reaction center (RC) but provide also an efficient cation sink to re-reduce chlorophyll a if the electron transfer (ET) from the Mn-cluster is inhibited. Using the newest PSII dimer crystal structure (3.0 Å resolution), in which 11 β-carotene molecules (Car) and 14 lipids are visible in the PSII monomer, we calculated the redox potentials (E_m) of one-electron oxidation for all Car (E_m(Car)) by solving the Poisson-Boltzmann equation. In each PSII monomer, the D1 protein harbors a previously unlocated Car (Car_D1) in van der Waals contact with the chlorin ring of Chl_Z(D1). Each Car_D1 in the PSII dimer complex is located in the interface between the D1 and CP47 subunits, together with another four Car of the other PSII monomer and several lipid molecules. The proximity of Car bridging between Car_D1 and plastoquinone/Q_A may imply a direct charge recombination of Car⁺Q_A⁻. The calculated E_m(Car_D1) and E_m(Chl_Z(D1)) are, respectively, 83 and 126 mV higher than E_m(Car_D2) and E_m(Chl_Z(D2)), which could explain why Car_D2⁺ and Chl_Z(D2)⁺ are observed rather than the corresponding Car_D1⁺ and Chl_Z(D1)⁺.     

Structurally conserved channels in cyanobacterial and plant photosystem II

In the cyanobacterial photosystem II (PSII), the O4-water chain in the D1 and CP43 proteins, a chain of water molecules that are directly H-bonded to O4 of the Mn₄Ca cluster, is linked with a channel that connects the protein bulk surface along with a membrane-extrinsic protein subunit, PsbU (O4-PsbU channel). The cyanobacterial PSII structure also shows that the O1 site of the Mn₄Ca cluster has a chain of H-bonded water molecules, which is linked with the channel that proceeds toward the bulk surface via PsbU and PsbV (O1-PsbU/V channel). Membrane-extrinsic protein subunits PsbU and PsbV in cyanobacterial PSII are replaced with PsbP and PsbQ in plant PSII. However, these four proteins have no structural similarity. It remains unknown whether the corresponding channels also exist in plant PSII, because water molecules are not identified in the plant PSII cryo-electron microscopy (cryo-EM) structure. Using the cyanobacterial and plant PSII structures, we analyzed the channels that proceed from the Mn₄Ca cluster. The cyanobacterial O4-PsbU and O1-PsbU/V channels were structurally conserved as the channel that proceeds along PsbP toward the protein bulk surface in the plant PSII (O4-PsbP and O1-PsbP channels, respectively). Calculated protonation states indicated that in contrast to the original geometry of the plant cryo-EM structure, protonated PsbP-Lys166 may form a salt-bridge with ionized D1-Glu329 and protonated PsbP-Lys173 may form a salt-bridge with ionized PsbQ-Asp28 near the O1-PsbP channel. The existence of these channels might explain the molecular mechanism of how PsbP can interact with the Mn₄Ca cluster.     

How does the QB site influence propagate to the QA site in photosystem II?

The redox potential of the primary quinone Q(A) [E(m)(Q(A))] in photosystem II (PSII) is lowered by replacement of the native plastoquinone (PQ) with bromoxynil (BR) at the secondary quinone Q(B) binding site. Using the BR-bound PSII structure presented in the previous Fourier transform infrared and docking calculation studies, we calculated E(m)(Q(A)) considering both the protein environment in atomic detail and the protonation pattern of the titratable residues. The calculated E(m)(Q(A)) shift in response to the replacement of PQ with deprotonated BR at the Q(B) binding site [ΔE(m)(Q(A))(PQ→BR)] was -55 mV when the three regions, Q(A), the non-heme iron complex, and Q(B) (Q(B) = PQ or BR), were treated as a conjugated supramolecule (Q(A)-Fe-Q(B)). The negative charge of BR apparently contributes to the downshift in ΔE(m)(Q(A))(PQ→BR). This downshift, however, is mostly offset by the influence of the residues near Q(B). The charge delocalization over the Q(A)-Fe-Q(B) complex and the resulting H-bond strength change between Q(A) and D2-His214 are crucial factors that yield a ΔE(m)(Q(A))(PQ→BR) of -55 mV by (i) altering the electrostatic influence of the H-bond donor D2-His214 on E(m)(Q(A)) and (ii) suppressing the proton uptake events of the titratable residues that could otherwise upshift ΔE(m)(Q(A))(PQ→BR) during replacement of PQ with BR at the Q(B) site.     

Function of redox-active tyrosine in photosystem II

Water oxidation at photosystem II Mn-cluster is mediated by the redox-active tyrosine Y(Z). We calculated the redox potential (E(m)) of Y(Z) and its symmetrical counterpart Y(D), by solving the linearized Poisson-Boltzmann equation. The calculated E(m)(Y( )/Y(-)) were +926 mV/+694 mV for Y(Z)/Y(D) with the Mn-cluster in S2 state. Together with the asymmetric position of the Mn-cluster relative to Y(Z/D), differences in H-bond network between Y(Z) (Y(Z)/D1-His(190)/D1-Asn(298)) and Y(D) (Y(D)/D2-His(189)/D2-Arg(294)/CP47-Glu(364)) are crucial for E(m)(Y(Z/D)). When D1-His(190) is protonated, corresponding to a thermally activated state, the calculated E(m)(Y(Z)) was +1216 mV, which is as high as the E(m) for P(D1/D2). We observed deprotonation at CP43-Arg(357) upon S-state transition, which may suggest its involvement in the proton exit pathway. E(m)(Y(D)) was affected by formation of P(D2)(+) (but not P(D1)(+)) and sensitive to the protonation state of D2-Arg(180). This points to an electrostatic link between Y(D) and P(D2).     

Function of two beta-carotenes near the D1 and D2 proteins in photosystem II dimers

The antenna proteins in photosystem II (PSII) not only promote energy transfer to the photosynthetic reaction center (RC) but provide also an efficient cation sink to re-reduce chlorophyll a if the electron transfer (ET) from the Mn-cluster is inhibited. Using the newest PSII dimer crystal structure (3.0 A resolution), in which 11 beta-carotene molecules (Car) and 14 lipids are visible in the PSII monomer, we calculated the redox potentials (Em) of one-electron oxidation for all Car (Em(Car)) by solving the Poisson-Boltzmann equation. In each PSII monomer, the D1 protein harbors a previously unlocated Car (CarD1) in van der Waals contact with the chlorin ring of ChlZ(D1). Each CarD1 in the PSII dimer complex is located in the interface between the D1 and CP47 subunits, together with another four Car of the other PSII monomer and several lipid molecules. The proximity of Car bridging between CarD1 and plastoquinone/Q(A) may imply a direct charge recombination of Car+Q(A)-. The calculated Em(CarD1) and Em(ChlZ(D1)) are, respectively, 83 and 126 mV higher than Em(CarD2) and Em(ChlZ(D2)), which could explain why CarD2+ and ChlZ(D2)+ are observed rather than the corresponding CarD1+ and ChlZ(D1)+.     

H atom positions and nuclear magnetic resonance chemical shifts of short H bonds in photoactive yellow protein

Recent neutron diffraction studies on photoactive yellow protein (PYP) proposed that the H bond between protonated Glu46 and the chromophore-ionized p-coumaric acid (pCA) is a low-barrier H bond (LBHB) mainly because the H atom position was assigned at the midpoint of the O(Glu46)-O(pCA) bond. However, the (1)H nuclear magnetic resonance (NMR) chemical shift (δ(H)) was 15.2 ppm, which is lower than the values of 17-19 ppm for typical LBHBs. We evaluated the dependence of δ(H) on an H atom position in the O(Glu46)-O(pCA) bond in the PYP ground state by using a quantum mechanical/molecular mechanical (QM/MM) approach. The calculated chemical shift unambiguously suggested that a δ(H) of 15.2 ppm for the O(Glu46)-O(pCA) bond in NMR studies should correspond to the QM/MM geometry (δ(H) = 14.5 ppm), where the H atom belongs to the Glu moiety, rather than the neutron diffraction geometry (δ(H) = 19.7 ppm), where the H atom is near the midpoint of the donor and acceptor atoms.