Binding Site Recognition and Docking Dynamics of a Single Electron Transport Protein: Cytochrome c2.

Small diffusible redox proteins facilitate electron transfer in respiration and photosynthesis by alternately binding to their redox partners and integral membrane proteins and exchanging electrons. Diffusive search, recognition, binding, and unbinding of these proteins often amount to kinetic bottlenecks in cellular energy conversion, but despite the availability of structures and intense study, the physical mechanisms controlling redox partner interactions remain largely unknown. The present molecular dynamics study provides an all-atom description of the cytochrome c2-docked bc1 complex in Rhodobacter sphaeroides in terms of an ensemble of favorable docking conformations and reveals an intricate series of conformational changes that allow cytochrome c2 to recognize the bc1 complex and bind or unbind in a redox state-dependent manner. In particular, the role of electron transfer in triggering a molecular switch and in altering water-mediated interface mobility, thereby strengthening and weakening complex formation, is described. The results resolve long-standing discrepancies between structural and functional data.

Molecular organization of cytochrome c2 near the binding domain of cytochrome bc1 studied by electron spin-lattice relaxation enhancement.

Measurements of specific interactions between proteins are challenging. In redox systems, interactions involve surfaces near the attachment sites of cofactors engaged in interprotein electron transfer (ET). Here we analyzed binding of cytochrome c2 to cytochrome bc1 by measuring paramagnetic relaxation enhancement (PRE) of spin label (SL) attached to cytochrome c2. PRE was exclusively induced by the iron atom of heme c1 of cytochrome bc1, which guaranteed that only the configurations with SL to heme c1 distances up to ∼30 Šwere detected. Changes in PRE were used to qualitatively and quantitatively characterize the binding. Our data suggest that at low ionic strength and under an excess of cytochrome c2 over cytochrome bc1, several cytochrome c2 molecules gather near the binding domain forming a "cloud" of molecules. When the cytochrome bc1 concentration increases, the cloud disperses to populate additional available binding domains. An increase in ionic strength weakens the attractive forces and the average distance between cytochrome c2 and cytochrome bc1 increases. The spatial arrangement of the protein complex at various ionic strengths is different. Above 150 mM NaCl the lifetime of the complexes becomes so short that they are undetectable. All together the results indicate that cytochrome c2 molecules, over the range of salt concentration encompassing physiological ionic strength, do not form stable, long-lived complexes but rather constantly collide with the surface of cytochrome bc1 and ET takes place coincidentally with one of these collisions.

Tyrosine triad at the interface between the Rieske iron-sulfur protein, cytochrome c1 and cytochrome c2 in the bc1 complex of Rhodobacter capsulatus.

A triad of tyrosine residues (Y152-154) in the cytochrome c(1) subunit (C1) of the Rhodobacter capsulatus cytochrome bc(1) complex (BC1) is ideally positioned to interact with cytochrome c(2) (C2). Mutational analysis of these three tyrosines showed that, of the three, Y154 is the most important, since its mutation to alanine resulted in significantly reduced levels, destabilization, and inactivation of BC1. A second-site revertant of this mutant that regained photosynthetic capacity was found to have acquired two further mutations-A181T and A200V. The Y152Q mutation did not change the spectral or electrochemical properties of C1, and showed wild-type enzymatic C2 reduction rates, indicating that this mutation did not introduce major structural changes in C1 nor affect overall activity. Mutations Y153Q and Y153A, on the other hand, clearly affect the redox properties of C1 (e.g. by lowering the midpoint potential as much as 117 mV in Y153Q) and the activity by 90% and 50%, respectively. A more conservative Y153F mutant on the other hand, behaves similarly to wild-type. This underscores the importance of an aromatic residue at position Y153, presumably to maintain close packing with P184, which modeling indicates is likely to stabilize the sixth heme ligand conformation.

Evidence from the structure and function of cytochromes c(2) that nonsulfur purple bacterial photosynthesis followed the evolution of oxygen respiration.

