Assignment of the Qy absorption spectrum of photosystem-I from Thermosynechococcus elongatus based on CAM-B3LYP calculations at the PW91-optimized protein structure

Structure of cryptophyte photosystem II-light-harvesting antennae supercomplex

Cryptophytes are ancestral photosynthetic organisms evolved from red algae through secondary endosymbiosis. They have developed alloxanthin-chlorophyll a/c2-binding proteins (ACPs) as light-harvesting complexes (LHCs). The distinctive properties of cryptophytes contribute to efficient oxygenic photosynthesis and underscore the evolutionary relationships of red-lineage plastids. Here we present the cryo-electron microscopy structure of the Photosystem II (PSII)-ACPII supercomplex from the cryptophyte Chroomonas placoidea. The structure includes a PSII dimer and twelve ACPII monomers forming four linear trimers. These trimers structurally resemble red algae LHCs and cryptophyte ACPI trimers that associate with Photosystem I (PSI), suggesting their close evolutionary links. We also determine a Chl a-binding subunit, Psb-γ, essential for stabilizing PSII-ACPII association. Furthermore, computational calculation provides insights into the excitation energy transfer pathways. Our study lays a solid structural foundation for understanding the light-energy capture and transfer in cryptophyte PSII-ACPII, evolutionary variations in PSII-LHCII, and the origin of red-lineage LHCIIs.

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.

Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach

In this paper, we present dyadic adaptive HOPS (DadHOPS), a new method for calculating linear absorption spectra for large molecular aggregates. This method combines the adaptive HOPS (adHOPS) framework, which uses locality to improve computational scaling, with the dyadic HOPS method previously developed to calculate linear and nonlinear spectroscopic signals. To construct a local representation of dyadic HOPS, we introduce an initial state decomposition that reconstructs the linear absorption spectra from a sum over locally excited initial conditions. We demonstrate the sum over initial conditions can be efficiently Monte Carlo sampled and that the corresponding calculations achieve size-invariant [i.e., O(1)] scaling for sufficiently large aggregates while trivially incorporating static disorder in the Hamiltonian. We present calculations on the photosystem I core complex to explore the behavior of the initial state decomposition in complex molecular aggregates as well as proof-of-concept DadHOPS calculations on an artificial molecular aggregate inspired by perylene bis-imide to demonstrate the size-invariance of the method.

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.

Quasiparticle Self-Consistent GW-Bethe-Salpeter Equation Calculations for Large Chromophoric Systems

The GW-Bethe–Salpeter equation (BSE) method is promising for calculating the low-lying excitonic states of molecular systems. However, so far it has only been applied to rather small molecules and in the commonly implemented diagonal approximations to the electronic self-energy, it depends on a mean-field starting point. We describe here an implementation of the self-consistent and starting-point-independent quasiparticle self-consistent (qsGW)-BSE approach, which is suitable for calculations on large molecules. We herein show that eigenvalue-only self-consistency can lead to an unfaithful description of some excitonic states for chlorophyll dimers while the qsGW-BSE vertical excitation energies (VEEs) are in excellent agreement with spectroscopic experiments for chlorophyll monomers and dimers measured in the gas phase. Furthermore, VEEs from time-dependent density functional theory calculations tend to disagree with experimental values and using different range-separated hybrid (RSH) kernels does change the VEEs by up to 0.5 eV. We use the new qsGW-BSE implementation to calculate the lowest excitation energies of the six chromophores of the photosystem II (PSII) reaction center (RC) with nearly 2000 correlated electrons. Using more than 11,000 (6000) basis functions, the calculation could be completed in less than 5 (2) days on a single modern compute node. In agreement with previous TD-DFT calculations using RSH kernels on models that also do not include environmental effects, our qsGW-BSE calculations only yield states with local characters in the low-energy spectrum of the hexameric complex. Earlier works with RSH kernels have demonstrated that the protein environment facilitates the experimentally observed interchromophoric charge transfer. Therefore, future research will need to combine correlation effects beyond TD-DFT with an explicit treatment of environmental electrostatics.

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.

Structure-based model of fucoxanthin-chlorophyll protein complex: Calculations of chlorophyll electronic couplings

