The role of charge-transfer states in the spectral tuning of antenna complexes of purple bacteria

Protein-Driven Electron-Transfer Process in a Fatty Acid Photodecarboxylase

Naturally occurring photoenzymes are rare in nature, but among them, fatty acid photodecarboxylases derived from Chlorella variabilis (CvFAPs) have emerged as promising photobiocatalysts capable of performing the redox-neutral, light-induced decarboxylation of free fatty acids (FAs) into C1-shortened n-alka(e)nes. Using a hybrid QM/MM approach combined with a polarizable embedding scheme, we identify the structural changes of the active site and determine the energetic landscape of the forward electron transfer (fET) from the FA substrate to the excited flavin adenine dinucleotide. We obtain a charge-transfer diradical structure where a water molecule rearranges spontaneously to form a H-bond interaction with the excited flavin, while the FA's carboxylate group twists and migrates away from it. Together, these structural modifications provide the driving force necessary for the fET to proceed in a downhill direction. Moreover, by examining the R451K mutant where the FA substrate is farther from the flavin core, we show that the marked reduction of the electronic coupling is counterbalanced by an increased driving force, resulting in a fET lifetime similar to the WT, thereby suggesting a resilience of the process to this mutation. Finally, through QM/MM molecular dynamic simulations, we reveal that, following fET, the decarboxylation of the FA radical occurs within tens of picoseconds, overcoming an energy barrier of ∼0.1 eV. Overall, by providing an atomistic characterization of the photoactivation of CvFAP, this work can be used for future protein engineering.

Spectral modulation of B850 bacteriochlorophyll a in light-harvesting complex 2 from purple photosynthetic bacterium Thermochromatium tepidum by detergents and calcium ions

Spectral variations of light-harvesting (LH) proteins of purple photosynthetic bacteria provide insight into the molecular mechanisms underlying spectral tuning of circular bacteriochlorophyll (BChl) arrays, which play crucial roles in photoenergy conversion in these organisms. Here we investigate spectral changes of the Qy band of B850 BChl a in LH2 protein from purple sulfur bacterium Thermochromatium tepidum (tepidum-LH2) by detergents and Ca2+. The tepidum-LH2 solubilized with lauryl dimethylamine N-oxide and n-octyl-β-D-glucoside (LH2LDAO and LH2OG, respectively) exhibited blue-shift of the B850 Qy band with hypochromism compared with the tepidum-LH2 solubilized with n-dodecyl-β-D-maltoside (LH2DDM), resulting in the LH3-like spectral features. Resonance Raman spectroscopy indicated that this blue-shift was ascribable to the loss of hydrogen-bonding between the C3-acetyl group in B850 BChl a and the LH2 polypeptides. Ca2+ produced red-shift of the B850 Qy band in LH2LDAO by forming hydrogen-bond for the C3-acetyl group in B850 BChl a, probably due to a change in the microenvironmental structure around B850. Ca2+-induced red-shift was also observed in LH2OG although the B850 acetyl group is still free from hydrogen-bonding. Therefore, the Ca2+-induced B850 red-shift in LH2OG would originate from an electrostatic effect of Ca2+. The current results suggest that the B850 Qy band in tepidum-LH2 is primarily tuned by two mechanisms, namely the hydrogen-bonding of the B850 acetyl group and the electrostatic effect.

Excitonic Approach for Nonadiabatic Dynamics: Extending Beyond the Frenkel Exciton Model

We report the formulation and implementation of an extended Frenkel exciton model (EFEM) designed for simulating the dynamics of multichromophoric systems, taking into account the possible presence of interchromophore charge transfer states, as well as other states in which two chromophores are simultaneously excited. Our approach involves constructing a Hamiltonian based on calculations performed on monomers and selected dimers within the multichromophoric aggregate. Nonadiabatic molecular dynamics is addressed using a surface hopping approach, while the electronic wave functions and energies required for constructing the EFEM are computed utilizing the semiempirical floating occupation molecular orbitals-configuration interaction (FOMO-CI) electronic structure method. To validate our approach, we simulate the singlet fission process in a trimer of 2,5-bis(fluorene-9-ylidene)-2,5-dihydrothiophene (ThBF) molecules, embedded in their crystal environment, comparing the results of the EFEM to the standard "supermolecule" approach.

Prominent Role of Charge Transfer in the Spectral Tuning of Photosynthetic Light-Harvesting I Complex

Purple bacteria possess two ring-shaped protein complexes, light-harvesting 1 (LH1) and 2 (LH2), both of which function as antennas for solar energy utilization for photosynthesis but exhibit distinct absorption properties. The two antennas have differing amounts of bacteriochlorophyll (BChl) a; however, their significance in spectral tuning remains elusive. Here, we report a high-precision evaluation of the physicochemical factors contributing to the variation in absorption maxima between LH1 and LH2, namely, BChl a structural distortion, protein electrostatic interaction, excitonic coupling, and charge transfer (CT) effects, as derived from detailed spectral calculations using an extended version of the exciton model, in the model purple bacterium Rhodospirillum rubrum. Spectral analysis confirmed that the electronic structure of the excited state in LH1 extended to the BChl a 16-mer. Further analysis revealed that the LH1-specific redshift (∼61% in energy) is predominantly accounted for by the CT effect resulting from the closer inter-BChl distance in LH1 than in LH2. Our analysis explains how LH1 and LH2, both with chemically identical BChl a chromophores, use distinct physicochemical effects to achieve a progressive redshift from LH2 to LH1, ensuring efficient energy transfer to the reaction center special pair.

Cavity-Mediated Enhancement of the Energy Transfer in the Reduced Fenna-Matthews-Olson Complex

