Color tuning in photofunctional proteins

Multi-scale simulation reveals that an amino acid substitution increases photosensitizing reaction inputs in Rhodopsins

Evaluating the availability of molecular oxygen (O₂ ) and energy of excited states in the retinal binding site of rhodopsin is a crucial challenging first step to understand photosensitizing reactions in wild-type (WT) and mutant rhodopsins by absorbing visible light. In the present work, energies of the ground and excited states related to 11-cis-retinal and the O₂ accessibility to the β-ionone ring are evaluated inside WT and human M207R mutant rhodopsins. Putative O₂ pathways within rhodopsins are identified by using molecular dynamics simulations, Voronoi-diagram analysis, and implicit ligand sampling while retinal energetic properties are investigated through density functional theory, and quantum mechanical/molecular mechanical methods. Here, the predictions reveal that an amino acid substitution can lead to enough energy and O₂ accessibility in the core hosting retinal of mutant rhodopsins to favor the photosensitized singlet oxygen generation, which can be useful in understanding retinal degeneration mechanisms and in designing blue-lighting-absorbing proteic photosensitizers.

Unified Model for Photophysical and Electro-Optical Properties of Green Fluorescent Proteins

Green fluorescent proteins (GFPs) have become indispensable imaging and optogenetic tools. Their absorption and emission properties can be optimized for specific applications. Currently, no unified framework exists to comprehensively describe these photophysical properties, namely the absorption maxima, emission maxima, Stokes shifts, vibronic progressions, extinction coefficients, Stark tuning rates, and spontaneous emission rates, especially one that includes the effects of the protein environment. In this work, we study the correlations among these properties from systematically tuned GFP environmental mutants and chromophore variants. Correlation plots reveal monotonic trends, suggesting that all these properties are governed by one underlying factor dependent on the chromophore's environment. By treating the anionic GFP chromophore as a mixed-valence compound existing as a superposition of two resonance forms, we argue that this underlying factor is defined as the difference in energy between the two forms, or the driving force, which is tuned by the environment. We then introduce a Marcus-Hush model with the bond length alternation vibrational mode, treating the GFP absorption band as an intervalence charge transfer band. This model explains all of the observed strong correlations among photophysical properties; related subtopics are extensively discussed in the Supporting Information. Finally, we demonstrate the model's predictive power by utilizing the additivity of the driving force. The model described here elucidates the role of the protein environment in modulating the photophysical properties of the chromophore, providing insights and limitations for designing new GFPs with desired phenotypes. We argue that this model should also be generally applicable to both biological and nonbiological polymethine dyes.

The Dynamic Origin of Color Tuning in Proteins Revealed by a Carotenoid Pigment

Understanding the microscopic origin of the color tuning in pigment-protein complexes is a challenging yet fundamental issue in photoactive biological systems. Here, we propose a possible interpretation by using a state-of-the-art multiscale strategy based on the integration of quantum chemistry and polarizable atomistic embeddings into a dynamic description. By means of such a strategy we are able to resolve the long-standing dispute over the coloration mechanism in the crustacyanin protein. It is shown that the combination of the dynamical flexibility of the carotenoid pigments (astaxanthin) with the responsive protein environment is essential to obtain quantitative predictions of the spectral tuning. The strong linear correlation between the excitation energies and the bond length alternation in the long-chain carotenoids modulated by the dynamical protein environment is a novel finding explaining the high color tunability in crustacyanin.

Electrostatic versus Resonance Interactions in Photoreceptor Proteins: The Case of Rhodopsin

Light sensing in photoreceptor proteins is subtly modulated by the multiple interactions between the chromophoric unit and its binding pocket. Many theoretical and experimental studies have tried to uncover the fundamental origin of these interactions but reached contradictory conclusions as to whether electrostatics, polarization, or intrinsically quantum effects prevail. Here, we select rhodopsin as a prototypical photoreceptor system to reveal the molecular mechanism underlying these interactions and regulating the spectral tuning. Combining a multireference perturbation method and density functional theory with a classical but atomistic and polarizable embedding scheme, we show that accounting for electrostatics only leads to a qualitatively wrong picture, while a responsive environment can successfully capture both the classical and quantum dominant effects. Several residues are found to tune the excitation by both differentially stabilizing ground and excited states and through nonclassical "inductive resonance" interactions. The results obtained with such a quantum-in-classical model are validated against both experimental data and fully quantum calculations.

