Nature Communications
Home Crystal structure of a photosynthetic LH1-RC in complex with its electron donor HiPIP
Crystal structure of a photosynthetic LH1-RC in complex with its electron donor HiPIP
Crystal structure of a photosynthetic LH1-RC in complex with its electron donor HiPIP

Article Type: research-article Article History
Abstract

Photosynthetic electron transfers occur through multiple components ranging from small soluble proteins to large integral membrane protein complexes. Co-crystallization of a bacterial photosynthetic electron transfer complex that employs weak hydrophobic interactions was achieved by using high-molar-ratio mixtures of a soluble donor protein (high-potential iron-sulfur protein, HiPIP) with a membrane-embedded acceptor protein (reaction center, RC) at acidic pH. The structure of the co-complex offers a snapshot of a transient bioenergetic event and revealed a molecular basis for thermodynamically unfavorable interprotein electron tunneling. HiPIP binds to the surface of the tetraheme cytochrome subunit in the light-harvesting (LH1) complex-associated RC in close proximity to the low-potential heme-1 group. The binding interface between the two proteins is primarily formed by uncharged residues and is characterized by hydrophobic features. This co-crystal structure provides a model for the detailed study of long-range trans-protein electron tunneling pathways in biological systems.

The high potential iron-sulfur (HiPIP) proteins are direct electron donors to the light-harvesting-reaction center complexes (LH1-RC) in photosynthetic β- and γ-Proteobacteria. Here, the authors present the 2.9 Å crystal structure of the HiPIP-bound LH1-RC complex from the thermophilic purple sulfur bacterium Thermochromatium tepidum and discuss mechanistic implications for the electron transfer pathway.

Keywords
Kawakami,Yu,Liang,Okazaki,Madigan,Kimura,and Wang-Otomo: Crystal structure of a photosynthetic LH1-RC in complex with its electron donor HiPIP

Introduction

Photosynthetic organisms oxidize either water or reduced organic or inorganic compounds as electron sources to produce new biomass. While noncyclic electron transport is the major pathway in oxygen-evolving photosynthetic organisms, electron transport in anoxygenic phototrophic bacteria is characterized by a light-driven cyclical process coupled to proton translocation1. In purple phototrophic bacteria, the cyclic electron transport chain is composed of three major components, two membrane-embedded protein complexes—reaction center (RC) and cytochrome (Cyt) bc1—and a group of soluble proteins that shuttle electrons between the RC and Cyt bc1 complexes. In most purple bacteria, the RC contains a tightly bound multi-heme Cyt subunit that functions as the immediate electron donor to the special pair bacteriochlorophylls (BChl) and as an acceptor from the soluble carriers2,3. A variety of soluble c-type cytochromes and high-potential iron-sulfur proteins (HiPIP) function as electron donors to the RC4,5. For many purple bacteria of the α-proteobacterial group6, Cyt c2 is the sole or major electron donor. This is the case for Rhodobacter (Rb.) sphaeroides, for example, where Cyt c2 transfers electrons from Cyt bc1 to the RC. The Rb. sphaeroides Cyt c2-bound RC complex has been crystallized and its structure determined7,8. The structure showed that Cyt c2 is positioned at the center of the periplasmic surface of the RC and that the binding site contains both long-range electrostatic interactions to facilitate rapid association and short-range hydrophobic interactions to facilitate electron tunneling.

By contrast to electron flow in photosynthetic α-Proteobacteria, the HiPIP functions as a direct electron donor to the RC in photosynthetic β- and γ-Proteobacteria9,10. Until now, no structure has been determined for an electron transfer complex formed between a soluble carrier and the RC bound with a multi-heme Cyt subunit. We present here the crystallographic structure of an HiPIP-bound RC in which the RC contains the tetraheme Cyt subunit and is associated with its core light-harvesting complex (LH1-RC, ~390 kDa). Both the HiPIP and LH1-RC were isolated from the thermophilic purple sulfur bacterium Thermochromatium (Tch.) tepidum (γ-Proteobacteria)11 and their individual structures have been determined to high resolutions12,13. Our structural results verify predictions of protein–protein interactions between the Tch. tepidum tetraheme Cyt subunit and HiPIP and go well beyond this to provide important new insight into the mechanism of electron transfer that occurs in this key step of photosynthetic energy generation.

Results

Both a high molar ratio of Tch. tepidum HiPIP to LH1-RC and acidic pH were required for co-crystal formation based on MALDI/TOF-MS measurements taken at the initial stage of screening (Supplementary Figs. 1 and 2). Further screenings of the two parameters by analyzing diffraction data over the ranges of HiPIP:LH1-RC = 1–30 and pH 4.5–7.5 revealed the optimum conditions to be near a ratio of 15 and at pH 5.0 (data not shown) as described in the “Methods”. Based on the absorption spectrum (Supplementary Fig. 1), Tch. tepidum HiPIP was purified in the reduced state.

Structural overview

The crystal structure of the HiPIP:LH1-RC co-complex was determined at 2.9 Å resolution. Each LH1-RC bound one HiPIP at the surface of the Cyt subunit, covering a range from Glu60 to Asp120 and in close proximity to the low-potential heme-1 (Fig. 1 and Supplementary Fig. 3). The amino acid sequences of these proteins in Tch. tepidum were compared with those of other species in Supplementary Fig. 3. The contact surface area between the two proteins was calculated to be 667 Å2, similar in order of magnitude to that of the Cyt c2:RC co-complex of the purple bacterium Rb. sphaeroides8. Upon binding to the Cyt subunit, the HiPIP exhibited an uneven B-factor distribution with smaller values for the amino acid residues on the binding interface where Gln62 and Leu63 showed the smallest B-factors (Fig. 1b and Supplementary Fig. 4). This indicates that the interface residues of HiPIP have significantly more restricted molecular motions compared with those on the opposite side exposed to the solvent.

