Structure-based calculations of optical spectra of photosystem I suggest an asymmetric light-harvesting process

Molecular Mechanisms behind Circular Dichroism Spectral Variations between Channelrhodopsin and Heliorhodopsin Dimers

Channelrhodopsin (ChR) and heliorhodopsin (HeR) are microbial rhodopsins with similar structures but different circular dichroism (CD) spectra: ChR shows biphasic negative and positive bands, whereas HeR shows a single positive band. We explored the physicochemical factors underlying these differences through computational methods. Using the exciton model based on first-principles computations, we obtained the CD spectra of ChR and HeR. The obtained spectra indicate that the protein dimer structures and the quantum mechanical treatment of the retinal chromophore and its interacting amino acids are crucial for accurately reproducing the experimental spectra. Further calculations revealed that the sign of the excitonic coupling was opposite between the ChR and HeR dimers, which was attributed to the contrasting second term of the orientation factor between the two retinal chromophores. These findings demonstrate that slight variations in the intermolecular orientation of the two chromophores can result in significant differences in the CD spectral shape.

Protein Effects on the Excitation Energies and Exciton Dynamics of the CP24 Antenna Complex

In this study, the site energy fluctuations, energy transfer dynamics, and some spectroscopic properties of the minor light-harvesting complex CP24 in a membrane environment were determined. For this purpose, a 3 μs-long classical molecular dynamics simulation was performed for the CP24 complex. Furthermore, using the density functional tight binding/molecular mechanics molecular dynamics (DFTB/MM MD) approach, we performed excited state calculations for the chlorophyll a and chlorophyll b molecules in the complex starting from five different positions of the MD trajectory. During the extended simulations, we observed variations in the site energies of the different sets as a result of the fluctuating protein environment. In particular, a water coordination to Chl-b 608 occurred only after about 1 μs in the simulations, demonstrating dynamic changes in the environment of this pigment. From the classical and the DFTB/MM MD simulations, spectral densities and the (time-dependent) Hamiltonian of the complex were determined. Based on these results, three independent strongly coupled chlorophyll clusters were revealed within the complex. In addition, absorption and fluorescence spectra were determined together with the exciton relaxation dynamics, which reasonably well agrees with experimental time scales.

High-Resolution Frequency-Domain Spectroscopic and Modeling Studies of Photosystem I (PSI), PSI Mutants and PSI Supercomplexes

Photosystem I (PSI) is one of the two main pigment-protein complexes where the primary steps of oxygenic photosynthesis take place. This review describes low-temperature frequency-domain experiments (absorption, emission, circular dichroism, resonant and non-resonant hole-burned spectra) and modeling efforts reported for PSI in recent years. In particular, we focus on the spectral hole-burning studies, which are not as common in photosynthesis research as the time-domain spectroscopies. Experimental and modeling data obtained for trimeric cyanobacterial Photosystem I (PSI3), PSI3 mutants, and PSI3-IsiA18 supercomplexes are analyzed to provide a more comprehensive understanding of their excitonic structure and excitation energy transfer (EET) processes. Detailed information on the excitonic structure of photosynthetic complexes is essential to determine the structure-function relationship. We will focus on the so-called "red antenna states" of cyanobacterial PSI, as these states play an important role in photochemical processes and EET pathways. The high-resolution data and modeling studies presented here provide additional information on the energetics of the lowest energy states and their chlorophyll (Chl) compositions, as well as the EET pathways and how they are altered by mutations. We present evidence that the low-energy traps observed in PSI are excitonically coupled states with significant charge-transfer (CT) character. The analysis presented for various optical spectra of PSI3 and PSI3-IsiA18 supercomplexes allowed us to make inferences about EET from the IsiA18 ring to the PSI3 core and demonstrate that the number of entry points varies between sample preparations studied by different groups. In our most recent samples, there most likely are three entry points for EET from the IsiA18 ring per the PSI core monomer, with two of these entry points likely being located next to each other. Therefore, there are nine entry points from the IsiA18 ring to the PSI3 trimer. We anticipate that the data discussed below will stimulate further research in this area, providing even more insight into the structure-based models of these important cyanobacterial photosystems.

Access to the Antenna System of Photosystem I via Single-Molecule Excitation-Emission Spectroscopy

In the development of single-molecule spectroscopy, the simultaneous detection of the excitation and emission spectra has been limited. The fluorescence excitation spectrum based on background-free signals is compatible with the fluorescence-emission-based detection of single molecules and can provide insight into the variations in the input energy of the different terminal emitters. Here, we implement single-molecule excitation-emission spectroscopy (SMEES) for photosystem I (PSI) via a cryogenic optical microscope. To this end, we extended our line-focus-based excitation-spectral microscope system to the cryogenic temperature-compatible version. PSI is one of the two photosystems embedded in the thylakoid membrane in oxygen-free photosynthetic organisms. PSI plays an essential role in electron transfer in the photosynthesis reaction. PSIs of many organisms contain a few red-shifted chlorophylls (Chls) with much lower excitation energies than ordinary antenna Chls. The fluorescence emission spectrum originates primarily from the red-shifted Chls, whereas the excitation spectrum is sensitive to the antenna Chls that are upstream of red-shifted Chls. Using SMEES, we obtained the inclining two-dimensional excitation-emission matrix (2D-EEM) of PSI particles isolated from a cyanobacterium, Thermosynechococcus vestitus (equivalent to elongatus), at about 80 K. Interestingly, by decomposing the inclining 2D-EEMs within time course observation, we found prominent variations in the excitation spectra of the red-shifted Chl pools with different emission wavelengths, strongly indicating the variable excitation energy transfer (EET) pathway from the antenna to the terminal emitting pools. SMEES helps us to directly gain information about the antenna system, which is fundamental to depicting the EET within pigment-protein complexes.

