Macrocycle ring deformation as the secondary design principle for light-harvesting complexes

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.

Theory of 2D electronic spectroscopy of water soluble chlorophyll-binding protein (WSCP): Signatures of Chl b derivate

We investigate how electronic excitations and subsequent dissipative dynamics in the water soluble chlorophyll-binding protein (WSCP) are connected to features in two-dimensional (2D) electronic spectra, thereby comparing results from our theoretical approach with experimental data from the literature. Our calculations rely on third-order response functions, which we derived from a second-order cumulant expansion of the dissipative dynamics involving the partial ordering prescription, assuming a fast vibrational relaxation in the potential energy surfaces of excitons. Depending on whether the WSCP complex containing a tetrameric arrangement of pigments composed of two dimers with weak excitonic coupling between them binds the chlorophyll variant Chl a or Chl b, the resulting linear absorption and circular dichroism spectra and particularly the 2D spectra exhibit substantial differences in line shapes. These differences between Chl a WSCP and Chl b WSCP cannot be explained by the slightly modified excitonic couplings within the two variants. In the case of Chl a WSCP, the assumption of equivalent dimer subunits facilitates a reproduction of substantial features from the experiment by the calculations. In contrast, for Chl b WSCP, we have to assume that the sample, in addition to Chl b dimers, contains a small but distinct fraction of chemically modified Chl b pigments. The existence of such Chl b derivates has been proposed by Pieper et al. [J. Phys. Chem. B 115, 4042 (2011)] based on low-temperature absorption and hole-burning spectroscopy. Here, we provide independent evidence.

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.

A Multiconfigurational Wave Function Implementation of the Frenkel Exciton Model for Molecular Aggregates

We present an implementation of the Frenkel exciton model into the OpenMolcas program package enabling calculations of collective electronic excited states of molecular aggregates based on a multiconfigurational wave function description of the individual monomers. The computational protocol avoids using diabatization schemes and, thus, supermolecule calculations. Additionally, the use of the Cholesky decomposition of the two-electron integrals entering pair interactions enhances the efficiency of the computational scheme. The application of the method is exemplified for two test systems, that is, a formaldehyde oxime and a bacteriochlorophyll-like dimer. For the sake of comparison with the dipole approximation, we restrict our considerations to situations where intermonomer exchange can be neglected. The protocol is expected to be beneficial for aggregates composed of molecules with extended π systems, unpaired electrons such as radicals or transition metal centers, where it should outperform widely used methods based on time-dependent density functional theory.

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.

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.

Hysteretic Pressure Dependence of Ca2+ Binding in LH1 Bacterial Membrane Chromoproteins

Much of the thermodynamic parameter values that support life are set by the properties of proteins. While the denaturing effects of pressure and temperature on proteins are well documented, their precise structural nature is rarely revealed. This work investigates the destabilization of multiple Ca²⁺ binding sites in the cyclic LH1 light-harvesting membrane chromoprotein complexes from two Ca-containing sulfur purple bacteria by hydrostatic high-pressure perturbation spectroscopy. The native (Ca-saturated) and denatured (Ca-depleted) phases of these complexes are well distinguishable by much-shifted bacteriochlorophyll a exciton absorption bands serving as innate optical probes in this study. The pressure-induced denaturation of the complexes related to the failure of the protein Ca-binding pockets and the concomitant breakage of hydrogen bonds between the pigment chromophores and protein environment were found cooperative, involving all or most of the Ca²⁺ binding sites, but irreversible. The strong hysteresis observed in the spectral and kinetic characteristics of phase transitions along the compression and decompression pathways implies asymmetry in the relevant free energy landscapes and activation free energy distributions. A phase transition pressure equal to about 1.9 kbar was evaluated for the complexes from Thiorhodovibrio strain 970 from the pressure dependence of biphasic kinetics observed in the minutes to 100 h time range.

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.

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.

Two-dimensional spectroscopy beyond the perturbative limit: The influence of finite pulses and detection modes

