Redox potential of quinones in photosynthetic reaction centers from Rhodobacter sphaeroides: dependence on protonation of Glu-L212 and Asp-L213

Absence of a link between stabilized charge-separated state and structural changes proposed from crystal structures of a photosynthetic reaction center

Structural differences between illuminated and unilluminated crystal structures led to the proposal that the charge-separated state was stabilized by structural changes in its membrane extrinsic protein subunit H in a bacterial photosynthetic reaction center [Katona, G. et al. Nat. Struct. Mol. Biol. 2005, 12, 630-631]. Here, we explored the proposal by titrating all titratable sites and calculating the redox potential (Em) values in these crystal structures. Contrary to the expected charge-separated states, Em for quinone, Em(QA/QA•-), is even lower in the proposed charge-separated structure than in the ground-state structure. The subunit-H residues, which were proposed to exhibit electron-density changes in the two crystal structures, contribute to an Em(QA/QA•-) difference of only <0.5 mV. Furthermore, the protonation states of the titratable residues in the entire reaction center are practically identical in the two structures. These findings indicate that the proposed structural differences are irrelevant to explaining the significant prolongation of the charge-separated-state lifetime.

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.

Electron Acceptor-Donor Iron Sites in the Iron-Sulfur Cluster of Photosynthetic Electron-Transfer Pathways

In photosystem I, two electron-transfer pathways via quinones (A_1A and A_1B) are merged at the iron-sulfur Fe₄S₄ cluster F_X into a single pathway toward the other two Fe₄S₄ clusters F_A and F_B. Using a quantum mechanical/molecular mechanical approach, we identify the redox-active Fe sites in the clusters. In F_A and F_B, the Fe site, which does not belong to the CxxCxxCxxxCP motif, serves as an electron acceptor/donor. F_X has two independent electron acceptor Fe sites for A- and B-branch electron transfers, depending on the Asp-B575 protonation state, which causes the A_1A-to-F_X electron transfer to be uphill and the A_1B-to-F_X electron transfer to be downhill. The two asymmetric electron-transfer pathways from A₁ to F_X and the separation of the electron acceptor and donor Fe sites are likely associated with the specific role of F_X in merging the two electron transfer pathways into the single pathway.

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.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

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.

Role of AP-endonuclease (Ape1) active site residues in stabilization of the reactant enzyme-DNA complex

Apurinic/apyrimidinic endonuclease 1 (Ape1) is an important metal-dependent enzyme in the base excision repair mechanism, responsible for the backbone cleavage of abasic DNA through a phosphate hydrolysis reaction. Molecular dynamics simulations of Ape1 complexed to its substrate DNA performed for models containing 1 or 2 Mg²⁺ -ions as cofactor located at different positions show a complex with 1 metal ion bound on the leaving group site of the scissile phosphate to be the most likely reaction-competent conformation. Active-site residue His309 is found to be protonated based on pKa calculations and the higher conformational stability of the Ape1-DNA substrate complex compared to scenarios with neutral His309. Simulations of the D210N mutant further support the prevalence of protonated His309 and strongly suggest Asp210 as the general base for proton acceptance by a nucleophilic water molecule.

Synthesis and characterization of chlorophyll a enol derivatives: Chlorophyll a tert-butyldimethylsilyl-enol ether and 131-deoxo-131, 132-didehydro-chlorophyll a