Exciton interactions of chlorophyll tetramer in water-soluble chlorophyll-binding protein BoWSCP

The exciton interaction of four chlorophyll a (Chl a) molecules in a symmetrical tetrameric complex of the water-soluble chlorophyll-binding protein BoWSCP was analyzed in the pH range of 3-11. Exciton splitting ΔE = 232 ± 2 cm-1 of the Qy band of Chl a into two subcomponents with relative intensities of 78.1 ± 0.7 % and 21.9 ± 0.7 % was determined by a joint decomposition of the absorption and circular dichroism spectra into Gaussian functions. The exciton coupling parameters were calculated based on the BoWSCP atomic structure in three approximations: the point dipole model, the distributed atomic monopoles, and direct ab initio calculations in the TDDFT/PCM approximation. The Coulomb interactions of monomers were calculated within the continuum model using three values of optical permittivity. The models based on the properties of free Chl a in solution suffer from significant errors both in estimating the absolute value of the exciton interaction and in the relative intensity of exciton transitions. Calculations within the TDDFT/PCM approximation reproduce the experimentally determined parameters of the exciton splitting and the relative intensities of the exciton bands. The following factors of pigment-protein and pigment-pigment interactions were examined: deviation of the macrocycle geometry from the planar conformation of free Chl; the formation of hydrogen bonds between the macrocycle and water molecules; the overlap of wave functions of monomers at close distances. The most significant factor is the geometrical deformation of the porphyrin macrocycle, which leads to an increase in the dipole moment of Chl monomer from 5.5 to 6.9 D and to a rotation of the dipole moment by 15° towards the cyclopentane ring. The contributions of resonant charge-transfer states to the wave functions of the Chl dimer were determined and the transition dipole moments of the symmetric and antisymmetric charge-transfer states were estimated.

Experimental and Theoretical Mutation of Exciton States on the Smallest Type-I Photosynthetic Reaction Center Complex of a Green Sulfur Bacterium Chlorobaclum tepidum

The exciton states on the smallest type-I photosynthetic reaction center complex of a green sulfur bacterium Chlorobaculum tepidum (GsbRC) consisting of 26 bacteriochlorophylls a (BChl a) and four chlorophylls a (Chl a) located on the homodimer of two PscA reaction center polypeptides were investigated. This analysis involved the study of exciton states through a combination of theoretical modeling and the genetic removal of BChl a pigments at eight sites. (1) A theoretical model of the pigment assembly exciton state on GsbRC was constructed using Poisson TrESP (P-TrESP) and charge density coupling (CDC) methods based on structural information. The model reproduced the experimentally obtained absorption spectrum, circular dichroism spectrum, and excitation transfer dynamics, as well as explained the effects of mutation. (2) Eight BChl a molecules at different locations on the GsbRC were selectively removed by genetic exchange of the His residue, which ligates the central Mg atom of BChl a, with the Leu residue on either one or two PscAs in the RC. His locations are conserved among all type-I RC plant polypeptide, cyanobacteria, and bacteria amino acid sequences. (3) Purified mutant-GsbRCs demonstrated distinct absorption and fluorescence spectra at 77 K, which were different from each other, suggesting successful pigment removal. (4) The same mutations were applied to the constructed theoretical model to analyze the outcomes of these mutations. (5) The combination of theoretical predictions and experimental mutations based on structural information is a new tool for studying the function and evolution of photosynthetic reaction centers.

The location of the low-energy states in Lhca1 favors excitation energy transfer to the core in the plant PSI-LHCI supercomplex

Lhca1 is one of the four pigment-protein complexes composing the outer antenna of plant Photosystem I-light-havesting I supercomplex (PSI-LHCI). It forms a functional dimer with Lhca4 but, differently from this complex, it does not contain 'red-forms,' i.e., pigments absorbing above 700 nm. Interestingly, the recent PSI-LHCI structures suggest that Lhca1 is the main point of delivering the energy harvested by the antenna to the core. To identify the excitation energy pathways in Lhca1, we developed a structure-based exciton model based on the simultaneous fit of the low-temperature absorption, linear dichroism, and fluorescence spectra of wild-type Lhca1 and two mutants, lacking chlorophylls contributing to the long-wavelength region of the absorption. The model enables us to define the locations of the lowest energy pigments in Lhca1 and estimate pathways and timescales of energy transfer within the complex and to the PSI core. We found that Lhca1 has a particular energy landscape with an unusual (compared to Lhca4, LHCII, and CP29) configuration of the low-energy states. Remarkably, these states are located near the core, facilitating direct energy transfer to it. Moreover, the low-energy states of Lhca1 are also coupled to the red-most state (red forms) of the neighboring Lhca4 antenna, providing a pathway for effective excitation energy transfer from Lhca4 to the core.

Thermal site energy fluctuations in photosystem I: new insights from MD/QM/MM calculations

Cyanobacterial photosystem I (PSI) is one of the most efficient photosynthetic machineries found in nature. Due to the large scale and complexity of the system, the energy transfer mechanism from the antenna complex to the reaction center is still not fully understood. A central element is the accurate evaluation of the individual chlorophyll excitation energies (site energies). Such an evaluation must include a detailed treatment of site specific environmental influences on structural and electrostatic properties, but also their evolution in the temporal domain, because of the dynamic nature of the energy transfer process. In this work, we calculate the site energies of all 96 chlorophylls in a membrane-embedded model of PSI. The employed hybrid QM/MM approach using the multireference DFT/MRCI method in the QM region allows to obtain accurate site energies under explicit consideration of the natural environment. We identify energy traps and barriers in the antenna complex and discuss their implications for energy transfer to the reaction center. Going beyond previous studies, our model also accounts for the molecular dynamics of the full trimeric PSI complex. Via statistical analysis we show that the thermal fluctuations of single chlorophylls prevent the formation of a single prominent energy funnel within the antenna complex. These findings are also supported by a dipole exciton model. We conclude that energy transfer pathways may form only transiently at physiological temperatures, as thermal fluctuations overcome energy barriers. The set of site energies provided in this work sets the stage for theoretical and experimental studies on the highly efficient energy transfer mechanisms in PSI.