Strong light-matter interaction leads to the formation of hybrid polariton states and can alter the light-harvesting properties of natural photosynthetic systems without modifying their chemical structure. In the present study, we computationally investigate the effect of the resonant cavity on the efficiency and the rate of the population transfer in a quantum system coupled to the cavity and the dissipative environment. The parameters of the model system were chosen to represent the Fenna-Matthews-Olson natural light-harvesting complex reduced to the three essential sites. The dynamics of the total system was propagated using the hierarchical equations of motion. Our results show that the strong light-matter interaction can accelerate the population transfer process compared to the cavity-free case but at the cost of lowering the transfer efficiency. The transition to the strong coupling regime was found to coincide with the degeneracy of polariton eigenvalues. Our findings indicate the potential and the limit of tuning the energy transfer in already efficient natural light-harvesting systems.

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.

Environment-dependent chlorophyll-chlorophyll charge transfer states in Lhca4 pigment-protein complex

Photosystem I (PSI) light-harvesting antenna complexes LHCI contain spectral forms that absorb and emit photons of lower energy than that of its primary electron donor, P700. The most red-shifted fluorescence is associated with the Lhca4 complex. It has been suggested that this red emission is related to the inter-chlorophyll charge transfer (CT) states. In this work we present a systematic quantum-chemical study of the CT states in Lhca4, accounting for the influence of the protein environment by estimating the electrostatic interactions. We show that significant energy shifts result from these interactions and propose that the emission of the Lhca4 complex is related not only to the previously proposed a603+-a608- state, but also to the a602+-a603- state. We also investigate how different protonation patterns of protein amino acids affect the energetics of the CT states.

Theoretical description of photoinduced electron transfer in donor-acceptor supramolecular complexes based on carbon buckybowls

Herein, we explore, from a theoretical perspective, the nonradiative photoinduced processes (charge separation and energy transfer) within a family of donor-acceptor supramolecular complexes based on the electron-donor truxene-tetrathiafulvalene (truxTTF) derivative and a series of curved fullerene fragments (buckybowls) of different shapes and sizes (C30H12, C32H12, and C38H14) as electron acceptors that successfully combine with truxTTF via non-covalent interactions. The resulting supramolecular complexes (truxTTF·C30H12, truxTTF·C32H12, and truxTTF·C38H14) undergo charge-separation processes upon photoexcitation through charge-transfer states involving the donor and acceptor units. Despite the not so different size of the buckybowls, they present noticeable differences in the charge-separation efficiency owing to a complex decay post-photoexcitation mechanism involving several low-lying excited states of different natures (local and charge-transfer excitations), all closely spaced in energy. In this intricate scenario, we have adopted a theoretical approach combining electronic structure calculations at (time-dependent) density functional theory, a multistate multifragment diabatization method, the Marcus-Levitch-Jortner semiclassical rate expression, and a kinetic model to estimate the charge separation rate constants of the supramolecular heterodimers. Our outcomes highlight that the efficiency of the photoinduced charge-separation process increases with the extension of the buckybowl backbone. The supramolecular heterodimer with the largest buckybowl (truxTTF·C38H14) displays multiple and efficient electron-transfer pathways, providing a global photoinduced charge separation in the ultrafast time scale in line with the experimental findings. The study reported indicates that modifications in the shape and size of buckybowl systems can give rise to attractive novel acceptors for potential photovoltaic applications.

Excitation landscape of the CP43 photosynthetic antenna complex from multiscale simulations

Photosystem II (PSII), the principal enzyme of oxygenic photosynthesis, contains two integral light harvesting proteins (CP43 and CP47) that bind chlorophylls and carotenoids. The two intrinsic antennae play crucial roles in excitation energy transfer and photoprotection. CP43 interacts most closely with the reaction center of PSII, specifically with the branch of the reaction center (D1) that is responsible for primary charge separation and electron transfer. Deciphering the function of CP43 requires detailed atomic-level insights into the properties of the embedded pigments. To advance this goal, we employ a range of multiscale computational approaches to determine the site energies and excitonic profile of CP43 chlorophylls, using large all-atom models of a membrane-bound PSII monomer. In addition to time-dependent density functional theory (TD-DFT) used in the context of a quantum-mechanics/molecular-mechanics setup (QM/MM), we present a thorough analysis using the perturbed matrix method (PMM), which enables us to utilize information from long-timescale molecular dynamics simulations of native PSII-complexed CP43. The excited state energetics and excitonic couplings have both similarities and differences compared with previous experimental fits and theoretical calculations. Both static TD-DFT and dynamic PMM results indicate a layered distribution of site energies and reveal specific groups of chlorophylls that have shared contributions to low-energy excitations. Importantly, the contribution to the lowest energy exciton does not arise from the same chlorophylls at each system configuration, but rather changes as a function of conformational dynamics. An unexpected finding is the identification of a low-energy charge-transfer excited state within CP43 that involves a lumenal (C2) and the central (C10) chlorophyll of the complex. The results provide a refined basis for structure-based interpretation of spectroscopic observations and for further deciphering excitation energy transfer in oxygenic photosynthesis.

Eclipsed and Twisted Excimers of Pyrene and 2-Azapyrene: How Nitrogen Substitution Impacts Excimer Emission

Due to their unique photophysical and electronic properties, pyrene and its analogues have been the subject of extensive research in recent decades. The propensity of pyrene and its derivatives to form excimers has found wide application in various fields. Nitrogen-substituted pyrene derivatives display similar photophysical properties, but for them, excimer emission has not been reported to date. Here, we use time-dependent density functional theory (TD-DFT) calculations to investigate the low-lying exciton states of dimers of pyrene and 2-azapyrene. The excimer equilibrium structures are determined and the contribution of charge transfer (CT) excitations and intermolecular interactions to the exciton states is disclosed using a diabatization procedure. The study reveals that the dimers formed by the two molecules have quite similar exciton-state patterns, in which the relevant CT contributions govern the formation of excimer states, along with the La/Lb state inversion. In contrast with pyrene, the dipole-dipole interactions in 2-azapyrene stabilize the dark eclipsed excimer structure and increase the barrier for conversion into a bright twisted excimer. It is suggested that these differences in the nitrogen-substituted derivative might influence the excimer emission properties.