Diatoms are a group of marine algae that are responsible for a significant part of global oxygen production. Adapted to life in an aqueous environment dominated by the blue-green light, their major light-harvesting antennae-fucoxanthin-chlorophyll protein complexes (FCPs)-exhibit different pigment compositions than of plants. Despite extensive experimental studies, until recently the theoretical description of excitation energy dynamics in these complexes was limited by the lack of high-resolution structural data. In this work, we use the recently resolved crystallographic information of the FCP complex from Phaeodactylum tricornutum diatom [Wang et al., Science 363, 6427 (2019)] and quantum chemistry-based calculations to evaluate the chlorophyll transition dipole moments, atomic transition charges from electrostatic potential, and the inter-chlorophyll couplings in this complex. The obtained structure-based excitonic couplings form the foundation for any modeling of stationary or time-resolved spectroscopic data. We also calculate the inter-pigment Förster energy transfer rates and identify two quickly equilibrating chlorophyll clusters.

Density functionals with asymptotic-potential corrections are required for the simulation of spectroscopic properties of materials

Five effects of correction of the asymptotic potential error in density functionals are identified that significantly improve calculated properties of molecular excited states involving charge-transfer character. Newly developed materials-science computational methods are used to demonstrate how these effects manifest in materials spectroscopy. Connection is made considering chlorophyll-a as a paradigm for molecular spectroscopy, 22 iconic materials as paradigms for 3D materials spectroscopy, and the VN - defect in hexagonal boron nitride as an example of the spectroscopy of defects in 2D materials pertaining to nanophotonics. Defects can equally be thought of as being "molecular" and "materials" in nature and hence bridge the relms of molecular and materials spectroscopies. It is concluded that the density functional HSE06, currently considered as the standard for accurate calculations of materials spectroscopy, should be replaced, in most instances, by the computationally similar but asymptotically corrected CAM-B3LYP functional, with some specific functionals for materials-use only providing further improvements.

Bridging Carotenoid-to-Bacteriochlorophyll Energy Transfer of Purple Bacteria LH2 With Temperature Variations: Insights From Conformational Changes

Photosynthesis is a key process for converting light energy into chemical energy and providing food for lives on Earth. Understanding the mechanism for the energy transfers could provide insights into regulating energy transfers in photosynthesis and designing artificial photosynthesis systems. Many efforts have been devoted to exploring the mechanism of temperature variations affecting the excitonic properties of LH2. In this study, we performed all-atom molecular dynamics (MD) simulations and quantum mechanics calculations for LH2 complex from purple bacteria along with its membrane environment under three typical temperatures: 270, 300, and 330 K. The structural analysis from validated MD simulations showed that the higher temperature impaired interactions at N-terminus of both α and β polypeptide helices and led to the dissociation of this hetero polypeptide dimer. Rhodopin-β-D-glucosides (RG1) moved centripetally with α polypeptide helices when temperature increased and enlarged their distances with bacteriochlorophylls molecules that have the absorption peak at 850 nm (B850), which resulted in reducing the coupling strengths between RG1 and B850 molecules. The present study reported a cascading mechanism for temperature regulating the energy transfers in LH2 of purple bacteria.

Two-Dimensional Electronic Spectroscopy of a Minimal Photosystem I Complex Reveals the Rate of Primary Charge Separation

Photosystem I (PSI), found in all oxygenic photosynthetic organisms, uses solar energy to drive electron transport with nearly 100% quantum efficiency, thanks to fast energy transfer among antenna chlorophylls and charge separation in the reaction center. There is no complete consensus regarding the kinetics of the elementary steps involved in the overall trapping, especially the rate of primary charge separation. In this work, we employed two-dimensional coherent electronic spectroscopy to follow the dynamics of energy and electron transfer in a monomeric PSI complex from Synechocystis PCC 6803, containing only subunits A-E, K, and M, at 77 K. We also determined the structure of the complex to 4.3 Å resolution by cryoelectron microscopy with refinements to 2.5 Å. We applied structure-based modeling using a combined Redfield-Förster theory to compute the excitation dynamics. The absorptive 2D electronic spectra revealed fast excitonic/vibronic relaxation on time scales of 50-100 fs from the high-energy side of the absorption spectrum. Antenna excitations were funneled within 1 ps to a small pool of chlorophylls absorbing around 687 nm, thereafter decaying with 4-20 ps lifetimes, independently of excitation wavelength. Redfield-Förster energy transfer computations showed that the kinetics is limited by transfer from these red-shifted pigments. The rate of primary charge separation, upon direct excitation of the reaction center, was determined to be 1.2-1.5 ps^-1. This result implies activationless electron transfer in PSI.

