Mechanism of the initial charge separation in bacterial photosynthetic reaction centers

Excited-state structural dynamics of nickel complexes probed by optical and X-ray transient absorption spectroscopies: insights and implications

Excited states of nickel complexes undergo a variety of photochemical processes, such as charge transfer, ligation/deligation, and redox reactions, relevant to solar energy conversion and photocatalysis. The efficiencies of the aforementioned processes are closely coupled to the molecular structures in the ground and excited states. The conventional optical transient absorption spectroscopy has revealed important excited-state pathways and kinetics, but information regarding the metal center, in particular transient structural and electronic properties, remains limited. These deficiencies are addressed by X-ray transient absorption (XTA) spectroscopy, a detailed probe of 3d orbital occupancy, oxidation state and coordination geometry. The examples of excited-state structural dynamics of nickel porphyrin and nickel phthalocyanine have been described from our previous studies with highlights on the unique structural information obtained by XTA spectroscopy. We close by surveying prospective applications of XTA spectroscopy to active areas of Ni-based photocatalysis based on the knowledge gained from our previous studies.

Dynamic band-shift signal in two-dimensional electronic spectroscopy: A case of bacterial reaction center

Optical nonlinear spectroscopies carry a high amount of information about the systems under investigation; however, as they report polarization signals, the resulting spectra are often congested and difficult to interpret. To recover the landscape of energy states and physical processes such as energy and electron transfer, a clear interpretation of the nonlinear signals is prerequisite. Here, we focus on the interpretation of the electrochromic band-shift signal, which is generated when an internal electric field is established in the system following optical excitation. Whereas the derivative shape of the band-shift signal is well understood in transient absorption spectroscopy, its emergence in two-dimensional electronic spectroscopy (2DES) has not been discussed. In this work, we employed 2DES to follow the dynamic band-shift signal in reaction centers of purple bacteria Rhodobacter sphaeroides at 77 K. The prominent two-dimensional derivative-shape signal appears with the characteristic formation time of the charge separated state. To explain and characterize the band-shift signal, we use expanded double-sided Feynman diagram formalism. We propose to distinguish two types of Feynman diagrams that lead to signals with negative amplitude: excited state absorption and re-excitation. The presented signal decomposition and modeling analysis allows us to recover precise electrochromic shifts of accessory bacteriochlorophylls, identify additional signals in the B band range, and gain a further insight into the electron transfer mechanism. In a broader perspective, expanded Feynman diagram formalism will allow for interpretation of all 2D signals in a clearer and more intuitive way and therefore facilitate studying the underlying photophysics.

Ultrafast structural changes within a photosynthetic reaction centre

Photosynthetic reaction centres harvest the energy content of sunlight by transporting electrons across an energy-transducing biological membrane. Here we use time-resolved serial femtosecond crystallography¹ using an X-ray free-electron laser² to observe light-induced structural changes in the photosynthetic reaction centre of Blastochloris viridis on a timescale of picoseconds. Structural perturbations first occur at the special pair of chlorophyll molecules of the photosynthetic reaction centre that are photo-oxidized by light. Electron transfer to the menaquinone acceptor on the opposite side of the membrane induces a movement of this cofactor together with lower amplitude protein rearrangements. These observations reveal how proteins use conformational dynamics to stabilize the charge-separation steps of electron-transfer reactions.

Coupled excitation energy and charge transfer dynamics in reaction centre inspired model systems

