The Chromophore Structure of the Cyanobacterial Phytochrome Cph1 As Predicted by Time-Dependent Density Functional Theory

Effect of the protein environment on the excited state phenomena in a bacteriophytochrome

The excited state processes of a bacteriophytochrome are studied using high-level multireference methods. The various non-radiative channels of deactivation are identified for the chromophore. The effects of the protein environment and substituents are elucidated for these excited state processes. It is observed that while the excited states are completely delocalized in the Franck-Condon (FC) region, they acquire significant charge transfer character near the conical intersections. Earlier studies have emphasized the delocalized nature of the excited states in the FC region, which leads to absorption spectra with minimal Stokes shift [Rumyantsev et al., Sci. Rep., 2015, 5, 18348]. The effect of the protein environment on the vertical excitation energies was minimal, while that on the conical intersection (CI) energetics was significant. This may lead one to believe that it is charge transfer driven. However, energy decomposition analysis shows that it is the effect of the dispersion of nearby residues and the steric effect on the rings and substituents that lead to the large effect of proteins on the energetics of the CIs.

Conformational Effects on the Absorption Spectra of Phytochromes

Bacteriophytochrome is a photoreceptor protein that contains the biliverdin (BV) chromophore as its active component. The spectra of BV upon mutation remain remarkably unchanged, as far as spectral positions are concerned. This points toward the minimal effect of electrostatic effects on the electronic structure of the chromophore. However, the relative intensities of the Q and Soret bands of the chromophore change dramatically upon mutation. In this work, we delve into the molecular origin of this unusual intensity modulation. Using extensive classical MD and QM/MM calculations, we show that due to mutation, the conformational population of the chromophore changes significantly. The noncovalent interactions, especially the stacking interactions, lead to extra stabilization of the cyclic form in the D207H mutated species as opposed to the open form in the wild-type BV. Thus, unlike the commonly observed direct electrostatic effect on the spectral shift, in the case of BV the difference observed is in varying intensities, and this in turn is driven by a conformational shift due to enhanced stacking interaction.

Estimation of the relative contributions to the electronic energy transfer rates based on Förster theory: The case of C-phycocyanin chromophores

In the present study, we provide a reformulation of the theory originally proposed by Förster which allows for simple and convenient formulas useful to estimate the relative contributions of transition dipole moments of a donor and acceptor (chemical factors), their orientation factors (intermolecular structural factors), intermolecular center-to-center distances (intermolecular structural factors), spectral overlaps of absorption and emission spectra (photophysical factors), and refractive index (material factor) to the excitation energy transfer (EET) rate constant. To benchmark their validity, we focused on the EET occurring in C-phycocyanin (C-PC) chromophores. To this aim, we resorted to quantum chemistry calculations to get optimized molecular structures of the C-PC chromophores within the density functional theory (DFT) framework. The absorption and emission spectra, as well as transition dipole moments, were computed by using the time-dependent DFT (TDDFT). Our method was applied to several types of C-PCs showing that the EET rates are determined by an interplay of their specific physical, chemical, and geometrical features. These results show that our formulas can become a useful tool for a reliable estimation of the relative contributions of the factors regulating the EET transfer rate.

Hydrogen-Bond Network Determines the Early Photoisomerization Processes of Cph1 and AnPixJ Phytochromes

Phytochrome proteins are light receptors that play a pivotal role in regulating the life cycles of plants and microorganisms. Intriguingly, while cyanobacterial phytochrome Cph1 and cyanobacteriochrome AnPixJ use the same phycocyanobilin (PCB) chromophore to absorb light, their excited‐state behavior is very different. We employ multiscale calculations to rationalize the different early photoisomerization mechanisms of PCB in Cph1 and AnPixJ. We found that their electronic S₁, T₁, and S₀ potential minima exhibit distinct geometric and electronic structures due to different hydrogen bond networks with the protein environment. These specific interactions influence the S₁ electronic structures along the photoisomerization paths, ultimately leading to internal conversion in Cph1 but intersystem crossing in AnPixJ. This explains why the excited‐state relaxation in AnPixJ is much slower (ca. 100 ns) than in Cph1 (ca. 30 ps). Further, we predict that efficient internal conversion in AnPixJ can be achieved upon protonating the carboxylic group that interacts with PCB.