Structure of the HiPIP:LH1-RC co-complex.
Fig. 1

Structure of the HiPIP:LH1-RC co-complex.

a Side view of the overall structure of the HiPIP:LH1-RC co-complex. Color scheme: HiPIP, red; LH1-α, green; LH1-β, slate blue; RC-Cyt subunit, cyan; RC-L subunit, violet; RC-M subunit, white tint; RC-H subunit, orange; BChl a, red sticks; Spirilloxanthin, yellow sticks; heme, magenta sticks; Fe-S, yellow-orange ball; Ca, gold ball. b Expanded view of the region marked in a where HiPIP is shown as sticks and Cyt subunit is shown in transparent surface representation. Hemes are shown in magenta sticks and 4Fe-4S cluster is shown as balls. The residues in HiPIP are color-coded by the values of their B-factor as indicated in the scale bar.

The Tch. tepidum HiPIP:LH1-RC co-complex had a more compact crystal packing than did the LH1-RC alone, as indicated by reductions of both the solvent content from 61.9 to 57.1% and the Matthew coefficient from 3.22 to 2.97 Å3/Da (Fig. 2a). The positions of the Fe atoms in both hemes and the 4Fe-4S cluster were clearly resolved from the strong anomalous difference Fourier map (Fig. 2b). Most sidechains of HiPIP residues could be modeled based on the electron density map with higher accuracies for the interface residues because of their lower B-factors (Fig. 2c). The interface between the two proteins displayed a zigzag shape in which the large sidechains of Tyr76 and Trp116 of the Cyt subunit and Phe48 of HiPIP likely function as shape-recognition residues. The nearest distance between the two prosthetic groups was 6.5 Å as measured from the sulfur S1 in HiPIP to the 21-methyl carbon of heme-1, and the distance between S1 and Fe (heme-1) was 12.2 Å (Fig. 2d). Coordinate error was estimated to be 0.50 Å based on the maximum-likelihood refinement14.

Structural arrangements of the HiPIP:LH1-RC co-complex at different scales.
Fig. 2

Structural arrangements of the HiPIP:LH1-RC co-complex at different scales.

a Crystal packing of the LH1-RC (blue) and HiPIP (red). The rectangular box (blue lines) shows a unit cell and each unit cell contains one HiPIP:LH1-RC. b Anomalous difference Fourier maps measured at 0.9 Å for the C-subunit and HiPIP. Red meshes represent the contour levels around Fe atoms at 4.2σ. c Electron density maps at a contour level of 1.4σ for the interface area between HiPIP (right) and Cyt subunit (left). The 4Fe-4S cluster is shown by a blue box. This figure was drawn using Coot software. d Distances (Å) between the 4Fe-4S cluster in HiPIP and the heme-1 in the Cyt subunit.

Interactions between Tch. tepidum HiPIP and Cyt subunit

Figure 3 shows the amino acid residues on the interface of the Tch. tepidum HiPIP and Cyt subunit. Most of the residues are distributed around the 4Fe-4S cluster and heme-1. The path between the two prosthetic groups is mainly composed of uncharged residues: Thr13, Leu17, Phe48, Leu63, and Ser77 in the HiPIP (Fig. 3c) and Tyr76, Thr91, Val95, Thr98, and Trp116 in the Cyt subunit (Fig. 3b). The closest distance between the two proteins was 2.3 Å as measured for the hydroxyl oxygen of Thr79 (HiPIP) and the aromatic nitrogen of Trp116 (Cyt subunit) which form a strong hydrogen bond (Fig. 3a and Supplementary Table 2). Other relatively weak hydrogen bonds may be formed between Arg33(HiPIP)/Asn108(Cyt subunit), Gln62(HiPIP)/Tyr76(Cyt subunit), and Gln62(HiPIP)/Thr91(Cyt subunit) (Supplementary Table 2).

Amino acid residues on the interface of HiPIP (orange) and the Cyt subunit (cyan).
Fig. 3

Amino acid residues on the interface of HiPIP (orange) and the Cyt subunit (cyan).

Side view (a) and front views of the Cyt subunit (b) and HiPIP (c). A hydrogen bond formed between Trp116 (Cyt subunit) and Thr79 (HiPIP) is shown by a dashed line in a. Sidechains of the residues are shown in sticks. Protein structures are drawn by transparent cartoons. Macrocycle of the heme-1 in the Cyt subunit is shown as red framework. Color scheme: oxygen, red; nitrogen, blue; sulfur, yellow.

It is notable that the binding surface of HiPIP is characterized by its hydrophobicity as indicated by the electrostatically neutral distribution formed by uncharged residues (Supplementary Fig. 5). As a result, the surface charge distribution is virtually unaffected by a change of pH (Supplementary Fig. 5a, b). In contrast, the opposite side of the binding surface of HiPIP is largely acidic at pH 7 (Supplementary Fig. 5c) and apparently shifts toward neutral at pH 5 (Supplementary Fig. 5d). By contrast, the binding surface of the Cyt subunit is relatively basic at pH 7 (Supplementary Fig. 5e), especially near heme-1, and becomes even more basic at pH 5 (Supplementary Fig. 5f). Considering that the HiPIP:LH1-RC co-crystal was preferentially formed at pH 5 and at high molar ratios of HiPIP to LH1-RC, the results imply that interactions between the acidic side of HiPIP (Supplementary Fig. 5c) and the basic surface of the Cyt subunit may weaken at acidic pH, thus increasing the probability of binding of the two proteins at the correct position required for efficient electron transfer.

