Ultrafast Excited-State Dynamics of Carotenoids and the Role of the SX State

N-Terminal domain homologs of the orange carotenoid protein increase quenching of cyanobacterial phycobilisomes

Stress exerted by excess captured light energy in cyanobacteria is prevented by the photoprotective activity of the orange carotenoid protein (OCP). Under high light, the OCP converts from an orange, inactive form (OCPO) into the red form (OCPR) that binds to and quenches the phycobilisome (PBS). Structurally, the OCP consists of 2 domains: the N-terminal effector domain and a C-terminal regulatory domain. Structural analysis of the OCP-PBS complex showed that the N-terminal domains of an OCP dimer interact with the PBS core. These N-terminal OCP domains have single-domain protein paralogs known as helical carotenoid proteins (HCPs). Using PBS quenching assays, we show that the HCP4 and HCP5 homologs efficiently quench PBS fluorescence in vitro, surpassing the quenching ability of the OCP. This is consistent with computational quantum mechanics/molecular mechanics results. Interestingly, when using a maximum quenching concentration of OCP with PBSs, HCP5 addition further increases PBS quenching. Our results provide mechanistic insight into the quenching capacity and roles of HCP4 and HCP5 in cyanobacteria, suggesting that they are more than simply functionally redundant to the OCP.

Initial Conditions for Excited-State Dynamics in Solvated Systems: A Case Study

Simulating excited-state dynamics or computing spectra for molecules in condensed phases requires sampling the ground state to generate initial conditions. Initial conditions (or snapshots for spectra) are typically produced by QM/MM Boltzmann sampling following MM equilibration or optimization. Given the switch from a MM to a QM/MM potential energy surface, one should discard a set period of time (which we call the "healing time") from the beginning of the QM/MM trajectory. Ideally, the healing time is as short as possible (to avoid unnecessary computational effort), but long enough to equilibrate to the QM/MM ground state distribution. Healing times in previous studies range from tens of femtoseconds to tens of picoseconds, suggesting the need for guidelines to choose a healing time. We examine the effect of healing time on the nonadiabatic dynamics and spectrum of a first-generation Donor-Acceptor Stenhouse Adduct in chloroform. Insufficient healing times skew the branching ratio of ground state products and alter the relaxation time for one pathway. The influence of the healing time on the absorption spectrum is less pronounced, warning that the spectrum is not a sensitive indicator for the quality of a set of initial conditions for dynamics. We demonstrate that a reasonable estimate for the healing time can be obtained by monitoring the solute temperature during the healing trajectory. We suggest that this procedure should become standard practice for determining healing times to generate initial conditions for nonadiabatic QM/MM simulations in large molecules and condensed phases.

Quantum-Classical Simulations Reveal the Photoisomerization Mechanism of a Prototypical First-Generation Molecular Motor

Light-driven molecular rotary motors convert the energy of absorbed light into unidirectional rotational motion and are key components in the design of molecular machines. The archetypal class of light-driven rotary motors is chiral overcrowded alkenes, where the rotational movement is achieved through consecutive cis-trans photoisomerization reactions and thermal helix inversion steps. While the thermal steps have been rather well understood by now, our understanding of the photoisomerization reactions of overcrowded alkene-based motors still misses key points that would explain the striking differences in operation efficiency of the known systems. Here, we employ quantum-chemical calculations and nonadiabatic molecular dynamics simulations to investigate the excited-state decay and photoisomerization mechanism in a prototypical alkene-based first-generation rotary motor. We show that the initially excited bright state undergoes an ultrafast relaxation to multiple excited-state minima separated by low energy barriers and reveal a slow picosecond-timescale decay to the ground state, which only occurs from a largely twisted dark excited-state minimum, far from any conical-intersection point. Additionally, we attribute the origin of the high yields of forward photoisomerization in our investigated motor to the favorable topography of the ground-state potential energy surface, which is controlled by the conformation of the central cyclopentene rings.

A Computational Approach to Modeling Excitation Energy Transfer and Quenching in Light-Harvesting Complexes