Symmetry breaking in photosystem I: ultrafast optical studies of variants near the accessory chlorophylls in the A- and B-branches of electron transfer cofactors

Femtosecond absorption spectroscopy of Photosystem I (PS I) complexes from the cyanobacterium Synechocystis sp. PCC 6803 was carried out on three pairs of complementary amino acid substitutions located near the second pair of chlorophyll molecules Chl_2A and Chl_2B (also termed A_-1A and A_-1B). The absorption dynamics at delays of 0.1-500 ps were analyzed by decomposition into discrete decay-associated spectra and continuously distributed exponential components. The multi-exponential deconvolution of the absorption changes revealed that the electron transfer reactions in the PsaA-N600M, PsaA-N600H, and PsaA-N600L variants near the B-branch of cofactors are similar to those of the wild type, while the PsaB-N582M, PsaB-N582H, and PsaB-N582L variants near the A-branch of cofactors cause significant alterations of the photochemical processes, making them heterogeneous and poorly described by a discrete exponential kinetic model. A redistribution of the unpaired electron between the second and the third monomers Chl_2A/Chl_2B and Chl_3A/Chl_3B was identified in the time range of 9-20 ps, and the subsequent reduction of A₁ was identified in the time range of 24-70 ps. In the PsaA-N600L and PsaB-N582H/L variants, the reduction of A₁ occurred with a decreased quantum yield of charge separation. The decreased quantum yield correlates with a slowing of the phylloquinone A₀ → A₁ reduction, but not with the initial transient spectra measured at the shortest time delay. The results support a branch competition model, where the electron is sheared between Chl_2A-Chl_3A and Chl_2B-Chl_3B cofactors before its transfer to phylloquinone in either A_1A or A_1B sites.

A dimeric chlorophyll electron acceptor differentiates type I from type II photosynthetic reaction centers

This research addresses one of the most compelling issues in the field of photosynthesis, namely, the role of the accessory chlorophyll molecules in primary charge separation. Using a combination of empirical and computational methods, we demonstrate that the primary acceptor of photosystem (PS) I is a dimer of accessory and secondary chlorophyll molecules, Chl_2A and Chl_3A, with an asymmetric electron charge density distribution. The incorporation of highly coupled donors and acceptors in PS I allows for extensive delocalization that prolongs the lifetime of the charge-separated state, providing for high quantum efficiency. The discovery of this motif has widespread implications ranging from the evolution of naturally occurring reaction centers to the development of a new generation of highly efficient artificial photosynthetic systems.

Video abstract

Accurate prediction of the properties of materials using the CAM-B3LYP density functional

Density functionals with asymptotic corrections to the long-range potential provide entry-level methods for calculations on molecules that can sustain charge transfer, but similar applications in materials science are rare. We describe an implementation of the CAM-B3LYP range-separated functional within the Vienna Ab-initio Simulation Package (VASP) framework, together with its analytical functional derivatives. Results obtained for eight representative materials: aluminum, diamond, graphene, silicon, NaCl, MgO, 2D h-BN, and 3D h-BN, indicate that CAM-B3LYP predictions embody mean-absolute deviations (MAD) compared to HSE06 that are reduced by a factor of six for lattice parameters, four for quasiparticle band gaps, three for the lowest optical excitation energies, and six for exciton binding energies. Further, CAM-B3LYP appears competitive compared to ab initio G₀ W₀ and Bethe-Salpeter equation approaches. The CAM-B3LYP implementation in VASP was verified by comparison of optimized geometries and reaction energies for isolated molecules taken from the ACCDB database, evaluated in large periodic unit cells, to analogous results obtained using Gaussian basis sets. Using standard GW pseudopotentials and energy cutoffs for the plane-wave calculations and the aug-cc-pV5Z basis set for the atomic-basis ones, the MAD in energy for 1738 chemical reactions was 0.34 kcal mol^-1 , while for 480 unique bond lengths this was 0.0036 Å; these values reduced to 0.28 kcal mol^-1 (largest error 0.94 kcal mol^-1 ) and 0.0009 Å by increasing the plane-wave cutoff energy to 850 eV.

