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

Theory of photosynthetic membrane influence on B800-B850 energy transfer in the LH2 complex

Photosynthetic organisms rely on a network of light-harvesting protein-pigment complexes to efficiently absorb sunlight and transfer excitation energy to reaction center proteins where charge separation occurs. In photosynthetic purple bacteria, these complexes are embedded within the cell membrane, with lipid composition affecting complex clustering, thereby impacting inter-complex energy transfer. However, the impact of the lipid bilayer on intra-complex excitation dynamics is less understood. Recent experiments have addressed this question by comparing photo-excitation dynamics in detergent-isolated light-harvesting complex 2 (LH2) to LH2 complexes embedded in membrane discs mimicking the biological environment, revealing differences in spectra and energy-transfer rates. In this paper, we use available quantum chemical and spectroscopy data to develop a complementary theoretical study on the excitonic structure and intra-complex energy-transfer kinetics of the LH2 of photosynthetic purple bacteria Rhodoblastus (Rbl.) acidophilus (formerly Rhodopseudomonas acidophila) in two different conditions: the LH2 in a membrane environment and detergent-isolated LH2. We find that dark excitonic states, crucial for B800-B850 energy transfer within LH2, are more delocalized in the membrane model. Using nonperturbative and generalized Förster calculations, we show that such increased quantum delocalization results in a 30% faster B800 to B850 transfer rate in the membrane model, in agreement with experimental results. We identify the dominant energy-transfer pathways in each environment and demonstrate how differences in the B800 to B850 transfer rate arise from changes in LH2's electronic properties when embedded in the membrane. Furthermore, by accounting for the quasi-static variations of electronic excitation energies in the LH2, we show that the broadening of the distribution of the B800-B850 transfer rates is affected by the lipid composition. We argue that such variation in broadening could be a signature of a speed-accuracy trade-off, commonly seen in biological process.

Computational Design of (B)Chl Models: Structural and Chemical Modifications toward Enriched Properties

The functional units of natural photosynthetic systems control the process of converting sunlight into chemical energy. In this article, we explore a series of chemically and structurally modified bacteriochlorophyll and chlorophyll pigments through computational chemistry to evaluate their electronic spectroscopy properties. More specifically, we use multiconfigurational and time-dependent density functional theory methods, along with molecular dynamics simulations, to compute the models' energetics both in an implicit and explicit solvent environment. Structural modifications aimed at reducing the planarity of the macrocycle through alkyl-bridge anchoring reveal the significant role of the curvature in fine-tuning spectral properties, which mimics protein scaffold effects on naturally occurring pigments. Furthermore, chemical substitutions with a carbonyl group show potential for expanding absorption spectra toward the blue region, while incorporating an additional double bond decreases absorption efficiency. These insights lay the groundwork to design novel synthetic pigments, with potential applications in artificial light-harvesting systems and more efficient photovoltaic devices.

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.

Tuning by Hydrogen Bonding in Photosynthesis

Hydrogen bonding plays a crucial role in stabilizing proteins throughout their folding process. In photosynthetic light-harvesting chromoproteins, enriched with pigment chromophores, hydrogen bonds also fine-tune optical absorption to align with the solar irradiation spectrum. Despite its significance for photosynthesis, the precise mechanism of spectral tuning through hydrogen bonding remains inadequately understood. This study investigates wild-type and genetically engineered LH2 and LH1 light-harvesting complexes from Rhodobacter sphaeroides using a unique set of advanced spectroscopic techniques combined with simple exciton modeling. Our findings reveal an intricate interplay between exciton and site energy shift mechanisms, challenging the prevailing belief that spectral changes observed in these complexes upon the modification of tertiary structure hydrogen bonds almost directly follow shifting site energies. These deeper insights into natural adaptation processes hold great promise for advancing sustainable solar energy conversion technologies.

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.

General framework for quantifying dissipation pathways in open quantum systems. II. Numerical validation and the role of non-Markovianity

In the previous paper [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111-1-214111-13 (2024)], we developed a theory called MQME-D, which allows us to decompose the overall energy dissipation process in open quantum system dynamics into contributions by individual components of the bath when the subsystem dynamics is governed by a Markovian quantum master equation (MQME). Here, we contrast the predictions of MQME-D against the numerically exact results obtained by combining hierarchical equations of motion (HEOM) with a recently reported protocol for monitoring the statistics of the bath. Overall, MQME-D accurately captures the contributions of specific bath components to the overall dissipation while greatly reducing the computational cost compared to exact computations using HEOM. The computations show that MQME-D exhibits errors originating from its inherent Markov approximation. We demonstrate that its accuracy can be significantly increased by incorporating non-Markovianity by exploiting time scale separations (TSS) in different components of the bath. Our work demonstrates that MQME-D combined with TSS can be reliably used to understand how energy is dissipated in realistic open quantum system dynamics.

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.

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.