Cytochromes c(2) are the nearest bacterial homologs of mitochondrial cytochrome c. The sequences of the known cytochromes c(2) can be placed in two subfamilies based upon insertions and deletions, one subfamily is most like mitochondrial cytochrome c (the small C2s, without significant insertions and deletions), and the other, designated large C2, shares 3- and 8-residue insertions as well as a single-residue deletion. C2s generally function between cytochrome bc(1) and cytochrome oxidase in respiration (ca 80 examples known to date) and between cytochrome bc(1) and the reaction center in nonsulfur purple bacterial photosynthesis (ca 21 examples). However, members of the large C2 subfamily are almost always involved in photosynthesis (12 of 14 examples). In addition, the gene for the large C2 (cycA) is associated with those for the photosynthetic reaction center (pufBALM). We hypothesize that the insertions in the large C2s, which were already functioning in photosynthesis, allowed them to replace the membrane-bound tetraheme cytochrome, PufC, that otherwise mediates between the small C2 or other redox proteins and photosynthetic reaction centers. Based upon our analysis, we propose that the involvement of C2 in nonsulfur purple bacterial photosynthesis was a metabolic feature subsequent to the evolution of oxygen respiration.

Electron transfer from cytochrome c(2) to the reaction center: a transition state model for ionic strength effects due to neutral mutations.

Interprotein electron transfer plays an important role in biological energy conversion. In this work, the electron transfer reaction between cytochrome c(2) (cyt) and the reaction center (RC) was studied to determine the mechanisms coupling association and electron transfer. Previous studies have shown that mutation of hydrophobic residues in the reaction interface, particularly Tyr L162, changes the binding affinity and rates of electron transfer at low ionic strengths. In this study, the effect of ionic strength on the second-order electron transfer rate constant, k(2), between cyt c(2) and native or mutant RCs was examined. Mutations of hydrophobic and hydrogen bonding residues caused k(2) to decrease more rapidly with an increase in ionic strength. This change is explained with a transition state model by a switch from a diffusion-limited reaction in native RCs, where electron transfer occurs upon each binding event, to a fast exchange reaction in the Tyr L162 mutant, where dissociation occurs before electron transfer and k(2) depends upon the equilibrium between bound and free protein complexes. The difference in ionic strength dependence is attributed to a smaller effect of ionic strength on the energy of the transition state compared to the bound state due to larger distances between charged residues in the transition state. This model explains the faster dissociation rate at higher ionic strengths that may assist rapid turnover that is important for biological function. These results provide a quantitative model for coupling protein association with electron transfer and elucidate the role of short-range interactions in determining the rate of electron transfer.

Interaction between cytochrome c2 and the photosynthetic reaction center from Rhodobacter sphaeroides: role of interprotein hydrogen bonds in binding and electron transfer.

The role of short-range hydrogen bond interactions at the interface between electron transfer proteins cytochrome c(2) (cyt) and the reaction center (RC) from Rhodobacter sphaeroides was studied by mutation (to Ala) of RC residues Asn M187, Asn M188, and Gln L258 which form interprotein hydrogen bonds to cyt in the cyt-RC complex. The largest decrease in binding constant K(A) (8-fold) for a single mutation was observed for Asn M187, which forms an intraprotein hydrogen bond to the key residue Tyr L162 in the center of the contact region with a low solvent accessibility. Interaction between Asn M187 and Tyr L162 was also implicated in binding by double mutation of the two residues. The hydrogen bond mutations did not significantly change the second-order rate constant, k(2), indicating the mutations did not change the association rate for formation of the cyt-RC complex but increased the dissociation rate. The first-order electron transfer rate, k(e), for the cyt-RC complex was reduced by a factor of up to 4 (for Asn M187). The changes in k(e) were correlated with the changes in binding affinity but were not accompanied by increases in activation energy. We conclude that short-range hydrogen bond interactions contribute to the close packing of residues in the central contact region between the cyt and RC near Asn M187 and Tyr L162. The close packing contributes to fast electron transfer by increasing the rate of electronic coupling and contributes to the binding energy holding the cyt in position for times sufficient for electron transfer to occur.

A mechanistic and electrochemical study of the interaction between dimethyl sulfide dehydrogenase and its electron transfer partner cytochrome c2.

