The Multiple Roles of the Protein in the Photoactivation of Orange Carotenoid Protein

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.

The Orange Carotenoid Protein Triggers Cyanobacterial Photoprotection by Quenching Bilins via a Structural Switch of Its Carotenoid

Cyanobacteria were the first microorganisms that released oxygen into the atmosphere billions of years ago. To do it safely under intense sunlight, they developed strategies that prevent photooxidation in the photosynthetic membrane, by regulating the light-harvesting activity of their antenna complexes-the phycobilisomes-via the orange-carotenoid protein (OCP). This water-soluble protein interacts with the phycobilisomes and triggers nonphotochemical quenching (NPQ), a mechanism that safely dissipates overexcitation in the membrane. To date, the mechanism of action of OCP in performing NPQ is unknown. In this work, we performed ultrafast spectroscopy on a minimal NPQ system composed of the active domain of OCP bound to the phycobilisome core. The use of this system allowed us to disentangle the signal of the carotenoid from that of the bilins. Our results demonstrate that the binding to the phycobilisomes modifies the structure of the ketocarotenoid associated with OCP. We show that this molecular switch activates NPQ, by enabling excitation-energy transfer from the antenna pigments to the ketocarotenoid.

Calculating absorption and fluorescence spectra for chromophores in solution with ensemble Franck-Condon methods

Accurately modeling absorption and fluorescence spectra for molecules in solution poses a challenge due to the need to incorporate both vibronic and environmental effects, as well as the necessity of accurate excited state electronic structure calculations. Nuclear ensemble approaches capture explicit environmental effects, Franck-Condon methods capture vibronic effects, and recently introduced ensemble-Franck-Condon approaches combine the advantages of both methods. In this study, we present and analyze simulated absorption and fluorescence spectra generated with combined ensemble-Franck-Condon approaches for three chromophore-solvent systems and compare them to standard ensemble and Franck-Condon spectra, as well as to the experiment. Employing configurations obtained from ground and excited state ab initio molecular dynamics, three combined ensemble-Franck-Condon approaches are directly compared to each other to assess the accuracy and relative computational time. We find that the approach employing an average finite-temperature Franck-Condon line shape generates spectra nearly identical to the direct summation of an ensemble of Franck-Condon spectra at one-fourth of the computational cost. We analyze how the spectral simulation method, as well as the level of electronic structure theory, affects spectral line shapes and associated Stokes shifts for 7-nitrobenz-2-oxa-1,3-diazol-4-yl and Nile red in dimethyl sulfoxide and 7-methoxy coumarin-4-acetic acid in methanol. For the first time, our studies show the capability of combined ensemble-Franck-Condon methods for both absorption and fluorescence spectroscopy and provide a powerful tool for simulating linear optical spectra.

Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol

The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. The strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.

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.

The Carbonyl Group in β2 of the Carotenoid Tunes the Photocycle Kinetics in Orange Carotenoid Protein

Adaptation to rapid environmental changes is crucial for maintaining optimal photosynthetic efficiency and is ultimately key to the survival of all photosynthetic organisms. Like most of them, cyanobacteria protect their photosynthetic apparatus against rapidly increasing light intensities by nonphotochemical quenching (NPQ). In cyanobacteria, NPQ is controlled by Orange Carotenoid Protein (OCP) photocycle. OCP is the only known photoreceptor that uses carotenoid for its light activation. How carotenoid drives and controls this unique photoactivation process is still unknown. However, understanding and potentially controlling the OCP photocycle may open up new possibilities for improving photosynthetic biomass. Here we investigate the effect of the carbonyl group in the β2 ring of the carotenoid on the OCP photocycle. We report microsecond to minute OCP light activation kinetics and Arrhenius plots of the two OCP forms: Canthaxanthin-bound OCP (OCPCAN) and echinenone-bound OCP (OCPECH). The difference between the two carotenoids is the presence of a carbonyl group in the β2-ring located in the N-terminal domain of the protein. A combination of temperature-dependent spectroscopy, flash photolysis, and pump-probe transient absorption allows us to report the previously unresolved OCP intermediate associated primarily with the absorption bleach (OCPB). OCPB dominates the photokinetics in the μs to subms time range for OCPCAN and in the μs to ms range for OCPECH. We show that in OCPCAN the OCP photocycle steps are always faster than in OCPECH: from 2 to almost 20 times depending on the step. These results suggest that the presence of the carbonyl group in the β2-ring of the carotenoid accelerates the OCP photocycle.

