Influence of Polar Mutations on the Electronic and Structural Properties of QA-· in Bacterial Reaction Centers.

Reaction centers from Rhodobacter sphaeroides with residue M265 mutated from isoleucine to threonine, serine, and asparagine (M265IT, M265IS, and M265IN, respectively) in the QA-· state are studied by high-resolution electron spin echo envelope modulation (ESEEM) and electron nuclear double resonance spectroscopy methods to investigate the structural characteristics of these mutants influencing the redox properties of the QA site. All three mutants decrease the redox midpoint potential (Em) of QA by ∼0.1 V, yet the mechanism for this drop in Em is unclear. In this work, we examine (i) the hydrogen bonding interactions between QA-· and residues histidine M219 and alanine M260, (ii) the electron spin density distribution of the semiquinone, and (iii) the orientations of the ubiquinone methoxy substituents. 13C measurements show no significant contribution of methoxy dihedral angles to the observed decrease in Em for the QA mutants. Instead, 14N three-pulse ESEEM data suggest that electrostatic or hydrogen bond formation between the mutated M265 side chain and His-M219 Nδ may be involved in the observed lowering of the QA midpoint potential. For mutant M265IN, analysis of the proton hyperfine couplings reveals a weakened hydrogen bond network, resulting in an altered QA-· spin density distribution. The magnetic resonance study presented here is most consistent with an electrostatic or structural perturbation of the His-M219 Nδ hydrogen bond in these mutants as a mechanism for the ∼0.1 V decrease in QA Em.

Determination of the Complete Spin Density Distribution in 13C-Labeled Protein-Bound Radical Intermediates Using Advanced 2D Electron Paramagnetic Resonance Spectroscopy and Density Functional Theory.

Determining the complete electron spin density distribution for protein-bound radicals, even with advanced pulsed electron paramagnetic resonance (EPR) methods, is a formidable task. Here we present a strategy to overcome this problem combining multifrequency HYSCORE and ENDOR measurements on site-specifically 13C-labeled samples with DFT calculations on model systems. As a demonstration of this approach, pulsed EPR experiments are performed on the primary QA and secondary QB ubisemiquinones of the photosynthetic reaction center from Rhodobacter sphaeroides 13C-labeled at the ring and tail positions. Despite the large number of nuclei interacting with the unpaired electron in these samples, two-dimensional X- and Q-band HYSCORE and orientation selective Q-band ENDOR resolve and allow for a characterization of the eight expected 13C resonances from significantly different hyperfine tensors for both semiquinones. From these results we construct, for the first time, the most complete experimentally determined maps of the s- and pπ-orbital spin density distributions for any protein organic cofactor radical to date. This work lays a foundation for understanding the relationship between the electronic structure of semiquinones and their functional properties, and introduces new techniques for mapping out the spin density distribution that are readily applicable to other systems.

Ultrafast Electron Transfer Kinetics in the LM Dimer of Bacterial Photosynthetic Reaction Center from Rhodobacter sphaeroides.

It has become increasingly clear that dynamics plays a major role in the function of many protein systems. One system that has proven particularly facile for studying the effects of dynamics on protein-mediated chemistry is the bacterial photosynthetic reaction center from Rhodobacter sphaeroides. Previous experimental and computational analysis have suggested that the dynamics of the protein matrix surrounding the primary quinone acceptor, QA, may be particularly important in electron transfer involving this cofactor. One can substantially increase the flexibility of this region by removing one of the reaction center subunits, the H-subunit. Even with this large change in structure, photoinduced electron transfer to the quinone still takes place. To evaluate the effect of H-subunit removal on electron transfer to QA, we have compared the kinetics of electron transfer and associated spectral evolution for the LM dimer with that of the intact reaction center complex on picosecond to millisecond time scales. The transient absorption spectra associated with all measured electron transfer reactions are similar, with the exception of a broadening in the QX transition and a blue-shift in the QY transition bands of the special pair of bacteriochlorophylls (P) in the LM dimer. The kinetics of the electron transfer reactions not involving quinones are unaffected. There is, however, a 4-fold decrease in the electron transfer rate from the reduced bacteriopheophytin to QA in the LM dimer compared to the intact reaction center and a similar decrease in the recombination rate of the resulting charge-separated state (P(+)QA(-)). These results are consistent with the concept that the removal of the H-subunit results in increased flexibility in the region around the quinone and an associated shift in the reorganization energy associated with charge separation and recombination.