Dimethyl sulfide dehydrogenase isolated from the photosynthetic bacterium Rhodovulum sulfidophilum is a heterotrimeric enzyme containing a molybdenum cofactor at its catalytic site, as well as five iron-sulfur clusters and a heme b cofactor. It oxidizes dimethyl sulfide (DMS) to dimethyl sulfoxide in its native role and transfers electrons to the photochemical reaction center. There is genetic evidence that cytochrome c2 mediates this process, and the steady state kinetics experiments reported here demonstrated that cytochrome c2 accepts electrons from DMS dehydrogenase. At saturating concentrations of both substrate (DMS) and cosubstrate (cytochrome c2), Michaelis constants, KM,DMS and KM,cyt of 53 and 21 microM, respectively, were determined at pH 8. Further kinetic analysis revealed a "ping-pong" enzyme reaction mechanism for DMS dehydrogenase with its two reactants. Direct cyclic voltammetry of cytochrome c2 immobilized within a polymer film cast on a glassy carbon electrode revealed a reversible FeIII/II couple at +328 mV versus the normal hydrogen electrode at pH 8. The FeIII/II redox potential exhibited only minor pH dependence. In the presence of DMS dehydrogenase and DMS, the peak-shaped voltammogram of cytochrome c2 is transformed into a sigmoidal curve consistent with a steady-state (catalytic) reaction. The cytochrome c2 effectively mediates electron transfer between the electrode and DMS dehydrogenase during turnover and a significantly lower apparent electrochemical Michaelis constant K'M,DMS of 13(+/-1) microM was obtained. The pH optimum for catalytic DMS oxidation by DMS dehydrogenase with cytochrome c2 as the electron acceptor was found to be approximately 8.3.

Conformational control of the Q(A) to Q(B) electron transfer in bacterial reaction centers: evidence for a frozen conformational landscape below -25 degrees C.

The competition between the P(+)Q(A)(-) --> PQ(A) charge recombination (P, bacteriochlorophyll pair acting as primary photochemical electron donor) and the electron transfer to the secondary quinone acceptor Q(A)(-)Q(B) --> Q(A)Q(B)(-) (Q(A) and Q(B), primary and secondary electron accepting quinones) was investigated in chromatophores of Rb. capsulatus, varying the temperature down to -65 degrees C. The analysis of the flash-induced pattern for the formation of P(+)Q(A)Q(B)(-) shows that the diminished yield, when lowering the temperature, is not due to a homogeneous slowing of the rate constant k(AB) of the Q(A)(-)Q(B) --> Q(A)Q(B)(-) electron transfer but to a distribution of conformations that modulate the electron transfer rate over more than 3 orders of magnitude. This distribution appears "frozen", as no dynamic redistribution was observed over time ranges > 10 s (below -25 degrees C). The kinetic pattern was analyzed to estimate the shape of the distribution of k(AB), showing a bell-shaped band on the high rate side and a fraction of "blocked" reaction centers (RCs) with very slow k(AB). When the temperature is lowered, the high rate band moves to slower rate regions and the fraction of blocked RCs increases at the expense of the high rate band. The RCs that recombine from the P(+)Q(A)Q(B)(-) state appear temporarily converted to a state with rapid k(AB), indicating that the stabilized state described by Kleinfeld et al. (Biochemistry 1984, 23, 5780-5786) is still accessible at -60 degrees C.

Optimisation of NMR dynamic models I. Minimisation algorithms and their performance within the model-free and Brownian rotational diffusion spaces.

