Delocalized electronic excitations and their role in directional charge transfer in the reaction center of Rhodobacter sphaeroides.

In purple bacteria, the fundamental charge-separation step that drives the conversion of radiation energy into chemical energy proceeds along one branch-the A branch-of a heterodimeric pigment-protein complex, the reaction center. Here, we use first principles time-dependent density functional theory (TDDFT) with an optimally-tuned range-separated hybrid functional to investigate the electronic and excited-state structure of the six primary pigments in the reaction center of Rhodobacter sphaeroides. By explicitly including amino-acid residues surrounding these six pigments in our TDDFT calculations, we systematically study the effect of the protein environment on energy and charge-transfer excitations. Our calculations show that a forward charge transfer into the A branch is significantly lower in energy than the first charge transfer into the B branch, in agreement with the unidirectional charge transfer observed experimentally. We further show that the inclusion of the protein environment redshifts this excitation significantly, allowing for energy transfer from the coupled Qx excitations. Through analysis of transition and difference densities, we demonstrate that most of the Q-band excitations are strongly delocalized over several pigments and that both their spatial delocalization and charge-transfer character determine how strongly affected they are by thermally-activated molecular vibrations. Our results suggest a mechanism for charge-transfer in this bacterial reaction center and pave the way for further first-principles investigations of the interplay between delocalized excited states, vibronic coupling, and the role of the protein environment in this and other complex light-harvesting systems.

The Charge Distribution on a Protein Surface Determines Whether Productive or Futile Encounter Complexes Are Formed.

Protein complex formation depends strongly on electrostatic interactions. The distribution of charges on the surface of redox proteins is often optimized by evolution to guide recognition and binding. To test the degree to which the electrostatic interactions between cytochrome c peroxidase (CcP) and cytochrome c (Cc) are optimized, we produced five CcP variants, each with a different charge distribution on the surface. Monte Carlo simulations show that the addition of negative charges attracts Cc to the new patches, and the neutralization of the charges in the regular, stereospecific binding site for Cc abolishes the electrostatic interactions in that region entirely. For CcP variants with the charges in the regular binding site intact, additional negative patches slightly enhance productive complex formation, despite disrupting the optimized charge distribution. Removal of the charges in the regular binding site results in a dramatic decrease in the complex formation rate, even in the presence of highly negative patches elsewhere on the surface. We conclude that additional charge patches can result in either productive or futile encounter complexes, depending on whether negative residues are located also in the regular binding site.

Self-interaction correction, electrostatic, and structural influences on time-dependent density functional theory excitations of bacteriochlorophylls from the light-harvesting complex 2.

First-principles calculations offer the chance to obtain a microscopic understanding of light-harvesting processes. Time-dependent density functional theory can have the computational efficiency to allow for such calculations. However, the (semi-)local exchange-correlation approximations that are computationally most efficient fail to describe charge-transfer excitations reliably. We here investigate whether the inexpensive average density self-interaction correction (ADSIC) remedies the problem. For the systems that we study, ADSIC is even more prone to the charge-transfer problem than the local density approximation. We further explore the recently reported finding that the electrostatic potential associated with the chromophores' protein environment in the light-harvesting complex 2 beneficially shifts spurious excitations. We find a great sensitivity on the chromophores' atomistic structure in this problem. Geometries obtained from classical molecular dynamics are more strongly affected by the spurious charge-transfer problem than the ones obtained from crystallography or density functional theory. For crystal structure geometries and density-functional theory optimized ones, our calculations confirm that the electrostatic potential shifts the spurious excitations out of the energetic range that is most relevant for electronic coupling.

Efficient Encounter Complex Formation and Electron Transfer to Cytochrome†c Peroxidase with an Additional, Distant Electrostatic Binding Site.

