On the primary event of phytochrome: quantum chemical comparison of photoreactions at C4, C10 and C15

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.

Photoisomerization of phytochrome chromophore models: an XMS-CASPT2 study

Phytochromes are a superfamily of photoreceptors that harbor linear tetrapyrroles as chromophores. Upon light illumination, the linear tetrapyrrole chromophore undergoes a double bond isomerization which starts a photocycle. In this work, we studied the photoisomerization of chromophore models designed based on the C- and D-rings of the phycocyanobilin (PCB) chromophore. In total, five different models with varying substitutions were investigated. Firstly, the vertical excitation energies were benchmarked using different computational methods to establish the relative order of the excited states. Based on these calculations, we computed the photoisomerization profiles using the extended multi-state (XMS) version of the CASPT2 method. The profiles were obtained for both the clockwise and counterclockwise rotations of the C₁₅C₁₆ bond in the Z and E isomers using a linear interpolation of internal coordinates between the Franck-Condon and MECI geometries. In the minimal chromophore model that lacks the substitutions at the pyrrole rings, the isomerization involves both C₁₄-C₁₅ and C₁₅C₁₆ bonds of the methine bridge between the C- and D-rings, resembling the hula-twist motion. The MECIs are characterized by a partial charge transfer between the two pyrrole rings pointing towards a twisted intramolecular charge transfer. Systematic introduction of substituents leads to an increase in the steric repulsion between the two pyrrole rings causing a pretwist of the dihedral around the C₁₅C₁₆ bond, which creates a preference for the counterclockwise isomerization. An introduction of the carbonyl group at the D-ring increases the extent of charge transfer which changes the isomerization mechanism from hula-twist to one-bond flip.

The impact of different geometrical restrictions on the nonadiabatic photoisomerization of biliverdin chromophores

The photoisomerization mechanism of the chromophore of bacterial biliverdin (BV) phytochromes is explored via nonadiabatic dynamics simulation by using the on-the-fly trajectory surface-hopping method at the semi-empirical OM2/MRCI level. Particularly, the current study focuses on the influence of geometrical constrains on the nonadiabatic photoisomerization dynamics of the BV chromophore. Here a rather simplified approach is employed in the nonadiabatic dynamics to capture the features of geometrical constrains, which adds mechanical restrictions to the specific moieties of the BV chromophore. This simplified method provides a rather quick approach to examine the influence of geometrical restrictions on photoisomerization. As expected, different constrains bring distinctive influences on the photoisomerization mechanism of the BV chromophore, giving either strong or minor modification of both involved reaction channels and excited-state lifetimes after the constrains are added in different ring moieties. These observations not only contribute to the primary understanding of the role of the spatial restriction caused by biological environments in photoinduced dynamics of the BV chromophore, but also provide useful ideas for the artificial regulation of the photoisomerization reaction channels of phytochrome proteins.

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.

The interplay between chromophore and protein determines the extended excited state dynamics in a single-domain phytochrome

Phytochromes are a diverse family of bilin-binding photoreceptors that regulate a wide range of physiological processes. Their photochemical properties make them attractive for applications in optogenetics and superresolution microscopy. Phytochromes undergo reversible photoconversion triggered by the Z ⇄ E photoisomerization about the double bond in the bilin chromophore. However, it is not fully understood at the molecular level how the protein framework facilitates the complex photoisomerization dynamics. We have studied a single-domain bilin-binding photoreceptor All2699g1 (Nostoc sp. PCC 7120) that exhibits photoconversion between the red light-absorbing (P_r) and far red-absorbing (P_fr) states just like canonical phytochromes. We present the crystal structure and examine the photoisomerization mechanism of the P_r form as well as the formation of the primary photoproduct Lumi-R using time-resolved spectroscopy and hybrid quantum mechanics/molecular mechanics simulations. We show that the unusually long excited state lifetime (broad lifetime distribution centered at ∼300 picoseconds) is due to the interactions between the isomerizing pyrrole ring D and an adjacent conserved Tyr142. The decay kinetics shows a strongly distributed character which is imposed by the nonexponential protein dynamics. Our findings offer a mechanistic insight into how the quantum efficiency of the bilin photoisomerization is tuned by the protein environment, thereby providing a structural framework for engineering bilin-based optical agents for imaging and optogenetics applications.

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.

Protonation of the Biliverdin IXα Chromophore in the Red and Far-Red Photoactive States of a Bacteriophytochrome

The tetrapyrrole chromophore biliverdin IXα (BV) in the bacteriophytochrome from Deinococcus radiodurans (DrBphP) is usually assumed to be fully protonated, but this assumption has not been systematically validated by experiments or extensive computations. Here, we use force field molecular dynamics simulations and quantum mechanics/molecular mechanics calculations with density functional theory and XMCQDPT2 methods to investigate the effect of the five most probable protonation forms of BV on structural stability, binding pocket interactions, and absorption spectra in the two photochromic states of DrBphP. While agreement with X-ray structural data and measured UV/vis spectra suggest that in both states the protonated form of the chromophore dominates, we also find that a minor population with a deprotonated D-ring could contribute to the red-shifted tail in the absorption spectra.

Role of the protein cavity in phytochrome chromoprotein assembly and double-bond isomerization: a comparison with model compounds

Difference patterns of (13)C NMR chemicals shifts for the protonation of a free model compound in organic solution, as reported in the literature (M. Stanek, K. Grubmayr [1998] Chem. Eur. J.4, 1653-1659), were compared with changes in the protonation state occurring during holophytochrome assembly from phycocyanobilin (PCB) and the apoprotein. Both processes induce identical changes in the NMR signals, indicating that the assembly process is linked to protonation of the chromophore, yielding a cationic cofactor in a heterogeneous, quasi-liquid protein environment. The identity of both difference patterns implies that the protonation of a model compound in solution causes a partial stretching of the geometry of the macrocycle as found in the protein. In fact, the similarity of the difference pattern within the bilin family for identical chemical transformations represents a basis for future theoretical analysis. On the other hand, the change of the (13)C NMR chemical shift pattern upon the Pr --> Pfr photoisomerization is very different to that of the free model compound upon ZZZ --> ZZE photoisomerization. Hence, the character of the double-bond isomerization in phytochrome is essentially different from that of a classical photoinduced double-bond isomerization, emphasizing the role of the protein environment in the modulation of this light-induced process.