The key to obtaining the model-free description of the dynamics of a macromolecule is the optimisation of the model-free and Brownian rotational diffusion parameters using the collected R (1), R (2) and steady-state NOE relaxation data. The problem of optimising the chi-squared value is often assumed to be trivial, however, the long chain of dependencies required for its calculation complicates the model-free chi-squared space. Convolutions are induced by the Lorentzian form of the spectral density functions, the linear recombinations of certain spectral density values to obtain the relaxation rates, the calculation of the NOE using the ratio of two of these rates, and finally the quadratic form of the chi-squared equation itself. Two major topological features of the model-free space complicate optimisation. The first is a long, shallow valley which commences at infinite correlation times and gradually approaches the minimum. The most severe convolution occurs for motions on two timescales in which the minimum is often located at the end of a long, deep, curved tunnel or multidimensional valley through the space. A large number of optimisation algorithms will be investigated and their performance compared to determine which techniques are suitable for use in model-free analysis. Local optimisation algorithms will be shown to be sufficient for minimisation not only within the model-free space but also for the minimisation of the Brownian rotational diffusion tensor. In addition the performance of the programs Modelfree and Dasha are investigated. A number of model-free optimisation failures were identified: the inability to slide along the limits, the singular matrix failure of the Levenberg-Marquardt minimisation algorithm, the low precision of both programs, and a bug in Modelfree. Significantly, the singular matrix failure of the Levenberg-Marquardt algorithm occurs when internal correlation times are undefined and is greatly amplified in model-free analysis by both the grid search and constraint algorithms. The program relax ( http://www.nmr-relax.com ) is also presented as a new software package designed for the analysis of macromolecular dynamics through the use of NMR relaxation data and which alleviates all of the problems inherent within model-free analysis.

A beta-galactosidase-based bacterial two-hybrid system to assess protein-protein interactions in the correct cellular environment.

The vast majority of proteins functions in complex with one or more of the same or other proteins, indicating that protein-protein interactions play crucial roles in biology. Here, we present a beta-galactosidase reconstitution-based bacterial two-hybrid system in which two proteins of interest are fused to two non-functional but complementing beta-galactosidase truncations (Delta alpha and Delta omega). The level of complemented beta-galactosidase activity, driven by the protein-protein recognition between both non-beta-galactosidase parts of the chimeras, reflects whether or not the proteins of interest interact. Our approach was validated by reconfirming some well-established Escherichia coli cytoplasmic and membranous interactions, including well-chosen mutants, and providing the first in vivo evidence for the transient periplasmic interaction between Rhodobacter capsulatus cytochrome c2 and cytochrome c peroxidase. We demonstrated the major advantages of this in vivo two-hybrid technique: i) analyses of interactions are not limited to particular cellular compartments, ii) the potential of using the system in mutation-driven structure-function studies, and iii) the possibility of its application to transiently interacting proteins. These benefits demonstrate the relevance of the method as a powerful new tool in the broad spectrum of interaction assessment methods.

Plasmon waveguide resonance spectroscopic evidence for differential binding of oxidized and reduced Rhodobacter capsulatus cytochrome c2 to the cytochrome bc1 complex mediated by the conformation of the Rieske iron-sulfur protein.

The dissociation constants for the binding of Rhodobacter capsulatus cytochrome c2 and its K93P mutant to the cytochrome bc1 complex embedded in a phospholipid bilayer were measured by plasmon waveguide resonance spectroscopy in the presence and absence of the inhibitor stigmatellin. The reduced form of cytochrome c2 strongly binds to reduced cytochrome bc1 (Kd = 0.02 microM) but binds much more weakly to the oxidized form (Kd = 3.1 microM). In contrast, oxidized cytochrome c2 binds to oxidized cytochrome bc1 in a biphasic fashion with Kd values of 0.11 and 0.58 microM. Such a biphasic interaction is consistent with binding to two separate sites or conformations of oxidized cytochrome c2 and/or cytochrome bc1. However, in the presence of stigmatellin, we find that oxidized cytochrome c2 binds to oxidized cytochrome bc1 in a monophasic fashion with high affinity (Kd = 0.06 microM) and reduced cytochrome c2 binds less strongly (Kd = 0.11 microM) but approximately 30-fold more tightly than in the absence of stigmatellin. Structural studies with cytochrome bc1, with and without the inhibitor stigmatellin, have led to the proposal that the Rieske protein is mobile, moving between the cytochrome b and cytochrome c1 components during turnover. In one conformation, the Rieske protein binds near the heme of cytochrome c1, while the cytochrome c2 binding site is also near the cytochrome c1 heme but on the opposite side from the Rieske site, where cytochrome c2 cannot directly interact with Rieske. However, the inhibitor, stigmatellin, freezes the Rieske protein iron-sulfur cluster in a conformation proximal to cytochrome b and distal to cytochrome c1. We conclude from this that the dual conformation of the Rieske protein is primarily responsible for biphasic binding of oxidized cytochrome c2 to cytochrome c1. This optimizes turnover by maximizing binding of the substrate, oxidized cytochrome c2, when the iron-sulfur cluster is proximal to cytochrome b and minimizing binding of the product, reduced cytochrome c2, when it is proximal to cytochrome c1.