Quantum mechanical analysis of excitation energy transfer couplings in photosystem II

We evaluated excitation energy transfer (EET) coupling (J) between all pairs of chlorophylls (Chls) and pheophytins (Pheos) in the protein environment of photosystem II based on the time-dependent density functional theory with a quantum mechanical/molecular mechanics approach. In the reaction center, the EET coupling between Chls PD1 and PD2 is weaker (|J(PD1/PD2)| = 79 cm-1), irrespective of a short edge-to-edge distance of 3.6 Å (Mg-to-Mg distance of 8.1 Å), than the couplings between PD1 and the accessory ChlD1 (|J(PD1/ChlD2)| = 104 cm-1) and between PD2 and ChlD2 (|J(PD2/ChlD1)| = 101 cm-1), suggesting that PD1 and PD2 are two monomeric Chls rather than a "special pair". There exist strongly coupled Chl pairs (|J| > ∼100 cm-1) in the CP47 and CP43 core antennas, which may be candidates for the red-shifted Chls observed in spectroscopic studies. In CP47 and CP43, Chls ligated to CP47-His26 and CP43-His56, which are located in the middle layer of the thylakoid membrane, play a role in the "hub" that mediates the EET from the lumenal to stromal layers. In the stromal layer, Chls ligated to CP47-His466, CP43-His441, and CP43-His444 mediate the EET from CP47 to ChlD2/PheoD2 and from CP43 to ChlD1/PheoD1 in the reaction center. Thus, the excitation energy from both CP47 and CP43 can always be utilized for the charge-separation reaction in the reaction center.

Frequency-Domain Spectroscopic Study of the Photosystem I Supercomplexes, Isolated IsiA Monomers, and the Intact IsiA Ring

The PSI₃-IsiA₁₈ supercomplex is one of the largest and most complicated assemblies in photosynthesis. The IsiA ring, composed of 18 IsiA monomers (IsiA₁₈) surrounding the PSI trimer (PSI₃), forms under iron-deficient conditions in cyanobacteria and acts as a peripheral antenna. Based on the supercomplex structure recently determined via cryo-EM imaging, we model various optical spectra of the IsiA monomers and IsiA₁₈ ring. Comparison of the absorption and emission spectra of the isolated IsiA monomers and the full ring reveals that about 2.7 chlorophylls (Chls) are lost in the isolated IsiA monomers. The best fits for isolated monomers spectra are obtained assuming the absence of Chl 508 and Chl 517 and 70% loss of Chl 511. The best model describing all three hexamers and the entire ring suggests that the lowest energy pigments are Chls 511, 514, and 517. Based on the modeling results presented in this work, we conclude that there are most likely three entry points for EET from the IsiA₆ hexamer to the PSI core monomer, with two of these entry points likely being located next to each other (i.e., nine entry points from IsiA₁₈ to the PSI₃ trimer). Finally, we show that excitation energy transfer inside individual monomers is fast (<2 ps at T = 5 K) and at least 20 times faster than intermonomer energy transfer.

Current state of the primary charge separation mechanism in photosystem I of cyanobacteria

This review analyzes new data on the mechanism of ultrafast reactions of primary charge separation in photosystem I (PS I) of cyanobacteria obtained in the last decade by methods of femtosecond absorption spectroscopy. Cyanobacterial PS I from many species harbours 96 chlorophyll a (Chl a) molecules, including six specialized Chls denoted Chl1A/Chl1B (dimer P700, or PAPB), Chl2A/Chl2B, and Chl3A/Chl3B arranged in two branches, which participate in electron transfer reactions. The current data indicate that the primary charge separation occurs in a symmetric exciplex, where the special pair P700 is electronically coupled to the symmetrically located monomers Chl2A and Chl2B, which can be considered together as a symmetric exciplex Chl2APAPBChl2B with the mixed excited (Chl2APAPBChl2B)* and two charge-transfer states P700 +Chl2A - and P700 +Chl2B -. The redistribution of electrons between the branches in favor of the A-branch occurs after reduction of the Chl2A and Chl2B monomers. The formation of charge-transfer states and the symmetry breaking mechanisms were clarified by measuring the electrochromic Stark shift of β-carotene and the absorption dynamics of PS I complexes with the genetically altered Chl 2B or Chl 2A monomers. The review gives a brief description of the main methods for analyzing data obtained using femtosecond absorption spectroscopy. The energy levels of excited and charge-transfer intermediates arising in the cyanobacterial PS I are critically analyzed.

Spectral densities and absorption spectra of the core antenna complex CP43 from photosystem II

Besides absorbing light, the core antenna complex CP43 of photosystem II is of great importance in transferring excitation energy from the antenna complexes to the reaction center. Excitation energies, spectral densities, and linear absorption spectra of the complex have been evaluated by a multiscale approach. In this scheme, quantum mechanics/molecular mechanics molecular dynamics simulations are performed employing the parameterized density functional tight binding (DFTB) while the time-dependent long-range-corrected DFTB scheme is applied for the excited state calculations. The obtained average spectral density of the CP43 complex shows a very good agreement with experimental results. Moreover, the excitonic Hamiltonian of the system along with the computed site-dependent spectral densities was used to determine the linear absorption. While a Redfield-like approximation has severe shortcomings in dealing with the CP43 complex due to quasi-degenerate states, the non-Markovian full second-order cumulant expansion formalism is able to overcome the drawbacks. Linear absorption spectra were obtained, which show a good agreement with the experimental counterparts at different temperatures. This study once more emphasizes that by combining diverse techniques from the areas of molecular dynamics simulations, quantum chemistry, and open quantum systems, it is possible to obtain first-principle results for photosynthetic complexes, which are in accord with experimental findings.