Regulation of the primary quinone binding conformation by the H subunit in reaction centers from Rhodobacter sphaeroides.

Unlike photosystem II (PSII) in higher plants, bacterial photosynthetic reaction centers (bRCs) from Proteobacteria have an additional peripheral membrane subunit "H". The H subunit is necessary for photosynthetic growth, but can be removed chemically in vitro. The remaining LM dimer retains its activity to perform light-induced charge separation. Here we investigate the influence of the H subunit on interactions between the primary semiquinone and the protein matrix, using a combination of site-specific isotope labeling, pulsed electron paramagnetic resonance (EPR), and density functional theory (DFT) calculations. The data reveal substantially weaker binding interactions between the primary semiquinone and the LM dimer than observed for the intact bRC; the amount of electron spin transferred to the nitrogen hydrogen bond donors is significantly reduced, the methoxy groups are more free to rotate, and the spectra indicate a heterogeneous mixture of bound semiquinone states. These results are consistent with a loosening of the primary quinone binding pocket in the absence of the H subunit.

Hydrogen bond network around the semiquinone of the secondary quinone acceptor Q(B) in bacterial photosynthetic reaction centers.

By utilizing a combined pulsed EPR and DFT approach, the high-resolution structure of the QB site semiquinone (SQB) was determined. The development of such a technique is crucial toward an understanding of protein-bound semiquinones on the structural level, as (i) membrane protein crystallography typically results in low resolution structures, and (ii) obtaining protein crystals in the semiquinone form is rarely feasible. The SQB hydrogen bond network was investigated with Q- (∼34 GHz) and X-band (∼9.7 GHz) pulsed EPR spectroscopy on fully deuterated reactions centers from Rhodobacter sphaeroides. Simulations in the SQB g-tensor reference frame provided the principal values and directions of the H-bond proton hyperfine tensors. Three protons were detected, one with an anisotropic tensor component, T = 4.6 MHz, assigned to the histidine NδH of His-L190, and two others with similar anisotropic constants T = 3.2 and 3.0 MHz assigned to the peptide NpH of Gly-L225 and Ile-L224, respectively. Despite the strong similarity in the peptide couplings, all hyperfine tensors were resolved in the Q-band ENDOR spectra. The Euler angles describing the series of rotations that bring the hyperfine tensors into the SQB g-tensor reference frame were obtained by least-squares fitting of the spectral simulations to the ENDOR data. These Euler angles show the locations of the hydrogen bonded protons with respect to the semiquinone. Our geometry optimized model of SQB used in previous DFT work is in strong agreement with the angular constraints from the spectral simulations, providing the foundation for future joint pulsed EPR and DFT semiquinone structural determinations in other proteins.

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

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

Protonated rhodosemiquinone at the Q(B) binding site of the M265IT mutant reaction center of photosynthetic bacterium Rhodobacter sphaeroides.

The second electron transfer from primary ubiquinone Q(A) to secondary ubiquinone Q(B) in the reaction center (RC) from Rhodobacter sphaeroides involves a protonated Q(B)(-) intermediate state whose low pK(a) makes direct observation impossible. Here, we replaced the native ubiquinone with low-potential rhodoquinone at the Q(B) binding site of the M265IT mutant RC. Because the in situ midpoint redox potential of Q(A) of this mutant was lowered approximately the same extent (≈100 mV) as that of Q(B) upon exchange of ubiquinone with low-potential rhodoquinone, the inter-quinone (Q(A) → Q(B)) electron transfer became energetically favorable. After subsequent saturating flash excitations, a period of two damped oscillations of the protonated rhodosemiquinone was observed. The Q(B)H(•) was identified by (1) the characteristic band at 420 nm of the absorption spectrum after the second flash and (2) weaker damping of the oscillation at 420 nm (due to the neutral form) than at 460 nm (attributed to the anionic form). The appearance of the neutral semiquinone was restricted to the acidic pH range, indicating a functional pK(a) of <5.5, slightly higher than that of the native ubisemiquinone (pK(a) < 4.5) at pH 7. The analysis of the pH and temperature dependencies of the rates of the second electron transfer supports the concept of the pH-dependent pK(a) of the semiquinone at the Q(B) binding site. The local electrostatic potential is severely modified by the strongly interacting neighboring acidic cluster, and the pK(a) of the semiquinone is in the middle of the pH range of the complex titration. The kinetic and thermodynamic data are discussed according to the proton-activated electron transfer mechanism combined with the pH-dependent functional pK(a) of the semiquinone at the Q(B) site of the RC.

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.

