Exciton properties and optical spectra of light harvesting complex II from a fully atomistic description

Combined Multireference-Multiscale Approach to the Description of Photosynthetic Reaction Centers

A first-principles description of the primary photochemical processes that drive photosynthesis and sustain life on our planet remains one of the grand challenges of modern science. Recent research established that explicit incorporation of protein electrostatics in excited-state calculations of photosynthetic pigments, achieved for example with quantum-mechanics/molecular-mechanics (QM/MM) approaches, is essential for a meaningful description of the properties and function of pigment-protein complexes. Although time-dependent density functional theory has been used productively so far in QM/MM approaches for the study of such systems, this methodology has limitations. Here we pursue for the first time a QM/MM description of the reaction center in the principal enzyme of oxygenic photosynthesis, Photosystem II, using multireference wave function theory for the high-level QM region. We identify best practices and establish guidelines regarding the rational choice of active space and appropriate state-averaging for the efficient and reliable use of complete active space self-consistent field (CASSCF) and the N-electron valence state perturbation theory (NEVPT2) in the prediction of low-lying excited states of chlorophyll and pheophytin pigments. Given that the Gouterman orbitals are inadequate as a minimal active space, we define specific minimal and extended active spaces for the NEVPT2 description of electronic states that fall within the Q and B bands. Subsequently, we apply our multireference-QM/MM protocol to the description of all pigments in the reaction center of Photosystem II. The calculations reproduce the electrochromic shifts induced by the protein matrix and the ordering of site energies consistent with the identity of the primary donor (ChlD1) and the experimentally known asymmetric and directional electron transfer. The optimized protocol sets the stage for future multireference treatments of multiple pigments, and hence for multireference studies of charge separation, while it is transferable to the study of any photoactive embedded tetrapyrrole system.

Insights into Energy Transfer in Light-Harvesting Complex II Through Machine-Learning Assisted Simulations

Light-harvesting complex II (LHCII) is the major antenna of higher plants. Energy transfer processes taking place inside its aggregate of chlorophylls have been experimentally investigated with time-resolved techniques, but a complete understanding of the most relevant energy transfer pathways and relative characteristic times remains elusive. Theoretical models to disentangle experimental data in LHCII have long been challenged by the large size and complex nature of the system. Here, we show that a fully first-principles approach combining molecular dynamics and machine learning can be successfully used to reproduce transient absorption spectra and characterize the EET pathways and the involved times.

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.

First-principles simulation of excitation energy transfer and transient absorption spectroscopy in the CP29 light-harvesting complex

The dynamics of delocalized excitons in light-harvesting complexes (LHCs) can be investigated using different experimental techniques, and transient absorption (TA) spectroscopy is one of the most valuable methods for this purpose. A careful interpretation of TA spectra is essential for the clarification of excitation energy transfer (EET) processes occurring during light-harvesting. However, even in the simplest LHCs, a physical model is needed to interpret transient spectra as the number of EET processes occurring at the same time is very large to be disentangled from measurements alone. Physical EET models are commonly built by fittings of the microscopic exciton Hamiltonians and exciton-vibrational parameters, an approach that can lead to biases. Here, we present a first-principles strategy to simulate EET and transient absorption spectra in LHCs, combining molecular dynamics and accurate multiscale quantum chemical calculations to obtain an independent estimate of the excitonic structure of the complex. The microscopic parameters thus obtained are then used in EET simulations to obtain the population dynamics and the related spectroscopic signature. We apply this approach to the CP29 minor antenna complex of plants for which we follow the EET dynamics and transient spectra after excitation in the chlorophyll b region. Our calculations reproduce all the main features observed in the transient absorption spectra and provide independent insight on the excited-state dynamics of CP29. The approach presented here lays the groundwork for the accurate simulation of EET and unbiased interpretation of transient spectra in multichromophoric systems.

Unraveling the electronic origin of a special feature in the triplet-minus-singlet spectra of carotenoids in natural photosystems