Theoretical Model of the Far-Red-Light-Adapted Photosystem I Reaction Center of Cyanobacterium Acaryochloris marina Using Chlorophyll d and the Effect of Chlorophyll Exchange

A theoretical model of the far-red-light-adapted photosystem I (PSI) reaction center (RC) complex of a cyanobacterium, Acaryochloris marina (AmPSI), was constructed based on the exciton theory and the recently identified molecular structure of AmPSI by Hamaguchi et al. (Nat. Commun., 2021, 12, 2333). A. marina performs photosynthesis under the visible to far-red light (400-750 nm), which is absorbed by chlorophyll d (Chl-d). It is in contrast to the situation of all the other oxygenic photosynthetic processes of cyanobacteria and plants, which contains chlorophyll a (Chl-a) that absorbs only 400-700 nm visible light. AmPSI contains 70 Chl-d, 1 Chl-d^', 2 pheophytin a (Pheo-a), and 12 carotenoids in the currently available structure. A special pair of Chl-d/Chl-d^' acts as the electron donor (P740) and two Pheo-a act as the primary electron acceptor A₀ as the counterparts of P700 and Chl-a, respectively, of Chl-a-type PSIs. The exciton Hamiltonian of AmPSI was constructed considering the excitonic coupling strength and site energy shift of individual pigments using the Poisson-TrESP (P-TrESP) and charge density coupling (CDC) methods. The model was constructed to fit the experimentally measured spectra of absorption and circular dichroism (CD) spectra during downhill/uphill excitation energy transfer processes. The constructed theoretical model of AmPSI was further compared with the Chl-a-type PSI of Thermosynechococcus elongatus (TePSI), which contains only Chl-a and Chl-a^'. The functional properties of AmPSI and TePSI were further examined by the in silico exchange of Chl-d by Chl-a in the models.

Structural basis for the absence of low-energy chlorophylls in a photosystem I trimer from Gloeobacter violaceus

Photosystem I (PSI) is a multi-subunit pigment-protein complex that functions in light-harvesting and photochemical charge-separation reactions, followed by reduction of NADP to NADPH required for CO₂ fixation in photosynthetic organisms. PSI from different photosynthetic organisms has a variety of chlorophylls (Chls), some of which are at lower-energy levels than its reaction center P700, a special pair of Chls, and are called low-energy Chls. However, the sites of low-energy Chls are still under debate. Here, we solved a 2.04-Å resolution structure of a PSI trimer by cryo-electron microscopy from a primordial cyanobacterium Gloeobacter violaceus PCC 7421, which has no low-energy Chls. The structure shows the absence of some subunits commonly found in other cyanobacteria, confirming the primordial nature of this cyanobacterium. Comparison with the known structures of PSI from other cyanobacteria and eukaryotic organisms reveals that one dimeric and one trimeric Chls are lacking in the Gloeobacter PSI. The dimeric and trimeric Chls are named Low1 and Low2, respectively. Low2 is missing in some cyanobacterial and eukaryotic PSIs, whereas Low1 is absent only in Gloeobacter. These findings provide insights into not only the identity of low-energy Chls in PSI, but also the evolutionary changes of low-energy Chls in oxyphototrophs.

Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches

A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0–0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang–Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang–Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0–0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling.

The importance of reaction field and screening effects on the excitonic couplings is demonstrated, and from quantum-chemical calculations a single scaling factor is derived that can be used to improve simple models based on the Poisson equation.

Recent advances in single-molecule spectroscopy studies on light-harvesting processes in oxygenic photosynthesis

Photosynthetic light-harvesting complexes (LHCs) play a crucial role in concentrating the photon energy from the sun that otherwise excites a typical pigment molecule, such as chlorophyll-a, only several times a second. Densely packed pigments in the complexes ensure efficient energy transfer to the reaction center. At the same time, LHCs have the ability to switch to an energy-quenching state and thus play a photoprotective role under excessive light conditions. Photoprotection is especially important for oxygenic photosynthetic organisms because toxic reactive oxygen species can be generated through photochemistry under aerobic conditions. Because of the extreme complexity of the systems in which various types of pigment molecules strongly interact with each other and with the surrounding protein matrixes, there has been long-standing difficulty in understanding the molecular mechanisms underlying the flexible switching between the light-harvesting and quenching states. Single-molecule spectroscopy studies are suitable to reveal the conformational dynamics of LHCs reflected in the fluorescence properties that are obscured in ordinary ensemble measurements. Recent advanced single-molecule spectroscopy studies have revealed the dynamical fluctuations of LHCs in their fluorescence peak position, intensity, and lifetime. The observed dynamics seem relevant to the conformational plasticity required for the flexible activations of photoprotective energy quenching. In this review, we survey recent advances in the single-molecule spectroscopy study of the light-harvesting systems of oxygenic photosynthesis.

Electronic Couplings and Electrostatic Interactions Behind the Light Absorption of Retinal Proteins

The photo-functional chromophore retinal exhibits a wide variety of optical absorption properties depending on its intermolecular interactions with surrounding proteins and other chromophores. By utilizing these properties, microbial and animal rhodopsins express biological functions such as ion-transport and signal transduction. In this review, we present the molecular mechanisms underlying light absorption in rhodopsins, as revealed by quantum chemical calculations. Here, symmetry-adapted cluster-configuration interaction (SAC-CI), combined quantum mechanical and molecular mechanical (QM/MM), and transition-density-fragment interaction (TDFI) methods are used to describe the electronic structure of the retinal, the surrounding protein environment, and the electronic coupling between chromophores, respectively. These computational approaches provide successful reproductions of experimentally observed absorption and circular dichroism (CD) spectra, as well as insights into the mechanisms of unique optical properties in terms of chromophore-protein electrostatic interactions and chromophore-chromophore electronic couplings. On the basis of the molecular mechanisms revealed in these studies, we also discuss strategies for artificial design of the optical absorption properties of rhodopsins.

