Theoretical analysis of the electronic asymmetry of the special pair in the photosynthetic reaction center: Effect of structural asymmetry and protein environment

Charge-Transfer Processes within the Isolated Tetramer Models of the Reaction Center Rhodobacter sphaeroides

Over the last two decades, advancements in structural resolution and spectral characterization have significantly enhanced our understanding of photosynthesis. However, the complexity of photosystem (PS) supercomplexes still presents challenges. In the Rhodobacter sphaeroides reaction center (RSRC), the charge separation process begins with a charge-transfer (CT) step at the special pair (P), a dimer of bacteriochlorophyll a (BChl), which acts as the donor, and continues with electron transport through the active pigments. Our computational study explores CT rectification in RSRC. We find that the CT rate is faster in the A branch compared with the B branch, which can be attributed to the orientation of the pigments near P and the influence of the surrounding protein complex on the dielectric constant. The calculated rate constants are derived using Fermi's golden rule, with a first-principles approach that employs an optimally tuned screened range-separated hybrid functional within a polarizable continuum model (SRSH-PCM).

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso-meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4-phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi's golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

On the Role of the Special Pair in Photosystems as a Charge Transfer Rectifier

The special pair, a bacteriochlorophyll a (BChl) dimer found at the core of bacterial reaction centers, is known to play a key role in the functionality of photosystems as a precursor to the photosynthesis process. In this paper, we analyze the inherent affinity of the special pair to rectify the intrapair photo-induced charge transfer (CT). In particular, we show that the molecular environment affects the nuclear geometry, resulting in symmetry breaking between the two possible intrapair CT processes. To this end, we study the relationships of the intrapair CT and the molecular geometry with respect to the effective dielectric constant provided by the molecular environment. We identify the special pair structural feature that breaks the symmetry between the two molecules, leading to CT rectification. Excited state energies, oscillator strengths, and electronic coupling values are obtained via time-dependent density functional theory, employing a recently developed framework based on a screened range-separated hybrid functional within a polarizable continuum model (SRSH-PCM). We analyze the rectification capability of the special pair by calculating the CT rates using a first-principles-based Fermi's golden rule approach.

Dimerization Behavior of Methyl Chlorophyllide a as the Model of Chlorophyll a in the Presence of Water Molecules-Theoretical Study

A dimerization of methyl chlorophyllide a molecules and a role of water in stabilization and properties of methyl chlorophyllide a dimers were studied by means of symmetry-adapted perturbation theory (SAPT), functional-group SAPT (F-SAPT), density-functional theory (DFT), and time-dependent DFT approaches. The quantification of various types of interactions, such as π-π stacking, coordinative, and hydrogen bonding by applying the F-SAPT energy decomposition scheme shows the major role of the magnesium atom and the pheophytin macrocycle in the stability of the complex. The examination of interaction energy components with respect to a mutual orientation of monomers and in the presence or absence of water molecules reveals that the dispersion energy is the main binding factor of the interaction, while water molecules tend to weaken the attraction between methyl chlorophyllide a species. The dimerization can be seen in computed UV-vis spectra, and results in a doubling of the lowest peaks, as compared to the monomer spectrum, and in an intensity rise of the lowest 1.8 and 2.4 eV peaks at a cost of the 3.5 eV peaks for the majority of dimer configurations. The complexation of water has little effect on the peaks' position; however, it affects the overall shape of simulated spectra through changes in peak intensities, which is strongly dependent on the structure of the complex. The VCD spectra for the dimers show several characteristic features attributed to the interaction of substituting groups and/or water ligand attached to macrocycle groups belonging to different monomers. VCD is sensitive to the type of the formed dimer, but not to the number of water molecules it contains. This and several other features, as well as the differential UV-vis spectra, may serve as the indicator of the presence of a given dimer structure in the experiment.

Electronic structures of heme a of cytochrome c oxidase in the redox states--charge density migration to the propionate groups of heme a