The rate of second electron transfer to QB(-) in bacterial reaction center of impaired proton delivery shows hydrogen-isotope effect.

The 2nd electron transfer in reaction center of photosynthetic bacterium Rhodobacter sphaeroides is a two step process in which protonation of QB(-) precedes interquinone electron transfer. The thermal activation and pH dependence of the overall rate constants of different RC variants were measured and compared in solvents of water (H2O) and heavy water (D2O). The electron transfer variants where the electron transfer is rate limiting (wild type and M17DN, L210DN and H173EQ mutants) do not show solvent isotope effect and the significant decrease of the rate constant of the second electron transfer in these mutants is due to lowering the operational pKa of QB(-)/QBH: 4.5 (native), 3.9 (L210DN), 3.7 (M17DN) and 3.1 (H173EQ) at pH7. On the other hand, the proton transfer variants where the proton transfer is rate limiting demonstrate solvent isotope effect of pH-independent moderate magnitude (2.11±0.26 (WT+Ni(2+)), 2.16±0.35 (WT+Cd(2+)) and 2.34±0.44 (L210DN/M17DN)) or pH-dependent large magnitude (5.7 at pH4 (L213DN)). Upon deuteration, the free energy and the enthalpy of activation increase in all proton transfer variants by about 1 kcal/mol and the entropy of activation becomes negligible in L210DN/M17DN mutant. The results are interpreted as manifestation of equilibrium and kinetic solvent isotope effects and the structural, energetic and kinetic possibility of alternate proton delivery pathways are discussed.

The 2-Methoxy Group Orientation Regulates the Redox Potential Difference between the Primary (QA) and Secondary (QB) Quinones of Type II Bacterial Photosynthetic Reaction Centers.

Recent studies have shown that only quinones with a 2-methoxy group can act simultaneously as the primary (QA) and secondary (QB) electron acceptors in photosynthetic reaction centers from purple bacteria such as Rb. sphaeroides. 13C HYSCORE measurements of the 2-methoxy group in the semiquinone states, SQA and SQB, were compared with DFT calculations of the 13C hyperfine couplings as a function of the 2-methoxy dihedral angle. X-ray structure comparisons support 2-methoxy dihedral angle assignments corresponding to a redox potential gap (ΔEm) between QA and QB of 175-193 mV. A model having a methyl group substituted for the 2-methoxy group exhibits no electron affinity difference. This is consistent with the failure of a 2-methyl ubiquinone analogue to function as QB in mutant reaction centers with a ΔEm of ∼160-195 mV. The conclusion reached is that the 2-methoxy group is the principal determinant of electron transfer from QA to QB in type II photosynthetic reaction centers with ubiquinone serving as both acceptor quinones.

Hyperfine and nuclear quadrupole tensors of nitrogen donors in the Q(A) site of bacterial reaction centers: correlation of the histidine N(δ) tensors with hydrogen bond strength.