Quantum Chemical Simulation of the Qy Absorption Spectrum of Zn Chlorin Aggregates for Artificial Photosynthesis

Zn chlorin (Znchl) is easy to synthesize and has similar optical properties to those of bacteriochlorophyll c in the nature, which is expected to be used as a light-harvesting antenna system in artificial photosynthesis. In order to further explore the optical characteristics of Znchl, various sizes of a parallel layered Znchl-aggregate model and the THF-Znchl explicit solvent monomer model were constructed in this study, and their Qy excited state properties were simulated by using time-dependent density functional theory (TDDFT) and exciton theory. For the Znchl monomer, with a combination of the explicit solvent model and the implicit solvation model based on density (SMD), the calculated Qy excitation energy agreed very well with the experimental one. The Znchl aggregates may be simplified to a Zn36 model to reproduce the experimental Qy absorption spectrum by the Förster coupling theory. The proposed Znchl aggregate model provides a good foundation for the future exploration of other properties of Znchl and simulations of artificial light-harvesting antennas. The results also indicate that J-aggregrates along z-direction, due to intermolecular coordination bonds, are the dominant factor in extending the Qy band of Znchl into the near infrared region.

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.

Generation of ion-radical chlorophyll states in the light-harvesting antenna and the reaction center of cyanobacterial photosystem I

The energy and charge-transfer processes in photosystem I (PS I) complexes isolated from cyanobacteria Thermosynechococcus elongatus and Synechocystis sp. PCC 6803 were investigated by pump-to-probe femtosecond spectroscopy. The formation of charge-transfer (CT) states in excitonically coupled chlorophyll a complexes (exciplexes) was monitored by measuring the electrochromic shift of β-carotene in the spectral range 500-510 nm. The excitation of high-energy chlorophyll in light-harvesting antenna of both species was not accompanied by immediate appearance of an electrochromic shift. In PS I from T. elongatus, the excitation of long-wavelength chlorophyll (LWC) caused a pronounced electrochromic effect at 502 nm assigned to the appearance of CT states of chlorophyll exciplexes. The formation of ion-radical pair P₇₀₀⁺A₁^- at 40 ps was limited by energy transfer from LWC to the primary donor P₇₀₀ and accompanied by carotenoid bleach at 498 nm. In PS I from Synechocystis 6803, the excitation at 720 nm produced an immediate bidentate bleach at 690/704 nm and synchronous carotenoid response at 508 nm. The bidentate bleach was assigned to the formation of primary ion-radical state P_B⁺Chl_2B^-, where negative charge is localized predominantly at the accessory chlorophyll molecule in the branch B, Chl_2B. The following decrease of carotenoid signal at ~ 5 ps was ascribed to electron transfer to the more distant molecule Chl_3B. The reduction of phylloquinone in the sites A_1A and A_1B was accompanied by a synchronous blue-shift of the carotenoid response to 498 nm, pointing to fast redistribution of unpaired electron between two branches in favor of the state P_B⁺A_1A^-.

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.

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.

Structural basis for the adaptation and function of chlorophyll f in photosystem I

Chlorophylls (Chl) play pivotal roles in energy capture, transfer and charge separation in photosynthesis. Among Chls functioning in oxygenic photosynthesis, Chl f is the most red-shifted type first found in a cyanobacterium Halomicronema hongdechloris. The location and function of Chl f in photosystems are not clear. Here we analyzed the high-resolution structures of photosystem I (PSI) core from H. hongdechloris grown under white or far-red light by cryo-electron microscopy. The structure showed that, far-red PSI binds 83 Chl a and 7 Chl f, and Chl f are associated at the periphery of PSI but not in the electron transfer chain. The appearance of Chl f is well correlated with the expression of PSI genes induced under far-red light. These results indicate that Chl f functions to harvest the far-red light and enhance uphill energy transfer, and changes in the gene sequences are essential for the binding of Chl f.