Spectral changes of light-harvesting complex 2 lacking B800 bacteriochlorophyll a under neutral pH conditions

Exchange of B800 bacteriochlorophyll (BChl) a in light-harvesting complex 2 (LH2) is promising for a better understanding of the mechanism on intracomplex excitation energy transfer of this protein. Structural and spectroscopic properties of LH2 lacking B800 BChl a (B800-depleted LH2), which is an important intermediate protein in the B800 exchange, will be useful to tackle the energy transfer mechanism in LH2 by the B800 exchange strategy. In this study, we report a unique spectral change of B800-depleted LH2, in which the Qy absorption band of B800 BChl a is automatically recovered under neutral pH conditions. This spectral change was facilitated by factors for destabilization of LH2, namely, a detergent, lauryl dimethylamine N-oxide, and an increase in temperature. Spectral analyses in the preparation of an LH2 variant denoted as B800-recovered LH2 indicated that most BChl a that was released by decomposition of part of B800-depleted LH2 was a source of the production of B800-recovered LH2. Characterization of purified B800-recovered LH2 demonstrated that its spectroscopic and structural features was quite similar to those of native LH2. The current results indicate that the recovery of the B800 Qy band of B800-depleted LH2 originates from the combination of decomposition of part of B800-depleted LH2 and in situ reconstitution of BChl a into the B800 binding pockets of residual B800-depleted LH2, resulting in the formation of stable B800-recovered LH2.

Mapping electronic decoherence pathways in molecules

Establishing the fundamental chemical principles that govern molecular electronic quantum decoherence has remained an outstanding challenge. Fundamental questions such as how solvent and intramolecular vibrations or chemical functionalization contribute to the decoherence remain unanswered and are beyond the reach of state-of-the-art theoretical and experimental approaches. Here we address this challenge by developing a strategy to isolate electronic decoherence pathways for molecular chromophores immersed in condensed phase environments that enables elucidating how electronic quantum coherence is lost. For this, we first identify resonance Raman spectroscopy as a general experimental method to reconstruct molecular spectral densities with full chemical complexity at room temperature, in solvent, and for fluorescent and non-fluorescent molecules. We then show how to quantitatively capture the decoherence dynamics from the spectral density and identify decoherence pathways by decomposing the overall coherence loss into contributions due to individual molecular vibrations and solvent modes. We illustrate the utility of the strategy by analyzing the electronic decoherence pathways of the DNA base thymine in water. Its electronic coherences decay in [Formula: see text]30 fs. The early-time decoherence is determined by intramolecular vibrations while the overall decay by solvent. Chemical substitution of thymine modulates the decoherence with hydrogen-bond interactions of the thymine ring with water leading to the fastest decoherence. Increasing temperature leads to faster decoherence as it enhances the importance of solvent contributions but leaves the early-time decoherence dynamics intact. The developed strategy opens key opportunities to establish the connection between molecular structure and quantum decoherence as needed to develop chemical strategies to rationally modulate it.

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.

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.

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

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

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.

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.

Charge transfer as a mechanism for chlorophyll fluorescence concentration quenching

Highly concentrated solutions of chlorophyll display rapid fluorescence quenching. The same devastating energy loss is not seen in photosynthetic light-harvesting antenna complexes, despite the need for chromophores to be in close proximity to facilitate energy transfer. A promising, though unconfirmed mechanism for the observed quenching is energy transfer from an excited chlorophyll monomer to a closely associated chlorophyll pair that subsequently undergoes rapid nonradiative decay to the ground state via a short-lived intermediate charge-transfer state. In this work, we make use of newly emerging fast methods in quantum chemistry to assess the feasibility of this proposed mechanism. We calculate rate constants for the initial charge separation, based on Marcus free-energy surfaces extracted from molecular dynamics simulations of solvated chlorophyll pairs, demonstrating that this pathway will compete with fluorescence (i.e., drive quenching) at experimentally measured quenching concentrations. We show that the rate of charge separation is highly sensitive to interchlorophyll distance and the relative orientations of chromophores within a quenching pair. We discuss possible solvent effects on the rate of charge separation (and consequently the degree of quenching), using the light-harvesting complex II (LH2) protein from rps. acidophila as a specific example of how this process might be controlled in a protein environment. Crucially, we reveal that the LH2 antenna protein prevents quenching, even at the high chlorophyll concentrations required for efficient energy transfer, by restricting the range of orientations that neighboring chlorophyll pairs can adopt.