Functional operation conditions of reaction centre core complexes require the tight coupling of exciton states to concomitant charge separation. Rigorous theoretical treatment of such integrated excitation energy transfer (EET) and charge transfer (CT) dynamics is particularly challenging due to (i) appreciable system sizes, (ii) inter-site and system-bath couplings of similar magnitude that render the Born-Markov approximation invalid, (iii) substantial reorganization energies of CT states, and (iv) the presence of complex structured spectral densities due to vibrational modes of the surroundings. We present numerical simulations on bacterial reaction centre (bRC) inspired model systems that utilize the recently developed MACGIC-iQUAPI method [Richter et al., J. Chem. Phys., 2017, 146, 214101]. The simulations demonstrate that the method provides a rigorous framework for the investigation of such integrated EET-CT dynamics. First, the applicability of the MACGIC-iQUAPI method is explored for a transition from monotonically decaying to oscillatory system-bath influence coefficients, a behavior inherently imposed by structured bath spectral densities. Tightly coupled EET and CT dynamics is further addressed for an excitonic subsystem that resembles strong coupling of special pair states and serves as donor towards a generic bridge-acceptor system. By solving the dissipative quantum dynamics of such bRC inspired model systems, the quenching of excitonic coherence on the hundreds of femtoseconds timescale is explored via a variation of the bridge state energetics, resembling a continuous transition from sequential to superexchange mediated CT regimes. Further, the simulations explore the influence of resonant vibrational modes on the quenching of excitonic coherence via CT. The results reveal a moderate influence of vibrational mode on charge separation dynamics in regimes of biologically relevant EET and CT dynamics.

Overlapping Electronic States with Nearly Parallel Transition Dipole Moments in Reduced Anionic Flavin Can Distort Photobiological Dynamics

Chromophoric biomolecules are exploited as reporters of a diverse set of phenomena, acting as internal distance monitors, environment and redox sensors, and endogenous imaging probes. The extent to which they can be exploited is dependent on an accurate knowledge of their fundamental electronic properties. Arguably of greatest importance is a precise knowledge of the direction(s) of the absorption transition dipole moment(s) (TDMs) in the molecular frame of reference. Such is the case for flavins, fluorescent redox cofactors utilized for ground- and excited-state redox and photochemical processes. The directions of the TDMs in oxidized and semiquinone flavins were characterized decades ago, and the details of charge redistribution in these forms have also been studied by Stark spectroscopy. The electronic structure of the fully reduced hydroquinone anionic state, FlH^-, however, has been the subject of unfounded assumptions and estimates about the number and direction of TDMs in FlH^-, as well the electronic structure changes that occur upon light absorption. Here we have used Stark spectroscopy to measure the magnitude and direction of charge redistribution in FlH^- upon optical excitation. These data were analyzed using TD-DFT calculations. The results show unequivocally that not one but two nearly orientation-degenerate electronic transitions are required to explain the 340-500 nm absorption spectral range, demolishing the commonly held assumption of a single transition. The difference dipole moments for these states show that electron density shifts toward the xylene ring for both transitions. These measurements force a reappraisal of previous studies that have used erroneous assumptions and unsubstantiated estimates of these quantities. The results put future optical studies of reduced flavins/flavoproteins on a firm photophysical footing.

Weak temperature dependence of P (+) H A (-) recombination in mutant Rhodobacter sphaeroides reaction centers

In contrast with findings on the wild-type Rhodobacter sphaeroides reaction center, biexponential P (+) H A (-)  → PH A charge recombination is shown to be weakly dependent on temperature between 78 and 298 K in three variants with single amino acids exchanged in the vicinity of primary electron acceptors. These mutated reaction centers have diverse overall kinetics of charge recombination, spanning an average lifetime from ~2 to ~20 ns. Despite these differences a protein relaxation model applied previously to wild-type reaction centers was successfully used to relate the observed kinetics to the temporal evolution of the free energy level of the state P (+) H A (-) relative to P (+) B A (-) . We conclude that the observed variety in the kinetics of charge recombination, together with their weak temperature dependence, is caused by a combination of factors that are each affected to a different extent by the point mutations in a particular mutant complex. These are as follows: (1) the initial free energy gap between the states P (+) B A (-) and P (+) H A (-) , (2) the intrinsic rate of P (+) B A (-)  → PB A charge recombination, and (3) the rate of protein relaxation in response to the appearance of the charge separated states. In the case of a mutant which displays rapid P (+) H A (-) recombination (ELL), most of this recombination occurs in an unrelaxed protein in which P (+) B A (-) and P (+) H A (-) are almost isoenergetic. In contrast, in a mutant in which P (+) H A (-) recombination is relatively slow (GML), most of the recombination occurs in a relaxed protein in which P (+) H A (-) is much lower in energy than P (+) H A (-) . The weak temperature dependence in the ELL reaction center and a YLH mutant was modeled in two ways: (1) by assuming that the initial P (+) B A (-) and P (+) H A (-) states in an unrelaxed protein are isoenergetic, whereas the final free energy gap between these states following the protein relaxation is large (~250 meV or more), independent of temperature and (2) by assuming that the initial and final free energy gaps between P (+) B A (-) and P (+) H A (-) are moderate and temperature dependent. In the case of the GML mutant, it was concluded that the free energy gap between P (+) B A (-) and P (+) H A (-) is large at all times.