Cytochrome c(2) Exit Strategy: Dissociation Studies and Evolutionary Implications.

Small, water-soluble, type c cytochromes form a transient network connecting major bioenergetic membrane protein complexes in both photosynthesis and respiration. In the photosynthesis cycle of Rhodobacter sphaeroides, cytochrome c2 (cyt c2) docks to the reaction center (RC), undergoes electron transfer, and exits for the cytochrome bc1 complex. Translations of cyt c2 about the RC-cyt c2 docking interface and surrounding membrane reveal possible exit pathways. A pathway at a minimal elevation allowed by the architecture of the RC is analyzed using both an all-atom steered molecular dynamics simulation of the RC-cyt c2 complex and a bioinformatic analysis of the structures and sequences of cyt c. The structure-based phylogenetic analysis allows for the identification of structural elements that have evolved to satisfy the requirements of having multiple functional partners. The patterns of evolutionary variation obtained from the phylogenetic analysis of both docking partners of cyt c2 reveal conservation of key residues involved in the interaction interfaces that would be candidates for further experimental studies. Additionally, using the molecular mechanics Poisson-Boltzmann surface area method we calculate that the binding free energy of reduced cyt c2 to the RC is nearly 6 kcal/mol more favorable than with oxidized cyt c2. The redox-dependent variations lead to changes in structural flexibility, behavior of the interfacial water molecules, and eventually changes in the binding free energy of the complex.

Structures of various cytochromes c evaluated from the redox behaviors using the optically active Co(III) complex-modified Au electrode.

Electrochemical studies of three c-type cytochromes (cyt c from horse heart, cyt c(2) from Rhodospirillum rubrum, and cyt c(553) from Alcaligenes xylosoxidans GIFU 1051) were performed by using the optically active Co(III) complex-modified Au electrode. Three cytochromes gave different redox behaviors reflected from the respective structural information. From relationship between their redox behaviors and their structural characteristics, we evaluated the solution structure around heme center of cyt c(553).

Local stability of Rhodobacter capsulatus cytochrome c2 probed by solution phase hydrogen/deuterium exchange and mass spectrometry.

The hydrogen/deuterium exchange kinetics of Rhodobacter capsulatus cytochrome c2 have been determined using mass spectrometry. As expected, the relative domain stability was generally similar to that of the cytochrome c2 structural homolog, horse heart cytochrome c, but we were able to find evidence to support the presence of a second, small beta-sheet not found in the horse cytochrome, which stabilizes a structural region dominated by Omega loops. Importantly, we find that the so-called hinge region, comprised of 15 amino acids, which include the methionine sixth heme ligand (M96), is destabilized on oxidation, and this destabilization is propagated to a portion of the second Omega loop, most likely through perturbation of two hydrogen bonds that couple these two domains in the three dimensional structure. The mutation of a lysine at position 93 to proline amplifies the destabilization observed on oxidation of the wild-type cytochrome c2 and results in further destabilization observed in regions 52-60, 75-82, and 83-97. This suggests that hydrogen bond interactions involving two bound waters, the T94 hydroxyl, the front heme propionate and the Y75 hydroxyl, are significantly compromised upon mutation. In summary, these observations are consistent with the approximately 20-fold increase in the movement of the hinge away from the heme face in the oxidized cytochrome c2 as determined by ligand binding kinetics. Thus, H/D exchange kinetics can be used to identify relatively subtle structural features and at least in some cases facilitate the understanding of the structural basis of the dynamic properties of proteins.

Reconstruction of a kinetic model of the chromatophore vesicles from Rhodobacter sphaeroides.

