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

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.

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.

Twisted Carotenoids Do Not Support Efficient Intramolecular Singlet Fission in the Orange Carotenoid Protein

Singlet exciton fission is the spin-allowed generation of two triplet electronic excited states from a singlet state. Intramolecular singlet fission has been suggested to occur on individual carotenoid molecules within protein complexes provided that the conjugated backbone is twisted out of plane. However, this hypothesis has been forwarded only in protein complexes containing multiple carotenoids and bacteriochlorophylls in close contact. To test the hypothesis on twisted carotenoids in a "minimal" one-carotenoid system, we study the orange carotenoid protein (OCP). OCP exists in two forms: in its orange form (OCPo), the single bound carotenoid is twisted, whereas in its red form (OCPr), the carotenoid is planar. To enable room-temperature spectroscopy on canthaxanthin-binding OCPo and OCPr without laser-induced photoconversion, we trap them in a trehalose glass. Using transient absorption spectroscopy, we show that there is no evidence of long-lived triplet generation through intramolecular singlet fission despite the canthaxanthin twist in OCPo.

Functional specialization of expanded orange carotenoid protein paralogs in subaerial Nostoc species

Orange carotenoid protein (OCP) is a photoactive protein that participates in the photoprotection of cyanobacteria. There are two full-length OCP proteins, four N-terminal paralogs (helical carotenoid protein, HCP), and one C-terminal domain-like carotenoid protein (CCP) found in Nostoc flagelliforme, a desert cyanobacterium. All HCPs (HCP1-3 and HCP6) from N. flagelliforme demonstrated their excellent singlet oxygen quenching activities, in which HCP2 was the strongest singlet oxygen quencher compared with others. Two OCPs, OCPx1 and OCPx2, were not involved in singlet oxygen scavenging; instead, they functioned as phycobilisome fluorescence quenchers. The fast-acting OCPx1 showed more effective photoactivation and stronger phycobilisome fluorescence quenching compared to OCPx2, which behaved differently from all reported OCP paralogs. The resolved crystal structure and mutant analysis revealed that Trp111 and Met125 play essential roles in OCPx2, which is dominant and long-acting. The resolved crystal structure of OCPx2 is maintained in a monomer state and showed more flexible regulation in energy quenching activities compared with the packed oligomer of OCPx1. The recombinant apo-CCP obtained the carotenoid pigment from holo-HCPs and holo-OCPx1 of N. flagelliforme. No such carotenoid transferring processes were observed between apo-CCP and holo-OCPx2. The close phylogenetic relationship of OCP paralogs from subaerial Nostoc species indicates an adaptive evolution toward development of photoprotection: protecting cellular metabolism against singlet oxygen damage using HCPs and against excess energy captured by active phycobilisomes using two different working modes of OCPx.

Reweighted Manifold Learning of Collective Variables from Enhanced Sampling Simulations

Enhanced sampling methods are indispensable in computational chemistry and physics, where atomistic simulations cannot exhaustively sample the high-dimensional configuration space of dynamical systems due to the sampling problem. A class of such enhanced sampling methods works by identifying a few slow degrees of freedom, termed collective variables (CVs), and enhancing the sampling along these CVs. Selecting CVs to analyze and drive the sampling is not trivial and often relies on chemical intuition. Despite routinely circumventing this issue using manifold learning to estimate CVs directly from standard simulations, such methods cannot provide mappings to a low-dimensional manifold from enhanced sampling simulations, as the geometry and density of the learned manifold are biased. Here, we address this crucial issue and provide a general reweighting framework based on anisotropic diffusion maps for manifold learning that takes into account that the learning data set is sampled from a biased probability distribution. We consider manifold learning methods based on constructing a Markov chain describing transition probabilities between high-dimensional samples. We show that our framework reverts the biasing effect, yielding CVs that correctly describe the equilibrium density. This advancement enables the construction of low-dimensional CVs using manifold learning directly from the data generated by enhanced sampling simulations. We call our framework reweighted manifold learning. We show that it can be used in many manifold learning techniques on data from both standard and enhanced sampling simulations.