Structural changes upon binding of Tch. tepidum HiPIP to LH1-RC

Superposition of the Cα carbons of the membrane-embedded RC L-subunit in the co-complex with those in free LH1-RC revealed that conformational changes mainly occurred on the binding surface of the Cyt subunit, with a maximum deviation of 1.2 Å (Fig. 4a, b). Residues closer to the surface displayed larger deviations from those of free LH1-RC, and large changes in position were observed for the residues of Thr68–Asn78 and Asn108–Ser118 around the heme-1-binding site in the Cyt subunit (Supplementary Fig. 3). No detectable conformational changes were found in the LH1 complex.

Structural changes in the HiPIP:LH1-RC co-complex.
Fig. 4

Structural changes in the HiPIP:LH1-RC co-complex.

a Partial structures for the RC portions obtained by superposition of the Cα carbons of the reaction center L-subunit in HiPIP:LH1-RC (colored) with those in free LH1-RC (gray, PDB ID: 5Y5S). Part of the HiPIP is shown in transparent orange cartoon, and the heme-1 in Cyt subunit is shown in transparent cyan sticks. b Expanded view of the region marked in a with several residues indicated that have deviations larger than 1.0 Å. The Cα traces of RC are color-coded by the values of root-mean-square deviation (RMSD) as indicated in the scale bar. c Side view of superposition of the carbon atoms of the pyrrole ring I in heme-1 (deepteal color) with those in free LH1-RC (gray, PDB ID: 5Y5S). d Structural changes of the 4Fe-4S cluster in the HiPIP:LH1-RC co-complex obtained by superposition of the Cα carbons of the HiPIP (colored) with those in free HiPIP (gray, PDB ID: 1EYT). Sulfur and Fe atoms are shown by yellow and orange balls, respectively.

A remarkable change was observed for heme-1 in the Cyt subunit of the co-complex (Fig. 4c). The heme-1 macrocycle displayed a significantly bent conformation along the diagonal of pyrrole rings I and III compared to that of heme-1 in the free LH1-RC (Supplementary Fig. 6a). Superposition of the pyrrole rings I of the heme-1 groups revealed a maximum deviation of 1.3 Å at the edge of pyrrole III (Fig. 4c). This feature may be related to the function of heme-1 because pyrrole ring I is nearest the 4Fe-4S cluster in HiPIP (Fig. 2d) and is predicted as the entry point for electrons donated by the 4Fe-4S cluster. By contrast, essentially no change was found for the conformation of HiPIP (Supplementary Fig. 6b). When superimposing the Cα carbons of the HiPIP in co-complex with those in free HiPIP, only a slight rotation was observed for the 4Fe-4S cluster (Fig. 4d), with a maximum deviation of 0.4 Å at sulfur S4. Moreover, no relative changes in distances and angles among the atoms in the cluster were found (Supplementary Fig. 6c).

Predicted electron transfer pathway

From the co-crystal structure, we could construct a five-step trans-protein electron tunneling chain from the HiPIP to the special pair of the RC via the tetraheme Cyt subunit—a distance that spans more than 70 Å (Fig. 5a and Supplementary Fig. 7). The center-to-center distances between each redox pair are in the range of 12.2–20.6 Å as given in Supplementary Table 3. Based on Pathways calculations15, the first step from 4Fe-4S to heme-1 is mediated by the sidechain of Leu63 in HiPIP, which is located in the middle, 3.5 Å away from both sulfur S1 in HiPIP and 21-methyl carbon in heme-1 (Fig. 5b). Based on its maximum coupling strength and close proximity to the Fe in heme-1, sulfur S1 was predicted as the atom in the 4Fe-4S cluster that donates electrons to heme-1 (Supplementary Table 3). The putative pathway within the Cyt subunit is an inter-heme electron tunneling process (center-to-center distances of 13.9–16.2 Å) with relatively greater coupling strengths compared to that of S1-to-heme-1 (Supplementary Table 3). The electrons in heme-1 were transferred along pyrrole rings I and II to the heme-2, with an edge-to-edge distance of 4.0 Å. While pyrrole rings I and II also participate in the electron paths in heme-2 and heme-3, the route in heme-4 was predicted as passing along the pyrrole rings II and IV. It is unclear whether the different routes are related to different coordinations of the central Fe atoms in which the Fe atoms in heme-1, heme-2, and heme-3 are ligated by His and Met residues whereas the Fe atom in heme-4 is ligated by two His residues. The last step of electron transfer to the special pair is relayed by the aromatic ring of Tyr171 in the L-subunit, which is located directly above the BChl a dimer (Fig. 5c). This Tyr residue is conserved in both Cyt subunit-possessing (Blastochloris viridis and Rubrivivax gelatinosus) and Cyt subunit-deficient (Rb. sphaeroides) species of purple bacteria. In the Cyt c2:RC co-crystal structure of Rb. sphaeroides, the corresponding Tyr(L162) was identified as the bridging residue for electron transfer from the heme group in Cyt c2 to the special pair8, and the heme plane is largely overlapped with that of heme-3 in the Tch. tepidum Cyt subunit. The final pathways to the central Mg atoms of the BChl a dimer in the L- and M-subunits were predicted to pass through pyrrole ring II in both BChl a molecules.