Short-Range Effects in the Special Pair of Photosystem II Reaction Centers: The Nonconservative Nature of Circular Dichroism

Photosystem II reaction centers extract electrons from water, providing the basis of oxygenic life on earth. Among the light-sensitive pigments of the reaction center, a central chlorophyll a dimer, known as the special pair, so far has escaped a complete theoretical characterization of its excited state properties. The close proximity of the special pair pigments gives rise to short-range effects that comprise a coupling between local and charge transfer (CT) excited states as well as other intermolecular quantum effects. Using a multiscale simulation and a diabatization technique, we show that the coupling to CT states is responsible for 45% of the excitonic coupling in the special pair. The other short-range effects cause a nonconservative nature of the circular dichroism spectrum of the reaction center by effectively rotating the electric transition dipole moments of the special pair pigments inverting and strongly enhancing their intrinsic rotational strength.

Molecular analyses of zebrafish V0v spinal interneurons and identification of transcriptional regulators downstream of Evx1 and Evx2 in these cells

BACKGROUND

V0v spinal interneurons are highly conserved, glutamatergic, commissural neurons that function in locomotor circuits. We have previously shown that Evx1 and Evx2 are required to specify the neurotransmitter phenotype of these cells. However, we still know very little about the gene regulatory networks that act downstream of these transcription factors in V0v cells.

METHODS

To identify candidate members of V0v gene regulatory networks, we FAC-sorted wild-type and evx1;evx2 double mutant zebrafish V0v spinal interneurons and expression-profiled them using microarrays and single cell RNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 double mutants and wild-type siblings.

RESULTS

Our data reveal two molecularly distinct subtypes of zebrafish V0v spinal interneurons at 48 h and suggest that, by this stage of development, evx1;evx2 double mutant cells transfate into either inhibitory spinal interneurons, or motoneurons. Our results also identify 25 transcriptional regulator genes that require Evx1/2 for their expression in V0v interneurons, plus a further 11 transcriptional regulator genes that are repressed in V0v interneurons by Evx1/2. Two of the latter genes are hmx2 and hmx3a. Intriguingly, we show that Hmx2/3a, repress dI2 interneuron expression of skor1a and nefma, two genes that require Evx1/2 for their expression in V0v interneurons. This suggests that Evx1/2 might regulate skor1a and nefma expression in V0v interneurons by repressing Hmx2/3a expression.

CONCLUSIONS

This study identifies two molecularly distinct subsets of zebrafish V0v spinal interneurons, as well as multiple transcriptional regulators that are strong candidates for acting downstream of Evx1/2 to specify the essential functional characteristics of these cells. Our data further suggest that in the absence of both Evx1 and Evx2, V0v spinal interneurons initially change their neurotransmitter phenotypes from excitatory to inhibitory and then, later, start to express markers of distinct types of inhibitory spinal interneurons, or motoneurons. Taken together, our findings significantly increase our knowledge of V0v and spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type.

Dominant role of excitons in photosynthetic color-tuning and light-harvesting

Photosynthesis is a vital process that converts sunlight into energy for the Earth's ecosystems. Color adaptation is crucial for different photosynthetic organisms to thrive in their ecological niches. Although the presence of collective excitons in light-harvesting complexes is well known, the role of delocalized excited states in color tuning and excitation energy transfer remains unclear. This study evaluates the characteristics of photosynthetic excitons in sulfur and non-sulfur purple bacteria using advanced optical spectroscopic techniques at reduced temperatures. The exciton effects in these bacteriochlorophyll a-containing species are generally much stronger than in plant systems that rely on chlorophylls. Their exciton bandwidth varies based on multiple factors such as chromoprotein structure, surroundings of the pigments, carotenoid content, hydrogen bonding, and metal ion inclusion. The study nevertheless establishes a linear relationship between the exciton bandwidth and Qy singlet exciton absorption peak, which in case of LH1 core complexes from different species covers almost 130 nm. These findings provide important insights into bacterial color tuning and light-harvesting, which can inspire sustainable energy strategies and devices.

Origin of Low-Lying Red States in the Lhca4 Light-Harvesting Complex of Photosystem I

The antenna complexes of Photosystem I present low-lying states visible as red-shifted and broadened absorption and fluorescence bands. Among these, Lhca4 has the most evident features of these "red" states, with a fluorescence band shifted by more than 25 nm from typical LHC emission. This signal arises from a mixing of exciton and charge-transfer (CT) states within the excitonically coupled a603-a609 chlorophyll (Chl) dimer. Here we combine molecular dynamics, multiscale quantum chemical calculations, and spectral simulations to uncover the molecular mechanism for the formation and tuning of exciton-CT interactions in Lhca4. We show that the coupling between exciton and CT states is extremely sensitive to tiny variations in the Chl dimer arrangement, explaining both the red-shifted bands and the switch between conformations with blue and red emission observed in single-molecule spectroscopy. Finally, we show that mutating the axial ligand of a603 diminishes the exciton-CT coupling, removing any red-state fingerprint.

Molecular Analyses of V0v Spinal Interneurons and Identification of Transcriptional Regulators Downstream of Evx1 and Evx2 in these Cells

Background

V0v spinal interneurons are highly conserved, glutamatergic, commissural neurons that function in locomotor circuits. We have previously shown that Evx1 and Evx2 are required to specify the neurotransmitter phenotype of these cells. However, we still know very little about the gene regulatory networks that act downstream of these transcription factors in V0v cells.

Methods