The structural changes in the signaling mechanism of bacteriophytochromes in solution revealed by a multiscale computational investigation

Phytochromes are red-light sensing proteins, with important light-regulatory roles in different organisms, which are capturing an increasing interest in bioimaging and optogenetics. Upon absorption of light by the embedded bilin chromophore, they undergo structural changes that extend from the chromophore to the protein and finally drive the biological function. Up to now, the underlying mechanism still has to be characterized fully. Here we investigate the Pfr activated form of a bacterial phytochrome, by combining extensive molecular dynamics simulations with a polarizable QM/MM description of the spectroscopic properties, revealing a large structure relaxation in solution, compared to the crystal structure, both in the chromophore-binding pocket and in the overall structure of the phytochrome. Our results indicate that the final opening of the dimeric structure is preceded by an important internal reorganization of the phytochrome specific (PHY) domain involving a bend of the helical spine connecting the PHY domain with the chromophore-binding domain, opening the way to a new understanding of the activation pathway.

Modeling photophysical properties of the bacteriophytochrome-based fluorescent protein IFP1.4

An enhanced interest in the phytochrome-based fluorescent proteins is explained by their ability to absorb and emit light in the far-red and infra-red regions particularly suitable for bioimaging. The fluorescent protein IFP1.4 was engineered from the chromophore-binding domain of a bacteriophytochrome in attempts to increase the fluorescence quantum yield. We report the results of simulations of structures in the ground S₀ and excited S₁ electronic states of IFP1.4 using the methods of quantum chemistry and quantum mechanics/molecular mechanics. We construct different protonation states of the biliverdin (BV) chromophore in the red-absorbing form of the protein by moving protons from the BV pyrrole rings to a suitable acceptor within the system and show that these structures are close in energy but differ by absorption bands. For the first time, we report structures of the minimum energy conical intersection points S₁/S₀ on the energy surfaces of BV in the protein environment and describe their connection to the local minima in the excited S₁ state. These simulations allow us to characterize the deactivation routes in IFP1.4.

Carbon Atoms Speaking Out: How the Geometric Sensitivity of 13C Chemical Shifts Leads to Understanding the Colour Tuning of Phycocyanobilin in Cph1 and AnPixJ

We present a combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics-statistical approach for the interpretation of nuclear magnetic resonance (NMR) chemical shift patterns in phycocyanobilin (PCB). These were originally associated with colour tuning upon photoproduct formation in red/green-absorbing cyanobacteriochrome AnPixJg2 and red/far-red-absorbing phytochrome Cph1Δ2. We pursue an indirect approach without computation of the absorption frequencies since the molecular geometry of cofactor and protein are not accurately known. Instead, we resort to a heuristic determination of the conjugation length in PCB through the experimental NMR chemical shift patterns, supported by quantum chemical calculations. We have found a characteristic correlation pattern of 13C chemical shifts to specific bond orders within the π-conjugated system, which rests on the relative position of carbon atoms with respect to electron-withdrawing groups and the polarisation of covalent bonds. We propose the inversion of this regioselective relationship using multivariate statistics and to apply it to the known experimental NMR chemical shifts in order to predict changes in the bond alternation pattern. Therefrom the extent of electronic conjugation, and eventually the change in absorption frequency, can be derived. In the process, the consultation of explicit mesomeric formulae plays an important role to qualitatively account for possible conjugation scenarios of the chromophore. While we are able to consistently associate the NMR chemical shifts with hypsochromic and bathochromic shifts in the Pg and Pfr, our approach represents an alternative method to increase the explanatory power of NMR spectroscopic data in proteins.

The Keto-Enol Tautomerism of Biliverdin in Bacteriophytochrome: Could it Explain the Bathochromic Shift in the Pfr Form?†