Predicted electron transfer pathways.
Fig. 5

Predicted electron transfer pathways.

a Trans-protein electron tunneling pathway from 4Fe-4S cluster to the special pair predicted by the Pathways plugin for VMD using standard parameters. Through-bond process and through-space hoppings are represented by thick magenta sticks and black arrows, respectively. b Expanded view of the portion of 4Fe-4S cluster in HiPIP and heme-1 in the Cyt subunit. Leu63 in HiPIP is shown by sticks. Reduction potentials are given in parentheses. c Expanded view of the portion of heme-3 in the Cyt subunit and the special pair BChl a dimer. The sidechain of the Tyr171 in the reaction center L-subunit is shown by sticks.

Discussion

Approximately equimolar mixtures of proteins have been used in co-crystallizations that mainly employ electrostatic interactions8,1620. In contrast, the high molar ratios (≥10) of HiPIP to LH1-RC were required for successful co-crystallization in our work. The requirement could be due to the weak interactions that allow for rapid equilibrium between the two proteins and exclude formation of a long-lived complex21. The high ratios needed likely increase the chance of forming an encounter complex because the HiPIP-to-RC electron transfer is a collisional reaction governed by hydrophobic forces21,22. Additionally, an acidic pH was also required for the formation of co-crystals; these conditions probably masked the negative charges on the protein surfaces (Supplementary Fig. 5), thus enhancing hydrophobic interactions. The hydrophobic nature of the interactions is evidenced by both the uncharged feature of the HiPIP-binding surface (Supplementary Fig. 5) and the ionic strength-dependence of electron transfer measurements21. This was confirmed by isothermal titration calorimetry (ITC) analysis that revealed positive peaks (endothermic heat change) for the binding of HiPIP to LH1-RC (Supplementary Fig. 8), indicating an entropically driven process. The entropy-driven signature is typically attributed to hydrophobic interactions accompanied by the release of bound water molecules or ions along with conformational changes23,24. Such has also been reported for the binding of ferredoxin, a water-soluble protein containing an iron-sulfur cluster, to photosystem I25 and ferredoxin-NADP+-reductase26,27 in oxygenic phototrophs.

The ultra-high-resolution structure (0.48 Å) of Tch. tepidum HiPIP revealed details for the charge distribution around its 4Fe-4S cluster at levels of iron 3d and sulfur 3p orbital electrons12. Charge-density analysis yielded a total atomic charge of –1.5 for the 4Fe-4S–4(Cys-Sγ) group, a value close to the total charge of this cluster in the reduced state. Based on an individual valence-charge analysis, the Fe1 (with its ligand Cys43-Sγ), Fe2, and S4 atoms were identified to be crucial for storing electron charges in the reduced free HiPIP12. In the Tch. tepidum HiPIP:LH1-RC structure, most of these atoms in HiPIP are located on the distal side of the binding surface (Fig. 4d), indicating a role as sources of electrons to the front-side (solvent-exposed) Fe and S atoms that have been proposed to form a delocalized 2Fe–2S subcluster in HiPIP28.

The reduction potential of the Tch. tepidum HiPIP was determined to be +340 mV29, and the two low-potential hemes (heme-1 and heme-4) in the Tch. tepidum Cyt subunit were estimated to have reduction potentials below +150 mV30. This points to a thermodynamically unfavorable electron transfer from HiPIP to heme-1, followed by a down–up–down–down flow to the special pair in the RC (Fig. 5a)31,32. This “uphill” nature of HiPIP oxidation has fueled debate over which heme group in the tetraheme Cyt subunit is actually reduced by HiPIP; for example, heme-2 is more electropositive (~+330 mV) than heme-1 and would therefore be a more thermodynamically favorable acceptor of electrons from HiPIP than would heme-1. However, despite its lower reduction potential, heme-1 has been convincingly shown to be the direct electron acceptor from HiPIP21,22,33 and various c-type cytochromes3436. Using an empirical formula for the endergonic tunneling proposed using fixed geometry by Page et al.37, the rate of electron transfer from the 4Fe-4S cluster of Tch. tepidum HiPIP to heme-1 was calculated to be 1.5 × 105 s−1 assuming a reorganization energy of 0.7 eV37, a volume fraction of 0.77, and a reduction potential of +10 mV for heme-1 derived from Allochromatium (Alc.) vinosum38, a mesophilic relative of Tch. tepidum. Despite the endergonic nature of this first step, the calculated rate is about four orders of magnitude faster than the measured value (~20 s1) for the Tch. tepidum HiPIP/RC system in the range of micromolar concentrations29. A similar pseudo-first-order rate constant was also observed from the Alc. vinosum HiPIP/RC system21. These results indicate that complex formation of HiPIP and LH1-RC is the rate-limiting step for the interprotein electron transfer and that the electron transfer from 4Fe-4S to heme-1 is the rate-limiting step in the entire pathway to the special pair in the RC.