Electrostatic interactions can strongly increase the efficiency of protein complex formation. The charge distribution in redox proteins is often optimized to steer a redox partner to the electron transfer active binding site. To test whether the optimized distribution is more important than the strength of the electrostatic interactions, an additional negative patch was introduced on the surface of cytochrome c peroxidase, away from the stereospecific binding site, and its effect on the encounter complex as well as the rate of complex formation was determined. Monte Carlo simulations and paramagnetic relaxation enhancement NMR experiments indicate that the partner, cytochrome c, interacts with the new patch. Unexpectedly, the rate of the active complex formation was not reduced, but rather slightly increased. The findings support the idea that for efficient protein complex formation the strength of the electrostatic interaction is more critical than an optimized charge distribution.

Assessing density functional theory in real-time and real-space as a tool for studying bacteriochlorophylls and the light-harvesting complex 2.

We use real-time density functional theory on a real-space grid to calculate electronic excitations of bacteriochlorophyll chromophores of the light-harvesting complex 2 (LH2). Comparison with Gaussian basis set calculations allows us to assess the numerical trust range for computing electron dynamics in coupled chromophores with both types of techniques. Tuned range-separated hybrid calculations for one bacteriochlorophyll as well as two coupled ones are used as a reference against which we compare results from the adiabatic time-dependent local density approximation (TDLDA). The tuned range-separated hybrid calculations lead to a qualitatively correct description of the electronic excitations and couplings. They allow us to identify spurious charge-transfer excitations that are obtained with the TDLDA. When we take into account the environment that the LH2 protein complex forms for the bacteriochlorophylls, we find that it substantially shifts the energy of the spurious charge-transfer excitations, restoring a qualitatively correct electronic coupling of the dominant excitations also for TDLDA.

The dynamic complex of cytochrome c6 and cytochrome f studied with paramagnetic NMR spectroscopy.

The rapid transfer of electrons in the photosynthetic redox chain is achieved by the formation of short-lived complexes of cytochrome b6f with the electron transfer proteins plastocyanin and cytochrome c6. A balance must exist between fast intermolecular electron transfer and rapid dissociation, which requires the formation of a complex that has limited specificity. The interaction of the soluble fragment of cytochrome f and cytochrome c6 from the cyanobacterium Nostoc sp. PCC 7119 was studied using NMR spectroscopy and X-ray diffraction. The crystal structures of wild type, M58H and M58C cytochrome c6 were determined. The M58C variant is an excellent low potential mimic of the wild type protein and was used in chemical shift perturbation and paramagnetic relaxation NMR experiments to characterize the complex with cytochrome f. The interaction is highly dynamic and can be described as a pure encounter complex, with no dominant stereospecific complex. Ensemble docking calculations and Monte-Carlo simulations suggest a model in which charge-charge interactions pre-orient cytochrome c6 with its haem edge toward cytochrome f to form an ensemble of orientations with extensive contacts between the hydrophobic patches on both cytochromes, bringing the two haem groups sufficiently close to allow for rapid electron transfer. This model of complex formation allows for a gradual increase and decrease of the hydrophobic interactions during association and dissociation, thus avoiding a high transition state barrier that would slow down the dissociation process.

An ensemble of rapidly interconverting orientations in electrostatic protein-peptide complexes characterized by NMR spectroscopy.

Protein complex formation involves an encounter state in which the proteins are associated in a nonspecific manner and often stabilized by interactions between charged surface patches. Such patches are thought to bind in many different orientations with similar affinity. To obtain experimental evidence for the dynamics in encounter complexes, a model was created using the electron transfer protein plastocyanin and short charged peptides. Three plastocyanins with distinct surface charge distributions were studied. The experimental results from chemical shift perturbations, paramagnetic relaxation enhancement (PRE) NMR, and theoretical results from Monte Carlo simulations indicate the presence of multiple binding orientations that interconvert quickly and are dominated by long-range charge interactions. The PRE data also suggest the presence of highly transient orientations stabilized by short-range interactions.

Loss of electrostatic interactions causes increase of dynamics within the plastocyanin-cytochrome f complex.