Relationship between non-photochemical quenching efficiency and the energy transfer rate from phycobilisomes to photosystem II

The chromophorylated PBLcm domain of the ApcE linker protein in the cyanobacterial phycobilisome (PBS) serves as a bottleneck for Förster resonance energy transfer (FRET) from the PBS to the antennal chlorophyll of photosystem II (PS II) and as a redirection point for energy distribution to the orange protein ketocarotenoid (OCP), which is excitonically coupled to the PBLcm chromophore in the process of non-photochemical quenching (NPQ) under high light conditions. The involvement of PBLcm in the quenching process was first directly demonstrated by measuring steady-state fluorescence spectra of cyanobacterial cells at different stages of NPQ development. The time required to transfer energy from the PBLcm to the OCP is several times shorter than the time it takes to transfer energy from the PBLcm to the PS II, ensuring quenching efficiency. The data obtained provide an explanation for the different rates of PBS quenching in vivo and in vitro according to the half ratio of OCP/PBS in the cyanobacterial cell, which is tens of times lower than that realized for an effective NPQ process in solution.

Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response

In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.

Intrinsic tryptophan fluorescence quenching by iodine in non-canonical amino acid reveals alteration of the hydrogen bond network in the photoactive orange carotenoid protein

The Orange Carotenoid Protein (OCP) regulates cyanobacterial photosynthetic activity through photoactivation in intense light. A hydrogen bonding network involving the keto-carotenoid oxygen and Y201 and W288 residues prevents the spontaneous activation of dark-adapted OCP. To investigate the role of the hydrogen bonds in OCP photocycling, we introduced non-canonical amino acids near the keto-carotenoid, particularly iodine at the meta-position of Y201. This modification significantly increased the yield of red OCP photoproducts, albeit with a shorter lifetime. Changes in tryptophan fluorescence during photocycling influenced by the presence of iodine near W288 revealed interactions between Y201 and W288 in the absence of the carotenoid in the C-domain. We propose that upon the relaxation of red states, a ternary complex with the carotenoid is formed. Analysis of spectral signatures and interaction energies indicates that the specific iodo-tyrosine configuration enhances interactions between the carotenoid and W288.

A Quasi Time-Reversible Scheme Based on Density Matrix Extrapolation on the Grassmann Manifold for Born-Oppenheimer Molecular Dynamics

This Letter introduces the so-called Quasi Time-Reversible scheme based on Grassmann extrapolation (QTR G-Ext) of density matrices for an accurate calculation of initial guesses in Born-Oppenheimer Molecular Dynamics (BOMD) simulations. The method shows excellent results on four large molecular systems that are representative of real-life production applications, ranging from 21 to 94 atoms simulated with Kohn-Sham (KS) density functional theory surrounded with a classical environment with 6k to 16k atoms. Namely, it clearly reduces the number of self-consistent field iterations while at the same time achieving energy-conserving simulations, resulting in a considerable speed-up of BOMD simulations even when tight convergence of the KS equations is required.

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.

Color-Tuning Mechanism of the Lit Form of Orange Carotenoid Protein

Orange carotenoid protein (OCP) of photosynthetic cyanobacteria binds to ketocarotenoids noncovalently and absorbs excess light to protect the host organism from light-induced oxidative damage. Herein, we found that mutating valine 40 in the α3 helix of Gloeocapsa sp. PCC 7513 (GlOCP1) resulted in blue- or red-shifts of 6-20 nm in the absorption maxima of the lit forms. We analyzed the origins of absorption maxima shifts by integrating X-ray crystallography, homology modeling, molecular dynamics simulations, and hybrid quantum mechanics/molecular mechanics calculations. Our analysis suggested that the single residue mutations alter the polar environment surrounding the bound canthaxanthin, thereby modulating the degree of charge transfer in the photoexcited state of the chromophore. Our integrated investigations reveal the mechanism of color adaptation specific to OCPs and suggest a design principle for color-specific photoswitches.

Computational infrared and Raman spectra by hybrid QM/MM techniques: a study on molecular and catalytic material systems

Vibrational spectroscopy is one of the most well-established and important techniques for characterizing chemical systems. To aid the interpretation of experimental infrared and Raman spectra, we report on recent theoretical developments in the ChemShell computational chemistry environment for modelling vibrational signatures. The hybrid quantum mechanical and molecular mechanical approach is employed, using density functional theory for the electronic structure calculations and classical forcefields for the environment. Computational vibrational intensities at chemical active sites are reported using electrostatic and fully polarizable embedding environments to achieve more realistic vibrational signatures for materials and molecular systems, including solvated molecules, proteins, zeolites and metal oxide surfaces, providing useful insight into the effect of the chemical environment on the signatures obtained from experiment. This work has been enabled by the efficient task-farming parallelism implemented in ChemShell for high-performance computing platforms.  This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