Using the sterically hindered base, 1,8-diazabicyclo[5.4.0]undec-7-ene, for enolization and tert-butyldimethylsilyl-trifluoromethanesulfonate for silylation, chlorophyll (Chl) [Formula: see text] produced after 15 min at 0 [Formula: see text]C in deaerated pyridine solution under argon, after work-up and chromatographic purification on a sucrose column, tert-butyldimethylsilyl-enol ether of Chl [Formula: see text] in a yield of 77%. The 131-deoxo-131,132-didehydro-chlorophyll [Formula: see text], was obtained in a yield of 23% through a reaction sequence, where Chl [Formula: see text] was first reduced with sodium borohydride to 13[Formula: see text]-hydroxy-Chl [Formula: see text], which via demetalation yielded 13[Formula: see text]-hydroxypheophytin [Formula: see text]. In the presence of the sterically hindered base, 1,8-bis(dimethylamino)naphthalene, trifluoroacetylimidazole dehydrated 13[Formula: see text]-hydroxypheophytin [Formula: see text] to 131-deoxo-131,132-didehydro-pheophytin [Formula: see text], which after metalation yielded 131-deoxo-131,132-didehydro-Chl [Formula: see text]. Using 1,8-bis(dimethylamino)naphthalene and trifluoroacetylimidazole, the straight conversion of 13[Formula: see text]-hydroxy-Chl [Formula: see text] to 131-deoxo-131,132-didehydro-Chl [Formula: see text] was found unsuccessful. The major products were characterized by electronic absorption spectra (UV-vis) and practically completely assigned 1H and [Formula: see text]C NMR spectra. Some intermediates of the syntheses were also characterized by ESI-TOF mass spectra. Compared with Chl [Formula: see text], the macrocyclic ring-current in the synthesized Chl [Formula: see text] enol derivatives was found weakened by the expansion of the [Formula: see text]-system to include the isocyclic ring E. Nevertheless, these enol derivatives were still considered to be diamagnetic and aromatic. The possibility of the functional role of the enol derivatives of chlorophyll in photosynthesis is discussed.

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.

Stigmatellin probes the electrostatic potential in the QB site of the photosynthetic reaction center

The electrostatic potential in the secondary quinone (QB) binding site of the reaction center (RC) of the photosynthetic bacterium Rhodobacter sphaeroides determines the rate and free energy change (driving force) of electron transfer to QB. It is controlled by the ionization states of residues in a strongly interacting cluster around the QB site. Reduction of the QB induces change of the ionization states of residues and binding of protons from the bulk. Stigmatellin, an inhibitor of the mitochondrial and photosynthetic respiratory chain, has been proven to be a unique voltage probe of the QB binding pocket. It binds to the QB site with high affinity, and the pK value of its phenolic group monitors the local electrostatic potential with high sensitivity. Investigations with different types of detergent as a model system of isolated RC revealed that the pK of stigmatellin was controlled overwhelmingly by electrostatic and slightly by hydrophobic interactions. Measurements showed a high pK value (>11) of stigmatellin in the QB pocket of the dark-state wild-type RC, indicating substantial negative potential. When the local electrostatics of the QB site was modulated by a single mutation, L213Asp → Ala, or double mutations, L213Asp-L212Glu → Ala-Ala (AA), the pK of stigmatellin dropped to 7.5 and 7.4, respectively, which corresponds to a >210 mV increase in the electrostatic potential relative to the wild-type RC. This significant pK drop (ΔpK > 3.5) decreased dramatically to (ΔpK > 0.75) in the RC of the compensatory mutant (AA+M44Asn → AA+M44Asp). Our results indicate that the L213Asp is the most important actor in the control of the electrostatic potential in the QB site of the dark-state wild-type RC, in good accordance with conclusions of former studies using theoretical calculations or light-induced charge recombination assay.

Electrochemical characterization of hydroquinone derivatives with different substituents in acetonitrile

The effect of carbonyl groups in the ortho position with respect to a hydroxyl group on the electrochemical oxidation of hydroquinones in acetonitrile is studied.

The nonheme iron in photosystem II