Light-harvesting complexes (LHCs) play a critical role in modulating energy flux within photosynthetic organisms in response to fluctuating light. Under high light conditions, they activate quenching mechanisms to mitigate photodamage. Despite their importance, the molecular mechanisms underlying these photoprotective processes remain incomplete. Herein, we present a computational protocol to model the energy pathways in the LHC, focusing specifically on the minor CP29 antenna complex of plants. We explore the factors that modulate the switch between the light-harvesting and quenched states. The protocol includes modeling the exciton Hamiltonian of the chlorophylls/lutein aggregate and calculating population dynamics using a kinetic model based on the Redfield-Förster approach. Our analysis reveals a highly tunable excited-state lifetime for the complex, that can switch between quenched and unquenched state depending on the excitation energy of the lutein, which acts as a final quencher, in accordance with recent experiments. Moreover, we observe that the s-trans lutein conformers are more likely to exhibit characteristics of the quencher.

Experimental evidence of the Sx state and fluorescence emission from the intramolecular charge transfer states in fucoxanthin

Fucoxanthin is a typical carotenoid that absorbs light in the blue region of the visible spectrum, and its detailed electronic structures remain to be clarified. It is well known that carotenoids harvest energy from sunlight and transfer it to chlorophylls (Chls) and/or bacteriochlorophylls (BChls) through its excited states as the intermediate states; however, some excited states still need evidence to be definitely confirmed. Through steady-state fluorescence emission spectroscopy and femtosecond time-resolved fluorescence up-conversion technique, we provide new evidence for the identification of the excited Sx state in fucoxanthin, a representative of carotenoids. The fluorescence emission from the intramolecular charge transfer (ICT) states was also observed and identified for the first time according to our limited survey. Our findings suggest that fucoxanthin absorbs the blue light and transfers most of the energy to BChls via Sx and ICT1 states for certain bacteria, while releasing them via the ICT1 state to protect against light-induced damage for algae.

Target Analysis Resolves the Ground and Excited State Properties from Femtosecond Stimulated Raman Spectra

Target analysis is employed to resolve the ground and excited state properties from simultaneously measured Femtosecond Stimulated Raman Spectra (FSRS) and Transient Absorption Spectra (TAS). FSRS is a three-pulse technique, involving picosecond Raman pump pulses and femtosecond visible pump and probe pulses. The TAS are needed to precisely estimate the properties of the Instrument Response Function. The prezero "coherent artifact" present during the overlap of the three pulses is described by a damped oscillation with a frequency (ω - ωn) where ωn is a ground state resonance Raman frequency. Simultaneous target analysis of the FSRS and TAS allows the complete excited state dynamics to be resolved with a time resolution better than 100 fs. The model system studied is the carotenoid lycopene in tetrahydrofuran. The lycopene dynamics show a spectral evolution with seven states, including a biphasic cooling process during the S2-S1 internal conversion, multiple S1 lifetimes, and an S* state decaying with a lifetime of 7 ps.

Raman Vibrational Signatures of Excited States of Echinenone in the Orange Carotenoid Protein (OCP) and Implications for its Photoactivation Mechanism

In this study, the vibrational characteristics of optically excited echinenone in various solvents and the Orange Carotenoid Protein (OCP) in red and orange states are systematically investigated through steady-state and time-resolved spectroscopy techniques. Time-resolved experiments, employing both Transient Absorption (TA) and Femtosecond Stimulated Raman Spectroscopy (FSRS), reveal different states in the OCP photoactivation process. The time-resolved studies indicate vibrational signatures of exited states positioned above the S1 state during the initial 140 fs of carotenoid evolution in OCP, an absence of a vibrational signature for the relaxed S1 state of echinenone in OCP, and more robust signatures of a highly excited ground state (GS) in OCP. Differences in S1 state vibration population signatures between OCP and solvents are attributed to distinct conformations of echinenone in OCP and hydrogen bonds at the keto group forming a short-lived intramolecular charge transfer (ICT) state. The vibrational dynamics of the hot GS in OCP show a more pronounced red shift of ground state CC vibration compared to echinenone in solvents, thus suggesting an unusually hot form of GS. The study proposes a hypothesis for the photoactivation mechanism of OCP, emphasizing the high level of vibrational excitation in longitudinal stretching modes as a driving force. In conclusion, the comparison of vibrational signatures reveals unique dynamics of energy dissipation in OCP, providing insights into the photoactivation mechanism and highlighting the impact of the protein environment on carotenoid behavior. The study underscores the importance of vibrational analysis in understanding the intricate processes involved in early phase OCP photoactivation.

