pK(a) of a Proton-Conducting Water Chain in Photosystem II

D1-Tyr246 and D2-Tyr244 in photosystem II: Insights into bicarbonate binding and electron transfer from QA•- to QB

In photosystem II (PSII), D1-Tyr246 and D2-Tyr244 are symmetrically located at the binding site of the bicarbonate ligand of the non-heme Fe complex. Here, we investigated the role of the symmetrically arranged tyrosine pair, D1-Tyr246 and D2-Tyr244, in the function of PSII, by generating four chloroplast mutants of PSII from Chlamydomonas reinhardtii: D1-Y246F, D1-Y246T, D2-Y244F, and D2-Y244T. The mutants exhibited altered photoautotrophic growth, reduced PSII protein accumulation, and impaired O2-evolving activity. Flash-induced fluorescence yield decay kinetics indicated a significant slowdown in electron transfer from QA•- to QB in all mutants. Bicarbonate reconstitution resulted in enhanced O2-evolving activity, suggesting destabilization of bicarbonate binding in the mutants. Structural analyses based on a quantum mechanical/molecular mechanical approach identified the existence of a water channel that leads to incorporation of bulk water molecules and destabilization of the bicarbonate binding site. The water intake channels, crucial for bicarbonate stability, exhibited distinct paths in the mutants. These findings shed light on the essential role of the tyrosine pair in maintaining bicarbonate stability and facilitating efficient electron transfer in native PSII.

The Effect of Removal of External Proteins PsbO, PsbP and PsbQ on Flash-Induced Molecular Oxygen Evolution and Its Biphasicity in Tobacco PSII

The oxygen evolution within photosystem II (PSII) is one of the most enigmatic processes occurring in nature. It is suggested that external proteins surrounding the oxygen-evolving complex (OEC) not only stabilize it and provide an appropriate ionic environment but also create water channels, which could be involved in triggering the ingress of water and the removal of O2 and protons outside the system. To investigate the influence of these proteins on the rate of oxygen release and the efficiency of OEC function, we developed a measurement protocol for the direct measurement of the kinetics of oxygen release from PSII using a Joliot-type electrode. PSII-enriched tobacco thylakoids were used in the experiments. The results revealed the existence of slow and fast modes of oxygen evolution. This observation is model-independent and requires no specific assumptions about the initial distribution of the OEC states. The gradual removal of exogenous proteins resulted in a slowdown of the rapid phase (~ms) of O2 release and its gradual disappearance while the slow phase (~tens of ms) accelerated. The role of external proteins in regulating the biphasicity and efficiency of oxygen release is discussed based on observed phenomena and current knowledge.

Structure Function Studies of Photosystem II Using X-Ray Free Electron Lasers

The structure and mechanism of the water-oxidation chemistry that occurs in photosystem II have been subjects of great interest. The advent of X-ray free electron lasers allowed the determination of structures of the stable intermediate states and of steps in the transitions between these intermediate states, bringing a new perspective to this field. The room-temperature structures collected as the photosynthetic water oxidation reaction proceeds in real time have provided important novel insights into the structural changes and the mechanism of the water oxidation reaction. The time-resolved measurements have also given us a view of how this reaction-which involves multielectron, multiproton processes-is facilitated by the interaction of the ligands and the protein residues in the oxygen-evolving complex. These structures have also provided a picture of the dynamics occurring in the channels within photosystem II that are involved in the transport of the substrate water to the catalytic center and protons to the bulk.

Energetics of the H-Bond Network in Exiguobacterium sibiricum Rhodopsin

Exiguobacterium sibiricum rhodopsin (ESR) functions as a light-driven proton pump utilizing Lys96 for proton uptake and maintaining its activity over a wide pH range. Using a combination of methodologies including the linear Poisson-Boltzmann equation and a quantum mechanical/molecular mechanical approach with a polarizable continuum model, we explore the microscopic mechanisms underlying its pumping activity. Lys96, the primary proton uptake site, remains deprotonated owing to the loss of solvation in the ESR protein environment. Asp85, serving as a proton acceptor group for Lys96, does not form a low-barrier H-bond with His57. Instead, deprotonated Asp85 forms a salt-bridge with protonated His57, and the proton is predominantly located at the His57 moiety. Glu214, the only acidic residue at the end of the H-bond network exhibits a pKa value of ∼6, slightly elevated due to solvation loss. It seems likely that the H-bond network [Asp85···His57···H2O···Glu214] serves as a proton-conducting pathway toward the protein bulk surface.

Exploring the Deprotonation Process during Incorporation of a Ligand Water Molecule at the Dangling Mn Site in Photosystem II