Analysis of the temperature-dependence of P(+)HA(-) charge recombination in the Rhodobacter sphaeroides reaction center suggests nanosecond temperature-independent protein relaxation

The temperature dependence of charge recombination of the pair P(+)HA(-) in isolated reaction centers from the purple bacterium Rhodobacter sphaeroides with prereduced quinone QA was studied by sub-nanosecond to microsecond time-scale transient absorption. Overall, the kinetics slowed down substantially upon cooling from room temperature to ∼200 K, and then remained virtually unchanged down to 77 K, indicating the coexistence of two competitive pathways of charge recombination, a thermally-activated pathway appearing only above ~200 K and a temperature-independent pathway. In our modelling, the thermally activated pathway includes an uphill electron transfer from HA(-) to BA(-) leading to transient formation of the state P(+)BA(-), whereas the temperature-independent pathway is due to direct downhill electron transfer from HA(-) to P(+). At all temperatures studied, the kinetics could be approximated by a four-component decay. Detailed analysis of the lifetimes and amplitudes of particular phases over the range of temperatures suggests that the kinetically resolved phases reveal the consecutive appearance of three conformational states characterized by an increasing free energy gap between the states P(+)BA(-) and P(+)HA(-). The initial gap between these states was estimated to be only ~8 meV, the intermediate gap being ~92 meV, and the final gap ~135 meV, with no dependence on temperature. It was also calculated through a very straightforward approach that the relaxation process from the initial to the intermediate state occurs within 0.6 ± 0.1 ns, whereas the second step of relaxation from the intermediate to the final state takes 11 ± 2 ns. Both phases of the protein relaxation process are essentially temperature-independent. Possible alternative models to describe the experimental data that cannot be definitely excluded are also discussed.

Analysis of the kinetics of P+ HA- recombination in membrane-embedded wild-type and mutant Rhodobacter sphaeroides reaction centers between 298 and 77 K indicates that the adjacent negatively charged QA ubiquinone modulates the free energy of P+ HA- and may influence the rate of the protein dielectric response

Time-resolved spectroscopic studies of recombination of the P(+)HA(-) radical pair in photosynthetic reaction centers (RCs) from Rhodobacter sphaeroides give an opportunity to study protein dynamics triggered by light and occurring over the lifetime of P(+)HA(-). The state P(+)HA(-) is formed after the ultrafast light-induced electron transfer from the primary donor pair of bacteriochlorophylls (P) to the acceptor bacteriopheophytin (HA). In order to increase the lifetime of this state, and thus increase the temporal window for the examination of protein dynamics, it is possible to block forward electron transfer from HA(-) to the secondary electron acceptor QA. In this contribution, the dynamics of P(+)HA(-) recombination were compared at a range of temperatures from 77 K to room temperature, electron transfer from HA(-) to QA being blocked either by prereduction of QA or by genetic removal of QA. The observed P(+)HA(-) charge recombination was significantly slower in the QA-deficient RCs, and in both types of complexes, lowering the temperature from RT to 77 K led to a slowing of charge recombination. The effects are explained in the frame of a model in which charge recombination occurs via competing pathways, one of which is thermally activated and includes transient formation of a higher-energy state, P(+)BA(-). An internal electrostatic field supplied by the negative charge on QA increases the free energy levels of the state P(+)HA(-), thus decreasing its energetic distance to the state P(+)BA(-). In addition, the dielectric response of the protein environment to the appearance of the state P(+)HA(-) is accelerated from ∼50-100 ns in the QA-deficient mutant RCs to ∼1-16 ns in WT RCs with a negatively charged QA(-). In both cases, the temperature dependence of the protein dynamics is weak.