An efficient protocol for excited states of large biochromophores

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Q_y transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory's accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.

Cryo-EM structures of light-harvesting 2 complexes from Rhodopseudomonas palustris reveal the molecular origin of absorption tuning

The light-harvesting (LH) complexes of phototrophic bacteria absorb solar energy for photosynthesis, and it is important to understand how the protein components influence the way bound pigments absorb light. We studied the LH2 complexes of Rhodopseudomonas palustris, which are encoded by a multigene family. Various combinations of LH2 genes were deleted, yielding strains that assemble only one of the four types of LH2. Following purification, the structures of four LH2 complexes were determined by cryogenic electron microscopy, revealing a basic nonameric ring structure comprising nine αβ-polypeptide pairs. An additional hitherto unknown polypeptide, γ, was found in each structure that binds six further bacteriochlorophylls. Comparison of these different structures shows how nature tunes their ability to absorb different wavelengths of light.

The genomes of some purple photosynthetic bacteria contain a multigene puc family encoding a series of α- and β-polypeptides that together form a heterogeneous antenna of light-harvesting 2 (LH2) complexes. To unravel this complexity, we generated four sets of puc deletion mutants in Rhodopseudomonas palustris, each encoding a single type of pucBA gene pair and enabling the purification of complexes designated as PucA-LH2, PucB-LH2, PucD-LH2, and PucE-LH2. The structures of all four purified LH2 complexes were determined by cryogenic electron microscopy (cryo-EM) at resolutions ranging from 2.7 to 3.6 Å. Uniquely, each of these complexes contains a hitherto unknown polypeptide, γ, that forms an extended undulating ribbon that lies in the plane of the membrane and that encloses six of the nine LH2 αβ-subunits. The γ-subunit, which is located near to the cytoplasmic side of the complex, breaks the C9 symmetry of the LH2 complex and binds six extra bacteriochlorophylls (BChls) that enhance the 800-nm absorption of each complex. The structures show that all four complexes have two complete rings of BChls, conferring absorption bands centered at 800 and 850 nm on the PucA-LH2, PucB-LH2, and PucE-LH2 complexes, but, unusually, the PucD-LH2 antenna has only a single strong near-infared (NIR) absorption peak at 803 nm. Comparison of the cryo-EM structures of these LH2 complexes reveals altered patterns of hydrogen bonds between LH2 αβ-side chains and the bacteriochlorin rings, further emphasizing the major role that H bonds play in spectral tuning of bacterial antenna complexes.

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.

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.

Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes

Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.

Electrostatic profiling of photosynthetic pigments: implications for directed spectral tuning

Photosynthetic pigment-protein complexes harvest solar energy with a high quantum efficiency. Protein scaffolds are known to tune the spectral properties of embedded pigments principally through structured electrostatic environments. Although the physical nature of electrostatic tuning is straightforward, the precise spatial principles of electrostatic preorganization remain poorly explored for different protein matrices and incompletely characterized with respect to the intrinsic properties of different photosynthetic pigments. In this work, we study the electronic structure features associated with the lowest excited state of a series of eight naturally occurring (bacterio)chlorophylls and pheophytins to describe the precise topological differences in electrostatic potentials and hence determine intrinsic differences in the expected mode and impact of electrostatic tuning. The difference electrostatic potentials between the ground and first excited states are used as fingerprints. Both the spatial profile and the propensity for spectral tuning are found to be unique for each pigment, indicating spatially and directionally distinct modes of electrostatic tuning. The results define a specific partitioning of the protein matrix around each pigment as an aid to identify regions with a maximal impact on spectral tuning and have direct implications for dimensionality reduction in protein design and engineering. Thus, a quantum mechanical basis is provided for understanding, predicting, and ultimately designing sequence-modified or pigment-exchanged biological systems, as suggested for selected examples of pigment-reconstituted proteins.

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.

Non-Uniform Excited State Electronic-Vibrational Coupling of Pigment-Protein Complexes