The role of the local environment on the structural heterogeneity of carotenoid β-ionone rings

Our analysis of the X-ray crystal structure of canthaxanthin (CAN) showed that its ketolated β-ionone rings can adopt two energetically equal, but structurally distinct puckers. Quantum chemistry calculations revealed that the potential energy surface of the β-ionone ring rotation over the plane of the conjugated π-system in carotenoids depends on the pucker state of the β-ring. Considering different pucker states and β-ionone ring rotation, we found six separate local minima on the potential energy surface defining the geometry of the keto-β-ionone ring-two cis and one trans orientation for each of two pucker states. We observed a small difference in energy and no difference in relative orientation for the cis-minima, but a pronounced difference for the position of trans-minimum in alternative pucker configurations. An energetic advantage of β-ionone ring rotation from a specific pucker type can reach up to 8 kJ/mol ([Formula: see text]). In addition, we performed the simulation of linear absorption of CAN in hexane and in a unit cell of the CAN crystal. The electronic energies of [Formula: see text] transition were estimated both for the CAN monomer and in the CAN crystal. The difference between them reached [Formula: see text], which roughly corresponds to the energy gap between A and B pucker states predicted by theoretical estimations. Finally, we have discussed the importance of such effects for biological systems whose local environment determines conformational mobility, and optical/functional characteristics of carotenoid.

Sustainable Production of Pigments from Cyanobacteria

Pigments are intensely coloured compounds used in many industries to colour other materials. The demand for naturally synthesised pigments is increasing and their production can be incorporated into circular bioeconomy approaches. Natural pigments are produced by bacteria, cyanobacteria, microalgae, macroalgae, plants and animals. There is a huge unexplored biodiversity of prokaryotic cyanobacteria which are microscopic phototrophic microorganisms that have the ability to capture solar energy and CO2 and use it to synthesise a diverse range of sugars, lipids, amino acids and biochemicals including pigments. This makes them attractive for the sustainable production of a wide range of high-value products including industrial chemicals, pharmaceuticals, nutraceuticals and animal-feed supplements. The advantages of cyanobacteria production platforms include comparatively high growth rates, their ability to use freshwater, seawater or brackish water and the ability to cultivate them on non-arable land. The pigments derived from cyanobacteria and microalgae include chlorophylls, carotenoids and phycobiliproteins that have useful properties for advanced technical and commercial products. Development and optimisation of strain-specific pigment-based cultivation strategies support the development of economically feasible pigment biorefinery scenarios with enhanced pigment yields, quality and price. Thus, this chapter discusses the origin, properties, strain selection, production techniques and market opportunities of cyanobacterial pigments.

Light activation of Orange Carotenoid Protein reveals bicycle-pedal single-bond isomerization

Orange Carotenoid protein (OCP) is the only known photoreceptor which uses carotenoid for its activation. It is found exclusively in cyanobacteria, where it functions to control light-harvesting of the photosynthetic machinery. However, the photochemical reactions and structural dynamics of this unique photosensing process are not yet resolved. We present time-resolved crystal structures at second-to-minute delays under bright illumination, capturing the early photoproduct and structures of the subsequent reaction intermediates. The first stable photoproduct shows concerted isomerization of C9'-C8' and C7'-C6' single bonds in the bicycle-pedal (s-BP) manner and structural changes in the N-terminal domain with minute timescale kinetics. These are followed by a thermally-driven recovery of the s-BP isomer to the dark state carotenoid configuration. Structural changes propagate to the C-terminal domain, resulting, at later time, in the H-bond rupture of the carotenoid keto group with protein residues. Solution FTIR and UV/Vis spectroscopy support the single bond isomerization of the carotenoid in the s-BP manner and subsequent thermal structural reactions as the basis of OCP photoreception.

Electronic Structure Effects Related to the Origin of the Remarkable Near-Infrared Absorption of Blastochloris viridis' Light Harvesting 1-Reaction Center Complex