The influence of carotenoid triplet states on the Qy electronic transitions of chlorophylls has been observed in experiments on light-harvesting complexes over the past three decades, but the interpretation of the resulting spectral feature in the triplet minus singlet (T-S) absorption spectra of photosystems is still debated, as the physical-chemical explanation of this feature has been elusive. Here, we resolve this debate, by explaining the T-S spectra of pigment complexes over the Qy-band spectral region through a comparative study of chlorophyll-carotenoid model dyads and larger pigment complexes from the main light harvesting complex of higher plants (LHCII). This goal is achieved by combining state-of-the-art time-dependent density functional theory with analysis of the relationship between electronic properties and nuclear structure, and by comparison to the experiment. We find that the special signature in the T-S spectra of both model and natural photosystems is determined by singlet-like triplet excitations that can be described as effective singlet excitations on chlorophylls influenced by a stable electronic triplet on the carotenoid. The comparison with earlier experiments on different light-harvesting complexes confirms our theoretical interpretation of the T-S spectra in the Qy spectral region. Our results indicate an important role for the chlorophyll-carotenoid electronic coupling, which is also responsible for the fast triplet-triplet energy transfer, suggesting a fast trapping of the triplet into the relaxed carotenoid structure. The gained understanding of the interplay between the electronic and nuclear structures is potentially informative for future studies of the mechanism of photoprotection by carotenoids.

Impact of Structure, Coupling Scheme, and State of Interest on the Energy Transfer in CP29

The Qy and Bx excitation energy transfer (EET) in the minor light-harvesting complex CP29 (LHCII B4.1) antenna complex of Pisum sativum was characterized using a computational approach. We applied Förster resonance energy transfer (FRET) and the transition density cube (TDC) method to estimate the Coulombic coupling, based on a combination of classical molecular dynamics and quantum mechanics/molecular mechanics calculations. Employing TDC instead of FRET mostly affects the EET between chlorophylls (Chls) and carotenoids (Crts), as expected due to the Crts being spatially more challenging for FRET. Only between Chls, effects are found to be small (about only 0.1 EET efficiency change when introducing TDC instead of FRET). Effects of structural sampling were found to be small, illustrated by a small average standard deviation for the Qy state coupling elements (FRET/TDC: 0.97/0.94 cm-1). Due to the higher flexibility of the Bx state, the corresponding deviations are larger (FRET/TDC between Chl-Chl pairs: 17.58/22.67 cm-1, between Crt-Chl pairs: 62.58/31.63 cm-1). In summary, it was found for the Q band that the coupling between Chls varies only slightly depending on FRET or TDC, resulting in a minute effect on EET acceptor preference. In contrast, the coupling in the B band spectral region is found to be more affected. Here, the S2 (1Bu) states of the spatially challenging Crts may act as acceptors in addition to the B states of the Chls. Depending on FRET or TDC, several Chls show different Chl-to-Crt couplings. Interestingly, the EET between Chls or Crts in the B band is found to often outcompete the corresponding decay processes. The individual efficiencies for B band EET to Crts vary however strongly with the chosen coupling scheme (e.g., up to 0.29/0.99 FRET/TDC efficiency for the Chl a604/neoxanthin pair). Thus, the choice of the coupling scheme must involve a consideration of the state of interest.

Tidying up the conformational ensemble of a disordered peptide by computational prediction of spectroscopic fingerprints