Recent studies on the electron transfer complex formed by cytochrome f and plastocyanin from Nostoc revealed that both hydrophobic and electrostatic interactions play a role in the process of complex formation. To study the balance between these two types of interactions in the encounter and the final state, the complex between plastocyanin from Phormidium laminosum and cytochrome f from Nostoc sp. PCC 7119 was investigated using NMR spectroscopy and Monte Carlo docking. Cytochrome f has a highly negative charge. Phormidium plastocyanin is similar to that from Nostoc, but the net charge of the protein is negative rather than positive. NMR titrations of Zn-substituted Phormidium plastocyanin and Nostoc cytochrome f indicated that a complex with an affinity intermediate between those of the Nostoc and Phormidium complexes is formed. Plastocyanin was found in a head-on orientation, as determined using pseudocontact shifts, similar to that in the Phormidium complex, in which the hydrophobic patch represents the main site of interaction on plastocyanin. However, the interaction in the cross-complex is dependent on electrostatics, similar to that in the Nostoc complex. The negative charge of plastocyanin decreases, but not abolishes, the attraction to cytochrome f, resulting in the formation of a more diffuse encounter complex than in the Nostoc case, as could be determined using paramagnetic relaxation spectroscopy. This work illustrates the subtle interplay of electrostatic and hydrophobic interactions in the formation of transient protein complexes. The results are discussed in the context of a model for association on the basis of hydrophobic contacts in the encounter state.

Role of hydrophobic interactions in the encounter complex formation of the plastocyanin and cytochrome f complex revealed by paramagnetic NMR spectroscopy.

Protein complex formation is thought to be at least a two-step process, in which the active complex is preceded by the formation of an encounter complex. The interactions in the encounter complex are usually dominated by electrostatic forces, whereas the active complex is also stabilized by noncovalent short-range forces. Here, the complex of cytochrome f and plastocyanin, electron-transfer proteins involved in photosynthesis, was studied using paramagnetic relaxation NMR spectroscopy. Spin labels were attached to cytochrome f, and the relaxation enhancements of plastocyanin nuclei were measured, demonstrating that a large part of the cytochrome f surface area is sampled by plastocyanin. In contrast, plastocyanin is always oriented with its hydrophobic patch toward cytochrome f. The complex was visualized using ensemble docking, showing that the encounter complex is stabilized by hydrophobic as well as electrostatic interactions. The results suggest a model of electrostatic preorientation before the proteins make contact, followed by the formation of an encounter complex that rapidly leads to electron-transfer active conformations by gradual increase of the overlap of nonpolar surface areas on cytochrome f and plastocyanin. In this model the distinction between the encounter and active complexes vanishes, at least in the case of electron-transfer complexes, which do not require a high degree of specificity.

Understanding Primary Charge Separation in the Heliobacterial Reaction Center.

The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.

Enhancing the population of the encounter complex affects protein complex formation efficiency.

Optimal charge distribution is considered to be important for efficient formation of protein complexes. Electrostatic interactions guide encounter complex formation that precedes the formation of an active protein complex. However, disturbing the optimized distribution by introduction of extra charged patches on cytochrome c peroxidase does not lead to a reduction in productive encounters with its partner cytochrome c. To test whether a complex with a high population of encounter complex is more easily affected by suboptimal charge distribution, the interactions of cytochrome c mutant R13A with wild-type cytochrome c peroxidase and a variant with an additional negative patch were studied. The complex of the peroxidase and cytochrome c R13A was reported to have an encounter state population of 80%, compared to 30% for the wild-type cytochrome c. NMR analysis confirms the dynamic nature of the interaction and demonstrates that the mutant cytochrome c samples the introduced negative patch. Kinetic experiments show that productive complex formation is fivefold to sevenfold slower at moderate and high ionic strength values for cytochrome c R13A but the association rate is not affected by the additional negative patch on cytochrome c peroxidase, showing that the total charge on the protein surface can compensate for less optimal charge distribution. At low ionic strength (44 mm), the association with the mutant cytochrome c reaches the same high rates as found for wild-type cytochrome c, approaching the diffusion limit.

Structural and Biophysical Analysis of the Phytochelatin-Synthase-Like Enzyme from Nostoc sp. Shows That Its Protease Activity is Sensitive to the Redox State of the Substrate.