Photosystem II (PSII), the light-driven water:plastoquinone (PQ) oxidoreductase of oxygenic photosynthesis, contains a nonheme iron (NHI) at its electron acceptor side. The NHI is situated between the two PQs QA and QB that serve as one-electron transmitter and substrate of the reductase part of PSII, respectively. Among the ligands of the NHI is a (bi)carbonate originating from CO2, the substrate of the dark reactions of oxygenic photosynthesis. Based on recent advances in the crystallography of PSII, we review the structure of the NHI in PSII and discuss ideas concerning its function and the role of bicarbonate along with a comparison to the reaction center of purple bacteria and other enzymes containing a mononuclear NHI site.

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.

Dynamical theory of primary processes of charge separation in the photosynthetic reaction center

A dynamical theory has been developed for primary separation of charges in the course of photosynthesis. The theory deals with both hopping and superexchange transfer mechanisms. Dynamics of electron transfer from dimeric bacteriochlorophyll to quinone has been calculated. The results obtained agree with experimental data and provide a unified explanation of both the hierarchy of the transfer time in the photosynthetic reaction center and the phenomenon of coherent oscillations accompanying the transfer process.

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(2)O was observed at ∼2200 cm(-1) in BR but was significantly higher in ASR (>2500 cm(-1)), 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.

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.

Light-induced quinone reduction in photosystem II

The photosystem II core complex is the water:plastoquinone oxidoreductase of oxygenic photosynthesis situated in the thylakoid membrane of cyanobacteria, algae and plants. It catalyzes the light-induced transfer of electrons from water to plastoquinone accompanied by the net transport of protons from the cytoplasm (stroma) to the lumen, the production of molecular oxygen and the release of plastoquinol into the membrane phase. In this review, we outline our present knowledge about the "acceptor side" of the photosystem II core complex covering the reaction center with focus on the primary (Q(A)) and secondary (Q(B)) quinones situated around the non-heme iron with bound (bi)carbonate and a comparison with the reaction center of purple bacteria. Related topics addressed are quinone diffusion channels for plastoquinone/plastoquinol exchange, the newly discovered third quinone Q(C), the relevance of lipids, the interactions of quinones with the still enigmatic cytochrome b559 and the role of Q(A) in photoinhibition and photoprotection mechanisms. This article is part of a Special Issue entitled: Photosystem II.

MCCE2: improving protein pKa calculations with extensive side chain rotamer sampling

Multiconformation continuum electrostatics (MCCE) explores different conformational degrees of freedom in Monte Carlo calculations of protein residue and ligand pK(a)s. Explicit changes in side chain conformations throughout a titration create a position dependent, heterogeneous dielectric response giving a more accurate picture of coupled ionization and position changes. The MCCE2 methods for choosing a group of input heavy atom and proton positions are described. The pK(a)s calculated with different isosteric conformers, heavy atom rotamers and proton positions, with different degrees of optimization are tested against a curated group of 305 experimental pK(a)s in 33 proteins. QUICK calculations, with rotation around Asn and Gln termini, sampling His tautomers and torsion minimum hydroxyls yield an RMSD of 1.34 with 84% of the errors being <1.5 pH units. FULL calculations adding heavy atom rotamers and side chain optimization yield an RMSD of 0.90 with 90% of the errors <1.5 pH unit. Good results are also found for pK(a)s in the membrane protein bacteriorhodopsin. The inclusion of extra side chain positions distorts the dielectric boundary and also biases the calculated pK(a)s by creating more neutral than ionized conformers. Methods for correcting these errors are introduced. Calculations are compared with multiple X-ray and NMR derived structures in 36 soluble proteins. Calculations with X-ray structures give significantly better pK(a)s. Results with the default protein dielectric constant of 4 are as good as those using a value of 8. The MCCE2 program can be downloaded from http://www.sci.ccny.cuny.edu/~mcce.

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.

Modulation of the protein environment in the hydrophilic pore of the ammonia transporter protein AmtB upon GlnK protein binding