The most advanced structure prediction methods are powerless in exploring the conformational ensemble of disordered peptides and proteins and for this reason the "protein folding problem" remains unsolved. We present a novel methodology that enables the accurate prediction of spectroscopic fingerprints (circular dichroism, infrared, Raman, and Raman optical activity), and by this allows for "tidying up" the conformational ensembles of disordered peptides and disordered regions in proteins. This concept is elaborated for and applied to a dodecapeptide, whose spectroscopic fingerprint is measured and theoretically predicted by means of enhanced-sampling molecular dynamics coupled with quantum mechanical calculations. Following this approach, we demonstrate that peptides lacking a clear propensity for ordered secondary-structure motifs are not randomly, but only conditionally disordered. This means that their conformational landscape, or phase-space, can be well represented by a basis-set of conformers including about 10 to 100 structures. The implications of this finding have profound consequences both for the interpretation of experimental electronic and vibrational spectral features of peptides in solution and for the theoretical prediction of these features using accurate and computationally expensive techniques. The here-derived methods and conclusions are expected to fundamentally impact the rationalization of so-far elusive structure-spectra relationships for disordered peptides and proteins, towards improved and versatile structure prediction methods.

Quantum chemical elucidation of a sevenfold symmetric bacterial antenna complex

The light-harvesting complex 2 (LH2) of purple bacteria is one of the most studied photosynthetic antenna complexes. Its symmetric structure and ring-like bacteriochlorophyll arrangement make it an ideal system for theoreticians and spectroscopists. LH2 complexes from most bacterial species are thought to have eightfold or ninefold symmetry, but recently a sevenfold symmetric LH2 structure from the bacterium Mch. purpuratum was solved by Cryo-Electron microscopy. This LH2 also possesses unique near-infrared absorption and circular dichroism (CD) spectral properties. Here we use an atomistic strategy to elucidate the spectral properties of Mch. purpuratum LH2 and understand the differences with the most commonly studied LH2 from Rbl. acidophilus. Our strategy exploits a combination of molecular dynamics simulations, multiscale polarizable quantum mechanics/molecular mechanics calculations, and lineshape simulations. Our calculations reveal that the spectral properties of LH2 complexes are tuned by site energies and exciton couplings, which in turn depend on the structural fluctuations of the bacteriochlorophylls. Our strategy proves effective in reproducing the absorption and CD spectra of the two LH2 complexes, and in uncovering the origin of their differences. This work proves that it is possible to obtain insight into the spectral tuning strategies of purple bacteria by quantitatively simulating the spectral properties of their antenna complexes.

Towards the description of charge transfer states in solubilised LHCII using subsystem DFT

Light harvesting complex II (LHCII) in plants and green algae have been shown to adapt their absorption properties, depending on the concentration of sunlight, switching between a light harvesting and a non-harvesting or quenched state. In a recent work, combining classical molecular dynamics (MD) simulations with quantum chemical calculations (Liguori et al. in Sci Rep 5:15661, 2015) on LHCII, it was shown that the Chl611-Chl612 cluster of the terminal emitter domain can play an important role in modifying the spectral properties of the complex. In that work the importance of charge transfer (CT) effects was highlighted, in re-shaping the absorption intensity of the chlorophyll dimer. Here in this work, we investigate the combined effect of the local excited (LE) and CT states in shaping the energy landscape of the chlorophyll dimer. Using subsystem Density Functional Theory over the classical [Formula: see text]s MD trajectory we look explicitly into the excitation energies of the LE and the CT states of the dimer and their corresponding couplings. Upon doing so, we observe a drop in the excitation energies of the CT states, accompanied by an increase in the couplings between the LE/LE and the LE/CT states facilitated by a shorter interchromophoric distance upon equilibration. Both these changes in conjunction, effectively produces a red-shift of the low-lying mixed exciton/CT states of the supramolecular chromophore pair.

Recent progress in atomistic modeling of light-harvesting complexes: a mini review

In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.

The origin of pigment-binding differences in CP29 and LHCII: the role of protein structure and dynamics