Dynamical theory of primary processes of charge separation in the photosynthetic reaction center

A dynamical theory has been developed for primary separation of charges in the course of photosynthesis. The theory deals with both hopping and superexchange transfer mechanisms. Dynamics of electron transfer from dimeric bacteriochlorophyll to quinone has been calculated. The results obtained agree with experimental data and provide a unified explanation of both the hierarchy of the transfer time in the photosynthetic reaction center and the phenomenon of coherent oscillations accompanying the transfer process.

Mechanism of recombination of the P+H(A)- radical pair in mutant Rhodobacter sphaeroides reaction centers with modified free energy gaps between P+B(A)- and P+H(A)-

The kinetics of recombination of the P(+)H(A)(-) radical pair were compared in wild-type reaction centers from Rhodobacter sphaeroides and in seven mutants in which the free energy gap, ΔG, between the charge separated states P(+)B(A)(-) and P(+)H(A)(-) was either increased or decreased. Five of the mutant RCs had been described previously, and X-ray crystal structures of two newly constructed complexes were determined by X-ray crystallography. The charge recombination reaction was accelerated in all mutants with a smaller ΔG than in the wild-type, and was slowed in a mutant having a larger ΔG. The free energy difference between the state P(+)H(A)(-) and the PH(A) ground state was unaffected by most of these mutations. These observations were consistent with a model in which the P(+)H(A)(-) → PH(A) charge recombination is thermally activated and occurs via the intermediate state P(+)B(A)(-), with a mean rate related to the size of the ΔG between the states P(+)B(A)(-) and P(+)H(A)(-) and not the ΔG between P(+)H(A)(-) and the ground state. A more detailed analysis of charge recombination in the mutants showed that the kinetics of the reaction were multiexponential, and characterized by ~0.5, ~1-3, and 7-17 ns lifetimes, similar to those measured for wild-type reaction centers. The exact lifetimes and relative amplitudes of the three components were strongly modulated by the mutations. Two models were considered in order to explain the observed multiexponentiality and modulation, involving heterogeneity or relaxation of P(+)H(A)(-) states, with the latter model giving a better fit to the experimental results.

Excitation and electron transfer in reaction centers from Rhodobacter sphaeroides probed and analyzed globally in the 1-nanosecond temporal window from 330 to 700 nm

Global analysis of a set of room temperature transient absorption spectra of Rhodobacter sphaeroides reaction centers, recorded in wide temporal and spectral ranges and triggered by femtosecond excitation of accessory bacteriochlorophylls at 800 nm, is presented. The data give a comprehensive review of all spectral dynamics features in the visible and near UV, from 330 to 700 nm, related to the primary events in the Rb. sphaeroides reaction center: excitation energy transfer from the accessory bacteriochlorophylls (B) to the primary donor (P), primary charge separation between the primary donor and primary acceptor (bacteriopheophytin, H), and electron transfer from the primary to the secondary electron acceptor (ubiquinone). In particular, engagement of the accessory bacteriochlorophyll in primary charge separation is shown as an intermediate electron acceptor, and the initial free energy gap of approximately 40 meV, between the states P(+)B(A)(-) and P(+)H(A)(-) is estimated. The size of this gap is shown to be constant for the whole 230 ps lifetime of the P(+)H(A)(-) state. The ultrafast spectral dynamics features recorded in the visible range are presented against a background of results from similar studies performed for the last two decades.

Measuring Electronic Coupling in the Reaction Center of Purple Photosynthetic Bacteria by Two-Color, Three-Pulse Photon Echo Peak Shift Spectroscopy