A significant distortion of heme-1 in the Cyt subunit of the co-complex was detected upon HiPIP binding (Fig. 4c and Supplementary Fig. 6a). Heme distortion has been suggested to modulate various properties of the heme group including its reduction potential. It was demonstrated that decreasing heme distortion resulted in a decrease in the reduction potential39. This infers that bending of heme-1 upon HiPIP binding should increase its reduction potential, making it energetically more favorable to accept electrons from HiPIP by reducing the energy gap. However, it should be pointed out that the effect of heme distortion on reduction potential may be less than that induced by electrostatic interactions and/or hydrophobicity around the heme pocket40. There is a strong correlation between the reduction potential and the surrounding hydrophobicity: reduction potential increases with decreasing heme exposure to solvent41. Based on this result, the reduction potential of heme-1 in the HiPIP:LH1-RC complex is expected to increase because the solvent-exposed edge of heme-1 as seen in Supplementary Fig. 5e, f is covered by the hydrophobic surface of HiPIP, leading to an increased hydrophobicity upon HiPIP binding. This effect again could reduce the energy gap, making the heme-1 more favorable for accepting electrons from HiPIP.

Based on the HiPIP:LH1-RC co-crystal structure determined in this work, Leu63 in Tch. tepidum HiPIP is predicted to be the bridging residue for the electron transfer from the 4Fe-4S cluster to heme-1 (Fig. 5b), consistent with the result of early molecular calculations using a docking model5. However, despite the high sequence identity (87%; Supplementary Fig. 3) and structural similarity between the Tch. tepidum and Alc. vinosum HiPIPs, the conserved Phe48 in Alc. vinosum HiPIP was predicted to be the mediating residue in the docking model5, indicating that the bridging residue could vary between species and that the distance between the donor and acceptor is the main determinant of the rate of electron tunneling37. The aromatic ring of Phe48 in Alc. vinosum HiPIP was almost parallel to the plane of the heme-1 group21; this is in contrast to that in the Tch. tepidum HiPIP:LH1-RC structure where the two planes are inclined by 55˚ toward each other (Fig. 2c). Moreover, the contact surface area of the Tch. tepidum HiPIP:LH1-RC in our study was found to be significantly smaller than that calculated for the HiPIP-docked RC of Alc. vinosum21. This points to difficulties in docking-model-based calculations that do not fully account for weak hydrophobic interactions.

In summary, our structural analysis has revealed a critical step of interprotein electron transfer that precedes the intra-protein steps in the purple bacterial photosynthetic RC, and a step that is the rate-determining one in the entire electron transport chain leading to the QB site. Moreover, the Tch. tepidum HiPIP:LH1-RC structure resolves uncertainty over which heme group in the tetraheme Cyt subunit is the direct electron acceptor from the HiPIP and provides a structural snapshot of how this thermodynamically unfavorable interprotein electron transfer occurs. The LH1-bound RC with HiPIP structure also unites light-harvesting, energy transfer, charge separation, and electron transport into one system and therefore provides a more robust model for studying LH-coupled long-range trans-protein electron tunneling processes in photosynthetic organisms.

Methods

Co-crystallization of the HiPIP:LH1-RC complex

Purifications of the HiPIP and LH1-RC from Tch. tepidum strain MC were conducted as previously described29,42. Soluble proteins were separated from membrane fractions by centrifugation of the sonicated cells. The HiPIP was purified by a DEAE column (Toyopearl 650S, TOSOH) with 200 mM NaCl in 20 mM Tris-HCl buffer (pH 8.5). LH1-RC was solubilized from the membrane fraction by two-step detergent treatments (0.3%w/v LDAO and 1%w/v DDM), followed by purification using the same DEAE column with 20 mM Tris-HCl buffer (pH 7.5) containing 0.05 w/v DDM and 50 mM CaCl2. Absorption spectra and a co-crystal of the purified proteins are shown in Supplementary Fig. 1. Based on the spectrum, HiPIP was purified in the reduced state. Co-crystallization was carried out using a micro-batch-under-oil method13 by mixing HiPIP and LH1-RC in a stoichiometric ratio of 15:1 in a buffer containing 30 mM sodium succinate (pH 5), 20 mM CaCl2, 3.4%(w/v) n-octyl-phosphocholine (OPC), and 30%(w/v) polyethylene glycol (PEG) 1500. Molar extinction coefficients used for calculations of the protein concentrations were 4322 cm−1 M−1 for LH1-RC at 915 nm43 and 17.3 cm1 M−1 for HiPIP at 389 nm obtained from inductively coupled plasma atomic emission spectroscopy44, which is consistent with the reported value21. Concentration of the LH1-RC used for co-crystallization was 0.0694 mM. Binding of the HiPIP to LH1-RC was confirmed by ITC of the individual proteins and by matrix-assisted laser desorption ionization time-of-flight mass spectroscopy (MALDI/TOF-MS) of the dissolved crystals using the methods described elsewhere44,45. ITC profile was obtained with titration of 0.134 mM HiPIP in a buffer containing 25 mM CaCl2 and 0.05% DDM to 0.0353 mM LH1-RC solution in the same buffer at 25 ˚C. Co-crystals were obtained at 20 °C with suitable sizes (~0.4 × 0.7 × 0.3 mm) after 5 days (Supplementary Fig. 1). Post-crystallization treatment followed the same method as that for the LH1-RC13. The co-crystals were transferred to into cryoprotectant solution containing 30 mM sodium succinate (pH 5), 20 mM CaCl2, 3.4%w/v OPC, 15%w/v glycerol, and 30%w/v PEG 1500, and flash-frozen immediately in a liquid-nitrogen container.