To identify candidate members of V0v gene regulatory networks, we FAC-sorted WT and evx1;evx2 double mutant zebrafish V0v spinal interneurons and expression-profiled them using microarrays and single cell RNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 double mutants and wild-type siblings.

Results

Our data reveal two molecularly distinct subtypes of V0v spinal interneurons at 48 h and suggest that, by this stage of development, evx1;evx2 double mutant cells transfate into either inhibitory spinal interneurons, or motoneurons. Our results also identify 25 transcriptional regulator genes that require Evx1/2 for their expression in V0v interneurons, plus a further 11 transcriptional regulator genes that are repressed in V0v interneurons by Evx1/2. Two of the latter genes are hmx2 and hmx3a. Intriguingly, we show that Hmx2/3a, repress dI2 interneuronal expression of skor1a and nefma, two genes that require Evx1/2 for their expression in V0v interneurons. This suggests that Evx1/2 might regulate skor1a and nefma expression in V0v interneurons by repressing Hmx2/3a expression.

Conclusions

This study identifies two molecularly distinct subsets of V0v spinal interneurons, as well as multiple transcriptional regulators that are strong candidates for acting downstream of Evx1/2 to specify the essential functional characteristics of these cells. Our data further suggest that in the absence of both Evx1 and Evx2, V0v spinal interneurons initially change their neurotransmitter phenotypes from excitatory to inhibitory and then, later, start to express markers of distinct types of inhibitory spinal interneurons, or motoneurons. Taken together, our findings significantly increase our knowledge of V0v and spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type.

Understanding Primary Charge Separation in the Heliobacterial Reaction Center

The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.

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.

Characterisation of the photosynthetic complexes from the marine gammaproteobacterium Congregibacter litoralis KT71

Possibly the most abundant group of anoxygenic phototrophs are marine photoheterotrophic Gammaproteobacteria belonging to the NOR5/OM60 clade. As little is known about their photosynthetic apparatus, the photosynthetic complexes from the marine phototrophic bacterium Congregibacter litoralis KT71 were purified and spectroscopically characterised. The intra-cytoplasmic membranes contain a smaller amount of photosynthetic complexes when compared with anaerobic purple bacteria. Moreover, the intra-cytoplasmic membranes contain only a minimum amount of peripheral LH2 complexes. The complexes are populated by bacteriochlorophyll a, spirilloxanthin and two novel ketocarotenoids, with biophysical and biochemical properties similar to previously characterised complexes from purple bacteria. The organization of the RC-LH1 complex has been further characterised using cryo-electron microscopy. The overall organisation is similar to the complex from the gammaproteobacterium Thermochromatium tepidum, with the type-II reaction centre surrounded by a slightly elliptical LH1 antenna ring composed of 16 αβ-subunits with no discernible gap or pore. The RC-LH1 and LH2 apoproteins are phylogenetically related to other halophilic species but LH2 also to some alphaproteobacterial species. It seems that the reduction of light-harvesting apparatus and acquisition of novel ketocarotenoids in Congregibacter litoralis KT71 represent specific adaptations for operating the anoxygenic photosynthesis under aerobic conditions at sea.

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.

Conversion of B800 Bacteriochlorophyll a to 3-Acetyl Chlorophyll a in the Light-Harvesting Complex 3 by In Situ Oxidation

The spectral features of energy donors and acceptors and the relationship between them in photosynthetic light-harvesting proteins are crucial for photofunctions of these proteins. Engineering energy donors and acceptors in light-harvesting proteins affords the means to increase our understanding of their photofunctional mechanisms. Herein, we demonstrate the conversion of energy-donating B800 bacteriochlorophyll (BChl) a to 3-acetyl chlorophyll (AcChl) a in light-harvesting complex 3 (LH3) from Rhodoblastus acidophilus by in situ oxidation with 2,3-dichloro-5,6-dicyano-1,4-benzoquinone. AcChl a in the B800 site exhibited a Q_y band that was 111 nm blue-shifted with respect to BChl a in oxidized LH3. The structure of LH3 was barely influenced by the oxidation process, based on circular dichroism spectroscopy and size-exclusion chromatography evidence. In oxidized LH3, AcChl a transferred excitation energy to B820 BChl a, but the rate of excitation energy transfer (EET) was lower than in native LH3. The intracomplex EET in oxidized LH3 was slightly faster than in oxidized light-harvesting complex 2 (LH2). This difference is rationalized by an increase in overlap of the luminescence band of AcChl a with the long tail of the B820 absorption band in oxidized LH3 compared with that of the B850 band in oxidized LH2.

Computing the electronic circular dichroism spectrum of DNA quadruple helices of different topology: A critical test for a generalized excitonic model based on a fragment diabatization

In this study, we exploit a recently developed fragment diabatization-based excitonic model, FrDEx, to simulate the electronic circular dichroism (ECD) spectra of three guanine-rich DNA sequences arranged in guanine quadruple helices with different topologies: thrombin binding aptamer (antiparallel), c-Myc promoter (parallel), and human telomeric sequence (3+1 hybrid). Starting from time-dependent density functional theory (TD-DFT) calculations with the M052X functional, we apply our protocol to parameterize the FrDEX Hamiltonian, which accounts for electron density overlap and includes both the coupling with charge transfer transitions and the effect of the surrounding bases on the local excitation of each chromophore. The TD-DFT/M052X spectral shapes are in good agreement with the experimental ones, the main source of discrepancy being related to the intrinsic error on the computed transition energies of guanine monomer. FrDEx spectra are fairly close to the reference TD-DFT ones, allowing a significant advance with respect to a more standard excitonic Hamiltonian. We also show that the ECD spectra are sensitive to the inclusion of the inner K + $$ {}^{+} $$ cation in the calculation.

Characterization of excited states in time-dependent density functional theory using localized molecular orbitals