The first step of photosynthesis in plants is performed by the light-harvesting complexes (LHC), a large family of pigment-binding proteins embedded in the photosynthetic membranes. These complexes are conserved across species, suggesting that each has a distinct role. However, they display a high degree of sequence homology and their static structures are almost identical. What are then the structural features that determine their different properties? In this work, we compared the two best-characterized LHCs of plants: LHCII and CP29. Using molecular dynamics simulations, we could rationalize the difference between them in terms of pigment-binding properties. The data also show that while the loops between the helices are very flexible, the structure of the transmembrane regions remains very similar in the crystal and the membranes. However, the small structural differences significantly affect the excitonic coupling between some pigment pairs. Finally, we analyzed in detail the structure of the long N-terminus of CP29, showing that it is structurally stable and it remains on top of the membrane even in the absence of other proteins. Although the structural changes upon phosphorylation are minor, they can explain the differences in the absorption properties of the pigments observed experimentally.

Machine Learning Exciton Hamiltonians in Light-Harvesting Complexes

We propose a machine learning (ML)-based strategy for an inexpensive calculation of excitonic properties of light-harvesting complexes (LHCs). The strategy uses classical molecular dynamics simulations of LHCs in their natural environment in combination with ML prediction of the excitonic Hamiltonian of the embedded aggregate of pigments. The proposed ML model can reproduce the effects of geometrical fluctuations together with those due to electrostatic and polarization interactions between the pigments and the protein. The training is performed on the chlorophylls of the major LHC of plants, but we demonstrate that the model is able to extrapolate well beyond the initial training set. Moreover, the accuracy in predicting the effects of the environment is tested on the simulation of the small changes observed in the absorption spectra of the wild-type and a mutant of a minor LHC.

Structural insights into blue-green light utilization by marine green algal light harvesting complex II at 2.78 Å

Light-harvesting complex II (LHCII) present in plants and green algae absorbs solar energy to promote photochemical reactions. A marine green macroalga, Codium fragile, exhibits the unique characteristic of absorbing blue-green light from the sun during photochemical reactions while being underwater owing to the presence of pigment-altered LHCII called siphonaxanthin-chlorophyll a/b-binding protein (SCP). In this study, we determined the structure of SCP at a resolution of 2.78 Å using cryogenic electron microscopy. SCP has a trimeric structure, wherein each monomer containing two lutein and two chlorophyll a molecules in the plant-type LHCII are replaced by siphonaxanthin and its ester and two chlorophyll b molecules, respectively. Siphonaxanthin occupies the binding site in SCP having a polarity in the trimeric inner core, and exhibits a distorted conjugated chain comprising a carbonyl group hydrogen bonded to a cysteine residue of apoprotein. These features suggest that the siphonaxanthin molecule is responsible for the characteristic green absorption of SCP. The replaced chlorophyll b molecules extend the region of the stromal side chlorophyll b cluster, spanning two adjacent monomers.

Bethe-Salpeter equation insights into the photo-absorption function and exciton structure of chlorophyll a and b in light-harvesting complex II

The photo-absorption process and the excitation of chlorophyll (Chl) is the primary and essential step of photosynthesis in green plants. By solving the Bethe-Salpeter equation (BSE) on top of the GW approximation within ab initio many-body perturbation theory, we calculate the photo-absorption function and the excitons structure of Chl a and b in their in vivo conformations as measured by X-ray diffraction in the light-harvesting complex (LHC) II. BSE optical absorption spectra are in good agreement with the experiment and we discuss residual discrepancies. The experimental evidence of multiple Chla forms in vivo is explained by BSE. The Chla and Chlb BSE exciton wavefunctions present important charge-transfer differences on the Soret band. Q excitons are almost identical, apart from charge (both electron and hole) localization on the Chlb C7 aldheide formyl group, absent on the Chla methyl C7, that is exactly the group where the two chlorophylls differ.

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.

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.

A fast method for electronic couplings in embedded multichromophoric systems