Diffraction data collection and structural analysis

X-ray diffraction data were collected at BL41XU and BL44XU of SPring-8 (Hyogo, Japan), and BL1A and BL17A of the Photon Factory (Ibaraki, Japan) at 100 K. The diffraction data were processed, integrated, and scaled using the XDS package46. Structure determination was performed by the molecular replacement method using the Tch. tepidum LH1-RC structure (PDB ID 5Y5S)13 as search model by the Phaser program47 in PHENIX48. The bound HiPIP structure was manually modeled based on the electron density map using the Tch. tepidum HiPIP structure (PDB ID 1EYT, 1.5-Å resolution)49 by Coot50. The initial HiPIP:LH1-RC model was then refined by Phenix.refine14. Positional and isotropic displacement parameters were refined in the resolution range of 30.0–2.9 Å. The final model was refined to 2.9 Å resolution with Rfree = 24.8% and Rwork = 21.9% (Supplementary Table 1). Accessible surface areas were estimated using AREAIMOL51 in CCP4 with a probe radius of 1.4 Å. Surface charge distributions were calculated using the Adaptive Poisson-Boltzmann Solver software suite52 and the PDB2PQR sever (ver. 2.1.1)53. The electron transfer pathways from the 4Fe-4S cluster in HiPIP to the special pair BChls in the RC through the tetraheme Cyt subunit were predicted using the Pathways plugin15 for VMD54. Figures were generated with Pymol55 unless otherwise stated. Coordinates of the HiPIP:LH1-RC co-crystal structure have been deposited in PDB under accession code 7C52.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Peer review information: Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
These authors contributed equally: Tomoaki Kawakami, Long-Jiang Yu.
These authors jointly supervised this work: Long-Jiang Yu, Yukihiro Kimura, Zheng-Yu Wang-Otomo.

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-021-21397-9.

Acknowledgements

We thank Dr. I. Balabin for instruction on the Pathways plugin for VMD, Kao Corporation for providing detergent LDAO. This work was performed in part using a synchrotron beamline BL44XU at SPring-8 under the Collaborative Research Program of Institute for Protein Research, Osaka University. X-ray diffraction data were collected at BL41XU (proposal numbers: 2017AB2512 and 2018AB2518) and at the Osaka University beamline BL44XU (proposal numbers: 2017AB6708, 2018AB6806, 2019AB6906, and 2020AB6506) at SPring-8 (Harima, Japan), and BL1A, BL17A, and NE3A at Photon Factory of KEK (proposal numbers: 2017G512 and 2019G513). We are grateful to E. Yamashita, N. Matsugaki, and Y. Yamada for their support in data collection. This work was supported in part by JSPS KAKENHI Grant Numbers JP16H04174, JP18H05153, JP20H05086, and JP20H02856, Takeda Science Foundation, the Kurata Memorial Hitachi Science and Technology Foundation, Japan, the grant of KINOU-KYOKA from Institute of Quantum Beam Science of Ibaraki University, and the National Key R&D Program of China (No. 2019YFA0904600).

Author contributions

Z.-Y.W.-O. designed the work; L.-J.Y. and Y.K. coordinated the experiments; M.T.M. provided materials; T.K., T.L., and K.O. performed the experiments; T.K., L.-J.Y., Y.K., and Z.-Y.W.-O. analyzed data; Z.-Y.W.-O. and M.T.M. wrote the manuscript.

Data availability

Coordinates and structure factors were deposited in the Protein Data Bank under accession code 7C52. All other data are available from the corresponding authors upon reasonable request.

Competing interests

The authors declare no competing interests.

References

1. 

Blankenship, R. E. Molecular Mechanisms of Photosynthesis (Blackwell Science, 2002).

2. 

Nitschke, W. & Dracheva, S. M. In Anoxygenic Photosynthetic Bacteria (eds. Blankenship, R. E., Madigan, M. T. & Bauer, C. E.) 775–805 (Kluwer Academic Publishers, 1995).

3. 

    Tsukatani Y, . Phylogenetic distribution of unusual triheme to tetraheme cytochrome subunit in the reaction center complex of purple bacteria. Photosynth. Res.2004. 79: 83-91 doi: 10.1023/B:PRES.0000011922.56394.92

4. 

    van Driessche G, . Amino acid sequences and distribution of high-potential iron-sulfur proteins that donate electrons to the photosynthetic reaction center in phototrophic proteobacteria. J. Mol. Evol.2003. 57: 181-199 doi: 10.1007/s00239-003-2465-y

5. 

    Ciurli S, Musiani F. High potential iron-sulfur proteins and their role as soluble electron carriers in bacterial photosynthesis: tale of a discovery. Photosynth. Res.2005. 85: 115-131 doi: 10.1007/s11120-004-6556-4

6. 

    Nagashima S, Nagashima KVP. Comparison of photosynthesis gene clusters retrieved from total genome sequences of purple bacteria. Adv. Bot. Res.2013. 66: 151-178 doi: 10.1016/B978-0-12-397923-0.00005-9

7. 

    Adir N, . Co-crystallization and characterization of the photosynthetic reaction center–cytochrome c2 complex from Rhodobacter sphaeroides. Biochemistry1996. 35: 2535-2547 doi: 10.1021/bi9522054

8. 

    Axelrod HL, . X-ray structure determination of the cytochrome c2:reaction center electron transfer complex from Rhodobacter sphaeroides. J. Mol. Biol.2002. 319: 501-515 doi: 10.1016/S0022-2836(02)00168-7