The Mn4CaO5 cluster, featuring four ligand water molecules (W1 to W4), serves as the water-splitting site in photosystem II (PSII). X-ray free electron laser (XFEL) structures exhibit an additional oxygen site (O6) adjacent to the O5 site in the fourth lowest oxidation state, S3, forming Mn4CaO6. Here, we investigate the mechanism of the second water ligand molecule at the dangling Mn (W2) as a potential incorporating species, using a quantum mechanical/molecular mechanical (QM/MM) approach. Previous QM/MM calculations demonstrated that W1 releases two protons through a low-barrier H-bond toward D1-Asp61 and subsequently releases an electron during the S2 to S3 transition, resulting in O•- at W1 and protonated D1-Asp61. During the process of Mn4CaO6 formation, O•-, rather than H2O or OH-, best reproduced the O5···O6 distance. Although the catalytic cluster with O•- at O6 is more stable than that with O•- at W1 in S3, it does not occur spontaneously due to the significantly uphill deprotonation process. Assuming O•- at W2 incorporates into the O6 site, an exergonic conversion from Mn1(III)Mn2(IV)Mn3(IV)Mn4(IV) (equivalent to the open-cubane S2 valence state) to Mn1(IV)Mn2(IV)Mn3(IV)Mn4(III) (equivalent to the closed-cubane S2 valence state) occurs. These findings provide energetic insights into the deprotonation and structural conversion events required for the Mn4CaO6 formation.

Superprotonic Conductivity in a Metalloporphyrin-Based SMOF (Supramolecular Metal-Organic Framework)

Metal-organic frameworks and supramolecular metal-organic frameworks (SMOFs) exhibit great potential for a broad range of applications taking advantage of the high surface area and pore sizes and tunable chemistry. In particular, metalloporphyrin-based MOFs and SMOFs are becoming of great importance in many fields due to the bioessential functions of these macrocycles that are being mimicked. On the other hand, during the last years, proton-conducting materials have aroused much interest, and those presenting high conductivity values are potential candidates to play a key role in some solid-state electrochemical devices such as batteries and fuel cells. In this way, using metalloporphyrins as building units we have obtained a new crystalline material with formula [H(bipy)]2[(MnTPPS)(H2O)2]·2bipy·14H2O, where bipy is 4,4'-bipyidine and TPPS4- is the meso-tetra(4-sulfonatephenyl) porphyrin. The crystal structure shows a zig-zag water chain along the [100] direction located between the sulfonate groups of the porphyrin. Taking into account those structural features, the compound was tested for proton conduction by complex electrochemical impedance spectroscopy (EIS). The as-obtained conductivity is 1 × 10-2 S·cm-1 at 40 °C and 98% relative humidity, which is a remarkably high value.

Mapping the Oxygens in the Oxygen-Evolving Complex of Photosystem II by Their Nucleophilicity Using Quantum Descriptors

The oxygen-evolving complex (OEC) of Photosystem II catalyzes the water-splitting reaction using solar energy. Thus, understanding the reaction mechanism will inspire the design of biomimetic artificial catalysts that convert solar energy to chemical energy. Conceptual Density Functional Theory (CDFT) focuses on understanding the reactivity of molecules and the atomic contribution to the overall nucleophilicity and electrophilicity of the molecule using quantum descriptors. However, this method has not been applied to the OEC before. Here, we use Fukui functions and the dual descriptor to provide quantitative measures of the nucleophilicity and electrophilicity of oxygens in the OEC for different models in different S states. Our results show that the μ-oxo bridges connected to terminal Mn4 are nucleophilic, and those in the cube formed by Mn1, Mn2, and Mn3 are mostly electrophilic. The dual descriptors of the bridging oxygens in the OEC showed a similar reactivity to that of bridging oxygens in Mn model compounds. However, the terminal water W1, which is bound to Mn4, showed very strong reactivity in some of the S3 models. Thus, our calculations support the model that proposes the formation of the O2 molecule through nucleophilic attack by a terminal water.

Interplay of two low-barrier hydrogen bonds in long-distance proton-coupled electron transfer for water oxidation

D1-Tyr161 (TyrZ) forms a low-barrier H-bond with D1-His190 and functions as a redox-active group in photosystem II. When oxidized to the radical form (TyrZ-O•), it accepts an electron from the oxygen-evolving Mn4CaO5 cluster, facilitating an increase in the oxidation state (Sn; n = 0-3). In this study, we investigated the mechanism of how TyrZ-O• drives proton-coupled electron transfer during the S2 to S3 transition using a quantum mechanical/molecular mechanical approach. In response to TyrZ-O• formation and subsequent loss of the low-barrier H-bond, the ligand water molecule at the Ca2+ site (W4) reorients away from TyrZ and donates an H-bond to D1-Glu189 at Mn4 of Mn4CaO5 together with an adjacent water molecule. The H-bond donation to the Mn4CaO5 cluster triggers the release of the proton from the lowest pKa site (W1 at Mn4) along the W1…D1-Asp61 low-barrier H-bond, leading to protonation of D1-Asp61. The interplay of the two low-barrier H-bonds, involving the Ca2+ interface and forming the extended Grotthuss-like network [TyrZ…D1-His190]-[Mn4CaO5]-[W1…D1-Asp61], rather than the direct electrostatic interaction, is likely a basis of the apparent long-distance interaction (11.4 Å) between TyrZ-O• formation and D1-Asp61 protonation.

How chloride functions to enable proton conduction in photosynthetic water oxidation: Time-resolved kinetics of intermediates (S-states) in vivo and bromide substitution