X- and Q-band pulsed EPR spectroscopy was applied to study the interaction of the QA site semiquinone (SQA) with nitrogens from the local protein environment in natural abundance (14)N and in (15)N uniformly labeled photosynthetic reaction centers of Rhodobacter sphaeroides. The hyperfine and nuclear quadrupole tensors for His-M219 Nδ and Ala-M260 peptide nitrogen (Np) were estimated through simultaneous simulation of the Q-band (15)N Davies ENDOR, X- and Q-band (14,15)N HYSCORE, and X-band (14)N three-pulse ESEEM spectra, with support from DFT calculations. The hyperfine coupling constants were found to be a((14)N) = 2.3 MHz, T = 0.3 MHz for His-M219 Nδ and a((14)N) = 2.6 MHz, T = 0.3 MHz for Ala-M260 Np. Despite that His-M219 Nδ is established as the stronger of the two H-bond donors, Ala-M260 Np is found to have the larger value of a((14)N). The nuclear quadrupole coupling constants were estimated as e(2)Qq/4h = 0.38 MHz, η = 0.97 and e(2)Qq/4h = 0.74 MHz, η = 0.59 for His-M219 Nδ and Ala-M260 Np, respectively. An analysis of the available data on nuclear quadrupole tensors for imidazole nitrogens found in semiquinone-binding proteins and copper complexes reveals these systems share similar electron occupancies of the protonated nitrogen orbitals. By applying the Townes-Dailey model, developed previously for copper complexes, to the semiquinones, we find the asymmetry parameter η to be a sensitive probe of the histidine Nδ-semiquinone hydrogen bond strength. This is supported by a strong correlation observed between η and the isotropic coupling constant a((14)N) and is consistent with previous computational works and our own semiquinone-histidine model calculations. The empirical relationship presented here for a((14)N) and η will provide an important structural characterization tool in future studies of semiquinone-binding proteins.

Nuclear hyperfine and quadrupole tensor characterization of the nitrogen hydrogen bond donors to the semiquinone of the QB site in bacterial reaction centers: a combined X- and S-band (14,15)N ESEEM and DFT study.

The secondary quinone anion radical QB(-) (SQB) in reaction centers of Rhodobacter sphaeroides interacts with Nδ of His-L190 and Np (peptide nitrogen) of Gly-L225 involved in hydrogen bonds to the QB carbonyls. In this work, S-band (∼3.6 GHz) ESEEM was used with the aim of obtaining a complete characterization of the nuclear quadrupole interaction (nqi) tensors for both nitrogens by approaching the cancelation condition between the isotropic hyperfine coupling and (14)N Zeeman frequency at lower microwave frequencies than traditional X-band (9.5 GHz). By performing measurements at S-band, we found a dominating contribution of Nδ in the form of a zero-field nqi triplet at 0.55, 0.92, and 1.47 MHz, defining the quadrupole coupling constant K = e(2)qQ/4h = 0.4 MHz and associated asymmetry parameter η = 0.69. Estimates of the hyperfine interaction (hfi) tensors for Nδ and Np were obtained from simulations of 1D and 2D (14,15)N X-band and three-pulse (14)N S-band spectra with all nuclear tensors defined in the SQB g-tensor coordinate system. From simulations, we conclude that the contribution of Np to the S-band spectrum is suppressed by its strong nqi and weak isotropic hfi comparable to the level of hyperfine anisotropy, despite the near-cancelation condition for Np at S-band. The excellent agreement between our EPR simulations and DFT calculations of the nitrogen hfi and nqi tensors to SQB is promising for the future application of powder ESEEM to full tensor characterizations.

Special issue of photosynthetic research.

This Special Issue of Photosynthesis Research honors Louis M. N. Duysens, Roderick K. Clayton, and George Feher, three pioneering researchers whose work on bacterial photosynthesis laid much of the groundwork for our understanding of the role of the reaction center in photosynthetic light energy conversion. Their key discoveries are briefly summarized and an overview of the special issue is presented.

Roderick K. Clayton: a life, and some personal recollections.

Roderick K. Clayton passed away on October 23, 2011, at the age of 89, shortly after the plan for this dedicatory issue of Photosynthesis Research had been hatched. I had just written a lengthy letter to him to re-establish contact after a hiatus of 2 or 3 years, and to suggest that I visit him to talk about his life. It isn't clear whether he saw the letter or not, but it was found at his home in Santa Rosa, California. Fortunately, Rod has written two memoirs for Photosynthesis Research that not only cover much of his research on reaction centers (Photosynth Res 73:63-71, 2002) but also provide a humorous and honest look at his personal life (Photosynth Res 19:207-224, 1988). I cannot hope to improve on these and will try, instead, to fill in some of the gaps that Rod's own writing has left, and offer some of my own personal recollections over the more recent years.