We present a molecular model of a chromatophore vesicle from Rhodobacter sphaeroides. These vesicles are ideal benchmark systems for molecular and systemic simulations, because they have been well studied, they are small, and they are naturally separated from their cellular environment. To set up a photosynthetic chain working under steady-state conditions, we compiled from the experimental literature the specific activities and geometries that have been determined for their constituents. This data then allowed defining the stoichiometries for all membrane proteins. This article contains the kinetic part of the reconstructed model, while the spatial reconstruction is presented in a companion article. By considering the transport properties of the Cytochrome c(2) and ubiquinone pools, we show that their size and oxidation states allow for an efficient buffering of the statistical fluctuations that arise from the small size of the vesicles. Stoichiometric and kinetic considerations indicate that a typical chromatophore vesicle of Rb. sphaeroides with a diameter of 45 nm should contain approximately five bc(1) monomers.

Electron transfer to nitrite reductase of Rhodobacter sphaeroides 2.4.3: examination of cytochromes c2 and cY.

The role of cytochrome c(2), encoded by cycA, and cytochrome c(Y), encoded by cycY, in electron transfer to the nitrite reductase of Rhodobacter sphaeroides 2.4.3 was investigated using both in vivo and in vitro approaches. Both cycA and cycY were isolated, sequenced and insertionally inactivated in strain 2.4.3. Deletion of either gene alone had no apparent effect on the ability of R. sphaeroides to reduce nitrite. In a cycA-cycY double mutant, nitrite reduction was largely inhibited. However, the expression of the nitrite reductase gene nirK from a heterologous promoter substantially restored nitrite reductase activity in the double mutant. Using purified protein, a turnover number of 5 s(-1) was observed for the oxidation of cytochrome c(2) by nitrite reductase. In contrast, oxidation of c(Y) only resulted in a turnover of approximately 0.1 s(-1). The turnover experiments indicate that c(2) is a major electron donor to nitrite reductase but c(Y) is probably not. Taken together, these results suggest that there is likely an unidentified electron donor, in addition to c(2), that transfers electrons to nitrite reductase, and that the decreased nitrite reductase activity observed in the cycA-cycY double mutant probably results from a change in nirK expression.

Interactions between cytochrome c2 and the photosynthetic reaction center from Rhodobacter sphaeroides: the cation-pi interaction.

The cation-pi interaction between positively charged and aromatic groups is a common feature of many proteins and protein complexes. The structure of the complex between cytochrome c(2) (cyt c(2)) and the photosynthetic reaction center (RC) from Rhodobacter sphaeroides exhibits a cation-pi complex formed between Arg-C32 on cyt c(2) and Tyr-M295 on the RC [Axelrod, H. L., et al. (2002) J. Mol. Biol. 319, 501-515]. The importance of the cation-pi interaction for binding and electron transfer was studied by mutating Tyr-M295 and Arg-C32. The first- and second-order rates for electron transfer were not affected by mutating Tyr-M295 to Ala, indicating that the cation-pi complex does not greatly affect the association process or structure of the state active in electron transfer. The dissociation constant K(D) showed a greater increase when Try-M295 was replaced with nonaromatic Ala (3-fold) as opposed to aromatic Phe (1.2-fold), which is characteristic of a cation-pi interaction. Replacement of Arg-C32 with Ala increased K(D) (80-fold) largely due to removal of electrostatic interactions with negatively charged residues on the RC. Replacement with Lys increased K(D) (6-fold), indicating that Lys does not form a cation-pi complex. This specificity for Arg may be due to a solvation effect. Double mutant analysis indicates an interaction energy between Tyr-M295 and Arg-C32 of approximately -24 meV (-0.6 kcal/mol). This energy is surprisingly small considering the widespread occurrence of cation-pi complexes and may be due to the tradeoff between the favorable cation-pi binding energy and the unfavorable desolvation energy needed to bury Arg-C32 in the short-range contact region between the two proteins.

Structural and functional studies on the tetraheme cytochrome subunit and its electron donor proteins: the possible docking mechanisms during the electron transfer reaction.