The electronic structures of heme a of cytochrome c oxidase in the redox states were studied, using hybrid density functional theory with a polarizable continuum model and a point charge model. We found that the most stable electronic configurations of the d electrons of the Fe ion are determined by the orbital interactions of the d orbitals of the Fe ion with the pi orbitals of the porphyrin ring and the His residues. The redox reaction of the Fe ion influences the charge density on the formyl group through the pi conjugation of the porphyrin ring. In addition, we found the charge transfer from the Fe ion to the propionate group of heme a in the redox change despite the lack of the pi-conjugation. We elucidated that the charge propagation originates from the heme a structure itself and that the origin of the charge delocalization to the heme propionate is the orbital interactions between the d orbital of the Fe ion and the p orbitals of the carboxylate part of the heme propionate via the pi conjugation of the porphyrin ring and the sigma* orbital of the C-C bond of the propionate group. The electrostatic effect by surrounding proteins enhances the charge transfer from the Fe ion to the propionate group. These results indicate that heme propionate groups serve electron mediators in electron transfer as well as electrostatic anchors, and that proteins surrounding the active site reinforce the congenital abilities of the cofactors.

Electronic pathway in reaction centers from Rhodobacter sphaeroides and Chloroflexus aurantiacus

The reaction centers (RC) of Chloroflexus aurantiacus and Rhodobacter sphaeroidesH(M182)L mutant were investigated. Prediction for electron transfer (ET) at very low temperatures was also performed. To describe the kinetics of the C. aurantiacus RCs, the incoherent model of electron transfer was used. It was shown that the asymmetry in electronic coupling parameters must be included to explain the experiments. For the description of R. sphaeroidesH(M182)L mutant RCs, the coherent and incoherent models of electron transfer were used. These two models are discussed with regard to the observed electron transfer kinetics. It seems likely that the electron transfer asymmetry in R. sphaeroides RCs is caused mainly by the asymmetry in the free energy levels of L- and M-side cofactors. In the case of C. aurantiacus RCs, the unidirectionality of the charge separation can be caused mainly by the difference in the electronic coupling parameters in two branches.

Effect of charge distribution over a chlorophyll dimer on the redox potential of P680 in photosystem II as studied by density functional theory calculations

The effect of charge distribution over a chlorophyll dimer on the redox potential of P680 in photosystem II was studied by density functional theory calculations using the P680 coordinates in the X-ray structure. From the calculated ionization potentials of the dimer and the monomeric constituents, the decrease in the redox potential by charge delocalization over the dimer was estimated to be approximately 140 mV. Such charge delocalization was previously observed in the isolated D1-D2-Cyt b 559 complexes, whereas the charge was primarily localized on P D1 in the core complexes. The calculated potential decrease of approximately 140 mV can explain the inhibition of Y Z oxidation in the former complexes and in turn implies that the charge localization on P D1 upon formation of the core complex increases the P680 potential to the level necessary for water oxidation.

Characterization of highly purified photosystem I complexes from the chlorophyll d-dominated cyanobacterium Acaryochloris marina MBIC 11017

Photochemically active photosystem (PS) I complexes were purified from the chlorophyll (Chl) d-dominated cyanobacterium Acaryochloris marina MBIC 11017, and several of their properties were characterized. PS I complexes consist of 11 subunits, including PsaK1 and PsaK2; a new small subunit was identified and named Psa27. The new subunit might replace the function of PsaI that is absent in A. marina. The amounts of pigments per one molecule of Chl d' were 97.0 +/- 11.0 Chl d, 1.9 +/- 0.5 Chl a, 25.2 +/- 2.4 alpha-carotene, and two phylloquinone molecules. The light-induced Fourier transform infrared difference spectroscopy and light-induced difference absorption spectra reconfirmed that the primary electron donor of PS I (P740) was the Chl d dimer. In addition to P740, the difference spectrum contained an additional band at 728 nm. The redox potentials of P740 were estimated to be 439 mV by spectroelectrochemistry; this value was comparable with the potential of P700 in other cyanobacteria and higher plants. This suggests that the overall energetics of the PS I reaction were adjusted to the electron acceptor side to utilize the lower light energy gained by P740. The distribution of charge in P740 was estimated by a density functional theory calculation, and a partial localization of charge was predicted to P1 Chl (special pair Chl on PsaA). Based on differences in the protein matrix and optical properties of P740, construction of the PS I core in A. marina was discussed.