9. 

    Verméglio A, Li J, Schoepp-Cothenet B, Pratt N, Knaff DB. The role of high-potential iron sulfur and cytochrome c8 as alternative electron donors to the reaction center of Chromatium vinosum. Biochemistry2002. 41: 8868-8875 doi: 10.1021/bi012037h

10. 

    Nagashima KVP, Matsuura K, Shimada K, Verméglio A. High-potential iron-sulfur protein (HiPIP) is the major electron donor to the reaction center complex in photosynthetically growing cells of the purple bacterium Rubrivivax gelatinosus. Biochemistry2002. 41: 14028-14032 doi: 10.1021/bi026511a

11. 

    Madigan MT. A novel photosynthetic bacterium isolated from a Yellowstone hot spring. Science1984. 225: 313-315 doi: 10.1126/science.225.4659.313

12. 

    Hirano Y, Takeda K, Miki K. Charge-density analysis of an iron-sulfur protein at an ultra-high resolution of 0.48 Å. Nature2016. 534: 281-284 doi: 10.1038/nature18001

13. 

    Yu L-J, Suga M, Wang-Otomo Z-Y, Shen J-R. Structure of photosynthetic LH1-RC supercomplex at 1.9 Å resolution. Nature2018. 556: 209-213 doi: 10.1038/s41586-018-0002-9

14. 

    Afonine PV, . Towards automated crystallographic structure refinement with phenix.refine. Acta Cryst.2012. D68: 352-367

15. 

    Balabin IA, Hu X, Beratan DN. Exploring biological electron transfer pathway dynamics with the Pathways plugin for VMD. J. Comput. Chem.2012. 33: 906-910 doi: 10.1002/jcc.22927

16. 

    Pelletier H, Kraut J. Crystal structure of a complex between electron transfer partners, cytochrome c peroxidase and cytochrome c. Science1992. 258: 1748-1755 doi: 10.1126/science.1334573

17. 

    Solmaz SRN, Hunte C. Structure of complex III with bound cytochrome c in reduced state and definition of minimal core interface for electron transfer. J. Biol. Chem.2008. 283: 17542-17549 doi: 10.1074/jbc.M710126200

18. 

    Díaz-Moreno I, . The dynamic complex of cytochrome c6 and cytochrome f studied with paramagnetic NMR spectroscopy. Biochim. Biophys. Acta Bioenerg.2014. 1837: 1305-1315 doi: 10.1016/j.bbabio.2014.03.009

19. 

    Terasaka E, . Dynamics of nitric oxide controlled by protein complex in bacterial system. Proc. Natl Acad. Sci. USA2017. 114: 9888-9893 doi: 10.1073/pnas.1621301114

20. 

    Shimada S, . Complex structure of cytochrome c–cytochrome c oxidase reveals a novel protein–protein mode. EMBO J.2017. 36: 291-300 doi: 10.15252/embj.201695021

21. 

    Venturoli G, . Electron transfer from HiPIP to the photooxidized tetraheme cytochrome subunit of Allochromatium vinosum reaction center: new insights from site-directed mutagenesis and computational studies. Biochemistry2004. 43: 437-445 doi: 10.1021/bi035384v

22. 

    Osyczka A, . Comparison of the binding sites for high-potential iron-sulfur protein and cytochrome c on the tetraheme cytochrome subunit bound to the bacterial photosynthetic reaction center. Biochemistry1999. 38: 15779-15790 doi: 10.1021/bi990907d

23. 

    Chodera JD, Mobley DL. Entropy–enthalpy compensation: role and ramifications in biomolecular ligand recognition and design.. Annu. Rev. Biophys.2013. 42: 121-142 doi: 10.1146/annurev-biophys-083012-130318

24. 

    Fox JM, Zhao M, Fink MJ, Kang K, Whitesides GM. The molecular origin of enthalpy/entropy compensation in biomolecular recognition. Annu. Rev. Biophys.2018. 47: 223-250 doi: 10.1146/annurev-biophys-070816-033743

25. 

    Marco P, . Binding of ferredoxin to algal photosystem I involves a single binding site and is composed of two thermodynamically distinct events. Biochim. Biophys. Acta Bioenerg.2018. 1859: 234-243 doi: 10.1016/j.bbabio.2018.01.001

26. 

    Jelesarov I, Bosshard HR. Thermodynamics of ferredoxin binding to ferredoxin:NADP+ reductase and the role of water at the complex interface. Biochemistry1994. 33: 13321-13328 doi: 10.1021/bi00249a019

27. 

    Martínez-Júlvez M, Medina M, Velázquez-Campoy A. Binding thermodynamics of ferredoxin:NADP+ reductase: two different protein substrates and one energetics. Biophys. J.2009. 96: 4966-4975 doi: 10.1016/j.bpj.2009.02.061

28. 

    Dey A, . Sulfur K-edge XAS and DFT calculations on [Fe4S4] clusters: effects of H-bonding and structural distortion on covalency and spin topology. Inorg. Chem.2005. 44: 8349-8354 doi: 10.1021/ic050981m

29. 

    Kobayashi M, Saito T, Takahashi K, Wang Z-Y, Nozawa T. Electronic properties and thermal stability of soluble redox proteins from a thermophilic purple sulfur photosynthetic bacterium, Thermochromatium tepidum. Bull. Chem. Soc. Jpn2005. 78: 2164-2170 doi: 10.1246/bcsj.78.2164