Localized molecular orbitals are often used for the analysis of chemical bonds, but they can also serve to efficiently and comprehensibly compute linear response properties. While conventional canonical molecular orbitals provide an adequate basis for the treatment of excited states, a chemically meaningful identification of the different excited-state processes is difficult within such a delocalized orbital basis. In this work, starting from an initial set of supermolecular canonical molecular orbitals, we provide a simple one-step top-down embedding procedure for generating a set of orbitals, which are localized in terms of the supermolecule but delocalized over each subsystem composing the supermolecule. Using an orbital partitioning scheme based on such sets of localized orbitals, we further present a procedure for the construction of local excitations and charge-transfer states within the linear response framework of time-dependent density functional theory (TDDFT). This procedure provides direct access to approximate diabatic excitation energies and, under the Tamm-Dancoff approximation, also their corresponding electronic couplings-quantities that are of primary importance in modeling energy transfer processes in complex biological systems. Our approach is compared with a recently developed diabatization procedure based on subsystem TDDFT using projection operators, which leads to a similar set of working equations. Although both of these methods differ in the general localization strategies adopted and the type of basis functions (Slaters vs Gaussians) employed, an overall decent agreement is obtained.

Intermolecular Interactions and Charge Resonance Contributions to Triplet and Singlet Exciton States of Oligoacene Aggregates

Intermolecular interactions modulate the electro-optical properties of molecular materials and the nature of low-lying exciton states. Molecular materials composed by oligoacenes are extensively investigated for their semiconducting and optoelectronic properties. Here, we analyze the exciton states derived from time-dependent density functional theory (TDDFT) calculations for two oligoacene model aggregates: naphthalene and anthracene dimers. To unravel the role of inter-molecular interactions, a set of diabatic states is selected, chosen to coincide with local (LE) and charge-transfer (CT) excitations within a restricted orbital space including two occupied and two unoccupied orbitals for each molecular monomer. We study energy profiles and disentangle inter-state couplings to disclose the (CT) character of singlet and triplet exciton states and assess the influence of inter-molecular orientation by displacing one molecule with respect to the other along the longitudinal translation coordinate. The analysis shows that (CT) contributions are relevant, although comparably less effective for triplet excitons, and induce a non-negligible mixed character to the low-lying exciton states for eclipsed monomers and for small translational displacements. Such (CT) contributions govern the L_a/L_b state inversion occurring for the low-lying singlet exciton states of naphthalene dimer and contribute to the switch from H- to J-aggregate type of the strongly allowed B_b transition of both oligoacene aggregates.

Phasor Analysis Reveals Multicomponent Fluorescence Kinetics in the LH2 Complex from Marichromatium purpuratum

We investigated the fluorescence kinetics of LH2 complexes from Marichromatium purpuratum, the cryo-EM structure of which has been recently elucidated with 2.4 Å resolution. The experiments have been carried out as a function of the excitation density by varying both the excitation fluence and the repetition rate of the laser excitation. Instead of the usual multiexponential fitting procedure, we applied the less common phasor formalism for evaluating the transients because this allows for a model-free analysis of the data without a priori knowledge about the number of processes that contribute to a particular decay. For the various excitation conditions, this analysis reproduces consistently three lifetime components with decay times below 100 ps, 500 ps, and 730 ps, which were associated with the quenched state, singlet-triplet annihilation, and fluorescence decay, respectively. Moreover, it reveals that the number of decay components that contribute to the transients depends on whether the excitation wavelength is in resonance with the B800 BChl a molecules or with the carotenoids. Based on the mutual arrangement of the chromophores in their binding pockets, this leads us to conclude that the energy transfer pathways within the LH2 complex of this species differ significantly from each other for exciting either the B800 BChl molecules or the carotenoids. Finally, we speculate whether the illumination with strong laser light converts the LH2 complexes studied here into a quenched conformation that might be related to the development of the non-photochemical quenching mechanism that occurs in higher plants.

Effects of Detergents on the Spectral Features of B820 Bacteriochlorophyll a in Light-Harvesting Complex 3

Excitonic coupling of bacteriochlorophyll (BChl) a in light-harvesting (LH) proteins of purple photosynthetic bacteria is key for efficient photon capture and energy transfer. Environmental factors can affect the spectral features of these BChl a pigments and investigating these effects can provide insight into the molecular mechanisms underlying the photosynthetic spectral tuning. The present study analyzes the spectral alterations of the Q_y band of B820 BChl a within the LH3 protein in relation to the type and concentration of detergents in the buffer. Changing the detergent from lauryl dimethylamine N-oxide (LDAO) to n-dodecyl-β-d-maltoside (DDM) caused a red shift in the B820 Q_y band accompanied by hyperchromism; these spectral alterations were completely reversed by exchanging back from DDM to LDAO. These results reflect the different effects of harsh vs mild detergents on the perturbation of LH3. The B820 Q_y band did not change when LDAO or NaCl concentration was altered, suggesting that electrostatic effects by external components have little influence on the spectral features of B820 BChl a in LH3.

Excitonic Nature of Carotenoid–Phthalocyanine Dyads and Its Role in Transient Absorption Spectra

Artificial carotenoid–tetrapyrrole dyads have been extensively used as model systems to understand the quenching mechanisms that occur in light-harvesting complexes during nonphotochemical quenching. In particular, dyads containing a carotenoid covalently linked to a zinc phthalocyanine have been studied by transient absorption spectroscopy, and the observed signals have been interpreted in terms of an excitonically coupled state involving the lowest excited states of the two fragments. If present, such excitonic delocalization would have significant implications on the mechanism of nonphotochemical quenching. Here, we use quantum chemical calculations to show that this delocalization is not needed to reproduce the transient absorption spectra. On the contrary, the observed signals can be explained through excitonic couplings in the higher-energy manifold of states. We also argue that the covalent linkage between the two fragments allows for electronic communications, which complicates the analysis of the spectra based on two independent but coupled moieties. These findings call for a more thorough reassessment of the photophysics in these dyads and its implications in the context of natural nonphotochemical quenching.