The conduction of ammonia/ammonium (NH3/NH4(+)) through the channel protein AmtB is inhibited by the binding of the signal transduction protein GlnK. In the AmtB-GlnK binding interface, there exists an NH3/NH4(+) binding site--Am6. The calculated pK(a) values at the Am6 sites in both the AmtB-GlnK complex and isolated AmtB implies the dominance of an uncharged NH3 state. The GlnK protein binding causes a significant downshift in the Am6 pK(a) value of the AmtB. However, this downshift is perfectly compensated by the reorientation of the protein backbone (carbonyl group of Cys312 from the AmtB part) upon AmtB-GlnK complex formation.

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.

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.

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.

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.

Modeling binding kinetics at the Q(A) site in bacterial reaction centers

Bacterial reaction centers (RCs) catalyze a series of electron-transfer reactions reducing a neutral quinone to a bound, anionic semiquinone. The dissociation constants and association rates of 13 tailless neutral and anionic benzo- and naphthoquinones for the Q(A) site were measured and compared. The K(d) values for these quinones range from 0.08 to 90 microM. For the eight neutral quinones, including duroquinone (DQ) and 2,3-dimethoxy-5-methyl-1,4-benzoquinone (UQ(0)), the quinone concentration and solvent viscosity dependence of the association rate indicate a second-order rate-determining step. The association rate constants (k(on)) range from 10(5) to 10(7) M(-)(1) s(-)(1). Association and dissociation rate constants were determined at pH values above the hydroxyl pK(a) for five hydroxyl naphthoquinones. These negatively charged compounds are competitive inhibitors for the Q(A) site. While the neutral quinones reach equilibrium in milliseconds, anionic hydroxyl quinones with similar K(d) values take minutes to bind or dissociate. These slow rates are independent of ionic strength, solvent viscosity, and quinone concentration, indicating a first-order rate-limiting step. The anionic semiquinone, formed by forward electron transfer at the Q(A) site, also dissociates slowly. It is not possible to measure the association rate of the unstable semiquinone. However, as the protein creates kinetic barriers for binding and releasing anionic hydroxyl quinones without greatly increasing the affinity relative to neutral quinones, it is suggested that the Q(A) site may do the same for anionic semiquinone. Thus, the slow semiquinone dissociation may not indicate significant thermodynamic stabilization of the reduced species in the Q(A) 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.

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).

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.

Accurate pKa Determination for a Heterogeneous Group of Organic Molecules

Single-molecule studies that allow to compute pKa values, proton affinities (gas-phase acidity/basicity) and the electrostatic energy of solvation have been performed for a heterogeneous set of 26 organic compounds. Quantum mechanical density functional theory (DFT) using the Becke-half&half and B3LYP functionals on optimized molecular geometries have been carried out to investigate the energetics of gas-phase protonation. The electrostatic contribution to the solvation energies of protonated and deprotonated compounds were calculated by solving the Poisson equation using atomic charges generated by fitting the electrostatic potential derived from the molecular wave functions in vacuum. The combination of gas-phase and electrostatic solvation energies by means of the thermodynamic cycle enabled us to compute pKa values for the 26 compounds, which cover six distinct chemical groups (carboxylic acids, benzoic acids, phenols, imides, pyridines and imidazoles). The computational procedure for determining pKa values is accurate and transferable with a root-mean-square deviation of 0.53 and 0.57 pKa units and a maximum error of 1.0 pKa and 1.3 pKa units for Becke-half&half and B3LYP DFT functionals, respectively.

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