Tuning cofactor redox potentials: the 2-methoxy dihedral angle generates a redox potential difference of >160 mV between the primary (Q(A)) and secondary (Q(B)) quinones of the bacterial photosynthetic reaction center.

Only quinones with a 2-methoxy group can act simultaneously as the primary (QA) and secondary (QB) electron acceptors in photosynthetic reaction centers from Rhodobacter sphaeroides. (13)C hyperfine sublevel correlation measurements of the 2-methoxy in the semiquinone states, SQA and SQB, were compared with quantum mechanics calculations of the (13)C couplings as a function of the dihedral angle. X-ray structures support dihedral angle assignments corresponding to a redox potential gap (ΔEm) between QA and QB of ~180 mV. This is consistent with the failure of a ubiquinone analogue lacking the 2-methoxy to function as QB in mutant reaction centers with a ΔEm of ≈160-195 mV.

Conformational differences between the methoxy groups of QA and QB site ubisemiquinones in bacterial reaction centers: a key role for methoxy group orientation in modulating ubiquinone redox potential.

Ubiquinone is an almost universal, membrane-associated redox mediator. Its ability to accept either one or two electrons allows it to function in critical roles in biological electron transport. The redox properties of ubiquinone in vivo are determined by its environment in the binding sites of proteins and by the dihedral angle of each methoxy group relative to the ring plane. This is an attribute unique to ubiquinone among natural quinones and could account for its widespread function with many different redox complexes. In this work, we use the photosynthetic reaction center as a model system for understanding the role of methoxy conformations in determining the redox potential of the ubiquinone/semiquinone couple. Despite the abundance of X-ray crystal structures for the reaction center, quinone site resolution has thus far been too low to provide a reliable measure of the methoxy dihedral angles of the primary and secondary quinones, QA and QB. We performed 2D ESEEM (HYSCORE) on isolated reaction centers with ubiquinones (13)C-labeled at the headgroup methyl and methoxy substituents, and have measured the (13)C isotropic and anisotropic components of the hyperfine tensors. Hyperfine couplings were compared to those derived by DFT calculations as a function of methoxy torsional angle allowing estimation of the methoxy dihedral angles for the semiquinones in the QA and QB sites. Based on this analysis, the orientation of the 2-methoxy groups are distinct in the two sites, with QB more out of plane by 20-25°. This corresponds to an ≈50 meV larger electron affinity for the QB quinone, indicating a substantial contribution to the experimental difference in redox potentials (60-75 mV) of the two quinones. The methods developed here can be readily extended to ubiquinone-binding sites in other protein complexes.

Hydrogen bonding between the Q(B) site ubisemiquinone and Ser-L223 in the bacterial reaction center: a combined spectroscopic and computational perspective.

In the Q(B) site of the Rhodobacter sphaeroides photosynthetic reaction center, the donation of a hydrogen bond from the hydroxyl group of Ser-L223 to the ubisemiquinone formed after the first flash is debatable. In this study, we use a combination of spectroscopy and quantum mechanics/molecular mechanics (QM/MM) calculations to comprehensively explore this topic. We show that ENDOR, ESEEM, and HYSCORE spectroscopic differences between mutant L223SA and the wild-type sample (WT) are negligible, indicating only minor perturbations in the ubisemiquinone spin density for the mutant sample. Qualitatively, this suggests that a strong hydrogen bond does not exist in the WT between the Ser-L223 hydroxyl group and the semiquinone O(1) atom, as removal of this hydrogen bond in the mutant should cause a significant redistribution of spin density in the semiquinone. We show quantitatively, using QM/MM calculations, that a WT model in which the Ser-L223 hydroxyl group is rotated to prevent hydrogen bond formation with the O(1) atom of the semiquinone predicts negligible change for the L223SA mutant. This, together with the better agreement between key QM/MM calculated and experimental hyperfine couplings for the non-hydrogen-bonded model, leads us to conclude that no strong hydrogen bond is formed between the Ser-L223 hydroxyl group and the semiquinone O(1) atom after the first flash. The implications of this finding for quinone reduction in photosynthetic reaction centers are discussed.

Hydrogen bonding and spin density distribution in the Qb semiquinone of bacterial reaction centers and comparison with the Qa site.