Structure-function-dynamics relationships in the peculiar Planktothrix PCC7805 OCP1: Impact of his-tagging and carotenoid type

The orange carotenoid protein (OCP) is a photoactive protein involved in cyanobacterial photoprotection. Here, we report on the functional, spectral and structural characteristics of the peculiar Planktothrix PCC7805 OCP (Plankto-OCP). We show that this OCP variant is characterized by higher photoactivation and recovery rates, and a stronger energy-quenching activity, compared to other OCP studied thus far. We characterize the effect of the functionalizing carotenoid and of his-tagging on these reactions, and identify the time scales on which these modifications affect photoactivation. The presence of a his-tag at the C-terminus has a large influence on photoactivation, thermal recovery and PBS-fluorescence quenching, and likewise for the nature of the carotenoid that additionally affects the yield and characteristics of excited states and the ns-s dynamics of photoactivated OCP. By solving the structures of Plankto-OCP in the ECN- and CAN-functionalized states, each in two closely-related crystal forms, we further unveil the molecular breathing motions that animate Plankto-OCP at the monomer and dimer levels. We finally discuss the structural changes that could explain the peculiar properties of Plankto-OCP.

A synthetic switch based on orange carotenoid protein to control blue-green light responses in chloroplasts

Synthetic biology approaches to engineer light-responsive systems are widely used, but their applications in plants are still limited due to the interference with endogenous photoreceptors and the intrinsic requirement of light for photosynthesis. Cyanobacteria possess a family of soluble carotenoid-associated proteins named orange carotenoid proteins (OCPs) that, when activated by blue-green light, undergo a reversible conformational change that enables the photoprotection mechanism that occurs on the phycobilisome. Exploiting this system, we developed a chloroplast-localized synthetic photoswitch based on a protein complementation assay where two nanoluciferase fragments were fused to separate polypeptides corresponding to the OCP2 domains. Since Arabidopsis (Arabidopsis thaliana) does not possess the prosthetic group needed for the assembly of the OCP2 complex, we first implemented the carotenoid biosynthetic pathway with a bacterial β-carotene ketolase enzyme (crtW) to generate keto-carotenoid-producing plants. The photoswitch was tested and characterized in Arabidopsis protoplasts and stably transformed plants with experiments aimed to uncover its regulation by a range of light intensities, wavelengths, and its conversion dynamics. Finally, we applied the OCP-based photoswitch to control transcriptional responses in chloroplasts in response to green light illumination by fusing the two OCP fragments with the plastidial SIGMA FACTOR 2 and bacteriophage T4 anti-sigma factor AsiA. This pioneering study establishes the basis for future implementation of plastid optogenetics to regulate organelle responses upon exposure to specific light spectra.

Integration of machine learning with computational structural biology of plants

Computational structural biology of proteins has developed rapidly in recent decades with the development of new computational tools and the advancement of computing hardware. However, while these techniques have widely been used to make advancements in human medicine, these methods have seen less utilization in the plant sciences. In the last several years, machine learning methods have gained popularity in computational structural biology. These methods have enabled the development of new tools which are able to address the major challenges that have hampered the wide adoption of the computational structural biology of plants. This perspective examines the remaining challenges in computational structural biology and how the development of machine learning techniques enables more in-depth computational structural biology of plants.

A favorable path to domain separation in the orange carotenoid protein