Spectral Tuning Mechanism of Photosynthetic Light-Harvesting Complex II Revealed by Ab Initio Dimer Exciton Model

Excited states of two kinds of bacteriochlorophyll (BChl) aggregates, B850 and B800, in photosynthetic light-harvesting complex II (LH2) are theoretically investigated by developing and using an extended exciton model considering efficiently evaluated excitonic coupling. Our exciton model based on dimer fragmentation is shown to reproduce the experimental absorption spectrum of LH2 with good accuracy, entailing their different redshifts originating from aggregations of B850 and B800. The systematic analysis has been performed on the spectra by quantitatively decomposing their spectral shift energies into the contributions of various effects: structural distortion, electrostatic, excitonic coupling, and charge-transfer (CT) effects. Our results show that the spectral redshift of B800 is mainly attributed to its electrostatic interaction with the protein environment, while that of B850 arises from the marked effect of the excitonic coupling between BChl units. The interchromophore CT excitation also plays a key role in the spectral redshift of B850. This CT effect can be effectively described using our dimer model. This suited characterization reveals that the pronounced CT effect originates from the characteristics of B850 that has closely spaced BChls as dimers. We highlight the importance of the refinement of the crystal structure with the use of quantum chemical methods for prediction of the spectrum.

Efficient propagation of the hierarchical equations of motion using the Tucker and hierarchical Tucker tensors

We develop new methods to efficiently propagate the hierarchical equations of motion (HEOM) by using the Tucker and hierarchical Tucker (HT) tensors to represent the reduced density operator and auxiliary density operators. We first show that by employing the split operator method, the specific structure of the HEOM allows a simple propagation scheme using the Tucker tensor. When the number of effective modes in the HEOM increases and the Tucker representation becomes intractable, the split operator method is extended to the binary tree structure of the HT representation. It is found that to update the binary tree nodes related to a specific effective mode, we only need to propagate a short matrix product state constructed from these nodes. Numerical results show that by further employing the mode combination technique commonly used in the multi-configuration time-dependent Hartree approaches, the binary tree representation can be applied to study excitation energy transfer dynamics in a fairly large system including over 10⁴ effective modes. The new methods may thus provide a promising tool in simulating quantum dynamics in condensed phases.

Structure of the far-red light utilizing photosystem I of Acaryochloris marina

Acaryochloris marina is one of the cyanobacterial species that can use far-red light to drive photochemical reactions for oxygenic photosynthesis. Here, we report the structure of A. marina photosystem I (PSI) reaction center, determined by cryo-electron microscopy at 2.58 Å resolution. The structure reveals an arrangement of electron carriers and light-harvesting pigments distinct from other type I reaction centers. The paired chlorophyll, or special pair (also referred to as P740 in this case), is a dimer of chlorophyll d and its epimer chlorophyll d'. The primary electron acceptor is pheophytin a, a metal-less chlorin. We show the architecture of this PSI reaction center is composed of 11 subunits and we identify key components that help explain how the low energy yield from far-red light is efficiently utilized for driving oxygenic photosynthesis.

Primary charge separation within the structurally symmetric tetrameric Chl2APAPBChl2B chlorophyll exciplex in photosystem I

In Photosystem I (PS I), the role of the accessory chlorophyll (Chl) molecules, Chl_2A and Chl_2B (also termed A_-1A and A_-1B), which are directly adjacent to the special pair P₇₀₀ and fork into the A- and B-branches of electron carriers, is incompletely understood. In this work, the Chl_2A and Chl_2B transient absorption ΔA₀(λ) at a time delay of 100 fs was identified by ultrafast pump-probe spectroscopy in three pairs of PS I complexes from Synechocystis sp. PCC 6803 with residues PsaA-N600 or PsaB-N582 (which ligate Chl_2B or Chl_2A through a H₂O molecule) substituted by Met, His, and Leu. The ΔA₀(λ) spectra were quantified using principal component analysis, the main component of which was interpreted as a mutation-induced shift of the equilibrium between the excited state of primary donor P₇₀₀^⁎ and the primary charge-separated state P₇₀₀⁺Chl₂^-. This equilibrium is shifted to the charge-separated state in wild-type PS I and to the excited P₇₀₀ in the PS I complexes with the substituted ligands to the Chl_2A and Chl_2B monomers. The results can be rationalized within the framework of an adiabatic model in which the P₇₀₀ is electronically coupled with the symmetrically arranged monomers Chl_2A and Chl_2B; such a structure can be considered a symmetric tetrameric exciplex Chl_2AP_AP_BChl_2B, in which the excited state (Chl_2AP_AP_BChl_2B)* is mixed with two charge-transfer states P₇₀₀⁺Chl_2A^- and P₇₀₀⁺Chl_2B^-. The electron redistribution between the two branches in favor of the A-branch apparently takes place in the picosecond time scale after reduction of the Chl_2A and Chl_2B monomers.

The structure of a red-shifted photosystem I reveals a red site in the core antenna

Photosystem I coordinates more than 90 chlorophylls in its core antenna while achieving near perfect quantum efficiency. Low energy chlorophylls (also known as red chlorophylls) residing in the antenna are important for energy transfer dynamics and yield, however, their precise location remained elusive. Here, we construct a chimeric Photosystem I complex in Synechocystis PCC 6803 that shows enhanced absorption in the red spectral region. We combine Cryo-EM and spectroscopy to determine the structure⁻function relationship in this red-shifted Photosystem I complex. Determining the structure of this complex reveals the precise architecture of the low energy site as well as large scale structural heterogeneity which is probably universal to all trimeric Photosystem I complexes. Identifying the structural elements that constitute red sites can expand the absorption spectrum of oxygenic photosynthetic and potentially modulate light harvesting efficiency.