Phytochromes are ubiquitous photoreceptors found in plants, eukaryotic algae, bacteria and fungi. Particularly, when bacteriophytochrome is irradiated with light, a Z-to-E (photo)isomerization takes place in the biliverdin chromophore as part of the Pr-to-Pfr conversion. This photoisomerization is concomitant with a bathochromic shift in the Q-band. Based on experimental evidence, we studied a possible keto-enol tautomerization of BV, as an alternative reaction channel after its photoisomerization. In this contribution, the noncatalyzed and water-assisted reaction pathways for the lactam-lactim interconversion through consecutive keto-enol tautomerization of a model BV species were studied deeply. It was found that in the absence of water molecules, the proton transfer reaction is unable to take place at normal conditions, due to large activation energies, and the endothermic formation of lactim derivatives prevents its occurrence. However, when a water molecule assists the process by catalyzing the proton transfer reaction, the activation free energy lowers considerably. The drastic lowering in the activation energy for the keto-enol tautomerism is due to the stabilization of the water moiety through hydrogen bonds along the reaction coordinate. The absorption spectra were computed for all tautomers. It was found that the UV-visible absorption bands are in reasonable agreement with the experimental data. Our results suggest that although the keto-enol equilibrium is likely favoring the lactam tautomer, the equilibrium could eventually be shifted in favor of the lactim, as it has been reported to occur in the dark reversion mechanism of bathy phytochromes.

Elucidating the role of structural fluctuations, and intermolecular and vibronic interactions in the spectroscopic response of a bacteriophytochrome

Molecular dynamics and a multiscale polarizable QM/MM strategy allow reproducing absorption, circular dichroism, and resonance Raman spectra of a bacteriophytochrome.

QM/MM Benchmarking of Cyanobacteriochrome Slr1393g3 Absorption Spectra

Cyanobacteriochromes are compact and spectrally diverse photoreceptor proteins that are promising candidates for biotechnological applications. Computational studies can contribute to an understanding at a molecular level of their wide spectral tuning and diversity. In this contribution, we benchmark methods to model a 110 nm shift in the UV/Vis absorption spectrum from a red- to a green-absorbing form of the cyanobacteriochrome Slr1393g3. Based on an assessment of semiempirical methods to describe the chromophore geometries of both forms in vacuo, we find that DFTB2+D leads to structures that are the closest to the reference method. The benchmark of the excited state calculations is based on snapshots from quantum mechanics/molecular mechanics molecular dynamics simulations. In our case, the methods RI-ADC(2) and sTD-DFT based on CAM-B3LYP ground state calculations perform the best, whereas no functional can be recommended to simulate the absorption spectra of both forms with time-dependent density functional theory. Furthermore, the difference in absorption for the lowest energy absorption maxima of both forms can already be modelled with optimized structures, but sampling is required to improve the shape of the absorption bands of both forms, in particular for the second band. This benchmark study can guide further computational studies, as it assesses essential components of a protocol to model the spectral tuning of both cyanobacteriochromes and the related phytochromes.

The Effective Conjugation Length Is Responsible for the Red/Green Spectral Tuning in the Cyanobacteriochrome Slr1393g3

The origin of the spectral shift from a red- to a green-absorbing form in a cyanobacteriochrome, Slr1393g3, was identified by combined quantum mechanics/molecular mechanics simulations. This protein, related to classical phytochromes, carries the open-chain tetrapyrrole chromophore phycocyanobilin. Our calculations reveal that the effective conjugation length in the chromophore becomes shorter upon conversion from the red to the green form. This is related to the planarity of the entire chromophore. A large distortion was found for the terminal pyrrole rings A and D; however, the D ring contributes more strongly to the photoproduct tuning, despite a larger change in the twist of the A ring. Our findings implicate that the D ring twist can be exploited to regulate the absorption of the photoproduct. Hence, mutations that affect the D ring twist can lead to rational tuning of the photoproduct absorption, allowing the tailoring of cyanobacteriochromes for biotechnological applications such as optogenetics and bioimaging.

Structural heterogeneity in a parent ground-state structure of AnPixJg2 revealed by theory and spectroscopy

We investigated the red absorbing, dark stable state (Pr state) of the second GAF domain of the cyanobacteriochrome AnPixJ (AnPixJg2) by a molecular dynamics simulation of 1 μs duration. Our results reveal two distinct conformational isoforms of the chromophore, from which only one was known from crystallographic experiments. The interconversion between both isoforms is accompanied by alterations in the hydrogen bond pattern between the chromophore and the protein and the solvation structure of the chromophore binding pocket. The existence of sub-states in the Pr form of AnPixJg2 is supported by the results from experimental ¹³C MAS NMR spectroscopy. Our finding is consistent with the observation of structural heterogeneity in other cyanobacteriochromes and phytochromes.