Ultra-fast and multi-dimensional spectroscopy gives a powerful looking glass into the dynamics of molecular systems. In particular, two-dimensional electronic spectroscopy (2DES) provides a probe of coherence and the flow of energy within quantum systems, which is not possible with more conventional techniques. While heterodyne-detected (HD) 2DES is increasingly common, more recently fluorescence-detected (FD) 2DES offers new opportunities, including single-molecule experiments. However, in both techniques, it can be difficult to unambiguously identify the pathways that dominate the signal. Therefore, the use of numerically modeling of 2DES is vitally important, which, in turn, requires approximating the pulsing scheme to some degree. Here, we employ non-perturbative time evolution to investigate the effects of finite pulse width and amplitude on 2DES signals. In doing so, we identify key differences in the response of HD and FD detection schemes, as well as the regions of parameter space where the signal is obscured by unwanted artifacts in either technique. Mapping out parameter space in this way provides a guide to choosing experimental conditions and also shows in which limits the usual theoretical approximations work well and in which limits more sophisticated approaches are required.

Spectral tuning of chlorophylls in proteins - electrostatics vs. ring deformation

In photosynthetic complexes, tuning of chlorophyll light-absorption spectra by the protein environment is crucial to their efficiency and robustness. Recombinant type II water soluble chlorophyll-binding proteins from Brassicaceae (WSCPs) are useful for studying spectral tuning mechanisms due to their symmetric homotetramer structure, and the ability to rigorously modify the chlorophyll's protein surroundings. Our previous comparison of the crystal structures of two WSCP homologues suggested that protein-induced chlorophyll ring deformation is the predominant spectral tuning mechanism. Here, we implement a more rigorous analysis based on hybrid quantum mechanics and molecular mechanics calculations to quantify the relative contributions of geometrical and electrostatic factors to the absorption spectra of WSCP-chlorophyll complexes. We show that when considering conformational dynamics, geometry distortions such as chlorophyll ring deformation accounts for about one-third of the spectral shift, whereas the direct polarization of the electron density accounts for the remaining two-thirds. From a practical perspective, protein electrostatics is easier to manipulate than chlorophyll conformations, thus, it may be more readily implemented in designing artificial protein-chlorophyll complexes.

Artificial Photosynthesis: Is Computation Ready for the Challenge Ahead?

A tremendous effort is currently devoted to the generation of novel hybrid materials with enhanced electronic properties for the creation of artificial photosynthetic systems. This compelling and challenging problem is well-defined from an experimental point of view, as the design of such materials relies on combining organic materials or metals with biological systems like light harvesting and redox-active proteins. Such hybrid systems can be used, e.g., as bio-sensors, bio-fuel cells, biohybrid photoelectrochemical cells, and nanostructured photoelectronic devices. Despite these efforts, the main bottleneck is the formation of efficient interfaces between the biological and the organic/metal counterparts for efficient electron transfer (ET). It is within this aspect that computation can make the difference and improve the current understanding of the mechanisms underneath the interface formation and the charge transfer efficiency. Yet, the systems considered (i.e., light harvesting protein, self-assembly monolayer and surface assembly) are more and more complex, reaching (and often passing) the limit of current computation power. In this review, recent developments in computational methods for studying complex interfaces for artificial photosynthesis will be provided and selected cases discussed, to assess the inherent ability of computation to leave a mark in this field of research.

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.

Color tuning of chlorophyll a and b pigments revealed from gas-phase spectroscopy