Various photosynthetic organisms have evolved to absorb light in different regions of the visible light spectrum, thus adapting to the various lighting conditions available on Earth. While most of these autotrophic organisms absorb wavelengths around the 700-800 nm region, some are capable of red-shifted absorptions above this range, but none as remarkably as Blastochloris viridis whose main absorption is observed at 1015 nm, approximately 220 nm (0.34 eV) lower in energy than their main constituent pigments, BChl-b, whose main absorption is observed at 795 nm. The structure of its light harvesting 1-reaction center was recently elucidated by cryo-EM; however, the electronic structure details behind this red-shifted absorption remain unattended. We used hybrid quantum mechanics/molecular mechanics (QM/MM) calculations to optimize one of the active centers and performed classical molecular dynamics (MD) simulations to sample conformations beyond the optimized structure. We did excited state calculations with the time-dependent density functional theory method at the CAM-B3LYP/cc-pVDZ level of theory. We reproduced the near IR absorption by sequentially modifying the number of components involved in our systems using representative structures from the calculated MD ensemble. Natural transition orbital analysis reveals the participation of the BChl-b fragments to the main transition in the native structure and the structures obtained from the QM/MM and MD simulations. H-bonding pigment-protein interactions play a role on the conformation stabilization and orientation; however, the bacteriochlorin ring conformations and the exciton delocalization are the most relevant factors to explain the red-shifting phenomenon.

Deciphering the Spectral Tuning Mechanism in Proteorhodopsin: The Dominant Role of Electrostatics Instead of Chromophore Geometry

Proteorhodopsin (PR) is a photoactive proton pump found in marine bacteria. There are two phenotypes of PR exhibiting an environmental adaptation to the ocean's depth which tunes their maximum absorption: blue‐absorbing proteorhodopsin (BPR) and green‐absorbing proteorhodopsin (GPR). This blue/green color‐shift is controlled by a glutamine to leucine substitution at position 105 which accounts for a 20 nm shift. Typically, spectral tuning in rhodopsins is rationalized by the external point charge model but the Q105L mutation is charge neutral. To study this tuning mechanism, we employed the hybrid QM/MM method with sampling from molecular dynamics. Our results reveal that the positive partial charge of glutamine near the C₁₄−C₁₅ bond of retinal shortens the effective conjugation length of the chromophore compared to the leucine residue. The derived mechanism can be applied to explain the color regulation in other retinal proteins and can serve as a guideline for rational design of spectral shifts.

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.

Grassmann Extrapolation of Density Matrices for Born-Oppenheimer Molecular Dynamics

Born-Oppenheimer molecular dynamics (BOMD) is a powerful but expensive technique. The main bottleneck in a density functional theory BOMD calculation is the solution to the Kohn-Sham (KS) equations that requires an iterative procedure that starts from a guess for the density matrix. Converged densities from previous points in the trajectory can be used to extrapolate a new guess; however, the nonlinear constraint that an idempotent density needs to satisfy makes the direct use of standard linear extrapolation techniques not possible. In this contribution, we introduce a locally bijective map between the manifold where the density is defined and its tangent space so that linear extrapolation can be performed in a vector space while, at the same time, retaining the correct physical properties of the extrapolated density using molecular descriptors. We apply the method to real-life, multiscale, polarizable QM/MM BOMD simulations, showing that sizeable performance gains can be achieved, especially when a tighter convergence to the KS equations is required.

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.

Role of hydrogen bond alternation and charge transfer states in photoactivation of the Orange Carotenoid Protein

Here, we propose a possible photoactivation mechanism of a 35-kDa blue light-triggered photoreceptor, the Orange Carotenoid Protein (OCP), suggesting that the reaction involves the transient formation of a protonated ketocarotenoid (oxocarbenium cation) state. Taking advantage of engineering an OCP variant carrying the Y201W mutation, which shows superior spectroscopic and structural properties, it is shown that the presence of Trp201 augments the impact of one critical H-bond between the ketocarotenoid and the protein. This confers an unprecedented homogeneity of the dark-adapted OCP state and substantially increases the yield of the excited photoproduct S*, which is important for the productive photocycle to proceed. A 1.37 Å crystal structure of OCP Y201W combined with femtosecond time-resolved absorption spectroscopy, kinetic analysis, and deconvolution of the spectral intermediates, as well as extensive quantum chemical calculations incorporating the effect of the local electric field, highlighted the role of charge-transfer states during OCP photoconversion.

Yaroshevich et al. present a chemical reaction mechanism of a 35-kDa blue light-triggered photoreceptor, the Orange Carotenoid Protein (OCP). They find that photoactivation critically involves the transient formation of a protonated ketocarotenoid (oxocarbenium cation) state. This study suggests the role of charge-transfer states during OCP photoconversion.