The orange carotenoid protein (OCP) is responsible for nonphotochemical quenching (NPQ) in cyanobacteria, a defense mechanism against potentially damaging effects of excess light conditions. This soluble two-domain protein undergoes profound conformational changes upon photoactivation, involving translocation of the ketocarotenoid inside the cavity followed by domain separation. Domain separation is a critical step in the photocycle of OCP because it exposes the N-terminal domain (NTD) to perform quenching of the phycobilisomes. Many details regarding the mechanism and energetics of OCP domain separation remain unknown. In this work, we apply metadynamics to elucidate the protein rearrangements that lead to the active, domain-separated, form of OCP. We find that translocation of the ketocarotenoid canthaxanthin has a profound effect on the energetic landscape and that domain separation only becomes favorable following translocation. We further explore, characterize, and validate the free energy surface (FES) using equilibrium simulations initiated from different states on the FES. Through pathway optimization methods, we characterize the most probable path to domain separation and reveal the barriers along that pathway. We find that the free energy barriers are relatively small (<5 kcal/mol), but the overall estimated kinetic rate is consistent with experimental measurements (>1 ms). Overall, our results provide detailed information on the requirement for canthaxanthin translocation to precede domain separation and an energetically feasible pathway to dissociation.

Enhancing the Inhomogeneous Photodynamics of Canonical Bacteriophytochrome

The ability of phytochromes to act as photoswitches in plants and microorganisms depends on interactions between a bilin-like chromophore and a host protein. The interconversion occurs between the spectrally distinct red (Pr) and far-red (Pfr) conformers. This conformational change is triggered by the photoisomerization of the chromophore D-ring pyrrole. In this study, as a representative example of a phytochrome-bilin system, we consider biliverdin IXα (BV) bound to bacteriophytochrome (BphP) from Deinococcus radiodurans. In the absence of light, we use an enhanced sampling molecular dynamics (MD) method to overcome the photoisomerization energy barrier. We find that the calculated free energy (FE) barriers between essential metastable states agree with spectroscopic results. We show that the enhanced dynamics of the BV chromophore in BphP contributes to triggering nanometer-scale conformational movements that propagate by two experimentally determined signal transduction pathways. Most importantly, we describe how the metastable states enable a thermal transition known as the dark reversion between Pfr and Pr, through a previously unknown intermediate state of Pfr. We present the heterogeneity of temperature-dependent Pfr states at the atomistic level. This work paves a way toward understanding the complete mechanism of the photoisomerization of a bilin-like chromophore in phytochromes.

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.

Excited States of Xanthophylls Revisited: Toward the Simulation of Biologically Relevant Systems

Xanthophylls are a class of oxygen-containing carotenoids, which play a fundamental role in light-harvesting pigment-protein complexes and in many photoresponsive proteins. The complexity of the manifold of the electronic states and the large sensitivity to the environment still prevent a clear and coherent interpretation of their photophysics and photochemistry. In this Letter, we compare cutting-edge ab initio methods (CC3 and DMRG/NEVPT2) with time-dependent DFT and semiempirical CI (SECI) on model keto-carotenoids and show that SECI represents the right compromise between accuracy and computational cost to be applied to real xanthophylls in their biological environment. As an example, we investigate canthaxanthin in the orange carotenoid protein and show that the conical intersections between excited states and excited-ground states are mostly determined by the effective bond length alternation coordinate, which is significantly tuned by the protein through geometrical constraints and electrostatic effects.

General method to stabilize mesophilic proteins in hyperthermal water

The stability of protein structures and biological functions at normal temperature is closely linked with the universal aqueous environment of organisms. Preserving bioactivities of proteins in hyperthermia water would expand their functional capabilities beyond those in native environments. However, only a limited number of proteins derived from hyperthermophiles are thermostable at elevated temperatures. Triggered by this, here we describe a general method to stabilize mesophilic proteins in hyperthermia water. The mesophilic proteins, protected by amphiphilic polymers with multiple binding sites, maintain their secondary and tertiary structures after incubation even in boiling water. This approach, outside the conventional environment for bioactivities of mesophilic proteins, provides a general strategy to dramatically increase the T_m (melting temperature) of mesophilic proteins without any changes to amino sequences of the native proteins. Current work offers a new insight with protein stability engineering for potential application, including vaccine storage and enzyme engineering.

•

Preserving bioactivities of proteins in hyperthermia water is promising.

•

Amphiphilic polymers could protect mesophilic proteins even in boiling water.

•

Mesophilic proteins protected by amphiphilic polymers show dramatically increased T_m.

•

The method offers application prospect for vaccine storage and enzyme engineering.