Exciton transport in molecular organic semiconductors boosted by transient quantum delocalization

Designing molecular materials with very large exciton diffusion lengths would remove some of the intrinsic limitations of present-day organic optoelectronic devices. Yet, the nature of excitons in these materials is still not sufficiently well understood. Here we present Frenkel exciton surface hopping, an efficient method to propagate excitons through truly nano-scale materials by solving the time-dependent Schrödinger equation coupled to nuclear motion. We find a clear correlation between diffusion constant and quantum delocalization of the exciton. In materials featuring some of the highest diffusion lengths to date, e.g. the non-fullerene acceptor Y6, the exciton propagates via a transient delocalization mechanism, reminiscent to what was recently proposed for charge transport. Yet, the extent of delocalization is rather modest, even in Y6, and found to be limited by the relatively large exciton reorganization energy. On this basis we chart out a path for rationally improving exciton transport in organic optoelectronic materials.

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.

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 Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory

The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.

Impact of Charge-Resonance Excitations on CT-Mediated J-Type Aggregation in Singlet and Triplet Exciton States of Perylene Di-Imide Aggregates: A TDDFT Investigation

The modulation of intermolecular interactions upon aggregation induces changes in excited state properties of organic molecules that can be detrimental for some optoelectronic applications but can be exploited for others. The time-dependent density functional theory (TDDFT) is a cost-effective approach to determining the exciton states of molecular aggregates, and it has been shown to provide reliable results when coupled with the appropriate choice of the functional. Here we apply a general procedure to analyze the aggregates’ exciton states derived from TDDFT calculations in terms of diabatic states chosen to coincide with local (LE) and charge-transfer (CT) excitations within a restricted orbital space. We apply the approach to study energy profiles, interstate couplings, and the charge-transfer character of singlet and triplet exciton states of perylene di-imide aggregates (PDI). We focus on the intermolecular displacement along the longitudinal translation coordinate, which mimics different amounts of slip-stacking observed in PDI crystals. The analysis, in terms of symmetry-adapted Frenkel excitations (FE) and charge-resonance (CR) states and their interactions, discloses how the interchange of the H/J character for small longitudinal shifts, previously reported for singlet exciton states, also occurs for triplet excitons.

Low-Scaling Excited State Calculation Using the Block Interaction Product State

We develop an automatic and efficient scheme for the accurate construction of the bases for excitonic models, which can enable "black-box" excited state structure calculations for large molecular systems. These new and optimized bases, which are named the block interaction product state (BIPS), can be expressed as the direct products of the local states for each chromophore. Each chromophore's local states are selected by diagonalization of its reduced density matrix, which is obtained by quantum chemical calculation of the small subsystem composed of the chromophore and its nearest neighbors. We implemented the BIPS framework with fragment-based calculations considering two- and three-body interactions. Test calculations for eight different molecular aggregates demonstrate that this framework provides an accurate description of not only the excitation energies but also the first-order wave function properties (dipole moment and transition dipole moment) of the low-lying excited states at a low-scaling computational cost.

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.

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.

Understanding the Relation between Structural and Spectral Properties of Light-Harvesting Complex II

Light-harvesting complex II (LHCII) is a pigment-protein complex present in higher plants and green algae. LHCII represents the main site of light absorption, and its role is to transfer the excitation energy toward the photosynthetic reaction centers, where primary energy conversion reactions take place. The optical properties of LHCII are known to depend on protein conformation. However, the relation between the structural and spectroscopic properties of the pigments is not fully understood yet. In this respect, previous classical molecular dynamics simulations of LHCII in a model membrane [Sci. Rep. 2015, 5, 1-10] have shown that the configuration and excitonic coupling of a chlorophyll (Chl) dimer functioning as the main terminal emitter of the complex are particularly sensitive to conformational changes. Here, we use quantum chemistry calculations to investigate in greater detail the effect of pigment-pigment interactions on the excited-state landscape. While most previous studies have used a local picture in which electrons are localized on single pigments, here we achieve a more accurate description of the Chl dimer by adopting a supramolecular picture where time-dependent density functional theory is applied to the whole system at once. Our results show that specific dimer configurations characterized by shorter inter-pigment distances can result in a sizable intensity decrease (up to 36%) of the Chl absorption bands in the visible spectral region. Such a decrease can be predicted only when accounting for Chl-Chl charge-transfer excitations, which is possible using the above-mentioned supramolecular approach. The charge-transfer character of the excitations is quantified by two types of analyses: one focusing on the composition of the excitations and the other directly on the observable total absorption intensities.

Addressing the Frenkel and charge transfer character of exciton states with a model Hamiltonian based on dimer calculations: Application to large aggregates of perylene bisimide

To understand the influence of interchromophoric arrangements on photo-induced processes and optical properties of aggregates, it is fundamental to assess the contribution of local excitations [charge transfer (CT) and Frenkel (FE)] to exciton states. Here, we apply a general procedure to analyze the adiabatic exciton states derived from time-dependent density functional theory calculations, in terms of diabatic states chosen to coincide with local excitations within a restricted orbital space. In parallel, motivated by the need of cost-effective approaches to afford the study of larger aggregates, we propose to build a model Hamiltonian based on calculations carried out on dimers composing the aggregate. Both approaches are applied to study excitation energy profiles and CT character modulation induced by interchromophore rearrangements in perylene bisimide aggregates up to a tetramer. The dimer-based approach closely reproduces the results of full-aggregate calculations, and an analysis in terms of symmetry-adapted diabatic states discloses the effects of CT/FE interactions on the interchange of the H-/J-character for small longitudinal shifts of the chromophores.

Electronic Couplings for Photoinduced Charge Transfer and Excitation Energy Transfer Based on Fragment Particle-Hole Densities