Cyanobacterial photosystem I has a highly conserved core antenna consisting of eleven subunits and more than 90 chlorophylls. Here via CryoEM and spectroscopy, the authors determine the location of a red-shifted low-energy chlorophyll that allows harvesting of longer wavelengths of light.

Cryogenic Single-Molecule Spectroscopy of the Primary Electron Acceptor in the Photosynthetic Reaction Center

The photosynthetic reaction center (RC) converts light energy into electrochemical energy. The RC of heliobacteria (hRC) is a primitive homodimeric RC containing 58 bacteriochlorophylls and 2 chlorophyll as. The chlorophyll serves as the primary electron acceptor (Chl a-A₀) responsible for light harvesting and charge separation. The single-molecule spectroscopy of Chl a-A₀ can be used to investigate heterogeneities of the RC photochemical function, though the low fluorescence quantum yield (0.1%) makes it difficult. Here, we show the fluorescence excitation spectroscopy of individual Chl a-A₀s in single hRCs at 6 K. The fluorescence quantum yield and absorption cross section of Chl a-A₀ increase 2- and 4-fold, respectively, compared to those at room temperature. The two Chl a-A₀s in single hRCs are identified as two distinct peaks in the fluorescence excitation spectrum, exhibiting different excitation polarization dependences. The spectral changes caused by photobleaching indicate the energy transfer across subunits in the hRC.

Evidence that chlorophyll f functions solely as an antenna pigment in far-red-light photosystem I from Fischerella thermalis PCC 7521

The Photosystem I (PSI) reaction center in cyanobacteria is comprised of ~96 chlorophyll (Chl) molecules, including six specialized Chl molecules denoted Chl1A/Chl1B (P₇₀₀), Chl2A/Chl2B, and Chl3A/Chl3B that are arranged in two branches and function in primary charge separation. It has recently been proposed that PSI from Chroococcidiopsis thermalis (Nürnberg et al. (2018) Science 360, 1210-1213) and Fischerella thermalis PCC 7521 (Hastings et al. (2019) Biochim. Biophys. Acta 1860, 452-460) contain Chl f in the positions Chl2A/Chl2B. We tested this proposal by exciting RCs from white-light grown (WL-PSI) and far-red light grown (FRL-PSI) F. thermalis PCC 7521 with femtosecond pulses and analyzing the optical dynamics. If Chl f were in the position Chl2A/Chl2B in FRL-PSI, excitation at 740 nm should have produced the charge-separated state P₇₀₀⁺A₀^- followed by electron transfer to A₁ with a τ of ≤25 ps. Instead, it takes ~230 ps for the charge-separated state to develop because the excitation migrates uphill from Chl f in the antenna to the trapping center. Further, we observe a strong electrochromic shift at 685 nm in the final P₇₀₀⁺A₁^- spectrum that can only be explained if Chl a is in the positions Chl2A/Chl2B. Similar arguments rule out the presence of Chl f in the positions Chl3A/Chl3B; hence, Chl f is likely to function solely as an antenna pigment in FRL-PSI. We additionally report the presence of an excitonically coupled homo- or heterodimer of Chl f absorbing around 790 nm that is kinetically independent of the Chl f population that absorbs around 740 nm.

PSI-SMALP, a Detergent-free Cyanobacterial Photosystem I, Reveals Faster Femtosecond Photochemistry

Cyanobacterial photosystem I (PSI) functions as a light-driven cyt c6-ferredoxin/oxidoreductase located in the thylakoid membrane. In this work, the energy and charge transfer processes in PSI complexes isolated from Thermosynechococcus elongatus via conventional n-dodecyl-β-D-maltoside solubilization (DM-PSI) and a, to our knowledge, new detergent-free method using styrene-maleic acid copolymers (SMA-PSI) have been investigated by pump-to-probe femtosecond laser spectroscopy. In DM-PSI preparations excited at 740 nm, the excitation remained localized on the long-wavelength chlorophyll forms within 0.1-20 ps and revealed little or no charge separation and oxidation of the special pair, P700. The formation of ion-radical pair P700+A1- occurred with a characteristic time of 36 ps, being kinetically controlled by energy transfer from the long-wavelength chlorophyll to P700. Quite surprisingly, the detergent-free SMA-PSI complexes upon excitation by these long-wave pulses undergo an ultrafast (<100 fs) charge separation in ∼45% of particles. In the remaining complexes (∼55%), the energy transfer to P700 occurred at ∼36 ps, similar to the DM-PSI. Both isolation methods result in a trimeric form of PSI, yet the SMA-PSI complexes display a heterogenous kinetic behavior. The much faster rate of charge separation suggests the existence of an ultrafast pathway for charge separation in the SMA-PSI that may be disrupted during detergent isolation.

A high-throughput computational approach to UV-Vis spectra in protein mutants

In this work we present a high-throughput approach to the computation of absorption UV-Vis spectra tailored to mutagenesis studies. The scheme makes use of a single molecular dynamics trajectory of a reference (non-mutated) species. The shifts in absorption energy caused by a residue mutation are evaluated by building an effective potential of the environment and computing a correction term based on perturbation theory. The sampling is only performed in the phase space of the initial protein. We analyze the robustness of the method by comparing different approximations for the effective potential, the sampling of mutant residue geometries and observing the impact in the prediction of both bathocromic and hypsochromic shifts. As a test subject, we consider a red fluorescent protein variant with potential biotechnological applications.