Phytochelatins (PCs) are nonribosomal thiol-rich oligopeptides synthetized from glutathione (GSH) in a γ-glutamylcysteinyl transpeptidation reaction catalyzed by PC synthases (PCSs). Ubiquitous in plant and present in some invertebrates, PCSs are involved in metal detoxification and homeostasis. The PCS-like enzyme from the cyanobacterium Nostoc sp. (NsPCS) is considered to be an evolutionary precursor enzyme of genuine PCSs because it shows sufficient sequence similarity for homology to the catalytic domain of the eukaryotic PCSs and shares the peptidase activity consisting in the deglycination of GSH. In this work, we investigate the catalytic mechanism of NsPCS by combining structural, spectroscopic, thermodynamic, and theoretical techniques. We report several crystal structures of NsPCS capturing different states of the catalyzed chemical reaction: (i) the structure of the wild-type enzyme (wt-NsPCS); (ii) the high-resolution structure of the γ-glutamyl-cysteine acyl-enzyme intermediate (acyl-NsPCS); and (iii) the structure of an inactive variant of NsPCS, with the catalytic cysteine mutated into serine (C70S-NsPCS). We characterize NsPCS as a relatively slow enzyme whose activity is sensitive to the redox state of the substrate. Namely, NsPCS is active with reduced glutathione (GSH), but is inhibited by oxidized glutathione (GSSG) because the cleavage product is not released from the enzyme. Our biophysical analysis led us to suggest that the biological function of NsPCS is being a part of a redox sensing system. In addition, we propose a mechanism how PCS-like enzymes may have evolved toward genuine PCS enzymes.

Investigating Primary Charge Separation in the Reaction Center of Heliobacterium modesticaldum.

We compute the primary charge separation step in the homodimeric reaction center (RC) of Heliobacterium modesticaldum from first principles. Using time-dependent density functional theory with the optimally tuned range-separated hybrid functional ωPBE, we calculate the excitations of a system comprising the special pair, the adjacent accessory bacteriochlorophylls, and the most relevant parts of the surrounding protein environment. The structure of the excitation spectrum can be rationalized from coupling of the individual bacteriochlorophyll pigments similar to molecular J- and H-aggregates. We find excited states corresponding to forward-charge transfer along the individual branches of the RC of H. modesticaldum. In the spectrum, these are located at an energy between the coupled Qy and Qx transitions. With ab initio Born-Oppenheimer molecular dynamics simulations, we reveal the influence of thermal vibrations on the excited states. The results show that the energy gap between the coupled Qy and the forward-charge transfer excitations is ∼0.4 eV, which we consider to conflict with the concept of a direct transfer mechanism. Our calculations, however, reveal a certain spectral overlap of the forward-charge transfer and the coupled Qx excitations. The reliability and robustness of the results are demonstrated by several numerical tests.

MCMap-A Computational Tool for Mapping Energy Landscapes of Transient Protein-Protein Interactions.

MCMap is a tool particularly well-suited for analyzing energy landscapes of transient macromolecular complexes. The program applies a Monte Carlo strategy, where the ligand moves randomly in the electrostatic field of the receptor. By applying importance sampling, the major interaction sites are mapped, resulting in a global distribution of ligand-receptor complexes. This approach displays the dynamic character of transiently interacting protein complexes where not a single complex but an ensemble of complexes better describes the protein interactions. The software provides a broad range of analysis options which allow for relating the simulations to experimental data and for interpreting them on a structural level. The application of MCMap is exemplified by the electron-transfer complex of cytochrome c peroxidase and cytochrome c from baker's yeast. The functionality of MCMap and the visualization of simulation data are in particular demonstrated by studying the dependence of the association on ionic strength and on the oxidation state of the binding partner. Furthermore, microscopically, a repulsion of a second ligand can be seen in the ternary complex upon the change of the oxidation state of the bound cytochrome c. The software is made available as open source software together with the example and can be downloaded free of charge from http://www.bisb.uni-bayreuth.de/index.php?page=downloads.