Ser-L223 is close to ubiquinone (Q(B)) in the B-branch of the bacterial photosynthetic reaction center (bRC) from Rhodobacter (Rb) sphaeroides. Therefore, the presence of a hydrogen bond (H bond) between the two was naturally proposed from the crystal structure. The hydrogen bonding pattern of Q(B) from the light-exposed structure was studied by generating hydrogen atom coordinates based on the CHARMM force field. In the Q(B) neutral charge state (Q(B)(0)), no H bond was found between the oxygen of the OH group from Ser-L223 and the carbonyl oxygen of Q(B) that is distal to the non-heme iron. In the reduced state (Q(B)(-)), however, Ser-L213 was found to form an H bond with Q(B) only when Asp-L213 is protonated by more than 0.75 H(+). This indicates the significance of the protonation of Asp-L213 in forming an H bond between Ser-L223 and Q(B). We found that the driving force to form the H bond between Ser-L223 and Q(B) is enhanced by the positively charged Arg-L217. The calculated Q(B) redox potentials with or without this H bond discriminated two ET rates, which are close to the faster and slower time phases observed in UV-Vis and FTIR studies. Together with the calculated redox potential of the quinones, this H-bond formation could play a key role in conformational gating for the ET process from Q(A) to Q(B).

Identification of a Novel Protonation Pattern for Carboxylic Acids upon QB Photoreduction in Rhodobacter sphaeroides Reaction Center Mutants at Asp-L213 and Glu-L212 Sites

In the reaction center from the photosynthetic purple bacterium Rhodobacter sphaeroides, light energy is rapidly converted to chemical energy through coupled electron-proton transfer to a buried quinone molecule Q(B). Involved in the proton uptake steps are carboxylic acids, which have characteristic infrared vibrations that are observable using light-induced Fourier transform infrared (FTIR) difference spectroscopy. Upon formation, Q(B)(-) induces protonation of Glu-L212, located within 5 A of Q(B), resulting in a IR signal at 1728 cm(-1). However, no other IR signal is observed within the classic absorption range of protonated carboxylic acids (1770-1700 cm(-1)). In particular, no signal for Asp-L213 is found despite its juxtaposition to Q(B) and importance for proton uptake on the second electron-transfer step. In an attempt to uncover the reason behind this lack of signal, the microscopic electrostatic environment in the vicinity of Q(B) was modified by interchanging Asp and Glu at the L213 and L212 positions. The Q(B)(-)/Q(B) FTIR spectrum of the Asp-L212/Glu-L213 swap mutant in the 1770-1700 cm(-1) range shows several distinct new signals, which are sensitive to (1)H/(2)H isotopic exchange, indicating that the reduction of Q(B) results in the change of the protonation state of several carboxylic acids. The new bands at 1752 and 1747 cm(-1) were assigned to an increase of protonation in response to Q(B) reduction of Glu-L213 and Asp-L212, respectively, based on the effect of replacing them with their amine analogues. Since other carboxylic acid signals were observed, it is concluded that the swap mutations at L212 and L213 affect a cluster of carboxylic acids larger than the L212/L213 acid pair. Implications for the native reaction center are discussed.

Redox potential of quinones in both electron transfer branches of photosystem I

The redox potentials of the two electron transfer (ET) active quinones in the central part of photosystem I (PSI) were determined by evaluating the electrostatic energies from the solution of the Poisson-Boltzmann equation based on the crystal structure. The calculated redox potentials are -531 mV for A1A and -686 mV for A1B. From these results we conclude the following. (i) Both branches are active with a much faster ET in the B-branch than in the A-branch. (ii) The measured lifetime of 200-290 ns of reduced quinones agrees with the estimate for the A-branch and corroborates with an uphill ET from this quinone to the iron-sulfur cluster as observed in recent kinetic measurements. (iii) The electron paramagnetic resonance spectroscopic data refer to the A-branch quinone where the corresponding ET is uphill in energy. The negative redox potential of A1 in PSI is primarily because of the influence from the negatively charged FX, in contrast to the positive shift on the quinone redox potential in bacterial reaction center and PSII that is attributed to the positively charged non-heme iron atom. The conserved residue Asp-B575 changes its protonation state after quinone reduction. The difference of 155 mV in the quinone redox potentials of the two branches were attributed to the conformation of the backbone with a large contribution from Ser-A692 and Ser-B672 and to the side chain of Asp-B575, whose protonation state couples differently with the formation of the quinone radicals.