Microenvironment control of porphyrin binding, organization, and function in peptide nanofiber assemblies

To take peptide materials from predominantly structural to functional assemblies, variations in cofactor binding sites must be engineered and controlled. Here, we have employed the peptide sequence c16-AHX3K3-CO2H where X3 represents the aliphatic structural component of the peptide design that dictates β-sheet formation and upon self-assembly yields a change in the overall microenvironment surrounding the Zn protoporphyrin IX ((PPIX)Zn) binding site. All peptides studied yield β-sheet rich nanofibers highlighting the materials' resiliency to amino acid substitution. We highlight that the (PPIX)Zn binding constants correlate strongly with amino acid side chain volume, where X = L or I yields the lowest dissociation constant values (KD). The resulting microenvironment highlights the materials' ability to control interchromophore electronic interactions such that slip-stacked cofacial arrangements are observed via exciton splitting in UV/visible and circular dichroism spectroscopy. Steady state and time-resolved photoluminescence suggests that greater interchromophore packing yields larger excimer populations and corresponding longer excimer association lifetimes (τA) which directly translates to shorter exciton diffusion lengths. In comparison to synthetic porphyrin molecular assemblies, this work demonstrates the ability to employ the peptide assembly to modulate the degree of cofactor arrangement, extent of excimer formation, and the exciton hopping rates all while in a platform amenable for producing polymer-like materials.

A theoretical study on the dynamics of light harvesting in the dimeric photosystem II core complex: regulation and robustness of energy transfer pathways

Here we present our theoretical investigations into the light reaction in the dimeric photosystem II (PSII) core complex. An effective model for excitation energy transfer (EET) and primary charge separation (CS) in the PSII core complex was developed, with model parameters constructed based on molecular dynamics (MD) simulation data. Compared to experimental results, we demonstrated that this model faithfully reproduces the absorption spectra of the RC and core light-harvesting complexes (CP43 and CP47) as well as the full EET dynamics among the chromophores in the PSII core complex. We then applied master equation simulations and network analysis to investigate detailed EET plus CS dynamics in the system, allowing us to identify key EET pathways and produce a coarse-grained cluster model for the light reaction in the dimeric PSII core complex. We show that non-equilibrium energy transfer channels play important roles in the efficient light harvesting process and that multiple EET pathways exist between subunits of PSII to ensure the robustness of light harvesting in the system. Furthermore, we revealed that inter-monomer energy transfer dominated by the coupling between the two CLA625 molecules enables efficient energy exchange between two CP47s in the dimeric PSII core complex, which leads to significant energy pooling in the CP47 domain during the light reaction. Our study provides a blueprint for the design of light harvesting in the PSII core and show that a structure-based approach using molecular dynamics simulations and quantum chemistry calculations can be effectively utilized to elucidate the dynamics of light harvesting in complex photosynthetic systems.

Redox-state dependent blinking of single photosystem I trimers at around liquid-nitrogen temperature

Efficient light harvesting in a photosynthetic antenna system is disturbed by a ragged and fluctuating energy landscape of the antenna pigments in response to the conformation dynamics of the protein. This situation is especially pronounced in Photosystem I (PSI) containing red shifted chlorophylls (red Chls) with the excitation energy much lower than the primary donor. The present study was conducted to clarify light-harvesting dynamics of PSI isolated from Synechocystis sp. PCC6803 by using single-molecule spectroscopy at liquid‑nitrogen temperatures. Fluorescence emission at around 720 nm from the red Chls in single PSI trimers was monitored at 80-100 K. Intermittent variations in the emission intensities, so-called blinking, were frequently observed. Its time scale lay in several tens of seconds. The blinking amplitude depended on the redox state of the phylloquinone (A₁). Electrochromic shifts of Chls induced by the negative charge on A₁ were calculated based on the X-ray crystallographic structure. A Chl molecule, Chl-A839 (numbering according to PDB 5OY0), bound near A₁ was found to have a large electrochromic shift. This Chl has strong exciton coupling with neighboring Chl (A838) whose site energy was predicted to be determined by interaction with an arginine residue (ArgF84) [Adolphs et al., 2010]. A possible scenario of the blinking was proposed. Conformational fluctuations of ArgF84 seesaw the excitation-energy of Chl-A838, which perturbs the branching ratio of excitation-energy between the red Chl and the cationic form of P700 as a quencher. The electrochromic shift of Chl-A839 enhances the effect of the conformation dynamics of ArgF84.

Ultrafast energy transfer within the photosystem II core complex

We report 2D electronic spectroscopy on the photosystem II core complex (PSII CC) at 77 K under different polarization conditions. A global analysis of the high time-resolution 2D data shows rapid, sub-100 fs energy transfer within the PSII CC. It also reveals the 2D spectral signatures of slower energy equilibration processes occurring on several to hundreds of picosecond time scales that are consistent with previous work. Using a recent structure-based model of the PSII CC [Y. Shibata, S. Nishi, K. Kawakami, J. R. Shen and T. Renger, J. Am. Chem. Soc., 2013, 135, 6903], we simulate the energy transfer in the PSII CC by calculating auxiliary time-resolved fluorescence spectra. We obtain the observed sub-100 fs evolution, even though the calculated electronic energy shows almost no dynamics at early times. On the other hand, the electronic-vibrational interaction energy increases considerably over the same time period. We conclude that interactions with vibrational degrees of freedom not only induce population transfer between the excitonic states in the PSII CC, but also reshape the energy landscape of the system. We suggest that the experimentally observed ultrafast energy transfer is a signature of excitonic-polaron formation.

Exciton states and optical properties of the CP26 photosynthetic protein