Chloride (Cl-) is essential for O2 evolution during photosynthetic water oxidation. Two chlorides near the water-oxidizing complex (WOC) in Photosystem II (PSII) structures from Thermosynechococcus elongatus (and T. vulcanus) have been postulated to transfer protons generated from water oxidation. We monitored four criteria: primary charge separation flash yield (P* → P+QA-), rates of water oxidation steps (S-states), rate of proton evolution, and flash O2 yield oscillations by measuring chlorophyll variable fluorescence (P* quenching), pH-sensitive dye changes, and oximetry. Br-substitution slows and destabilizes cellular growth, resulting from lower light-saturated O2 evolution rate (-20 %) and proton release (-36 % ΔpH gradient). The latter implies less ATP production. In Br- cultures, protonogenic S-state transitions (S2 → S3 → S0') slow with increasing light intensity and during O2/water exchange (S0' → S0 → S1), while the non-protonogenic S1 → S2 transition is kinetically unaffected. As flash rate increases in Cl- cultures, both rate and extent of acidification of the lumen increase, while charge recombination is suppressed relative to Br-. The Cl- advantage in rapid proton escape from the WOC to lumen is attributed to correlated ion-pair movement of H3O+Cl- in dry water channels vs. separated Br- and H+ ion movement through different regions (>200-fold difference in Bronsted acidities). By contrast, at low flash rates a previously unreported reversal occurs that favors Br- cultures for both proton evolution and less PSII charge recombination. In Br- cultures, slower proton transfer rate is attributed to stronger ion-pairing of Br- with AA residues lining the water channels. Both anions charge-neutralize protons and shepherd them to the lumen using dry aqueous channels.

Evolutionary diversity of proton and water channels on the oxidizing side of photosystem II and their relevance to function

One of the reasons for the high efficiency and selectivity of biological catalysts arise from their ability to control the pathways of substrates and products using protein channels, and by modulating the transport in the channels using the interaction with the protein residues and the water/hydrogen-bonding network. This process is clearly demonstrated in Photosystem II (PS II), where its light-driven water oxidation reaction catalyzed by the Mn4CaO5 cluster occurs deep inside the protein complex and thus requires the transport of two water molecules to and four protons from the metal center to the bulk water. Based on the recent advances in structural studies of PS II from X-ray crystallography and cryo-electron microscopy, in this review we compare the channels that have been proposed to facilitate this mass transport in cyanobacteria, red and green algae, diatoms, and higher plants. The three major channels (O1, O4, and Cl1 channels) are present in all species investigated; however, some differences exist in the reported structures that arise from the different composition and arrangement of membrane extrinsic subunits between the species. Among the three channels, the Cl1 channel, including the proton gate, is the most conserved among all photosynthetic species. We also found at least one branch for the O1 channel in all organisms, extending all the way from Ca/O1 via the 'water wheel' to the lumen. However, the extending path after the water wheel varies between most species. The O4 channel is, like the Cl1 channel, highly conserved among all species while having different orientations at the end of the path near the bulk. The comparison suggests that the previously proposed functionality of the channels in T. vestitus (Ibrahim et al., Proc Natl Acad Sci USA 117:12624-12635, 2020; Hussein et al., Nat Commun 12:6531, 2021) is conserved through the species, i.e. the O1-like channel is used for substrate water intake, and the tighter Cl1 and O4 channels for proton release. The comparison does not eliminate the potential role of O4 channel as a water intake channel. However, the highly ordered hydrogen-bonded water wire connected to the Mn4CaO5 cluster via the O4 may strongly suggest that it functions in proton release, especially during the S0 → S1 transition (Saito et al., Nat Commun 6:8488, 2015; Kern et al., Nature 563:421-425, 2018; Ibrahim et al., Proc Natl Acad Sci USA 117:12624-12635, 2020; Sakashita et al., Phys Chem Chem Phys 22:15831-15841, 2020; Hussein et al., Nat Commun 12:6531, 2021).

New insights into the proton pumping mechanism of ba3 cytochrome c oxidase: the functions of key residues and water

As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.

Structural and energetic insights into Mn-to-Fe substitution in the oxygen-evolving complex

Manganese (Mn) serves as the catalytic center for water splitting in photosystem II (PSII), despite the abundance of iron (Fe) on earth. As a first step toward why Mn and not Fe is employed by Nature in the water oxidation catalyst, we investigated the Fe₄CaO₅ cluster in the PSII protein environment using a quantum mechanical/molecular mechanical (QM/MM) approach, assuming an equivalence between Mn(III/IV) and Fe(II/III). Substituting Mn with Fe resulted in the protonation of μ-oxo bridges at sites O2 and O3 by Arg357 and D1-His337, respectively. While the Mn₄CaO₅ cluster exhibits distinct open- and closed-cubane S₂ conformations, the Fe₄CaO₅ cluster lacks this variability due to an equal spin distribution over sites Fe1 and Fe4. The absence of a low-barrier H-bond between a ligand water molecule (W1) and D1-Asp61 in the Fe₄CaO₅ cluster may underlie its incapability for ligand water deprotonation, highlighting the relevance of Mn in natural water splitting.

Dioxygen and glucose force motion of the electron-transfer switch in the iron(III) flavohemoglobin-type nitric oxide dioxygenase