Electronic couplings are key to understanding exciton delocalization and transport in natural and artificial light harvesting processes. We develop a method to compute couplings in multichromophoric aggregates embedded in complex environments without running expensive quantum chemical calculations. We use a transition charge approximation to represent the quantum mechanical transition densities of the chromophores and an atomistic and polarizable classical model to describe the environment atoms. We extend our framework to estimate transition charges directly from the chromophore geometry, i.e., bypassing completely the quantum mechanical calculations using a regression approach. The method allows to rapidly compute accurate couplings for a large number of geometries along molecular dynamics trajectories.

Violaxanthin and Zeaxanthin May Replace Lutein at the L1 Site of LHCII, Conserving the Interactions with Surrounding Chlorophylls and the Capability of Triplet-Triplet Energy Transfer

Carotenoids represent the first line of defence of photosystems against singlet oxygen (¹O₂) toxicity, because of their capacity to quench the chlorophyll triplet state (³Chl) through a physical mechanism based on the transfer of triplet excitation (triplet-triplet energy transfer, TTET). In previous works, we showed that the antenna LHCII is characterised by a robust photoprotective mechanism, able to adapt to the removal of individual chlorophylls while maintaining a remarkable capacity for ³Chl quenching. In this work, we investigated the effects on this quenching induced in LHCII by the replacement of the lutein bound at the L1 site with violaxanthin and zeaxanthin. We studied LHCII isolated from the Arabidopsis thaliana mutants lut2-in which lutein is replaced by violaxanthin-and lut2 npq2, in which all xanthophylls are replaced constitutively by zeaxanthin. We characterised the photophysics of these systems via optically detected magnetic resonance (ODMR) and time-resolved electron paramagnetic resonance (TR-EPR). We concluded that, in LHCII, lutein-binding sites have conserved characteristics, and ensure efficient TTET regardless of the identity of the carotenoid accommodated.

The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol

The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.

Structure of the stress-related LHCSR1 complex determined by an integrated computational strategy

Light-harvesting complexes (LHCs) are pigment-protein complexes whose main function is to capture sunlight and transfer the energy to reaction centers of photosystems. In response to varying light conditions, LH complexes also play photoregulation and photoprotection roles. In algae and mosses, a sub-family of LHCs, light-harvesting complex stress-related (LHCSR), is responsible for photoprotective quenching. Despite their functional and evolutionary importance, no direct structural information on LHCSRs is available that can explain their unique properties. In this work, we propose a structural model of LHCSR1 from the moss P. patens, obtained through an integrated computational strategy that combines homology modeling, molecular dynamics, and multiscale quantum chemical calculations. The model is validated by reproducing the spectral properties of LHCSR1. Our model reveals the structural specificity of LHCSR1, as compared with the CP29 LH complex, and poses the basis for understanding photoprotective quenching in mosses.

The structure of the moss P. patens light-harvesting complex stress-related 1 (LHCSR1) is determined using a multi-scale computational approach for investigations of its photoprotective function.

Effect of protein matrix on CP29 spectra and energy transfer pathways

We investigate energy transfer pathways between strongly coupled chlorophylls (Chls) in the CP29 (LHCII B4.1) antenna complex of Pisum sativum, including the possibility of higher energy states. We test for the environmental effects caused by the protein, membrane and solvent using a hybrid QM/MM approach. Classical molecular dynamics simulations of the full CP29 complex embedded in a DOPC membrane have been performed, followed by calculations of the time dependent DFT spectra of all Chls at several timesteps. The relative orientations of transition dipole moments (TDMs) were specifically analyzed, including and excluding the point charge field (PCF) of the surrounding environment. The PCF is found to drastically shift the spectra of specific Chls, while the majority of Chls is mostly unaffected. The net effect on the sum spectrum is however found to be negligible: The few strong changes in Chl spectra cancel each other due to being opposite in sign. We further find that the spectra of the Chls coordinating to water show a blue shift upon introduction of the environment. Conversely, the spectra of the Chls coordinating to glutamine show a red shift upon activation of the PCF. As the main influence of the PCF for tuning the couplings, we identify the energetic position of the individual chromophores. The fine-tuning, especially for states energetically above the Qy state, is however controlled by the changes in the TDM orientations. We also find an indication for the PCF to steer potentially harmful high energy excitations away from the PSII core complex.