Chlorophyll f (Chl f) is the most red-shifted Chl in oxygenic photosynthesis but its localization in photosystem I (PSI) has been unknown so far. Here the authors determine the cryo-EM structures of PSI complexes from a Chl f-containing cyanobacterium grown either under white light or far-red light conditions and identify seven Chls f in the far-red light PSI structure, whereas PSI from cells grown under white light contains only Chl a.

Single Molecule Fluorescence Spectroscopy of PSI Trimers from Arthrospira platensis: A Computational Approach

Based on single molecule spectroscopy analysis and our preliminary theoretical studies, the linear and fluorescence spectra of the PSI trimer from Arthrospira platensis with different realizations of the static disorder were modeled at cryogenic temperature. Considering the previously calculated spectral density of chlorophyll, an exciton model for the PSI monomer and trimer including the red antenna states was developed taking into account the supposed similarity of PSI antenna structures from Thermosynechococcus e., Synechocystis sp. PCC6803, and Arthrospira platensis. The red Chls in the PSI monomer were assumed to be in the nearest proximity of the reaction center. The PSI trimer model allowed the simulation of experimentally measured zero phonon line distribution of the red states considering a weak electron-phonon coupling for the antenna exciton states. However, the broad absorption and fluorescence spectra of an individual emitter at 760 nm were calculated by adjusting the Huang-Rhys factors of the chlorophyll lower phonon modes assuming strong electron-phonon coupling.

Unipolar to ambipolar semiconductivity switching in charge transfer cocrystals of 2,7-di-tert-butylpyrene

Charge transfer cocrystals of 2,7-di-tert-butylpyrene donor and tetracyanoquinodimethane, tetracyanobenzene and 1,3-dinitrobenzene acceptor exhibited switchable semi-conductivity.

Screening for the selective inhibitors of MMP-9 from natural products based on pharmacophore modeling and molecular docking in combination with bioassay experiment, hybrid QM/MM calculation, and MD simulation

Matrix metalloproteinase-9 (MMP-9) has been considered as an attractive target involving cancer therapy. In this study, the 3D QSAR pharmacophore model of MMP-9 inhibitors is built, and its reliability is subsequently validated based on different methods. The built pharmacophore model consists of the four chemical features, including two hydrogen bond acceptors (HBA), one hydrophobic (HY), and one ring aromatic (RA). Among them, both HY and RA are found to be especially important features because they involve the interactions of inhibitors with the S1' pocket of MMP-9, which determines the selectivity of MMP-9 inhibitors. By combining pharmacophore model with molecular docking, the virtual screening is carried out to identify the selective MMP-9 inhibitors from natural products. The four potential selective MMP-9 inhibitors of natural products are found. One of them was used to carry out the bioassay experiment inhibiting MMP-9, and the estimated IC₅₀ value of only 26.94 µM clearly shows its strongly inhibitory activity; besides, both the hybrid quantum mechanics/molecular mechanics (QM/MM) calculation and the molecular dynamics simulation are performed to examine the reliability regarding the binding mode of this inhibitor with MMP-9 active sites predicted by molecular docking. All the screened four natural products are found to well bind with the MMP-9 active sites by different kinds of interactions. Finally, the ADMET properties of screened four natural products are assessed. These screened MMP-9 inhibitors of natural products could be used as the lead compounds to perform structural modifications and optimizations in the future work. Communicated by Ramaswamy H. Sarma.

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.

Towards an ab initio description of the optical spectra of light-harvesting antennae: application to the CP29 complex of photosystem II

Light-harvesting pigment-protein complexes (PPC) represent the fundamental units through which the photosynthetic organisms absorb sunlight and funnel the energy to the reaction centre for carrying out the primary energy conversion reactions of photosynthesis. Here we apply a multiscale computational strategy to a specific PPC present in the photosystem II of plants and algae (CP29) to investigate in what detail should the environment effects due to protein and membrane/solvent be included for an accurate description of optical spectra. We find that a refinement of the crystal structure is needed before any meaningful quantum chemical calculations of pigment transition energies can be performed. For this purpose we apply classical molecular dynamics simulations of the PPC within its natural environment and we perform ab initio computations of the exciton Hamiltonian of the complex, including the environment either implicitly by the polarizable continuum model (PCM) or explicitly using the polarizable QM/MM methodology (MMPol). However, PCM essentially leads to an unspecific redshift of all transition energies, and MMPol is able to reveal site-specific changes in the optical properties of the pigments. Based on the latter and the excitonic couplings obtained within a polarizable QM/MM methodology, optical spectra are calculated, which are in good qualitative agreement with experimental data. A weakness of the approach is however found in the overestimation of the fluctuations of the excitonic parameters of the pigments along the MD trajectory. An explanation for such a finding in terms of the limits of the force fields commonly used for protein cofactors is presented and discussed.