Multiscale Simulations on Spectral Tuning and the Photoisomerization Mechanism in Fluorescent RNA Spinach

Fluorescent RNA aptamer Spinach can bind and activate a green fluorescent protein (GFP)-like chromophore (an anionic DFHBDI chromophore) displaying green fluorescence. Spectroscopic properties, spectral tuning, and the photoisomerization mechanism in the Spinach-DFHBDI complex have been investigated by high-level QM and hybrid ONIOM(QM:AMBER) methods (QM method: (TD)DFT, SF-BHHLYP, SAC-CI, LT-DF-LCC2, CASSCF, or MS-CASPT2), as well as classical molecular dynamics (MD) simulations. First, our benchmark calculations have shown that TD-DFT and spin-flip (SF) TD-DFT (SF-BHHLYP) failed to give a satisfactory description of absorption and emission of the anionic DFHBDI chromophore. Comparatively, SAC-CI, LT-DF-LCC2, and MS-CASPT2 can give more reliable transition energies and are mainly used to further study the spectra of the anionic DFHBDI chromophore in Spinach. The RNA environmental effects on the spectral tuning and the photoisomerization mechanism have been elucidated. Our simulations show that interactions of the anionic cis-DFHBDI chromophore with two G-quadruplexes as well as a UAU base triple suppress photoisomerization of DFHBDI. In addition, strong hydrogen bonds between the anionic cis-DFHBDI chromophore and nearby nucleotides facilitate its binding to Spinach and further inhibit the cis-trans photoisomerization of DFHBDI. Solvent molecules, ions, and loss of key hydrogen bonds with nearby nucleotides could induce dissociation of the anionic trans-DFHBDI chromophore from the binding site. These results provide new insights into fluorescent RNA Spinach and may help rational design of other fluorescent RNAs.

Role of the Molecular Environment in Flavoprotein Color and Redox Tuning: QM Cluster versus QM/MM Modeling

We investigate the origin of the excitation energy shifts induced by the apoprotein in the active site of the bacterial photoreceptor BLUF (Blue Light sensor Using Flavin adenine dinucleotide). In order to compute the vertical excitation energies of three low-lying electronic states, including two π-π* states of flavin (S1 and S2) and a π-π* tyrosine-flavin electron-transfer state (ET), with respect to the energy of the closed-shell ground state (S0), we prepared alternative quantum mechanical (QM) cluster and quantum mechanics/molecular mechanics (QM/MM) models. We found that the excitation energies computed with both types of models correlate with the magnitude of the charge transfer character of the excitation. Accordingly, we conclude that the small charge transfer character of the light absorbing S0-S1 transition and the substantial charge transfer character of the nonabsorbing but redox active S0-ET transition explain the small color changes but substantial redox tuning in BLUF and also in other flavoproteins. Further analysis showed that redox tuning is governed by the electrostatic interaction in the QM/MM model and transfer of charge between the active site and its environment in the QM cluster. Moreover, the wave function polarization of the QM subsystem by the MM subsystem influences the magnitude of the charge transfer, resulting in the QM/MM and QM excitation energies that are not entirely consistent.

Local electric fields and molecular properties in heterogeneous environments through polarizable embedding

In spectroscopies, the local field experienced by a molecule embedded in an environment will be different from the externally applied electromagnetic field, and this difference may significantly alter the response and transition properties of the molecule. The polarizable embedding (PE) model has previously been developed to model the local field contribution stemming from the direct molecule-environment coupling of the electromagnetic response properties of molecules in solution as well as in heterogeneous environments, such as proteins. Here we present an extension of this approach to address the additional effective external field effect, i.e., the manifestations of the environment polarization induced by the external field, which allows for the calculation of properties defined in terms of the external field. Within a response framework, we report calculations of the one- and two-photon absorption (1PA and 2PA, respectively) properties of PRODAN-methanol clusters as well as the fluorescent protein DsRed. Our results demonstrate the necessity of accounting for both the dynamical reaction field and effective external field contributions to the local field in order to reproduce full quantum chemical reference calculations. For the lowest π→π* transition in DsRed, inclusion of effective external field effects gives rise to a 1.9- and 3.5-fold reduction in the 1PA and 2PA cross-sections, respectively. The effective external field is, however, strongly influenced by the heterogeneity of the protein matrix, and the resulting effect can lead to either screening or enhancement depending on the nature of the transition under consideration.