Kinetic and structural investigations of the flavohemoglobin-type NO dioxygenase have suggested critical roles for transient Fe(III)O2 complex formation and O2-forced movements affecting hydride transfer to the FAD cofactor and electron-transfer to the Fe(III)O2 complex. Stark-effect theory together with structural models and dipole and internal electrostatic field determinations provided a semi-quantitative spectroscopic method for investigating the proposed Fe(III)O2 complex and O2-forced movements. Deoxygenation of the enzyme causes Stark effects on the ferric heme Soret and charge-transfer bands revealing the Fe(III)O2 complex. Deoxygenation also elicits Stark effects on the FAD that expose forces and motions that create a more restricted NADH access to FAD for hydride transfer and switch electron-transfer off. Glucose also forces the enzyme toward an off state. Amino acid substitutions at the B10, E7, E11, G8, D5, and F7 positions influence the Stark effects of O2 on resting heme spin states and FAD consistent with the proposed roles of the side chains in the enzyme mechanism. Deoxygenation of ferric myoglobin and hemoglobin A also induces Stark effects on the hemes suggesting a common 'oxy-met' state. The ferric myoglobin and hemoglobin heme spectra are also glucose-responsive. A conserved glucose or glucose-6-phosphate binding site is found bridging the BC-corner and G-helix in flavohemoglobin and myoglobin suggesting novel allosteric effector roles for glucose or glucose-6-phosphate in the NO dioxygenase and O2 storage functions. The results support the proposed roles of a ferric O2 intermediate and protein motions in regulating electron-transfer during NO dioxygenase turnover.

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.

Energetics and proton release in photosystem II from Thermosynechococcus elongatus with a D1 protein encoded by either the psbA2 or psbA3 gene

In the cyanobacterium Thermosynechococcus elongatus, there are three psbA genes coding for the Photosystem II (PSII) D1 subunit that interacts with most of the main cofactors involved in the electron transfers. Recently, the 3D crystal structures of both PsbA2-PSII and PsbA3-PSII have been solved [Nakajima et al., J. Biol. Chem. 298 (2022) 102668.]. It was proposed that the loss of one hydrogen bond of Phe_D1 due to the D1-Y147F exchange in PsbA2-PSII resulted in a more negative E_m of Phe_D1 in PsbA2-PSII when compared to PsbA3-PSII. In addition, the loss of two water molecules in the Cl-1 channel was attributed to the D1-P173M substitution in PsbA2-PSII. This exchange, by narrowing the Cl-1 proton channel, could be at the origin of a slowing down of the proton release. Here, we have continued the characterization of PsbA2-PSII by measuring the thermoluminescence from the S₂Q_A^-/DCMU charge recombination and by measuring proton release kinetics using time-resolved absorption changes of the dye bromocresol purple. It was found that i) the E_m of Phe_D1^-•/Phe_D1 was decreased by ~30 mV in PsbA2-PSII when compared to PsbA3-PSII and ii) the kinetics of the proton release into the bulk was significantly slowed down in PsbA2-PSII in the S₂Tyr_Z^• to S₃Tyr_Z and S₃Tyr_Z^• → (S₃Tyr_Z^•)' transitions. This slowing down was partially reversed by the PsbA2/M173P mutation and induced by the PsbA3/P173M mutation thus confirming a role of the D1-173 residue in the egress of protons trough the Cl-1 channel.

Identification of the protonation and oxidation states of the oxygen-evolving complex in the low-dose X-ray crystal structure of photosystem II

In photosystem II (PSII), the O3 and O4 sites of the Mn₄CaO₅ cluster form hydrogen bonds with D1-His337 and a water molecule (W539), respectively. The low-dose X-ray structure shows that these hydrogen bond distances differ between the two homogeneous monomer units (A and B) [Tanaka et al., J. Am Chem. Soc. 2017, 139, 1718]. We investigated the origin of the differences using a quantum mechanical/molecular mechanical (QM/MM) approach. QM/MM calculations show that the short O4-O_W539 hydrogen bond (~2.5 Å) of the B monomer is reproduced when O4 is protonated in the S₁ state. The short O3-Nε_His337 hydrogen bond of the A monomer is due to the formation of a low-barrier hydrogen bond between O3 and doubly-protonated D1-His337 in the overreduced states (S_-1 or S_-2). It seems plausible that the oxidation state differs between the two monomer units in the crystal.

Substitution of Ca2+ and changes in the H-bond network near the oxygen-evolving complex of photosystem II

Ca²⁺, which provides binding sites for ligand water molecules W3 and W4 in the Mn₄CaO₅ cluster, is a prerequisite for O₂ evolution in photosystem II (PSII). We report structural changes in the H-bond network and the catalytic cluster itself upon the replacement of Ca²⁺ with other alkaline earth metals, using a quantum mechanical/molecular mechanical approach. The small radius of Mg²⁺ makes W3 donate an H-bond to D1-Glu189 in Mg²⁺-PSII. If an additional water molecule binds at the large surface of Ba²⁺, it donates H-bonds to D1-Glu189 and the ligand water molecule at the dangling Mn, altering the H-bond network. The potential energy profiles of the H-bond between D1-Tyr161 (TyrZ) and D1-His190 and the interconversion between the open- and closed-cubane S₂ conformations remain substantially unaltered upon the replacement of Ca²⁺. Remarkably, the O5⋯Ca²⁺ distance is shortest among all O5⋯metal distances irrespective of the radius being larger than that of Mg²⁺. Furthermore, Ca²⁺ is the only alkaline earth metal that equalizes the O5⋯metal and O2⋯metal distances and facilitates the formation of the symmetric cubane structure.

Electron-Transfer Route in the Early Oxidation States of the Mn4CaO5 Cluster in Photosystem II