Putting David Craig’s Legacy to Work in Nanotechnology and Biotechnology

David Craig (1919–2015) left us with a lasting legacy concerning basic understanding of chemical spectroscopy and bonding. This is expressed in terms of some of the recent achievements of my own research career, with a focus on integration of Craig’s theories with those of Noel Hush to solve fundamental problems in photosynthesis, molecular electronics (particularly in regard to the molecules synthesized by Maxwell Crossley), and self-assembled monolayer structure and function. Reviewed in particular is the relation of Craig’s legacy to: the 50-year struggle to assign the visible absorption spectrum of arguably the world’s most significant chromophore, chlorophyll; general theories for chemical bonding and structure extending Hush’s adiabatic theory of electron-transfer processes; inelastic electron-tunnelling spectroscopy (IETS); chemical quantum entanglement and the Penrose–Hameroff model for quantum consciousness; synthetic design strategies for NMR quantum computing; Gibbs free-energy measurements and calculations for formation and polymorphism of organic self-assembled monolayers on graphite surfaces from organic solution; and understanding the basic chemical processes involved in the formation of gold surfaces and nanoparticles protected by sulfur-bound ligands, ligands whose form is that of Au0-thiyl rather than its commonly believed AuI-thiolate tautomer.

Explicit calculation of the excited electronic states of the photosystem II reaction centre

The excited states of the photosystem II reaction centre cofactors have been calculated as a single “supermolecule”. Charge transfer states are shown to be dependent on electrostatic environment.

Investigations on the mechanisms of interactions between matrix metalloproteinase 9 and its flavonoid inhibitors using a combination of molecular docking, hybrid quantum mechanical/molecular mechanical calculations, and molecular dynamics simulations

Many experimental studies have found that flavonoids can inhibit the activities of matrix metalloproteinases (MMPs), but the relevant mechanisms are still unclear. In this paper, the interaction mechanisms of MMP-9 with its five flavonoid inhibitors are investigated using a combination of molecular docking, hybrid quantum mechanical and molecular mechanical (QM/MM) calculations, and molecular dynamics simulations. The molecular dynamics simulation results show a good linear correlation between the calculated binding free energies of QM/MM−Poisson–Boltzmann surface area (PBSA) and the experimental −log(EC50) regarding the studied five flavonoids on MMP-9 inhibition in explicit solvent. It is found that compared with the MM−PBSA method, the QM/MM−PBSA method can obviously improve the accuracy for the calculated binding free energies. The predicted binding modes of the five flavonoid−MMP-9 complexes reveal that the different hydrogen bond networks can form besides producing the Zn−O coordination bonds, which can reasonably explain previous experimental results. The agreement between our calculated results and the previous experimental facts indicates that the force field parameters used here are effective and reliable for investigating the systems of flavonoid−MMP-9 interactions, and thus, these simulations and analyses could be reproduced for the other related systems involving protein−ligand interactions. This paper may be helpful for designing the new MMP-9 inhibitors having higher biological activities by carrying out the structural modifications of flavonoid molecules.

Assignment of the Q-bands of the chlorophylls: coherence loss via Qx - Qy mixing

We provide a new and definitive spectral assignment for the absorption, emission, high-resolution fluorescence excitation, linear dichroism, and/or magnetic circular dichroism spectra of 32 chlorophyllides in various environments. This encompases all data used to justify previous assignments and provides a simple interpretation of unexplained complex decoherence phenomena associated with Q_x → Q_y relaxation. Whilst most chlorophylls conform to the Gouterman model and display two independent transitions Q_x (S₂) and Q_y (S₁), strong vibronic coupling inseparably mixes these states in chlorophyll-a. This spreads x-polarized absorption intensity over the entireQ-band system to influence all exciton-transport, relaxation and coherence properties of chlorophyll-based photosystems. The fraction of the total absorption intensity attributed to Q_x ranges between 7% and 33%, depending on chlorophyllide and coordination, and is between 10% and 25% for chlorophyll-a. CAM-B3LYP density-functional-theory calculations of the band origins, relative intensities, vibrational Huang-Rhys factors, and vibronic coupling strengths fully support this new assignment.