Emission shaping in fluorescent proteins: role of electrostatics and π-stacking

We obtained the fluorescence spectrum of the GFP with trajectory simulations, and revealed the role of the protein sidechains in emission shifts.

Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone

We present an implementation of analytical quantum mechanical molecular gradients within the polarizable embedding (PE) model to allow for efficient geometry optimizations and vibrational analysis of molecules embedded in large, geometrically frozen environments. We consider a variational ansatz for the quantum region, covering (multiconfigurational) self-consistent-field and Kohn-Sham density functional theory. As the first application of the implementation, we consider the internal vibrational Stark effect of the C=O group of acetophenone in different solvents and derive its vibrational linear Stark tuning rate using harmonic frequencies calculated from analytical gradients and computed local electric fields. Comparisons to PE calculations employing an enlarged quantum region as well as to a non-polarizable embedding scheme show that the inclusion of mutual polarization between acetophenone and water is essential in order to capture the structural modifications and the associated frequency shifts observed in water. For more apolar solvents, a proper description of dispersion and exchange-repulsion becomes increasingly important, and the quality of the optimized structures relies to a larger extent on the quality of the Lennard-Jones parameters.

Analysis of computational models for an accurate study of electronic excitations in GFP

The accuracy of PERI-CC2 is successfully assessed against RVS-CC2 for GFP model clusters and subsequently applied to a whole protein model.

A theoretical analysis of the potential role of π-π stacking interactions in the photoprotolytic cycle of firefly luciferin

Firefly oxyluciferin is a photoacid that presents a pH-sensitive fluorescence, which results from pH-dependent changes on the conformation of self-aggregated π-π stacking complexes. Luciferin is a derivative of oxyluciferin with very similar fluorescence and photoacidic properties. This similarity indicates that luciferin is also expected to be able to form π-π stacking complexes, but no pH-sensitive fluorescence is found for this compound. Here, a theoretical approach is used to rationalize this finding. We have found that luciferin only forms π-π stacking complexes in the ground state at acidic pH. At basic pH and in the excited state, luciferin is present as a dianion. This species is not able to self-aggregate, owing to repulsive electrostatic interactions. Thus, this emissive species is not subject to π-π stacking interactions; this explains its pH-insensitive fluorescence.

Quantum mechanical molecular interactions for calculating the excitation energy in molecular environments: a first-order interacting space approach

Intermolecular interactions regulate the molecular properties in proteins and solutions such as solvatochromic systems. Some of the interactions have to be described at an electronic-structure level. In this study, a commutator for calculating the excitation energy is used for deriving a first-order interacting space (FOIS) to describe the environmental response to solute excitation. The FOIS wave function for a solute-in-solvent cluster is solved by second-order perturbation theory. The contributions to the excitation energy are decomposed into each interaction and for each solvent.

Oxyluciferin photoacidity: the missing element for solving the keto-enol mystery?

The oxyluciferin family of fluorophores has been receiving much attention from the research community and several systematic studies have been performed in order to gain more insight regarding their photophysical properties and photoprotolytic cycles. In this minireview, we summarize the knowledge obtained so far and define several possible lines for future research. More importantly, we analyze the impact of the discoveries on the firefly bioluminescence phenomenon made so far and explain how they re-open again the discussion regarding the identity (keto or enol species) of the bioluminophore.

Rhodopsin Absorption from First Principles: Bypassing Common Pitfalls

Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.

Are the bio- and chemiluminescence states of the firefly oxyluciferin the same as the fluorescence state?

A usual strategy in both experimental and theoretical studies on bio- and chemiluminescence is to analyze the fluorescent properties of the bio- and chemiluminescence reaction product. Recent findings in a coelenteramide and Cypridina oxyluciferin model raise a concern on the validity of this procedure, showing that the light emitters in each of these luminescent processes might differ. Here, the thermal decomposition path of the firefly dioxetanone and the light emission states of the Firefly oxyluciferin responsible for the bio-, chemiluminescence, and fluorescence of the molecule are characterized using ab initio quantum chemistry and hybrid quantum chemistry/molecular mechanics methods to determine if the scenario found in the coelenteramide and Cypridina oxyluciferin study does also apply to the Firefly bioluminescent systems. The results point out to a unique emission state in the bio-, chemiluminescence, and fluorescence phenomena of the Firefly oxyluciferin and, therefore, using fluorescence properties of this system is reasonable.