A new scheme is proposed to calculate the electronic couplings for photoinduced charge transfer and excitation energy transfer for both singlet and triplet states. In this scheme, the locally excited and charge-transfer states are constructed from the adiabatic ones by maximally localizing the particle (i.e., electron) and hole densities in terms of predefined molecular fragments. The construction process, after which the electronic couplings are directly obtained, is highly efficient and can be combined with various kinds of preliminary electronic structure calculations as long as the adiabatic excitation energies and transition densities are available. The method also applies to the systems with multiple charge or excitation centers. Its validity is demonstrated by the applications to the 6,13-dichloropentacene dimer and tetramer and the C₆₀-Zn porphyrin dyad. The results reveal that the environment has a strong impact on the electronic couplings and can even enlarge those for long-range charge transfer by several orders of magnitude.

Excitonic Model for Strongly Coupled Multichromophoric Systems: The Electronic Circular Dichroism Spectra of Guanine Quadruplexes as Test Cases

We here propose a general and flexible approach, based on fragment diabatization, which incorporates charge transfer states and significantly increases the reliability of excitonic Hamiltonians for systems where the chromophores are very close. This model (FrDEx) is used to compute the electronic circular dichroism and absorption spectra of two prototype guanine-rich DNA sequences folded in quadruple helices (GQs), i.e., a fragment of the human telomeric sequence (Tel21, antiparallel), and (TGGGGT)₄ (TG4T, parallel). Calculations on different subsets of Tel21 and TG4T, from dimers to tetramers, show that FrDEx provides spectra close to the reference full quantum mechanical (QM) ones (obtained with time-dependent density functional theory), with significant improvements with respect to "standard" excitonic Hamiltonians. Furthermore, these tests enable the most cost-effective procedure for the whole GQ to be determined. FrDEx spectra of Tel21 and TG4T are also in good agreement with the QM and experimental ones and give access to interesting insights into the chemical-physical effects modulating the spectral signals. FrDEx could be profitably used to investigate many other biological and nanotechnological materials, from DNA to (opto)electronic polymers.

Anharmonic Molecular Motion Drives Resonance Energy Transfer in peri-Arylene Dyads

Spectral and dynamical properties of molecular donor-acceptor systems strongly depend on the steric arrangement of the constituents with exciton coupling J as a key control parameter. In the present work we study two peri-arylene based dyads with orthogonal and parallel transition dipoles for donor and acceptor moieties, respectively. We show that the anharmonic multi-well character of the orthogonal dyad's intramolecular potential explains findings from both stationary and time-resolved absorption experiments. While for a parallel dyad, standard quantum chemical estimates of J at 0 K are in good agreement with experimental observations, J becomes vanishingly small for the orthogonal dyad, in contrast to its ultrafast experimental transfer times. This discrepancy is not resolved even by accounting for harmonic fluctuations along normal coordinates. We resolve this problem by supplementing quantum chemical approaches with dynamical sampling of fluctuating geometries. In contrast to the moderate Gaussian fluctuations of J for the parallel dyad, fluctuations for the orthogonal dyad are found to follow non-Gaussian statistics leading to significantly higher effective J in good agreement with experimental observations. In effort to apply a unified framework for treating the dynamics of optical coherence and excitonic populations of both dyads, we employ a vibronic approach treating electronic and selected vibrational degrees on an equal footing. This vibronic model is used to model absorption and fluorescence spectra as well as donor-acceptor transport dynamics and covers the more traditional categories of Förster and Redfield transport as limiting cases.

Electronic couplings for photo-induced processes from subsystem time-dependent density-functional theory: The role of the diabatization

Subsystem time-dependent density-functional theory (sTDDFT) making use of approximate non-additive kinetic energy (NAKE) functionals is known to be capable of describing excitation energy transfer processes in a variety of applications. Here, we show that sTDDFT, especially when combined with projection-based embedding (PbE), can be employed for the entire range of photo-induced electronic couplings essential for modeling photophysical properties of complex chemical and biological systems and therefore represents a complete toolbox for this class of problems. This means that it is capable of capturing the interaction/coupling associated with local- and charge-transfer (CT) excitons. However, this requires the choice of a reasonable diabatic basis. We therefore propose different diabatization strategies of the virtual orbital space in PbE-sTDDFT and show how CT excitations can be included in sTDDFT using NAKE functionals via a phenomenological approach. Finally, these electronic couplings are compared to couplings from a multistate fragment excitation difference (FED)-fragment charge difference (FCD) diabatization procedure. We show that both procedures, multistate FED-FCD and sTDDFT (with the right diabatization procedure chosen), lead to an overall good agreement for the electronic couplings, despite differences in their general diabatization strategy. We conclude that the entire range of photo-induced electronic couplings can be obtained using sTDDFT (with the right diabatization procedure chosen) in a black-box manner.

DFTB/MM Molecular Dynamics Simulations of the FMO Light-Harvesting Complex

Because of the size of light-harvesting complexes and the involvement of electronic degrees of freedom, computationally these systems need to be treated with a combined quantum-classical description. To this end, Born-Oppenheimer molecular dynamics simulations have been employed in a quantum mechanics/molecular mechanics (QM/MM) fashion for the ground state followed by excitation energy calculations again in a QM/MM scheme for the Fenna-Matthews-Olson (FMO) complex. The self-consistent-charge density functional tight-binding (DFTB) method electrostatically coupled to a classical description of the environment was applied to perform the ground-state dynamics. Subsequently, long-range-corrected time-dependent DFTB calculations were performed to determine the excitation energy fluctuations of the individual bacteriochlorophyll a molecules. The spectral densities obtained using this approach show an excellent agreement with experimental findings. In addition, the fluctuating site energies and couplings were used to estimate the exciton transfer dynamics.