Understanding photosynthetic light-harvesting: a bottom up theoretical approach

We discuss a bottom up approach for modeling photosynthetic light-harvesting. Methods are reviewed for a full structure-based parameterization of the Hamiltonian of pigment-protein complexes (PPCs). These parameters comprise (i) the local transition energies of the pigments in their binding sites in the protein, the site energies; (ii) the couplings between optical transitions of the pigments, the excitonic couplings; and (iii) the spectral density characterizing the dynamic modulation of pigment transition energies and excitonic couplings by protein vibrations. Starting with quantum mechanics perturbation theory, we provide a microscopic foundation for the standard PPC Hamiltonian and relate the expressions obtained for its matrix elements to quantities that can be calculated with classical molecular mechanics/electrostatics approaches including the whole PPC in atomic detail and using charge and transition densities obtained with quantum chemical calculations on the isolated building blocks of the PPC. In the second part of this perspective, the Hamiltonian is utilized to describe the quantum dynamics of excitons. Situations are discussed that differ in the relative strength of excitonic and exciton-vibrational coupling. The predictive power of the approaches is demonstrated in application to different PPCs, and challenges for future work are outlined.

Quantum chemical description of absorption properties and excited-state processes in photosynthetic systems

The theoretical description of the initial steps in photosynthesis has gained increasing importance over the past few years. This is caused by more and more structural data becoming available for light-harvesting complexes and reaction centers which form the basis for atomistic calculations and by the progress made in the development of first-principles methods for excited electronic states of large molecules. In this Review, we discuss the advantages and pitfalls of theoretical methods applicable to photosynthetic pigments. Besides methodological aspects of excited-state electronic-structure methods, studies on chlorophyll-type and carotenoid-like molecules are discussed. We also address the concepts of exciton coupling and excitation-energy transfer (EET) and compare the different theoretical methods for the calculation of EET coupling constants. Applications to photosynthetic light-harvesting complexes and reaction centers based on such models are also analyzed.

Studies of the ground and excited-state surfaces of the retinal chromophore using CAM-B3LYP

The isomerization of the 11-cis isomer (PSB11) of the retinal chromophore to its all-trans isomer (PSBT) is examined. Optimized structures on both the ground state and the excited state are calculated, and the dependence on torsional angles in the carbon chain is investigated. Time-dependent density functional theory is used to produce excitation energies and the excited-state surface. To avoid problems with the description of excited states that can arise with standard DFT methods, the CAM-B3LYP functional was used. Comparing CAM-B3LYP with B3LYP results indicates that the former is significantly more accurate, as a consequence of which detailed cross sections of the retinal excited-state surface are obtained.

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

Optical line shape theory is combined with a quantum-chemical/electrostatic calculation of the site energies of the 96 chlorophyll a pigments and their excitonic couplings to simulate optical spectra of photosystem I core complexes from Thermosynechococcus elongatus. The absorbance, linear dichroism and circular dichroism spectra, calculated on the basis of the 2.5 A crystal structure, match the experimental data semiquantitatively allowing for a detailed analysis of the pigment-protein interaction. The majority of site energies are determined by multiple interactions with a large number (>20) of amino acid residues, a result which demonstrates the importance of long-range electrostatic interactions. The low-energy exciton states of the antenna are found to be located at a nearest distance of about 25 A from the special pair of the reaction center. The intermediate pigments form a high-energy bridge, the site energies of which are stabilized by a particularly large number (>100) of amino acid residues. The concentration of low energy exciton states in the antenna is larger on the side of the A-branch of the reaction center, implying an asymmetric delivery of excitation energy to the latter. This asymmetry in light-harvesting may provide the key for understanding the asymmetric use of the two branches in primary electron transfer reactions. Experiments are suggested to check for this possibility.