Spectral Features of Canthaxanthin in HCP2. A QM/MM Approach

The increased interest in sequencing cyanobacterial genomes has allowed the identification of new homologs to both the N-terminal domain (NTD) and C-terminal domain (CTD) of the Orange Carotenoid Protein (OCP). The N-terminal domain homologs are known as Helical Carotenoid Proteins (HCPs). Although some of these paralogs have been reported to act as singlet oxygen quenchers, their distinct functional roles remain unclear. One of these paralogs (HCP2) exclusively binds canthaxanthin (CAN) and its crystal structure has been recently characterized. Its absorption spectrum is significantly red-shifted, in comparison to the protein in solution, due to a dimerization where the two carotenoids are closely placed, favoring an electronic coupling interaction. Both the crystal and solution spectra are red-shifted by more than 50 nm when compared to canthaxanthin in solution. Using molecular dynamics (MD) and quantum mechanical/molecular mechanical (QM/MM) studies of HCP2, we aim to simulate these shifts as well as obtain insight into the environmental and coupling effects of carotenoid-protein interactions.

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

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

Uncovering the interactions driving carotenoid binding in light-harvesting complexes

Carotenoids are essential constituents of plant light-harvesting complexes (LHCs), being involved in protein stability, light harvesting, and photoprotection. Unlike chlorophylls, whose binding to LHCs is known to require coordination of the central magnesium, carotenoid binding relies on weaker intermolecular interactions (such as hydrogen bonds and van der Waals forces), whose character is far more elusive. Here we addressed the key interactions responsible for carotenoid binding to LHCs by combining molecular dynamics simulations and polarizable quantum mechanics/molecular mechanics calculations on the major LHC, LHCII. We found that carotenoid binding is mainly stabilized by van der Waals interactions with the surrounding chlorophyll macrocycles rather than by hydrogen bonds to the protein, the latter being more labile than predicted from structural data. Furthermore, the interaction network in the binding pockets is relatively insensitive to the chemical structure of the embedded carotenoid. Our results are consistent with a number of experimental data and challenge the role played by specific interactions in the assembly of pigment-protein complexes.

Multilevel Density Functional Theory

Following recent developments in multilevel embedding methods, we introduce a novel density matrix-based multilevel approach within the framework of density functional theory (DFT). In this multilevel DFT, the system is partitioned in an active and an inactive fragment, and all interactions are retained between the two parts. The decomposition of the total system is performed upon the density matrix. The orthogonality between the two parts is maintained by solving the Kohn–Sham equations in the MO basis for the active part only, while keeping the inactive density matrix frozen. This results in the reduction of computational cost. We outline the theory and implementation and discuss the differences and similarities with state-of-the-art DFT embedding methods. We present applications to aqueous solutions of methyloxirane and glycidol.

Molecular Mechanisms of Activation in the Orange Carotenoid Protein Revealed by Molecular Dynamics

Light-harvesting in photosynthesis is accompanied by photoprotective processes. In cyanobacteria, the photoprotective role is played by a specialized complex, the orange carotenoid protein, which is activated by strong blue-green light. This photoactivation involves a unique series of structural changes which terminate with an opening of the complex into two separate domains, one of which acts as a quencher for the light-harvesting complexes. Many experimental studies have tried to reveal the molecular mechanisms through which the energy absorbed by the carotenoid finally leads to the large conformational change of the complex. Here, for the first time, these mechanisms are revealed by simulating at the atomistic level the whole dynamics of the complex through an effective combination of enhanced sampling techniques. On the basis of our findings, we can conclude that the carotenoid does not act as a spring that, releasing its internal strain, induces the dissociation, as was previously proposed, but as a “latch” locking together the two domains. The photochemically triggered displacement of the carotenoid breaks this balance, allowing the complex to dissociate.

Spectral Signatures of Canthaxanthin Translocation in the Orange Carotenoid Protein

The orange carotenoid protein (OCP) is involved in the photoprotective processes in cyanobacteria via nonphotochemical quenching. Triggered by blue-green light absorption, the carotenoid chromophore undergoes translocation, displacing around 12 Å from the C-terminal domain (CTD) to the N-terminal domain (NTD). The detailed molecular rearrangements that occur within the carotenoid and the protein during this process remain largely elusive. By using a combination of molecular dynamics, well-tempered metadynamics, and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations, we were able to mimic the translocation of the carotenoid from the inactive OCP^O and obtain metastable red-shifted states in the photoactivation mechanism, replicating the λ_max values of reference experimental spectra. In addition, our simulations give insight into the structure of the red-shifted form of the inactive state of OCP.