A comparative look at structural variation among RC-LH1 'Core' complexes present in anoxygenic phototrophic bacteria

All purple photosynthetic bacteria contain RC-LH1 'Core' complexes. The structure of this complex from Rhodobacter sphaeroides, Rhodopseudomonas palustris and Thermochromatium tepidum has been solved using X-ray crystallography. Recently, the application of single particle cryo-EM has revolutionised structural biology and the structure of the RC-LH1 'Core' complex from Blastochloris viridis has been solved using this technique, as well as the complex from the non-purple Chloroflexi species, Roseiflexus castenholzii. It is apparent that these structures are variations on a theme, although with a greater degree of structural diversity within them than previously thought. Furthermore, it has recently been discovered that the only phototrophic representative from the phylum Gemmatimonadetes, Gemmatimonas phototrophica, also contains a RC-LH1 'Core' complex. At present only a low-resolution EM-projection map exists but this shows that the Gemmatimonas phototrophica complex contains a double LH1 ring. This short review compares these different structures and looks at the functional significance of these variations from two main standpoints: energy transfer and quinone exchange.

Comparison of the Energy-Transfer Rates in Structural and Spectral Variants of the B800-850 Complex from Purple Bacteria

Photosynthetic light harvesting can occur with a remarkable near-unity quantum efficiency. The B800-850 complex, also known as light-harvesting complex 2 (LH2), is the primary light-harvesting complex in purple bacteria and has been extensively studied as a model system. The bacteriochlorophylls of the B800-850 complex are organized into two concentric rings, known as the B800 and B850 rings. However, depending on the species and growth conditions, the number of constituent subunits, the pigment geometry, and the absorption energies vary. While the dynamics of some B800-850 variants have been exhaustively characterized, others have not been measured. Furthermore, a direct and simultaneous comparison of how both structural and spectral differences between variants affect these dynamics has not been performed. In this work, we utilize ultrafast transient absorption measurements to compare the B800 to B850 energy-transfer rates in the B800-850 complex as a function of the number of subunits, geometry, and absorption energies. The nonameric B800-850 complex from Rhodobacter (Rb.) sphaeroides is 40% faster than the octameric B800-850 complex from Rhodospirillum (Rs.) molischianum, consistent with structure-based predictions. In contrast, the blue-shifted B800-820 complex from Rs. molischianum is only 20% faster than the B800-850 complex from Rs. molischianum despite an increase in the spectral overlap between the rings that would be expected to produce a larger increase in the energy-transfer rate. These measurements support current models that contain dark, higher-lying excitonic states to bridge the energy gap between rings, thereby maintaining similar energy-transfer dynamics. Overall, these results demonstrate that energy-transfer dynamics in the B800-850 complex are robust to the spectral and structural variations between species used to optimize energy capture and flow in purple bacteria.

Charge transfer from the carotenoid can quench chlorophyll excitation in antenna complexes of plants

The photosynthetic apparatus of higher plants can dissipate excess excitation energy during high light exposure, by deactivating excited chlorophylls through a mechanism called nonphotochemical quenching (NPQ). However, the precise molecular details of quenching and the mechanism regulating the quenching level are still not completely understood. Focusing on the major light-harvesting complex LHCII of Photosystem II, we show that a charge transfer state involving Lutein can efficiently quench chlorophyll excitation, and reduce the excitation lifetime of LHCII to the levels measured in the deeply quenched LHCII aggregates. Through a combination of molecular dynamics simulations, multiscale quantum chemical calculations, and kinetic modeling, we demonstrate that the quenching level can be finely tuned by the protein, by regulating the energy of the charge transfer state. Our results suggest that a limited conformational rearrangement of the protein scaffold could act as a molecular switch to activate or deactivate the quenching mechanism.

On the arrangement of chromophores in light harvesting complexes: chance versus design

We used a homogeneous computational approach to derive the excitonic Hamiltonian for five light harvesting complexes containing only one type of chromophore and compare them in terms of statistical descriptors. We then studied the approximate exciton dynamics for the five complexes introducing a measure, the (averaged and time-dependent) inverse participation ratio, that enables the comparison between very diverse complexes on the same ground. We find that the global dynamics are very similar across the set of systems despite the variety of geometric structures of the complexes. In particular, the dynamics of four out of five light harvesting complexes are barely distinguishable with a small variation from the norm seen only for the Fenna-Matthews-Olson complex. We use the information from the realistic Hamiltonians to build a reduced model system that shows how the global dynamics are ultimately dominated by a single parameter, the degree of localization of the excitonic Hamiltonian eigenstates. Considering the physically plausible range of system parameters, the reduced model explains why the dynamics are so similar across most light harvesting complexes containing a single type of chromophore regardless of the detailed pattern of the inter-chromophore excitonic coupling.

The molecular mechanisms of light adaption in light-harvesting complexes of purple bacteria revealed by a multiscale modeling

The light-harvesting in photosynthetic purple bacteria can be tuned in response to the light conditions during cell growth. One of the used strategies is to change the energy of the excitons in the major fight-harvesting complex, commonly known as LH2. In the present study we report the first systematic investigation of the microscopic origin of the exciton tuning using three complexes, namely the common (high-light) and the low-light forms of LH2 from Rps. acidophila plus a third complex analogous to the PucD complex from Rps. palustris. The study is based on the combination of classical molecular dynamics of each complex in a lipid membrane and excitonic calculations based on a multiscale quantum mechanics/molecular mechanics approach including a polarizable embedding. From the comparative analysis, it comes out that the mechanisms that govern the adaptation of the complex to different light conditions use the different H-bonding environment around the bacteriochlorophyll pigments to dynamically control both internal and inter-pigment degrees of freedom. While the former have a large effect on the site energies, the latter significantly change the electronic couplings, but only the combination of the two effects can fully reproduce the tuning of the final excitons and explain the observed spectroscopic differences.