Altering the exciton landscape by removal of specific chlorophylls in monomeric LHCII provides information on the sites of triplet formation and quenching by means of ODMR and EPR spectroscopies

The triplet states populated under illumination in the monomeric light-harvesting complex II (LHCII) were analyzed by EPR and Optically Detected Magnetic Resonance (ODMR) in order to fully characterize the perturbations introduced by site-directed mutations leading to the removal of key chlorophylls. We considered the A2 and A5 mutants, lacking Chls a612(a611) and Chl a603 respectively, since these Chls have been proposed as the sites of formation of triplet states which are subsequently quenched by the luteins. Chls a612 and Chl a603 belong to the two clusters determining the low energy exciton states in the complex. Their removal is expected to significantly alter the excitation energy transfer pathways. On the basis of the TR- and pulse EPR triplet spectra, the two symmetrically related pairs constituted by Chl a612/Lut620 and Chl a603/Lut621 were both possible candidate for triplet-triplet energy transfer (TTET). However, the ODMR results clearly show that only Lut620 is involved in triplet quenching. In the A5 mutant, the Chl a612/Lut620 pair retains this pivotal photoprotective role, while the A2 mutant was found to activate an alternative pathway involving the Chl a603/Lut621pair. These results shows that LHCII is characterized by a robust photoprotective mechanism, able to adapt to the removal of individual chromophores while maintaining a remarkable degree of Chl triplet quenching. Small amounts of unquenched Chl triplet states were also detected. The analysis of the results allowed us to assign the sites of "unquenched" chlorophyll triplets to Chl a610 and Chl a602.

Photosynthetic Light-Harvesting (Antenna) Complexes-Structures and Functions

Chlorophylls and bacteriochlorophylls, together with carotenoids, serve, noncovalently bound to specific apoproteins, as principal light-harvesting and energy-transforming pigments in photosynthetic organisms. In recent years, enormous progress has been achieved in the elucidation of structures and functions of light-harvesting (antenna) complexes, photosynthetic reaction centers and even entire photosystems. It is becoming increasingly clear that light-harvesting complexes not only serve to enlarge the absorption cross sections of the respective reaction centers but are vitally important in short- and long-term adaptation of the photosynthetic apparatus and regulation of the energy-transforming processes in response to external and internal conditions. Thus, the wide variety of structural diversity in photosynthetic antenna “designs” becomes conceivable. It is, however, common for LHCs to form trimeric (or multiples thereof) structures. We propose a simple, tentative explanation of the trimer issue, based on the 2D world created by photosynthetic membrane systems.

Multiscale QM/MM molecular dynamics simulations of the trimeric major light-harvesting complex II

Photosynthetic processes are driven by sunlight. Too little of it and the photosynthetic machinery cannot produce the reductive power to drive the anabolic pathways. Too much sunlight and the machinery can get damaged. In higher plants, the major Light-Harvesting Complex (LHCII) efficiently absorbs the light energy, but can also dissipate it when in excess (quenching). In order to study the dynamics related to the quenching process but also the exciton dynamics in general, one needs to accurately determine the so-called spectral density which describes the coupling between the relevant pigment modes and the environmental degrees of freedom. To this end, Born-Oppenheimer molecular dynamics simulations in a quantum mechanics/molecular mechanics (QM/MM) fashion utilizing the density functional based tight binding (DFTB) method have been performed for the ground state dynamics. Subsequently, the time-dependent extension of the long-range-corrected DFTB scheme has been employed for the excited state calculations of the individual chlorophyll-a molecules in the LHCII complex. The analysis of this data resulted in spectral densities showing an astonishing agreement with the experimental counterpart in this rather large system. This consistency with an experimental observable also supports the accuracy, robustness, and reliability of the present multi-scale scheme. To the best of our knowledge, this is the first theoretical attempt on this large complex system is ever made to accurately simulate the spectral density. In addition, the resulting spectral densities and site energies were used to determine the exciton transfer rate within a special pigment pair consisting of a chlorophyll-a and a carotenoid molecule which is assumed to play a role in the balance between the light harvesting and quenching modes.