The electron transfer from the oxygen-evolving Mn4CaO5 cluster to the electron acceptor D1-Tyr161 (TyrZ) is a prerequisite for water oxidation and O2 evolution. Here, we analyzed the electronic coupling in the rate-limiting electron-transfer transitions using a combined quantum mechanical/molecular mechanical/polarizable continuum model approach. In the S0 to S1 transition, the electronic coupling between the electron-donor Mn3(III) and TyrZ is small (2 meV). In contrast, the electronic coupling between the dangling Mn4(III) and TyrZ is significantly large (172 meV), which suggests that the electron transfer proceeds from Mn3(III) to TyrZ via Mn4(III). In the S1 to S2 transition, the electronic coupling between Mn4(III) and TyrZ is also larger (124 meV) than that between Mn1(III) and TyrZ (1 meV), which favors the formation of the open-cubane S2 conformation with Mn4(IV) over the formation of the closed-cubane S2 conformation with Mn1(IV). In the S0 to S1 and S1 to S2 transitions, the Mn4 d-orbital and the TyrZ π-orbital are hybridized via D1-Asp170, which suggests that D1-Asp170 commonly provides a dominant electron-transfer route.

Functional Water Networks in Fully Hydrated Photosystem II

Water channels and networks within photosystem II (PSII) of oxygenic photosynthesis are critical for enzyme structure and function. They control substrate delivery to the oxygen-evolving center and mediate proton transfer at both the oxidative and reductive endpoints. Current views on PSII hydration are derived from protein crystallography, but structural information may be compromised by sample dehydration and technical limitations. Here, we simulate the physiological hydration structure of a cyanobacterial PSII model following a thorough hydration procedure and large-scale unconstrained all-atom molecular dynamics enabled by massively parallel simulations. We show that crystallographic models of PSII are moderately to severely dehydrated and that this problem is particularly acute for models derived from X-ray free electron laser (XFEL) serial femtosecond crystallography. We present a fully hydrated representation of cyanobacterial PSII and map all water channels, both static and dynamic, associated with the electron donor and acceptor sides. Among them, we describe a series of transient channels and the attendant conformational gating role of protein components. On the acceptor side, we characterize a channel system that is absent from existing crystallographic models but is likely functionally important for the reduction of the terminal electron acceptor plastoquinone Q_B. The results of the present work build a foundation for properly (re)evaluating crystallographic models and for eliciting new insights into PSII structure and function.

Protonation structure of the closed-cubane conformation of the O2-evolving complex in photosystem II

In photosystem II (PSII), one-electron oxidation of the most stable state of the oxygen-evolving Mn4CaO5 cluster (S1) leads to the S2 state formation, Mn1(III)Mn2(IV)Mn3(IV)Mn4(IV) (open-cubane S2) or Mn1(IV)Mn2(IV)Mn3(IV)Mn4(III) (closed-cubane S2). In electron paramagnetic resonance (EPR) spectroscopy, the g = 4.1 signal is not observed in cyanobacterial PSII but in plant PSII, whereas the g = 4.8 signal is observed in cyanobacterial PSII and extrinsic-subunit-depleted plant PSII. Here, we investigated the closed-cubane S2 conformation, a candidate for a higher spin configuration that accounts for g > 4.1 EPR signal, considering all pairwise exchange couplings in the PSII protein environment (i.e. instead of considering only a single exchange coupling between the [Mn3(CaO4)] cubane region and the dangling Mn4 site). Only when a ligand water molecule that forms an H-bond with D1-Asp61 (W1) is deprotonated at dangling Mn4(IV), the g = 4.1 EPR spectra can be reproduced using the cyanobacterial PSII crystal structure. The closed-cubane S2 is less stable than the open-cubane S2 in cyanobacterial PSII, which may explain why the g = 4.1 EPR signal is absent in cyanobacterial PSII.

Role of D1-Glu65 in Proton Transfer during Photosynthetic Water Oxidation in Photosystem II

Photosynthetic water oxidation takes place at the Mn₄CaO₅ cluster in photosystem II (PSII) through a light-driven cycle of five intermediates called S states (S₀-S₄). Although the PSII structures have shown the presence of several channels around the Mn₄CaO₅ cluster leading to the lumen, the pathways for proton release in the individual S-state transitions remain unidentified. Here, we studied the involvement of the so-called Cl channel in proton transfer during water oxidation by examining the effect of the mutation of D1-Glu65, a key residue in this channel, to Ala using Fourier transform infrared difference and time-resolved infrared spectroscopies together with thermoluminescence and delayed luminescence measurements. It was shown that the structure and the redox property of the catalytic site were little affected by the D1-Glu65Ala mutation. In the S₂ → S₃ transition, the efficiency was still high and the transition rate was only moderately retarded in the D1-Glu65Ala mutant. In contrast, the S₃ → S₀ transition was significantly inhibited by this mutation. These results suggest that proton transfer in the S₂ → S₃ transition occurs through multiple pathways including the Cl channel, whereas this channel likely serves as a single pathway for proton exit in the S₃ → S₀ transition.

Binding and functions of the two chloride ions in the oxygen-evolving center of photosystem II