Dynamic and Energetic Aspects of Carotenoids In-and-Around Model Lipid Membranes Revealed in Molecular Modelling

In contrast to plants, humans are unable to synthesise carotenoids and have to obtain them from diet. Carotenoids fulfil several crucial biological functions in the organism; however, due to poor solubility in water, their bioavailability from plant-based food is low. The processes of carotenoid absorption and availability in the human body have been intensively studied. The recent experimental findings concerning these processes are briefly presented in the introductory part of this review, together with a summary of such topics as carotenoid carriers, body transport and tissue delivery, to finally report on molecular-level studies of carotenoid binding by membrane receptors. The main message of the review is contained in the section describing computational investigations of carotenoid intercalation and dynamic behaviour in lipid bilayers. The relevance of these computational studies lies in showing the direct link between the microscopic behaviour of molecules and the characteristics of their macroscopic ensembles. Furthermore, studying the interactions between carotenoids and lipid bilayers, and certainly proteins, on the molecular- and atomic-level using computational methods facilitates the interpretation and explanation of their macroscopic properties and, hopefully, helps to better understand the biological functions of carotenoids.

Predicting Solvatochromism of Chromophores in Proteins through QM/MM and Machine Learning

Solvatochromism occurs in both homogeneous solvents and more complex biological environments, such as proteins. While in both cases the solvatochromic effects report on the surroundings of the chromophore, their interpretation in proteins becomes more complicated not only because of structural effects induced by the protein pocket but also because the protein environment is highly anisotropic. This is particularly evident for highly conjugated and flexible molecules such as carotenoids, whose excitation energy is strongly dependent on both the geometry and the electrostatics of the environment. Here, we introduce a machine learning (ML) strategy trained on quantum mechanics/molecular mechanics calculations of geometrical and electrochromic contributions to carotenoids' excitation energies. We employ this strategy to compare solvatochromism in protein and solvent environments. Despite the important specifities of the protein, ML models trained on solvents can faithfully predict excitation energies in the protein environment, demonstrating the robustness of the chosen descriptors.

Structural and quantum chemical basis for OCP-mediated quenching of phycobilisomes

Cyanobacteria use large antenna complexes called phycobilisomes (PBSs) for light harvesting. However, intense light triggers non-photochemical quenching, where the orange carotenoid protein (OCP) binds to PBS, dissipating excess energy as heat. The mechanism of efficiently transferring energy from phycocyanobilins in PBS to canthaxanthin in OCP remains insufficiently understood. Using cryo-electron microscopy, we unveiled the OCP-PBS complex structure at 1.6- to 2.1-angstrom resolution, showcasing its inherent flexibility. Using multiscale quantum chemistry, we disclosed the quenching mechanism. Identifying key protein residues, we clarified how canthaxanthin's transition dipole moment in its lowest-energy dark state becomes large enough for efficient energy transfer from phycocyanobilins. Our energy transfer model offers a detailed understanding of the atomic determinants of light harvesting regulation and antenna architecture in cyanobacteria.

Ultrafast excited-state dynamics of Luteins in the major light-harvesting complex LHCII

Carotenoid pigments are known to present a functional versatility when bound to light-harvesting complexes. This versatility originates from a strong correlation between a complex electronic structure and a flexible geometry that is easily tunable by the surrounding protein environment. Here, we investigated how the different L1 and L2 sites of the major trimeric light-harvesting complex (LHCII) of green plants tune the electronic structure of the two embedded luteins, and how this reflects on their ultrafast dynamics upon excitation. By combining molecular dynamics and quantum mechanics/molecular mechanics calculations, we found that the two luteins feature a different conformation around the second dihedral angle in the lumenal side. The s-cis preference of the lutein in site L2 allows for a more planar geometry of the π -conjugated backbone, which results in an increased degree of delocalization and a reduced excitation energy, explaining the experimentally observed red shift. Despite these remarkable differences, according to surface hopping simulations the two luteins present analogous ultrafast dynamics upon excitation: the bright S 2 state quickly decays (in ∼ 50 fs) to the dark intermediate S x , eventually ending up in the S 1 state. Furthermore, by employing two different theoretical approaches (i.e., Förster theory and an excitonic version of surface hopping), we investigated the experimentally debated energy transfer between the two luteins. With both approaches, no evident energy transfer was observed in the ultrafast timescale.