Photosynthetic organisms exploit interacting quantum degrees of freedom, namely intrapigment electron-vibrational (vibronic) and interpigment dipolar couplings (J-coupling), to rapidly and efficiently convert light into chemical energy. These interactions result in wave function configurations that delocalize excitation between pigments and pigment vibrations. Our study uses multidimensional spectroscopy to compare two model photosynthetic proteins, the Fenna-Matthews Olson (FMO) complex and light harvesting 2 (LH2), and confirm that long-lived excited state coherences originate from the vibrational modes of the pigment. Within this framework, the J-coupling of vibronic pigments should have a cascading effect in modifying the structured spectral density of excitonic states. We show that FMO effectively couples all of its excitations to a uniform set of vibrations while in LH2, its two chromophore rings each couple to a unique vibrational environment. We simulate energy transfer in a simple model system with non-uniform vibrational coupling to demonstrate how modification of the vibronic coupling strength can modulate energy transfer. Because increasing vibronic coupling increases internal relaxation, strongly coupled vibronic states can act as an energy funnel, which can potentially benefit energy transport.

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.

Protein Matrix Control of Reaction Center Excitation in Photosystem II

Photosystem II (PSII) is a multisubunit pigment-protein complex that uses light-induced charge separation to power oxygenic photosynthesis. Its reaction center chromophores, where the charge transfer cascade is initiated, are arranged symmetrically along the D1 and D2 core polypeptides and comprise four chlorophyll (PD1, PD2, ChlD1, ChlD2) and two pheophytin molecules (PheoD1 and PheoD2). Evolution favored productive electron transfer only via the D1 branch, with the precise nature of primary excitation and the factors that control asymmetric charge transfer remaining under investigation. Here we present a detailed atomistic description for both. We combine large-scale simulations of membrane-embedded PSII with high-level quantum-mechanics/molecular-mechanics (QM/MM) calculations of individual and coupled reaction center chromophores to describe reaction center excited states. We employ both range-separated time-dependent density functional theory and the recently developed domain based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD), the first coupled cluster QM/MM calculations of the reaction center. We find that the protein matrix is exclusively responsible for both transverse (chlorophylls versus pheophytins) and lateral (D1 versus D2 branch) excitation asymmetry, making ChlD1 the chromophore with the lowest site energy. Multipigment calculations show that the protein matrix renders the ChlD1 → PheoD1 charge-transfer the lowest energy excitation globally within the reaction center, lower than any pigment-centered local excitation. Remarkably, no low-energy charge transfer states are located within the "special pair" PD1-PD2, which is therefore excluded as the site of initial charge separation in PSII. Finally, molecular dynamics simulations suggest that modulation of the electrostatic environment due to protein conformational flexibility enables direct excitation of low-lying charge transfer states by far-red light.

Absorption and Circular Dichroism Spectra of Molecular Aggregates With the Full Cumulant Expansion

The exciton Hamiltonian of multichromophoric aggregates can be probed by spectroscopic techniques such as linear absorption and circular dichroism. To compare calculated Hamiltonians to experiments, a lineshape theory is needed, which takes into account the coupling of the excitons with inter- and intramolecular vibrations. This coupling is normally introduced in a perturbative way through the cumulant expansion formalism and further approximated by assuming a Markovian exciton dynamics, for example with the modified Redfield theory. Here, we present the implementation of the full cumulant expansion (FCE) formalism (J. Chem. Phys.142, 2015, 09410625747060) to efficiently compute absorption and circular dichroism spectra of molecular aggregates beyond the Markov approximation, without restrictions on the form of exciton–phonon coupling. By employing the LH2 system of purple bacteria as a challenging test case, we compare the FCE lineshapes with the Markovian lineshapes obtained with the modified Redfield theory, showing that the latter presents a less satisfying agreement with experiments. The FCE approach instead accurately describes the lineshapes, especially in the vibronic sideband of the B800 peak. We envision that the FCE approach will become a valuable tool for accurately comparing model exciton Hamiltonians with optical spectroscopy experiments.

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

We present a fully atomistic simulation of linear optical spectra (absorption, fluorescence and circular dichroism) of the Light Harvesting Complex II (LHCII) trimer using a hybrid approach, which couples a quantum chemical description of the chlorophylls with a classical model for the protein and the external environment (membrane and water). The classical model uses a polarizable Molecular Mechanics force field, thus allowing mutual polarization effects in the calculations of the excitonic properties. The investigation is performed both on the crystal structure and on structures generated by a μs long classical molecular dynamics simulation of the complex within a solvated membrane. The results show that this integrated approach not only provides a good description of the excitonic properties and optical spectra without the need for additional refinements of the excitonic parameters, but it also allows an atomistic investigation of the relative importance of electronic, structural and environment effects in determining the optical spectra.

The conformational phase space of the photoprotective switch in the major light harvesting complex II

Extensive conformational sampling of the major light harvesting complex II defines the collective variables of the photoprotective switch.