Light-driven water oxidation in photosynthesis occurs at the oxygen-evolving center (OEC) of photosystem II (PSII). Chloride ions (Cl^-) are essential for oxygen evolution by PSII, and two Cl^- ions have been found to specifically bind near the Mn₄CaO₅ cluster in the OEC. The retention of these Cl^- ions within the OEC is critically supported by some of the membrane-extrinsic subunits of PSII. The functions of these two Cl^- ions and the mechanisms of their retention both remain to be fully elucidated. However, intensive studies performed recently have advanced our understanding of the functions of these Cl^- ions, and PSII structures from various species have been reported, aiding the interpretation of previous findings regarding Cl^- retention by extrinsic subunits. In this review, we summarize the findings to date on the roles of the two Cl^- ions bound within the OEC. Additionally, together with a short summary of the functions of PSII membrane-extrinsic subunits, we discuss the mechanisms of Cl^- retention by these extrinsic subunits.

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.

D139N mutation of PsbP enhances the oxygen-evolving activity of photosystem II through stabilized binding of a chloride ion

Photosystem II (PSII) is a multisubunit membrane protein complex that catalyzes light-driven oxidation of water to molecular oxygen. The chloride ion (Cl-) has long been known as an essential cofactor for oxygen evolution by PSII, and two Cl- ions (Cl-1 and Cl-2) have been found to specifically bind near the Mn4CaO5 cluster within the oxygen-evolving center (OEC). However, despite intensive studies on these Cl- ions, little is known about the function of Cl-2, the Cl- ion that is associated with the backbone nitrogens of D1-Asn338, D1-Phe339, and CP43-Glu354. In green plant PSII, the membrane extrinsic subunits-PsbP and PsbQ-are responsible for Cl- retention within the OEC. The Loop 4 region of PsbP, consisting of highly conserved residues Thr135-Gly142, is inserted close to Cl-2, but its importance has not been examined to date. Here, we investigated the importance of PsbP-Loop 4 using spinach PSII membranes reconstituted with spinach PsbP proteins harboring mutations in this region. Mutations in PsbP-Loop 4 had remarkable effects on the rate of oxygen evolution by PSII. Moreover, we found that a specific mutation, PsbP-D139N, significantly enhances the oxygen-evolving activity in the absence of PsbQ, but not significantly in its presence. The D139N mutation increased the Cl- retention ability of PsbP and induced a unique structural change in the OEC, as indicated by light-induced Fourier transform infrared (FTIR) difference spectroscopy and theoretical calculations. Our findings provide insight into the functional significance of Cl-2 in the water-oxidizing reaction of PSII.

Release of a Proton and Formation of a Low-Barrier Hydrogen Bond between Tyrosine D and D2-His189 in Photosystem II

In photosystem II (PSII), the second-lowest oxidation state (S₁) of the oxygen-evolving Mn₄CaO₅ cluster is the most stable, as the radical form of the redox-active D2-Tyr160 is considered to be a candidate that accepts an electron from the lowest oxidation state (S₀) in the dark. Using quantum mechanical/molecular mechanical calculations, we investigated the redox potential (E _m) of TyrD and its H-bond partner, D2-His189. The potential energy profile indicates that the release of a proton from the TyrD...D2-His189 pair leads to the formation of a low-barrier H-bond. The E _m depends on the H⁺ position along the low-barrier H-bond, e.g., 680 mV when the H⁺ is at the D2-His189 moiety and 800 mV when the H⁺ is at the TyrD moiety, which can explain why TyrD mediates both the S₀ to S₁ oxidation and the S₂ to S₁ reduction.

Glycerol binding at the narrow channel of photosystem II stabilizes the low-spin S2 state of the oxygen-evolving complex

The oxygen-evolving complex (OEC) of photosystem II (PSII) cycles through redox intermediate states Si (i = 0-4) during the photochemical oxidation of water. The S2 state involves an equilibrium of two isomers including the low-spin S2 (LS-S2) state with its characteristic electron paramagnetic resonance (EPR) multiline signal centered at g = 2.0, and a high-spin S2 (HS-S2) state with its g = 4.1 EPR signal. The relative intensities of the two EPR signals change under experimental conditions that shift the HS-S2/LS-S2 state equilibrium. Here, we analyze the effect of glycerol on the relative stability of the LS-S2 and HS-S2 states when bound at the narrow channel of PSII, as reported in an X-ray crystal structure of cyanobacterial PSII. Our quantum mechanics/molecular mechanics (QM/MM) hybrid models of cyanobacterial PSII show that the glycerol molecule perturbs the hydrogen-bond network in the narrow channel, increasing the pKa of D1-Asp61 and stabilizing the LS-S2 state relative to the HS-S2 state. The reported results are consistent with the absence of the HS-S2 state EPR signal in native cyanobacterial PSII EPR spectra and suggest that the narrow water channel hydrogen-bond network regulates the relative stability of OEC catalytic intermediates during water oxidation.

Requirement of Chloride for the Downhill Electron Transfer Pathway from the Water-Splitting Center in Natural Photosynthesis