The nature of carotenoid S* state and its role in the nonphotochemical quenching of plants

In plants, light-harvesting complexes serve as antennas to collect and transfer the absorbed energy to reaction centers, but also regulate energy transport by dissipating the excitation energy of chlorophylls. This process, known as nonphotochemical quenching, seems to be activated by conformational changes within the light-harvesting complex, but the quenching mechanisms remain elusive. Recent spectroscopic measurements suggest the carotenoid S* dark state as the quencher of chlorophylls' excitation. By investigating lutein embedded in different conformations of CP29 (a minor antenna in plants) via nonadiabatic excited state dynamics simulations, we reveal that different conformations of the complex differently stabilize the lutein s-trans conformer with respect to the dominant s-cis one. We show that the s-trans conformer presents the spectroscopic signatures of the S* state and rationalize its ability to accept energy from the closest excited chlorophylls, providing thus a relationship between the complex's conformation and the nonphotochemical quenching.

Optical Control of Crossing the Conical Intersection in β-Carotene

Optical control of dynamic processes has been challenging yet has only been demonstrated in several chemical and biological systems. The control of a reaction passing the widely present conical intersection has not been realized. Here, we modulated the phase of the excitation pulse to control the dynamics of β-carotene through accessing the conical intersection (CI). We observed different dynamics in 110-220 fs into the CI and the consecutive process in 400-600 fs through another CI by various chirped excitation pulses. We successfully controlled those ultrafast wavepacket dynamics passing the CIs on the femtosecond time scales. The method developed here can be used to control a various of ultrafast chemical and biological reactions through the CI(s).

How orange carotenoid protein controls the excited state dynamics of canthaxanthin

Orange Carotenoid Protein (OCP) is a ketocarotenoid-binding protein essential for photoprotection in cyanobacteria. The main steps of the photoactivated conversion which converts OCP from its resting state to the active one have been extensively investigated. However, the initial photochemical event in the ketocarotenoid which triggers the large structural changes finally leading to the active state is still not understood. Here we employ QM/MM surface hopping nonadiabatic dynamics to investigate the excited-state decay of canthaxanthin in OCP, both in the ultrafast S2 to S1 internal conversion and the slower decay leading back to the ground state. For the former step we show the involvement of an additional excited state, which in the literature has been often named the SX state, and we characterize its nature. For the latter step, we reveal an excited state decay characterized by multiple timescales, which are related to the ground-state conformational heterogeneity of the ketocarotenoid. We assigned the slowly decaying population to the so-called S* state. Finally, we identify a minor decay pathway involving double-bond photoisomerization, which could be the initial trigger to photoactivation of OCP.

Understanding the Excited-State Relaxation Mechanisms of Xanthophyll Lutein by Multi-configurational Electronic Structure Calculations

The contradictory behaviors in light harvesting and non-photochemical quenching make xanthophyll lutein the most attractive functional molecule in photosynthesis. Despite several theoretical simulations on the spectral properties and excited-state dynamics, the atomic-level photophysical mechanisms need to be further studied and established, especially for an accurate description of geometric and electronic structures of conical intersections for the lowest several electronic states of lutein. In the present work, semiempirical OM2/MRCI and multi-configurational restricted active space self-consistent field methods were performed to optimize the minima and conical intersections in and between the 1Ag-, 2Ag-, 1Bu+, and 1Bu- states. Meanwhile, the relative energies were refined by MS-CASPT2(10,8)/6-31G*, which can reproduce correct electronic state properties as those in the spectroscopic experiments. Based on the above calculation results, we proposed a possible excited-state relaxation mechanism for lutein from its initially populated 1Bu+ state. Once excited to the optically bright 1Bu+ state, the system will propagate along the key reaction coordinate, i.e., the stretching vibration of the conjugated carbon chain. During this period of time, the 1Bu- state will participate in and forms a resonance state between the 1Bu- and 1Bu+ states. Later, the system will rapidly hop to the 2Ag- state via the 1Bu+/2Ag- conical intersection. Finally, the lutein molecule will survive in the 2Ag- state for a relatively long time before it internally converts to the ground state directly or via a twisted S1/S0 conical intersection. Notably, though the photophysical picture may be very different in solvents and proteins, the current theoretical study proposed a promising calculation protocol and also provided many valuable mechanistic insights for lutein and similar carotenoids.