30. 

    Drepper F, Saito T, Kobayashi M, Nozawa T, Mathis P. Electron transfer reactions of high-potential cytochromes in the reaction center of Chromatium tepidum. Photosynth. Res.1998. 55: 325-330 doi: 10.1023/A:1005903802573

31. 

    Chen I-P, Mathis P, Koepke J, Michel H. Uphill electron transfer in the tetraheme cytochrome subunit of the Rhodopseudomonas viridis photosynthetic reaction center: evidence from site-directed mutagenesis. Biochemistry2000. 39: 3592-3602 doi: 10.1021/bi992443p

32. 

    Alric J, . Kinetic performance and energy profile in a roller coaster electron transfer chain: a study of modified tetraheme-reaction center constructs. J. Am. Chem. Soc.2006. 128: 4136-4145 doi: 10.1021/ja058131t

33. 

    Osyczka A, Nagashima KVP, Shimada K, Matsuura K. Interaction site for high-potential iron-sulfur protein on the tetraheme cytochrome subunit bound to the photosynthetic reaction center of Rubrivivax gelatinosus. Biochemistry1999. 38: 2861-2865 doi: 10.1021/bi982747w

34. 

    Knaff DB, . Reaction of cytochrome c2 with photosynthetic reaction centers from Rhodopseudomonas viridis. Biochemistry1991. 30: 1303-1310 doi: 10.1021/bi00219a021

35. 

    Osyczka A, . Interaction site for soluble cytochromes on the tetraheme cytochrome subunit bound to the bacterial photosynthetic reaction center mapped by site-directed mutagenesis. Biochemistry1998. 37: 11732-11744 doi: 10.1021/bi980910h

36. 

    Ortega JM, Drepper F, Mathis P. Electron transfer between cytochrome c2 and the tetraheme cytochrome c in Rhodopseudomonas viridis. Photosynth. Res.1999. 59: 147-157 doi: 10.1023/A:1006149621029

37. 

    Page CC, Moser CC, Chen X, Dutton PL. Natural engineering principles of electron tunneling in biological oxidation-reduction. Nature1999. 402: 47-52 doi: 10.1038/46972

38. 

    Nitschke W, Jubault-Bregler M, Rutherford AW. The reaction center associated tetraheme cytochrome subunit from Chromatium vinosum revisited: a reexamination of its EPS properties. Biochemistry1993. 32: 8871-8879 doi: 10.1021/bi00085a019

39. 

    Olea J, . Probing the function of heme distortion in the H-NOX family. ACS Chem. Biol.2008. 3: 703-710 doi: 10.1021/cb800185h

40. 

    Saito K, . Deformation of chlorin rings in the photosystem II crystal structure. Biochemistry2012. 51: 4290-4299 doi: 10.1021/bi300428s

41. 

    Tezcan FA, Winkler JR, Gray HB. Effects of ligation and folding on reduction potentials of heme proteins. J. Am. Chem. Soc.1998. 120: 13383-13388 doi: 10.1021/ja982536e

42. 

    Suzuki H, . Purification, characterization and crystallization of the core complex from thermophilic purple sulfur bacterium Thermochromatium tepidum. Biochim. Biophs. Acta Bioenerg.2007. 1767: 1057-1063 doi: 10.1016/j.bbabio.2007.06.002

43. 

    Nagatsuma S, . Phospholipid distributions in purple phototrophic bacteria and LH1-RC core complexes. Biochim. Biophys. Acta Bioenerg.2019. 1860: 461-468 doi: 10.1016/j.bbabio.2019.04.001

44. 

    Kimura Y, . Effects of calcium ions on the thermostability and spectroscopic properties of the LH1-RC complex from a new thermophilic purple bacterium Allochromatium tepidum. J. Phys. Chem. B2017. 121: 5025-5032 doi: 10.1021/acs.jpcb.7b03341

45. 

    Kimura Y, . C-terminal cleavage of the LH1 α-polypeptide in the Sr2+-cultured Thermochromatium tepidum. Photosynth. Res.2018. 135: 23-31 doi: 10.1007/s11120-017-0393-8

46. 

    Kabsch W. XDS. Acta Cryst.2010. D66: 125-132

47. 

48. 

    Adams PD, . PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Cryst.2010. D66: 213-221

49. 

    Nogi T, Fathir I, Kobayashi M, Nozawa T, Miki K. Crystal structures of photosynthetic reaction center and high-potential iron-sulfur protein from Thermochromatium tepidum: thermostability and electron transfer. Proc. Natl Acad. Sci. USA2000. 97: 13561-13566 doi: 10.1073/pnas.240224997

50. 

    Emsley P, Lohkamp B, Scott WG, Cowtan K. Features and development of Coot. Acta Cryst.2010. D66: 486-501

51. 

    Lee B, Richards FM. The interpretation of protein structures: Estimation of static accessibility. J. Mol. Biol.1971. 55: 379-400 doi: 10.1016/0022-2836(71)90324-X

52. 

    Jurrus E, . Improvements to the APBS biomolecular solvation software suite. Protein Sci.2018. 27: 112-128 doi: 10.1002/pro.3280

53. 

    Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA. PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatic calculations. Nucleic Acids Res.2004. 32: W665-W667 doi: 10.1093/nar/gkh381

54. 

55. 

DeLano, W. L. The PyMOL Molecular Graphics System (DeLano Scientific, LCC, 2004).