In photosystem II (PSII), Cl^- is a prerequisite for the second flash-induced oxidation of the Mn₄CaO₅ cluster (the S₂ to S₃ transition). We report proton transfer from the substrate water molecule via D1-Asp61 and electron transfer via redox-active D1-Tyr161 (TyrZ) to the chlorophyll pair in Cl^--depleted PSII using a quantum mechanical/molecular mechanical approach. The low-barrier H-bond formation between the substrate water molecule and D1-Asp61 remained unaffected upon the depletion of Cl^-. However, the binding site, D2-Lys317, formed a salt bridge with D1-Asp61, leading to the inhibition of the subsequent proton transfer. Remarkably, the redox potential (E_m) of S₂/S₃ increased significantly, making electron transfer from S₂ to TyrZ energetically uphill, as observed in Ca²⁺-depleted PSII. The uphill electron transfer pathway was induced by the significant increase in E_m(S₂/S₃) caused by the loss of charge compensation for D2-Lys317 upon the depletion of Cl^-, whereas it was induced by the significant decrease in E_m(TyrZ) caused by the rearrangement of the water molecules at the Ca²⁺ binding moiety upon the depletion of Ca²⁺.

C-Graphs Tool with Graphical User Interface to Dissect Conserved Hydrogen-Bond Networks: Applications to Visual Rhodopsins

Dynamic hydrogen-bond networks provide proteins with structural plasticity required to translate signals such as ligand binding into a cellular response or to transport ions and larger solutes across membranes and, thus, are of central interest to understand protein reaction mechanisms. Here, we present C-Graphs, an efficient tool with graphical user interface that analyzes data sets of static protein structures or of independent numerical simulations to identify conserved, vs unique, hydrogen bonds and hydrogen-bond networks. For static structures, which may belong to the same protein or to proteins with different sequences, C-Graphs uses a clustering algorithm to identify sites of the hydrogen-bond network where waters are conserved among the structures. Using C-Graphs, we identify an internal protein-water hydrogen-bond network common to static structures of visual rhodopsins and adenosine A2A G protein-coupled receptors (GPCRs). Molecular dynamics simulations of a visual rhodopsin indicate that the conserved hydrogen-bond network from static structure can recruit dynamic hydrogen bonds and extend throughout most of the receptor. We release with this work the code for C-Graphs and its graphical user interface.

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

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

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 exit pathways surrounding the oxygen evolving complex of photosystem II

Photosystem II allows water to be the primary electron source for the photosynthetic electron transfer chain. Water is oxidized to dioxygen at the Oxygen Evolving Complex (OEC), a Mn₄CaO₅ inorganic core embedded on the lumenal side of PSII. Water-filled channels surrounding the OEC must bring in substrate water molecules, remove the product protons to the lumen, and may transport the product oxygen. Three water-filled channels, denoted large, narrow, and broad, extend from the OEC towards the aqueous surface more than 15 Å away. However, the role of each pathway in the transport in and out of the OEC is yet to be established. Here, we combine Molecular Dynamics (MD), Multi Conformation Continuum Electrostatics (MCCE) and Network Analysis to compare and contrast the three potential proton transfer paths. Hydrogen bond network analysis shows that near the OEC the waters are highly interconnected with similar free energy for hydronium at all locations. The paths diverge as they move towards the lumen. The water chain in the broad channel is better connected than in the narrow and large channels, where disruptions in the network are observed approximately 10 Å from the OEC. In addition, the barrier for hydronium translocation is lower in the broad channel. Thus, a proton released from any location on the OEC can access all paths, but the likely exit to the lumen passes through PsbO via the broad channel.

Two Distinct Oxygen-Radical Conformations in the X-ray Free Electron Laser Structures of Photosystem II

We report the existence of two distinct oxygen-radical-containing Mn₄CaO_5/6 conformations with short O···O bonds in the crystal structures of the oxygen-evolving enzyme photosystem II (PSII), obtained using an X-ray free electron laser (XFEL). A short O···O distance of <2.3 Å between the O4 site of the Mn₄CaO₅ complex and the adjacent water molecule (W539) in the proton-conducting O4-water chain was observed in the second flash-induced (2F) XFEL structure (2F-XFEL), which may correspond to S₃. By use of a quantum mechanical/molecular mechanical approach, the OH^• formation at W539 and the short O4···O_W539 distance (<2.3 Å) were reproduced in S₂ and S₃ with reduced Mn1(III), which lacks the additional sixth water molecule O6. As the O^•- formation at O6 and the short O5···O6 distance (1.9 Å) have been reported in another 2F-XFEL structure with reduced Mn4(III), two distinct oxygen-radical conformations exist in the 2F-XFEL crystals.

Rigidly hydrogen-bonded water molecules facilitate proton transfer in photosystem II

In the water-splitting enzyme photosystem II (PSII), the proton is released from the catalytic site and transferred to the protein bulk surface via the proton-relay mechanism. Proton transfer occurs in a proton-conducting channel consisting of a series of water molecules connected by hydrogen-bonded (H-bonded) chains. The water-transport protein aquaporin (AQP) also contains a water chain with structure similar to that of the PSII proton channel, although the water chain does not transport protons. We compared the PSII proton channel with the AQP water channel from the following standpoints: (1) the energetics of proton transfer based on crystal structures obtained from quantum mechanical/molecular mechanical calculations, and (2) fluctuations in water molecules obtained from molecular dynamics simulations. The results showed that residues facing the channel and acting as H-bonded partners of water molecules predominantly determined the proton-transfer ability. In PSII, the water chain is surrounded by H-bond acceptors (e.g., carbonyl groups), and the water chain transports protons where the water molecules are rigidly fixed. In AQP, the water chain is surrounded by hydrophobic sidechains or H-bond donors (e.g., NH2 groups), and it does not transport protons where the water molecules are flexible and fluctuating.

