Estimating the effect of protein dynamics on electron transfer to the special pair in the photosynthetic reaction center

Structural and spectropotentiometric analysis of Blastochloris viridis heterodimer mutant reaction center

Heterodimer mutant reaction centers (RCs) of Blastochloris viridis were crystallized using microfluidic technology. In this mutant, a leucine residue replaced the histidine residue which had acted as a fifth ligand to the bacteriochlorophyll (BChl) of the primary electron donor dimer M site (HisM200). With the loss of the histidine-coordinated Mg, one bacteriochlorophyll of the special pair was converted into a bacteriopheophytin (BPhe), and the primary donor became a heterodimer supermolecule. The crystals had dimensions 400 x 100 x 100 microm, belonged to space group P4(3)2(1)2, and were isomorphous to the ones reported earlier for the wild type (WT) strain. The structure was solved to a 2.5 A resolution limit. Electron-density maps confirmed the replacement of the histidine residue and the absence of Mg. Structural changes in the heterodimer mutant RC relative to the WT included the absence of the water molecule that is typically positioned between the M side of the primary donor and the accessory BChl, a slight shift in the position of amino acids surrounding the site of the mutation, and the rotation of the M194 phenylalanine. The cytochrome subunit was anchored similarly as in the WT and had no detectable changes in its overall position. The highly conserved tyrosine L162, located between the primary donor and the highest potential heme C(380), revealed only a minor deviation of its hydroxyl group. Concomitantly to modification of the BChl molecule, the redox potential of the heterodimer primary donor increased relative to that of the WT organism (772 mV vs. 517 mV). The availability of this heterodimer mutant and its crystal structure provides opportunities for investigating changes in light-induced electron transfer that reflect differences in redox cascades.

Conformationally averaged score functions for electronic propagation in proteins

We explore the influence of conformational dynamics on protein-mediated electron donor-acceptor interactions. We introduce a thermally averaged score function to characterize electronic propagation from redox cofactors into the protein and solvent. The score function is explored for myoglobin at the extended-Hückel level, and the results are compared with those of simpler models. The conformationally averaged quantum results are consistent with the empirical analysis of the Pathways model. Notably, subtle effects of quantum interference among multiple coupling pathways that arise in static structures are largely averaged out when protein thermal motion is included. Propagation through bulk water near the single-protein interface decays rapidly with distance.

Simulation of Electron Transfer between Cytochrome c2 and the Bacterial Photosynthetic Reaction Center: Brownian Dynamics Analysis of the Native Proteins and Double Mutants

Electron transfer is essential for bacterial photosynthesis which converts light energy into chemical energy. This paper theoretically studies the interprotein electron transfer from cytochrome c(2) of Rhodobacter capsulatus to the photosynthetic reaction center of Rhodobacter sphaeroides in native and mutated systems. Brownian dynamics is used with an exponential distance-dependent electron-transfer rate model to compute bimolecular rate constants, which are consistent with experimental data when reasonable prefactors and decay constants are used. Interestingly, switching of the reaction mechanism from the diffusion-controlled limit in the native proteins to the activation-controlled limit in one of the mutants (DK(L261)/KE(C99)) was found. We also predict that the second-order rate for the native reaction center/cytochrome c(2) system will decrease with increasing ionic strength, a characteristic of electrostatically controlled docking.

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.

Interference, Fluctuation, and Alternation of Electron Tunneling in Protein Media. 1. Two Tunneling Routes in Photosynthetic Reaction Center Alternate Due to Thermal Fluctuation of Protein Conformation

Electron tunneling routes for the electron transfer from the bacteriopheophytin anion to the primary quinone in the bacterial photosynthetic reaction center of Rhodobactor sphaeroides are investigated by a combined method of molecular dynamics simulations for the protein conformation fluctuation and quantum chemical calculations for the electronic states of the donor, acceptor, and protein medium. The analysis of the tunneling route is made by mapping interatomic electron tunneling currents for each protein conformation. We found that there are two dominant routes mainly passing through Trp(M252) (Trp route) or mainly passing through Met(M218) (Met route). Actual electron tunneling pathways alternate between the two routes, depending on the protein conformation which varies with time. When either the Trp route or the Met route dominates, the electron tunneling matrix element /T(DA)/ becomes large. When both the Trp route and the Met route dominate, /T(DA)/ becomes very small due to the destructive interference of the electron tunneling currents between the two routes. We found that a linear relationship exists between the value of /T(DA)/ and the inverse of the degree of destructive interference Q for a wide range of values (ca. 3-10(3) for Q). A similar relationship was also found previously for electron transfer in ruthenium-modified azurins, suggesting that this relationship holds true in general. From these results, we are led to the conclusion that /T(DA)/ cannot exceed a maximum value at Q = 1, even if much variation of /T(DA)/ happens due to the fluctuation of protein conformation. We also conclude that the property of the electron transfer alternates between constructive and destructive interference, due to the fluctuation of protein conformation. It is impossible to keep a system in either constructive or destructive interference because thermal fluctuation of protein conformation takes place.

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.

Transition state and encounter complex for fast association of cytochrome c2 with bacterial reaction center