Chlorophyll (Chl) pigments are responsible for vital mechanisms in photosynthetic proteins: light harvesting, energy transfer and charge separation. A complex interplay between the Chl molecule and its microenvironment determines its transition energy. Interactions such as excitonic coupling with one or more pigments (Chls or carotenoids), axial ligation to the magnesium center, or electrostatic interactions between Chl and nearby amino-acid residues all influence the photophysical properties. Here we use time-resolved photodissociation action spectroscopy to determine transition energies of Chla/b complexes in vacuo to directly compare the impact of a negatively charged axial ligand (formate) to that of exciton coupling between two Chls. Experiments carried out at the electrostatic ion storage ring ELISA allow dissociation to be sampled on hundreds of milliseconds time scale. Absorption-band maxima of Chla-formate complexes are found at 433 ± 4 nm/2.86 ± 0.03 eV (Soret band) and in the region 654-675 nm/1.84-1.90 eV (Q band) and those of Chla dimers tagged by a quaternary ammonium ion at 419 ± 5 nm/2.96 ± 0.04 eV (Soret band) and 647 nm/1.92 eV (Q band). The axial ligand strongly affects the Chla transition energies causing redshifts of 0.21 eV of the Soret band and 0.04-0.1 eV of the Q band compared to Chla tagged by a quaternary ammonium. Slightly smaller shifts were found in case of Chlb. The redshifts are approximately twice that induced by excitonic coupling between two Chlas, also tagged by a quaternary ammonium ion. Axial ligation brings the absorption by isolated Chls very close to that of photosynthetic proteins.

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

How the Protein Environment Can Tune the Energy, the Coupling, and the Ultrafast Dynamics of Interacting Chlorophylls: The Example of the Water-Soluble Chlorophyll Protein

The interplay between active molecules and the protein environment in light-harvesting complexes tunes the photophysics and the dynamical properties of pigment-protein complexes in a subtle way, which is not fully understood. Here we characterized the photophysics and the ultrafast dynamics of four variants of the water-soluble chlorophyll protein (WSCP) as an ideal model system to study the behavior of strongly interacting chlorophylls. We found that when coordinated by the WSCP protein, the presence of the formyl group in chlorophyll b replacing the methyl group in chlorophyll a strongly affects the exciton energy and the dynamics of the system, opening up the possibility of tuning the photophysics and the transport properties of multichromophores by engineering specific interactions with the surroundings.

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.

Optical excitations of chlorophyll a and b monomers and dimers

A necessary first step in the development of technologies such as artificial photosynthesis is understanding the photoexcitation process within the basic building blocks of naturally occurring light harvesting complexes (LHCs). The most important of these building blocks in biological LHCs such as LHC II from green plants are the chlorophyll a (Chl a) and chlorophyll b (Chl b) chromophores dispersed throughout the protein matrix. However, efforts to describe such systems are still hampered by the lack of computationally efficient and accurate methods that are able to describe optical absorption in large biomolecules. In this work, we employ a highly efficient linear combination of atomic orbitals (LCAOs) to represent the Kohn-Sham (KS) wave functions at the density functional theory (DFT) level and perform time-dependent density functional theory (TDDFT) calculations in either the reciprocal space and frequency domain (LCAO-TDDFT-k-ω) or real space and time domain (LCAO-TDDFT-r-t) of the optical absorption spectra of Chl a and b monomers and dimers. We find that our LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t calculations reproduce results obtained with a plane-wave (PW) representation of the KS wave functions (PW-TDDFT-k-ω) but with a significant reduction in computational effort. Moreover, by applying the Gritsenko, van Leeuwen, van Lenthe, and Baerends solid and correlation derivative discontinuity correction Δ_x to the KS eigenenergies, with both LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t methods, we are able to semiquantitatively reproduce the experimentally measured photoinduced dissociation results. This work opens the path to first principles calculations of optical excitations in macromolecular systems.

Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation

Light-harvesting is a crucial step of photosynthesis. Its mechanisms and related energetics have been revealed by a combination of experimental investigations and theoretical modeling. The success of theoretical modeling is largely due to the application of atomistic descriptions combining quantum chemistry, classical models and molecular dynamics techniques. Besides the important achievements obtained so far, a complete and quantitative understanding of how the many different light-harvesting complexes exploit their structural specificity is still missing. Moreover, many questions remain unanswered regarding the mechanisms through which light-harvesting is regulated in response to variable light conditions. Here we show that, in both fields, a major role will be played once more by atomistic descriptions, possibly generalized to tackle the numerous time and space scales on which the regulation takes place: going from the ultrafast electronic excitation of the multichromophoric aggregate, through the subsequent conformational changes in the embedding protein, up to the interaction between proteins.