A Many-Body, Fully Polarizable Approach to QM/MM Simulations

We present a new development in quantum mechanics/molecular mechanics (QM/MM) methods by replacing conventional MM models with data-driven many-body (MB) representations rigorously derived from high-level QM calculations. The new QM/MM approach builds on top of mutually polarizable QM/MM schemes developed for polarizable force fields with inducible dipoles and uses permutationally invariant polynomials to effectively account for quantum-mechanical contributions (e.g., exchange-repulsion and charge transfer and penetration) that are difficult to describe by classical expressions adopted by conventional MM models. Using the many-body MB-pol and MB-DFT potential energy functions for water, which include explicit two-body and three-body terms fitted to reproduce the corresponding CCSD(T) and PBE0 two-body and three-body energies for water, we demonstrate a smooth energetic transition as molecules are transferred between QM and MM regions, without the need of a transition layer. By effectively elevating the accuracy of both the MM region and the QM/MM interface to that of the QM region, the new QM/MB-MM approach achieves an accuracy comparable to that obtained with a fully QM treatment of the entire system.

Structural analysis of a new carotenoid-binding protein: the C-terminal domain homolog of the OCP

The Orange Carotenoid Protein (OCP) is a water-soluble protein that governs photoprotection in many cyanobacteria. The 35 kDa OCP is structurally and functionally modular, consisting of an N-terminal effector domain (NTD) and a C-terminal regulatory domain (CTD); a carotenoid spans the two domains. The CTD is a member of the ubiquitous Nuclear Transport Factor-2 (NTF2) superfamily (pfam02136). With the increasing availability of cyanobacterial genomes, bioinformatic analysis has revealed the existence of a new family of proteins, homologs to the CTD, the C-terminal domain-like carotenoid proteins (CCPs). Here we purify holo-CCP2 directly from cyanobacteria and establish that it natively binds canthaxanthin (CAN). We use small-angle X-ray scattering (SAXS) to characterize the structure of this carotenoprotein in two distinct oligomeric states. A single carotenoid molecule spans the two CCPs in the dimer. Our analysis with X-ray footprinting-mass spectrometry (XFMS) identifies critical residues for carotenoid binding that likely contribute to the extreme red shift (ca. 80 nm) of the absorption maximum of the carotenoid bound by the CCP2 dimer and a further 10 nm shift in the tetramer form. These data provide the first structural description of carotenoid binding by a protein consisting of only an NTF2 domain.

Probing of carotenoid-tryptophan hydrogen bonding dynamics in the single-tryptophan photoactive Orange Carotenoid Protein

The photoactive Orange Carotenoid Protein (OCP) plays a key role in cyanobacterial photoprotection. In OCP, a single non-covalently bound keto-carotenoid molecule acts as a light intensity sensor, while the protein is responsible for forming molecular contacts with the light-harvesting antenna, the fluorescence of which is quenched by OCP. Activation of this physiological interaction requires signal transduction from the photoexcited carotenoid to the protein matrix. Recent works revealed an asynchrony between conformational transitions of the carotenoid and the protein. Intrinsic tryptophan (Trp) fluorescence has provided valuable information about the protein part of OCP during its photocycle. However, wild-type OCP contains five Trp residues, which makes extraction of site-specific information impossible. In this work, we overcame this problem by characterizing the photocycle of a fully photoactive OCP variant (OCP-3FH) with only the most critical tryptophan residue (Trp-288) in place. Trp-288 is of special interest because it forms a hydrogen bond to the carotenoid's keto-oxygen to keep OCP in its dark-adapted state. Using femtosecond pump-probe fluorescence spectroscopy we analyzed the photocycle of OCP-3FH and determined the formation rate of the very first intermediate suggesting that generation of the recently discovered S* state of the carotenoid in OCP precedes the breakage of the hydrogen bonds. Therefore, following Trp fluorescence of the unique photoactive OCP-3FH variant, we identified the rate of the H-bond breakage and provided novel insights into early events accompanying photoactivation of wild-type OCP.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.

Photophysics of deinoxanthin, the keto-carotenoid bound to the main S-layer unit of Deinococcus radiodurans

An ultrafast transient absorption experiment on the SDBC, which binds the carotenoid deinoxanthin, reveals a non-specific binding site that loosely binds the carotenoid, but protects the carotenoid from the outer environment.