Active and silent chromophore isoforms for phytochrome Pr photoisomerization: An alternative evolutionary strategy to optimize photoreaction quantum yields

Photoisomerization of a protein bound chromophore is the basis of light sensing of many photoreceptors. We tracked Z-to-E photoisomerization of Cph1 phytochrome chromophore PCB in the Pr form in real-time. Two different phycocyanobilin (PCB) ground state geometries with different ring D orientations have been identified. The pre-twisted and hydrogen bonded PCB(a) geometry exhibits a time constant of 30 ps and a quantum yield of photoproduct formation of 29%, about six times slower and ten times higher than that for the non-hydrogen bonded PCB(b) geometry. This new mechanism of pre-twisting the chromophore by protein-cofactor interaction optimizes yields of slow photoreactions and provides a scaffold for photoreceptor engineering.

Phycocyanobilin in solution – a solvent triggered molecular switch

The chromophore phycocyanobilin changes its spectroscopic behaviour upon solvent change. Our calculations trace this effect back to conformational switching, induced by changes in the hydrogen bonding pattern.

Conformational homogeneity and excited-state isomerization dynamics of the bilin chromophore in phytochrome Cph1 from resonance Raman intensities

The ground-state structure and excited-state isomerization dynamics of the P(r) and P(fr) forms of phytochrome Cph1 are investigated using resonance Raman intensity analysis. Electronic absorption and stimulated resonance Raman spectra of P(r) and P(fr) are presented; vibronic analysis of the Raman intensities and absorption spectra reveals that both conformers exist as a single, homogeneous population of molecules in the ground state. The homogeneous and inhomogeneous contributions to the overall electronic broadening are determined, and it is found that the broadening is largely homogeneous in nature, pointing to fast excited-state decay. Franck-Condon displacements derived from the Raman intensity analysis reveal the initial atomic motions in the excited state, including the highly displaced, nontotally symmetric torsional and C(15)-H HOOP modes that appear because of symmetry-reducing distortions about the C(14)-C(15) and C(15)=C(16) bonds. P(fr) is especially well primed for ultrafast isomerization and torsional Franck-Condon analysis predicts a <200 fs P(fr) → P(r) isomerization. This time is significantly faster than the observed 700 fs reaction time, indicating that the P(fr) S(1) surface has a D-ring rotational barrier caused by steric interactions with the protein.

Homogeneity of phytochrome Cph1 vibronic absorption revealed by resonance Raman intensity analysis

Phytochromes are an important class of red/far-red responsive photoreceptors that act as light-activated biological switches, ultimately driving growth and development in plants, bacteria, and fungi. The composition of the red-absorbing ground-state has been widely debated due to the presence of a shoulder feature on the blue edge of electronic absorption spectra, which many have attributed to the presence of multiple ground-state conformers. Here we use resonance Raman intensity analysis to calculate the vibronic absorption profile of cyanobacterial phytochrome Cph1 and show that this shoulder feature is due simply to vibronic transitions from a single species, thus reflecting a homogeneous ground-state population.

A switchable light-input, light-output system modelled and constructed in yeast

Background

Advances in synthetic biology will require spatio-temporal regulation of biological processes in heterologous host cells. We develop a light-switchable, two-hybrid interaction in yeast, based upon the Arabidopsis proteins PHYTOCHROME A and FAR-RED ELONGATED HYPOCOTYL 1-LIKE. Light input to this regulatory module allows dynamic control of a light-emitting LUCIFERASE reporter gene, which we detect by real-time imaging of yeast colonies on solid media.

Results

The reversible activation of the phytochrome by red light, and its inactivation by far-red light, is retained. We use this quantitative readout to construct a mathematical model that matches the system's behaviour and predicts the molecular targets for future manipulation.

Conclusion

Our model, methods and materials together constitute a novel system for a eukaryotic host with the potential to convert a dynamic pattern of light input into a predictable gene expression response. This system could be applied for the regulation of genetic networks - both known and synthetic.