In the photosynthetic reaction center from Rhodobacter sphaeroides, the primary (Q(A)) and secondary (Q(B)) electron acceptors are both ubiquinone-10, but with very different properties and functions. To investigate the protein environment that imparts these functional differences, we have applied X-band HYSCORE, a 2D pulsed EPR technique, to characterize the exchangeable protons around the semiquinone (SQ) in the Q(A) and Q(B) sites, using samples of (15)N-labeled reaction centers, with the native high spin Fe(2+) exchanged for diamagnetic Zn(2+), prepared in (1)H(2)O and (2)H(2)O solvent. The powder HYSCORE method is first validated against the orientation-selected Q-band ENDOR study of the Q(A) SQ by Flores et al. (Biophys. J.2007, 92, 671-682), with good agreement for two exchangeable protons with anisotropic hyperfine tensor components, T, both in the range 4.6-5.4 MHz. HYSCORE was then applied to the Q(B) SQ where we found proton lines corresponding to T ≈ 5.2, 3.7 MHz and T ≈ 1.9 MHz. Density functional-based quantum mechanics/molecular mechanics (QM/MM) calculations, employing a model of the Q(B) site, were used to assign the observed couplings to specific hydrogen bonding interactions with the Q(B) SQ. These calculations allow us to assign the T = 5.2 MHz proton to the His-L190 N(δ)H···O(4) (carbonyl) hydrogen bonding interaction. The T = 3.7 MHz spectral feature most likely results from hydrogen bonding interactions of O1 (carbonyl) with both Gly-L225 peptide NH and Ser-L223 hydroxyl OH, which possess calculated couplings very close to this value. The smaller 1.9 MHz coupling is assigned to a weakly bound peptide NH proton of Ile-L224. The calculations performed with this structural model of the Q(B) site show less asymmetric distribution of unpaired spin density over the SQ than seen for the Q(A) site, consistent with available experimental data for (13)C and (17)O carbonyl hyperfine couplings. The implications of these interactions for Q(B) function and comparisons with the Q(A) site are discussed.

The binding interface of cytochrome c and cytochrome c1 in the bc1 complex: rationalizing the role of key residues.

The interaction of cytochrome c with ubiquinol-cytochrome c oxidoreductase (bc1 complex) has been studied for >30 years, yet many aspects remain unclear or controversial. We report the first molecular dynamic simulations of the cyt c-bc1 complex interaction. Contrary to the results of crystallographic studies, our results show that there are multiple dynamic hydrogen bonds and salt bridges in the cyt c-c1 interface. These include most of the basic cyt c residues previously implicated in chemical modification studies. We suggest that the static nature of x-ray structures can obscure the quantitative significance of electrostatic interactions between highly mobile residues. This provides a clear resolution of the discrepancy between the structural data and functional studies. It also suggests a general need to consider dynamic interactions of charged residues in protein-protein interfaces. In addition, a novel structural change in cyt c is reported, involving residues 21-25, which may be responsible for cyt c destabilization upon binding. We also propose a mechanism of interaction between cyt c1 monomers responsible for limiting the binding of cyt c to only one molecule per bc1 dimer by altering the affinity of the cytochrome c binding site on the second cyt c1 monomer.

Photoelectrochemical complexes for solar energy conversion that chemically and autonomously regenerate.

Naturally occurring photosynthetic systems use elaborate pathways of self-repair to limit the impact of photo-damage. Here, we demonstrate a complex consisting of two recombinant proteins, phospholipids and a carbon nanotube that mimics this process. The components self-assemble into a configuration in which an array of lipid bilayers aggregate on the surface of the carbon nanotube, creating a platform for the attachment of light-converting proteins. The system can disassemble upon the addition of a surfactant and reassemble upon its removal over an indefinite number of cycles. The assembly is thermodynamically metastable and can only transition reversibly if the rate of surfactant removal exceeds a threshold value. Only in the assembled state do the complexes exhibit photoelectrochemical activity. We demonstrate a regeneration cycle that uses surfactant to switch between assembled and disassembled states, resulting in an increased photoconversion efficiency of more than 300% over 168 hours and an indefinite extension of the system lifetime.