One- and two-color, three-pulse photon echo peak shift spectroscopy (1C and 2C3PEPS) was used to estimate the electronic coupling between the accessory bacteriochlorophyll (B) and the bacteriopheophytin (H) in the reaction center of the purple photosynthetic bacterium Rhodobacter sphaeroides as approximately 170 +/- 30 cm-1. This is the first direct experimental determination of this parameter; it is within the range of values found in previously published calculations. The 1C3PEPS signal of the Qy band of the bacteriochlorophyll B shows that it is weakly coupled to nuclear motions of the bath, whereas the 1C3PEPS signal of the Qy band of the bacteriopheophytin, H, shows that it is more strongly coupled to the bath, but has minimal inhomogeneous broadening. Our simulations capture the major features of the data with the theoretical framework developed in our group to separately calculate the response functions and population dynamics.

Dynamics of Photoinduced Electron Transfer in a Molecular Donor−Acceptor Quartet

The electronic structures and dynamics of photoinduced charge separation and recombination in a new donor/acceptor quartet molecule with bis-oligothiophene (BOTH) and bis-perylenediimide (BPDI) blocks attached to a benzene ring were described. Detailed transient spectroscopic studies were carried out on this compound and reference compounds at isolated molecular levels in solution. Two different dynamics of charge separation and recombination associated with two types of donor/acceptor pair conformations in solution were observed. These results were discussed based on Marcus theory and ascribed to both through-bond and through-space electron-transfer processes associated with two different orientations of the acceptors relative to the donor group. This molecular system exhibits a more efficient charge separation than charge recombination processes in both polar and nonpolar organic solvents, indicating that the material is an interesting candidate for photovoltaic studies in solid state.

A nonadiabatic theory for ultrafast catalytic electron transfer: a model for the photosynthetic reaction center

A non-adiabatic theory of Electron Transfer (ET), which improves the standard theory near the inversion point and becomes equivalent to it far from the inversion point, is presented. The complex amplitudes of the electronic wavefunctions at different sites are used as Kramers variables for describing the quantum tunneling of the electron in the deformable potential generated by its environment (nonadiabaticity) which is modeled as a harmonic classical thermal bath. After exact elimination of the bath, the effective electron dynamics is described by a discrete nonlinear Schrödinger equation with norm preserving dissipative terms and a Langevin random force, with a frequency cut-off, due to the thermalized phonons.This theory reveals the existence of a specially interesting marginal case when the linear and nonlinear coefficients of a two electronic states system are appropriately tuned for forming a Coherent Electron-Phonon Oscillator (CEPO). An electron injected on one of the electronic states of a CEPO generates large amplitude charge oscillations (even at zero temperature) associated with coherent phonon oscillations and electronic level oscillations. This fluctuating electronic level may resonate with a third site which captures the electron so that Ultrafast Electron Transfer (UFET) becomes possible. Numerical results are shown where two weakly interacting sites, a donor and a catalyst, form a CEPO that triggers an UFET to an acceptor. Without a catalytic site, a very large energy barrier prevents any direct ET. This UFET is shown to have many qualitative features similar to those observed in the primary charge separation in photosynthetic reaction centers. We suggest that more generally, CEPO could be a paradigm for understanding many selective chemical reactions involving electron transfer in biosystems.

Energy trapping and detrapping in reaction center mutants from Rhodobacter sphaeroides

Time-resolved fluorescence of chromatophores isolated from strains of Rhodobacter sphaeroides containing light harvesting complex I (LHI) and reaction center (RC) (no light harvesting complex II) was measured at several temperatures between 295 K and 10 K. Measurements were performed to investigate energy trapping from LHI to the RC in RC mutants that have a P/P(+) midpoint potential either above or below wild-type (WT). Six different strains were investigated: WT + LHI, four mutants with altered RC P/P(+) midpoint potentials, and an LHI-only strain. In the mutants with the highest P/P(+) midpoint potentials, the electron transfer rate decreases significantly, and at low temperatures it is possible to directly observe energy transfer from LHI to the RC by detecting the fluorescence kinetics from both complexes. In all mutants, fluorescence kinetics are multiexponential. To explain this, RC + LHI fluorescence kinetics were analyzed using target analysis in which specific kinetic models were compared. The kinetics at all temperatures can be well described with a model which accounts for the energy transfer between LHI and the RC and also includes the relaxation of the charge separated state P(+)H(A)(-), created in the RC as a result of the primary charge separation.