The photosynthetic complex CP26, one of the minor antennae of the photosystem II, plays an important role in regulation of the excitation energy transfer in the PSII. Due to instability during isolation and purification, it remained poorly studied from the viewpoint of theoretical chemistry because of the absence of X-ray crystallography data. In this work, using the recently determined three-dimensional structure of the complex we apply the quantum chemical approach to study the properties of exciton states in it. Spectral properties, structure of exciton states and roles of the pigments in the complex and photosystem II are discussed.

Oxidative species-induced excitonic transport in tubulin aromatic networks: Potential implications for neurodegenerative disease

Oxidative stress is a pathological hallmark of neurodegenerative tauopathic disorders such as Alzheimer's disease and Parkinson's disease-related dementia, which are characterized by altered forms of the microtubule-associated protein (MAP) tau. MAP tau is a key protein in stabilizing the microtubule architecture that regulates neuron morphology and synaptic strength. When MAP tau is degraded in tauopathic disorders, neuron dysfunction results. The precise role of reactive oxygen species (ROS) in the tauopathic disease process, however, is poorly understood. Classically, mitochondrial dysfunction has been viewed as the major source of oxidative stress and has been shown to precede tau and amyloid pathology in various dementias, but the exact mechanisms are not clear. It is known that the production of ROS by mitochondria can result in ultraweak photon emission (UPE) within cells. While of low intensity, surrounding proteins within the cytosol can still absorb these energetic photons via aromatic amino acids (e.g., tryptophan and tyrosine). One likely absorber of these photons is the microtubule cytoskeleton, as it forms a vast network spanning neurons, is highly co-localized with mitochondria, and shows a high density of aromatic amino acids. Functional microtubule networks may traffic this ROS-generated endogenous photon energy for cellular signaling, or they may serve as dissipaters/conduits of such energy to protect the cell from potentially harmful effects. Experimentally, after in vitro exposure to exogenous photons, microtubules have been shown to reorient and reorganize in a dose-dependent manner with the greatest effect being observed around 280nm, in the tryptophan and tyrosine absorption range. In this paper, recent modeling efforts based on ambient temperature experiment are presented, showing that tubulin polymers can feasibly absorb and channel these photoexcitations via resonance energy transfer, on the order of dendritic length scales and neuronal fine structure. Since microtubule networks are compromised in tauopathic diseases such as Alzheimer's and Parkinson's dementias, patients with these illnesses would be unable to support effective channeling of these photons for signaling or dissipation. Consequent emission surplus due to increased UPE production or decreased ability to absorb and transfer may lead to increased cellular oxidative damage, thus hastening the neurodegenerative process.

Toward an Understanding of the Excitonic Structure of the CP47 Antenna Protein Complex of Photosystem II Revealed via Circularly Polarized Luminescence

Identification of the lowest energy pigments in the photosynthetic CP47 antenna protein complex of Photosystem II (PSII) is essential for understanding its excitonic structure, as well as excitation energy pathways in the PSII core complex. Unfortunately, there is no consensus concerning the nature of the low-energy state(s), nor chlorophyll (Chl) site energies in this important photosynthetic antenna. Although we raised concerns regarding the estimations of Chl site energies obtained from modeling studies of various types of CP47 optical spectra [Reinot, T; et al., Anal. Chem. Insights 2016, 11, 35-48] recent new assignments imposed by the shape of the circularly polarized luminescence (CPL) spectrum [Hall, J.; et al., Biochim. Biophys. Acta 2016, 1857, 1580-1593] necessitate our comments. We demonstrate that other combinations of low-energy Chls provide equally good or improved simultaneous fits of various optical spectra (absorption, emission, CPL, circular dichroism, and nonresonant hole-burned spectra), but more importantly, we expose the heterogeneous nature of the recently studied complexes and argue that the published composite nature of the CPL (contributed to by CPL₆₈₅, CPL₆₉₁, and CPL₆₉₅) does not represent an intact CP47 protein. A positive CPL₆₉₅ is extracted for the intact protein, which, when simultaneously fitted with multiple other optical spectra, provides new information on the excitonic structure of intact and destabilized CP47 complexes and their lowest energy state(s).

Applications of large-scale density functional theory in biology

Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality.

Light Absorption and Energy Transfer in the Antenna Complexes of Photosynthetic Organisms

The process of photosynthesis is initiated by the capture of sunlight by a network of light-absorbing molecules (chromophores), which are also responsible for the subsequent funneling of the excitation energy to the reaction centers. Through evolution, genetic drift, and speciation, photosynthetic organisms have discovered many solutions for light harvesting. In this review, we describe the underlying photophysical principles by which this energy is absorbed, as well as the mechanisms of electronic excitation energy transfer (EET). First, optical properties of the individual pigment chromophores present in light-harvesting antenna complexes are introduced, and then we examine the collective behavior of pigment-pigment and pigment-protein interactions. The description of energy transfer, in particular multichromophoric antenna structures, is shown to vary depending on the spatial and energetic landscape, which dictates the relative coupling strength between constituent pigment molecules. In the latter half of the article, we focus on the light-harvesting complexes of purple bacteria as a model to illustrate the present understanding of the synergetic effects leading to EET optimization of light-harvesting antenna systems while exploring the structure and function of the integral chromophores. We end this review with a brief overview of the energy-transfer dynamics and pathways in the light-harvesting antennas of various photosynthetic organisms.

Challenges facing an understanding of the nature of low-energy excited states in photosynthesis

While the majority of the photochemical states and pathways related to the biological capture of solar energy are now well understood and provide paradigms for artificial device design, additional low-energy states have been discovered in many systems with obscure origins and significance. However, as low-energy states are naively expected to be critical to function, these observations pose important challenges. A review of known properties of low energy states covering eight photochemical systems, and options for their interpretation, are presented. A concerted experimental and theoretical research strategy is suggested and outlined, this being aimed at providing a fully comprehensive understanding.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.