Electrostatic interactions strongly enhance the electron transfer reaction between cytochrome (Cyt) c(2) and reaction center (RC) from photosynthetic bacteria, yielding a second-order rate constant, k(2) approximately 10(9) s(-1).M(-1), close to the diffusion limit. The proposed mechanism involves an encounter complex (EC) stabilized by electrostatic interactions, followed by a transition state (TS), leading to the bound complex active in electron transfer. The effect of electrostatic interactions was previously studied by Tetreault et al. [Tetreault, M., Cusanovich, M., Meyer, T., Axelrod, H. & Okamura, M. Y. (2002) Biochemistry 41, 5807-5815] by measuring k(2) for RC and Cyt molecules with modified charged residues at the binding interface. The present work is a computational analysis of this kinetic study to determine the ensemble of configurations of the TS and EC. Changes in the TS energies due to different mutations were compared with differences in the calculated electrostatic energies for a wide range of Cyt/RC configurations. The TS ensemble, obtained from structures having the highest correlation coefficients in the comparison with experimental data, has the Cyt displaced by approximately 10 A from its position in x-ray crystal structure, close to the average position of the EC ensemble, with strong electrostatic interactions between Cyt on the M subunit side of the RC surface. The heme of the Cyt is oriented toward Tyr L162 on the RC, the tunneling contact in the bound final state on the RC. The similarity between the structures of the EC, TS, and bound state can account for the rapid rate of association responsible for fast diffusion-controlled electron transfer.

Interactions between cytochrome c2 and photosynthetic reaction center from Rhodobacter sphaeroides: changes in binding affinity and electron transfer rate due to mutation of interfacial hydrophobic residues are strongly correlated

The structure of the complex between cytochrome c(2) (cyt) and the photosynthetic reaction center (RC) from Rhodobacter sphaeroides shows contacts between hydrophobic residues Tyr L162, Leu M191, and Val M192 on the RC and the surface of the cyt [Axelrod et al. (2002) J. Mol. Biol. 319, 501-515]. The role of these hydrophobic residues in binding and electron transfer was investigated by replacing them with Ala and other residues. Mutations of the hydrophobic residues generally resulted in relatively small changes in the second-order electron-transfer rate k(2) (Brönsted coefficient, alpha( )()= 0.15 +/- 0.05) indicating that the transition state for association occurs before short-range hydrophobic contacts are established. Larger changes in k(2), found in some cases, were attributed to a change in the second-order mechanism from a diffusion controlled regime to a rapidly reversible binding regime. The association constant, K(A), of the cyt and the rate of electron transfer from the bound cyt, k(e), were both decreased by mutation. Replacement of Tyr L162, Leu M191, or Val M192 by Ala decreased K(A) and k(e) by factors of 130, 10, 0.6, and 120, 9, 0.6, respectively. The largest changes were obtained by mutation of Tyr L162, showing that this residue plays a key role in both binding and electron transfer. The binding affinity, K(A), and electron-transfer rate, k(e) were strongly correlated, showing that changes of hydrophobic residues affect both binding and electron transfer. This correlation suggests that changes in distance across hydrophobic interprotein contacts have similar effects on both electron tunneling and binding interactions.

X-ray structure determination of the cytochrome c2: reaction center electron transfer complex from Rhodobacter sphaeroides

In the photosynthetic bacterium Rhodobacter sphaeroides, a water soluble cytochrome c2 (cyt c2) is the electron donor to the reaction center (RC), the membrane-bound pigment-protein complex that is the site of the primary light-induced electron transfer. To determine the interactions important for docking and electron transfer within the transiently bound complex of the two proteins, RC and cyt c2 were co-crystallized in two monoclinic crystal forms. Cyt c2 reduces the photo-oxidized RC donor (D+), a bacteriochlorophyll dimer, in the co-crystals in approximately 0.9 micros, which is the same time as measured in solution. This provides strong evidence that the structure of the complex in the region of electron transfer is the same in the crystal and in solution. X-ray diffraction data were collected from co-crystals to a maximum resolution of 2.40 A and refined to an R-factor of 22% (R(free)=26%). The structure shows the cyt c2 to be positioned at the center of the periplasmic surface of the RC, with the heme edge located above the bacteriochlorophyll dimer. The distance between the closest atoms of the two cofactors is 8.4 A. The side-chain of Tyr L162 makes van der Waals contacts with both cofactors along the shortest intermolecular electron transfer pathway. The binding interface can be divided into two domains: (i) A short-range interaction domain that includes Tyr L162, and groups exhibiting non-polar interactions, hydrogen bonding, and a cation-pi interaction. This domain contributes to the strength and specificity of cyt c2 binding. (ii) A long-range, electrostatic interaction domain that contains solvated complementary charges on the RC and cyt c2. This domain, in addition to contributing to the binding, may help steer the unbound proteins toward the right conformation.

Dynamically controlled protein tunneling paths in photosynthetic reaction centers

Marcus theory has explained how thermal nuclear motions modulate the energy gap between donor and acceptor sites in protein electron transfer reactions. Thermal motions, however, may also modulate electron tunneling between these reactions. Here we identify a new mechanism of nuclear dynamics amplification that plays a central role when interference among the dominant tunneling pathway tubes is destructive. In these cases, tunneling takes place in protein conformations far from equilibrium that minimize destructive interference. As an example, we demonstrate how this dynamical amplification mechanism affects certain reaction rates in the photosynthetic reaction center and therefore may be critical for biological function.

Electron transfer mechanisms

The tunneling pathway framework description of protein electron transfer reactions has prompted a lively discussion of how structure and evolution influence electron transfer rates. Recent protein and model system experiments, performed in solution and in organized media, are providing answers. The molecular mechanisms of DNA electron transfer reactions are being probed as well with new theoretical and experimental strategies.