Selective perturbation of the second electron transfer step in mutant bacterial reaction centers

In order to specifically perturb the primary electron acceptor B(A) -- a monomeric bacteriochlorophyll (BChl) a -- involved in bacterial photosynthetic charge separation (CS), the protein environment of B(A) in the reaction center (RC) of Rhodobacter sphaeroides was modified by site-directed mutagenesis. Isolated RCs were characterized by redox titrations, low temperature optical spectroscopy, ENDOR/TRIPLE resonance spectroscopy and femtosecond time-resolved spectroscopy. Two mutations were studied: In the GS(M203) mutant a serine is introduced near the ring E keto group of B(A), while in FY(L146) a phenylalanine near the ring A acetyl group of B(A) is replaced by tyrosine. In all mutations the oxidation potential of the primary electron donor P as well as the electronic structure of both the P(*+) radical cation and the radical anion of the secondary electron acceptor, H(A)(*-), are not significantly altered compared to the wild type (WT), while changes of the optical absorption spectra at 77 K in the BChl Q(X) and Q(Y) regions are observed. The GS(M203) mutation only leads to a minor retardation of the CS reactions at room temperature, whereas for FY(L146) significant deviations from the native electron transfer (ET) rates could be detected: In addition to a faster first (2.9 ps) and a slower second (1 ps) ET step, a new 8-ps time constant was found in the FY(L146) mutant, which can be ascribed to a fraction of RCs with slowed down secondary ET. The results allow us to address the functional role of the acetyl group of B(A) and question the role of the free energy changes as the main determining factor of ET rates in RCs. It is concluded that structural rearrangements alter the electronic coupling between the pigments and thereby influence the rate of fast CS.

Nuclear wavepacket motion producing a reversible charge separation in bacterial reaction centers

The excitation of bacterial reaction centers (RCs) at 870 nm by 30 fs pulses induces the nuclear wavepacket motions on the potential energy surface of the primary electron donor excited state P*, which lead to the fs oscillations in stimulated emission from P* [M.H. Vos, M.R. Jones, C.N. Hunter, J. Breton, J.-C. Lambry and J.-L. Martin (1994) Biochemistry 33, 6750-6757] and in Qy absorption band of the primary electron acceptor, bacteriochlorophyll monomer B(A) [A.M. Streltsov, S.I.E. Vulto, A.Y. Shkuropatov, A.J. Hoff, T.J. Aartsma and V.A. Shuvalov (1998) J. Phys. Chem. B 102, 7293-7298] with a set of fundamental frequencies in the range of 10-300 cm(-1). We have found that in pheophytin-modified RCs, the fs oscillations with frequency around 130 cm(-1) observed in the P*-stimulated emission as well as in the B(A) absorption band at 800 nm are accompanied by remarkable and reversible formation of the 1020 nm absorption band which is characteristic of the radical anion band of bacteriochlorophyll monomer B(A)-. These results are discussed in terms of a reversible electron transfer between P* and B(A) induced by a motion of the wavepacket near the intersection of potential energy surfaces of P* and P+B(A)-, when a maximal value of the Franck-Condon factor is created.

Ab initio calculation of electronic coupling in the photosynthetic reaction center

We have carried out an ab initio electronic structure calculations of electron transfer couplings between chromophores in the bacterial photosynthetic reaction center. The couplings agree remarkably well with parameters obtained from recent quantum dynamical modeling of experimental data assuming an explicit intermediate mechanism. We also have computed couplings on the M-side of the reaction center and have found that the interaction of the primary donor to the M-side intermediate bacteriochlorophyll is quite small because of destructive interference of the two localized coupling matrix elements. This may explain the slow rate of electron transfer down the M-side of the reaction center.