Confronting FCP structure with ultrafast spectroscopy data: evidence for structural variations

Diatoms are a major group of algae, responsible for a quarter of the global primary production on our planet. Their adaptation to marine environments is ensured by their light-harvesting antenna - the fucoxanthin-chlorophyll protein (FCP) complex, which absorbs strongly in the blue-green spectral region. Although these essential proteins have been the subject of many studies, for a long time their comprehensive description was not possible in the absence of structural data. Last year, the 3D structures of several FCP complexes were revealed. The structure of an FCP dimer was resolved by crystallography for the pennate diatom Phaeodactylum tricornutum [W. Wang et al., Science, 2019, 363, 6427] and the structure of the PSII supercomplex from the centric diatom Chaetoceros gracilis, containing several FCPs, was obtained by electron microscopy [X. Pi et al., Science, 2019, 365, 6452; R. Nagao et al., Nat. Plants, 2019, 5, 890]. In this Perspective article, we evaluate how precisely these structures may account for previously published ultrafast spectroscopy results, describing the excitation energy transfer in the FCP from another centric diatom Cyclotella meneghiniana. Surprisingly, we find that the published FCP structures cannot explain several observations obtained from ultrafast spectroscopy. Using the available structures, and results from electron microscopy, we construct a trimer-based FCP model for Cyclotella meneghiniana, consistent with ultrafast experimental data. As a whole, our observations suggest that the structures from the proteins belonging to the FCP family display larger variations than the equivalent LHC proteins in plants, which may reflect species-specific adaptations or original strategies for adapting to rapidly changing marine environments.

Static Disorder in Excitation Energies of the Fenna-Matthews-Olson Protein: Structure-Based Theory Meets Experiment

Inhomogeneous broadening of optical lines of the Fenna-Matthews-Olson (FMO) light-harvesting protein is investigated by combining a Monte Carlo sampling of low-energy conformational substates of the protein with a quantum chemical/electrostatic calculation of local transition energies (site energies) of the pigments. The good agreement between the optical spectra calculated for the inhomogeneous ensemble and the experimental data demonstrates that electrostatics is the dominant contributor to static disorder in site energies. Rotamers of polar amino acid side chains are found to cause bimodal distribution functions of site energy shifts, which can be probed by hole burning and single-molecule spectroscopy. When summing over the large number of contributions, the resulting distribution functions of the site energies become Gaussians, and the correlations in site energy fluctuations at different sites practically average to zero. These results demonstrate that static disorder in the FMO protein is in the realm of the central limit theorem of statistics.

Multiscale modeling of molecular structure and optical properties of complex supramolecular aggregates

Supramolecular aggregates of synthetic dye molecules offer great perspectives to prepare biomimetic functional materials for light-harvesting and energy transport. The design is complicated by the fact that structure–property relationships are hard to establish, because the molecular packing results from a delicate balance of interactions and the excitonic properties that dictate the optics and excited state dynamics, in turn sensitively depend on this packing. Here we show how an iterative multiscale approach combining molecular dynamics and quantum mechanical exciton modeling can be used to obtain accurate insight into the packing of thousands of cyanine dye molecules in a complex double-walled tubular aggregate in close interaction with its solvent environment. Our approach allows us to answer open questions not only on the structure of these prototypical aggregates, but also about their molecular-scale structural and energetic heterogeneity, as well as on the microscopic origin of their photophysical properties. This opens the route to accurate predictions of energy transport and other functional properties.

Multiscale modeling resolves the molecular structure of a synthetic light-harvesting complex, unraveling the microscopic origin of its photophysical properties.