The photosynthetic reaction centers (RCs) classified as the group II possess a peripheral cytochrome (Cyt) subunit, which serves as the electron mediator to the special-pair. In the cycle of the photosynthetic electron transfer reactions, the Cyt subunit accepts electrons from soluble electron carrier proteins, and re-reduces the photo-oxidized special-pair of the bacteriochlorophyll. Physiologically, high-potential cytochromes such as the cytochrome c2 and the high-potential iron-sulfur protein (HiPIP) function as the electron donors to the Cyt subunit. Most of the Cyt subunits possess four heme c groups, and it was unclear which heme group first accepts the electron from the electron donor. The most distal heme to the special-pair, the heme-1, has a lower redox potential than the electron donors, which makes it difficult to understand the electron transfer mechanism mediated by the Cyt subunit. Extensive mutagenesis combined with kinetic studies has made a great contribution to our understanding of the molecular interaction mechanisms, and has demonstrated the importance of the region close to the heme-1 in the electron transfer. Moreover, crystallographic studies have elucidated two high-resolution three-dimensional structures for the RCs containing the Cyt subunit, the Blastochloris viridis and Thermochromatium tepidum RCs, as well as the structures of their electron donors. An examination of the structural data also suggested that the binding sites for both the cytochrome c2 and the HiPIP are located adjacent to the solvent-accessible edge of the heme-1. In addition, it is also indicated by the structural and biochemical data that the cytochrome c2 and the HiPIP dock with the Cyt subunit by c2 is recognized through electrostatic interactions while hydrophobic interactions are important in the HiPIP docking.

The structure and function of the cytochrome c2: reaction center electron transfer complex from Rhodobacter sphaeroides.

In the photosynthetic bacterium, Rhodobacter sphaeroides, the mobile electron carrier, cytochrome c2 (cyt c2) transfers an electron from reduced heme to the photooxidized bacteriochlorophyll dimer in the membrane bound reaction center (RC) as part of the light induced cyclic electron transfer chain. A complex between these two proteins that is active in electron transfer has been crystallized and its structure determined by X-ray diffraction. The structure of the cyt:RC complex shows the cyt c2 (cyt c2) positioned at the center of the periplasmic surface of the RC. The exposed heme edge from cyt c2 is in close tunneling contact with the electron acceptor through an intervening bridging residue, Tyr L162 located on the RC surface directly above the bacteriochlorophyll dimer. The binding interface between the two proteins can be divided into two regions: a short-range interaction domain and a long-range interaction domain. The short-range domain includes residues immediately surrounding the tunneling contact region around the heme and Tyr L162 that display close intermolecular contacts optimized for electron transfer. These include a small number of hydrophobic interactions, hydrogen bonds and a pi-cation interaction. The long-range interaction domain consists of solvated complementary charged residues; positively charged residues from the cyt and negatively charged residues from the RC that provide long range electrostatic interactions that can steer the two proteins into position for rapid association.

Interprotein electron transfer from cytochrome c2 to photosynthetic reaction center: tunneling across an aqueous interface.

Interprotein electron transfer (ET) reactions play an important role in biological energy conversion processes. One of these reactions, the ET between cytochrome c(2) (cyt) and reaction center from photosynthetic bacteria, is the focus of this theoretical study. The changes in the ET rate constant at fixed distances during the association process were calculated as the cyt moved from the electrostatically stabilized encounter complex to the bound state having short range van der Waals contacts in the tunneling region. Multiple conformations of the protein were generated by molecular dynamics simulations including explicit water molecules. For each of these conformations, the ET rate was calculated by using the Pathways model. The ET rate increased smoothly as the cyt approached from the encounter complex to the bound state, with a tunneling decay factor beta = 1.1 A(-1). This relatively efficient coupling between redox centers is due to the ability of interfacial water molecules to form multiple strong hydrogen bonding pathways connecting tunneling pathways on the surfaces of the two proteins. The ET rate determined for the encounter complex ensemble of states is only about a factor of 100 slower than that of the bound state (tau = 100 micros, compared with 1 micros), because of fluctuations of the cyt within the encounter complex ensemble through configurations having strong tunneling pathways. The ET rate for the encounter complex is in agreement with rates observed in mutant reaction centers modified to remove shortrange hydrophobic interactions, suggesting that in this case, ET occurs within the solvent-separated, electrostatically stabilized encounter complex.