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.

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.

Ultrafast spectroscopy of biliverdin dimethyl ester in solution: pathways of excited-state depopulation

Biliverdin is a bile pigment that has a very low fluorescence quantum yield in solution, but serves as a chromophore in far-red fluorescent proteins being developed for bio-imaging. In this work, excited-state dynamics of biliverdin dimethyl ether (BVE) in solvents were investigated using femtosecond (fs) and picosecond (ps) time-resolved absorption and fluorescence spectroscopy. This study is the first fs timescale investigation of BVE in solvents, and therefore revealed numerous dynamics that were not resolved in previous, 200 ps time resolution measurements. Viscosity- and isotope-dependent experiments were performed to identify the contributions of isomerization and proton transfer to the excited-state dynamics. In aprotic solvents, a ∼2 ps non-radiative decay accounts for 95% of the excited-state population loss. In addition, a minor ∼30 ps emissive decay pathway is likely associated with an incomplete isomerization process around the C15[double bond, length as m-dash]C16 double bond that results in a flip of the D-ring. In protic solvents, the dynamics are more complex due to hydrogen bond interactions between solute and solvent. In this case, the ∼2 ps decay pathway is a minor channel (15%), whereas ∼70% of the excited-state population decays through an 800 fs emissive pathway. The ∼30 ps timescale associated with isomerization is also observed in protic solvents. The most significant difference in protic solvents is the presence of a >300 ps timescale in which BVE can decay through an emissive state, in parallel with excited-state proton transfer to the solvent. Interestingly, a small fraction of a luminous species, which we designate lumin-BVE (LBVE), is present in protic solvents.

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.

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.