The Nature of the Short Oxygen-Oxygen Distance in the Mn4CaO6 Complex of Photosystem II Crystals

The O···O distance for a typical H-bond is ∼2.8 Å, whereas the radiation-damage-free structures of photosystem II (PSII), obtained using the X-ray free electron laser (XFEL), shows remarkably short O···O distances of ∼2 Å in the oxygen-evolving Mn₄CaO_5/6 complex. Herein, we report the protonation/oxidation states of the short O···O atoms in the XFEL structures using a quantum mechanical/molecular mechanical approach. The O5···O6 distance of 1.9 Å is reproduced only when O6 is an unprotonated O radical (O^•-) with Mn(IV)₃Mn(III), i.e., the S₃ state. The potential energy profile shows a barrier-less energy minimum region when O5···O6 = 1.90-2.05 Å (O^•- ↓) or 2.05-2.20 Å (O^•- ↑). Formation of such a short O5···O6 distance is not possible when O6 is OH^- with Mn(IV)₄. In the case in which the O5···O6 distance is 1.9 Å, it seems likely that the O radical species exists in the oxygen-evolving complex of the XFEL-S₃ crystals.

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.

pKa of the ligand water molecules in the oxygen-evolving Mn4CaO5 cluster in photosystem II

Release of the protons from the substrate water molecules is prerequisite for O₂ evolution in photosystem II (PSII). Proton-releasing water molecules with low pK_a values at the catalytic moiety can be the substrate water molecules. In some studies, one of the ligand water molecules, W2, is regarded as OH^-. However, the PSII crystal structure shows neither proton acceptor nor proton-transfer pathway for W2, which is not consistent with the assumption of W2 = OH^-. Here we report the pK_a values of the four ligand water molecules, W1 and W2 at Mn4 and W3 and W4 at Ca²⁺, of the Mn₄CaO₅ cluster. pK_a(W1) ≈ pK_a(W2) << pK_a(W3) ≈ pK_a(W4) in the Mn₄CaO₅ cluster in water. However, pK_a(W1) ≈ pK_a(D1-Asp61) << pK_a(W2) in the PSII protein environment. These results suggest that in PSII, deprotonation of W2 is energetically disfavored as far as W1 exists.

Redox Potential of the Oxygen-Evolving Complex in the Electron Transfer Cascade of Photosystem II

In photosystem II (PSII), water oxidation occurs in the Mn₄CaO₅ cluster with the release of electrons via the redox-active tyrosine (TyrZ) to the reaction-center chlorophylls (P_D1/P_D2). Using a quantum mechanical/molecular mechanical approach, we report the redox potentials (E_m) of these cofactors in the PSII protein environment. The E_m values suggest that the Mn₄CaO₅ cluster, TyrZ, and P_D1/P_D2 form a downhill electron transfer pathway. E_m for the first oxidation step, E_m(S₀/S₁), is uniquely low (730 mV) and is ∼100 mV lower than that for the second oxidation step, E_m(S₁/S₂) (830 mV) only when the O4 site of the Mn₄CaO₅ cluster is protonated in S₀. The O4-water chain, which directly forms a low-barrier H-bond with the Mn₄CaO₅ cluster and mediates proton-coupled electron transfer in the S₀ to S₁ transition, explains why the second lowest oxidation state, S₁, is the most stable and S₀ is converted to S₁ even in the dark.

Prediction of Passive Membrane Permeability by Semi-Empirical Method Considering Viscous and Inertial Resistances and Different Rates of Conformational Change and Diffusion

Membrane permeability is an important property of drugs in adsorption. Many prediction methods work well for small molecules, but the prediction of middle-molecule permeability is still difficult. In the present study, we modified a classical permeability model based on Fick's law to study passive membrane permeability. The model consisted of the distribution of solute from water to membrane and the diffusion of solute in each solvent. The diffusion coefficient is the inverse of the resistance, and we examined the inertial resistance in addition to the viscous resistance, the latter of which has been widely used in permeability prediction. Also, we examined three models changing the balance between the diffusion of solute in membrane and the conformational change of solute. The inertial resistance improved the prediction results in addition to the viscous resistance. The models worked well not only for small molecules but also for middle molecules, whose structures have more conformational freedom.

Tracing the Pathways of Waters and Protons in Photosystem II and Cytochrome c Oxidase

Photosystem II (PSII) uses water as the terminal electron donor, producing oxygen in the Mn4CaO5 oxygen evolving complex (OEC), while cytochrome c oxidase (CcO) reduces O2 to water in its heme–Cu binuclear center (BNC). Each protein is oriented in the membrane to add to the proton gradient. The OEC, which releases protons, is located near the P-side (positive, at low-pH) of the membrane. In contrast, the BNC is in the middle of CcO, so the protons needed for O2 reduction must be transferred from the N-side (negative, at high pH). In addition, CcO pumps protons from N- to P-side, coupled to the O2 reduction chemistry, to store additional energy. Thus, proton transfers are directly coupled to the OEC and BNC redox chemistry, as well as needed for CcO proton pumping. The simulations that study the changes in proton affinity of the redox active sites and the surrounding protein at different states of the reaction cycle, as well as the changes in hydration that modulate proton transfer paths, are described.

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.