All-Optical Sensing of the Components of the Internal Local Electric Field in Proteins

Here, we present a new all-optical method of interrogation of the internal electric field vector inside proteins. The method is based on experimental evaluation of the permanent dipole moment change upon excitation and the pure electronic transition frequency of a fluorophore embedded in a protein matrix. The permanent dipole moment change can be obtained from two-photon absorption measurements. In addition, permanent dipole moment change, tensor of polarizability change, and transition frequency for the free chromophore should be calculated quantum-mechanically. This allows obtaining the components of the electric field by considering the second-order Stark shift. We use the fluorescent protein mCherry as an example to demonstrate the applicability of the method.

Advances in the knowledge of light emission by firefly luciferin and oxyluciferin

Firefly luciferase is the most important and studied bioluminescence system. Due to very interesting characteristics, this system has gained numerous biomedical, pharmaceutical and bioanalytical applications, among others. In order to improve the use of this system, various researchers have tried to understand experimentally the colour of bioluminescence, and to create ways of tuning the colour emitted. The objective of this manuscript is to review the experimental studies of firefly luciferin and oxyluciferin, and related analogues, fluorescence and bioluminescence.

Interpolated mechanics-molecular mechanics study of internal rotation dynamics of the chromophore unit in blue fluorescent protein and its variants

The excited state dynamics of the blue fluorescent protein (BFP) and its variants, azurite, EBFP1.2, and EBFP2.0, are studied using molecular dynamics (MD) simulations on potential energy surfaces (PESs) generated with the interpolated mechanics-molecular mechanics (IM/MM) scheme. This IM/MM strategy adopts the interpolated PES for an important area of the complex and the conventional force field for the remaining part. We focus on the internal rotation dynamics of the chromophore unit, which is directly related to its fluorescence property, and analyze the time evolutions of the nonrotated chromophore fractions based on trajectories over 10 μs of aggregate simulation time. The characteristics obtained from the calculated time progresses of the nonrotated chromophore fractions in BFP and other variants agree well with experimentally observed properties. The results show that the MD simulation with an IM/MM potential is an attractive approach for studying excited state dynamics of fluorescent proteins in consideration of its efficiency and reliability. We also attempt to investigate the detailed roles that the mutated residues play in delaying the excited state chromophore twisting and thus improving the fluorescence property, and discuss the contributions by the Coulombic and the steric interactions between the chromophore and the mutated residues.

Firefly chemiluminescence and bioluminescence: efficient generation of excited states

Firefly luciferase catalyzes a light-emitting reaction in which an excited-state product is formed. Both experimental and theoretical methodologies are used to study this system, and the reactions catalyzed by luciferase are relatively well characterized. However, the mechanism by which an excited-state product is formed is still unknown. This Minireview deals with the current understanding of firefly bioluminescence and chemiluminescence. Thermal decomposition of simple 1,2-dioxetanes is also discussed, due to their role in formation of the excited-state bioluminophore.

Color Tuning in rhodopsins: the origin of the spectral shift between the chloride-bound and anion-free forms of halorhodopsin

Detailed knowledge of the molecular mechanisms that control the spectral properties in the rhodopsin protein family is important for understanding the functions of these photoreceptors and for the rational design of artificial photosensitive proteins. Here we used a high-level ab initio QM/MM method to investigate the mechanism of spectral tuning in the chloride-bound and anion-free forms of halorhodopsin from Natronobacterium pharaonis (phR) and the interprotein spectral shift between them. We demonstrate that the chloride ion tunes the spectral properties of phR via two distinct mechanisms: (i) electrostatic interaction with the chromophore, which results in a 95 nm difference between the absorption maxima of the two forms, and (ii) induction of a structural reorganization in the protein, which changes the positions of charged and polar residues and reduces this difference to 29 nm. The present study expands our knowledge concerning the role of the reorganization of the internal H-bond network for color tuning in general and provides a detailed investigation of the tuning mechanism in phR in particular.