Computational evidence in favor of a protonated chromophore in the photoactivation of phytochrome

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.

E to Z Photoisomerization of Phytochrome Cph1Δ Exceeds the Born-Oppenheimer Adiabatic Limit

The Born-Oppenheimer adiabatic limit applies broadly in chemistry because most reactions occur on the ground electronic state. Photochemical reactions involve two or more electronic states and need not be subject to this adiabatic limit. The spectroscopic signatures of nonadiabatic processes are subtle, and therefore, experimental investigations have been limited to the few systems dominated by single photochemical outcomes. Systems with branched excited-state pathways have been neglected, despite their potential to reveal insights into photochemical reactivity. Here we present experimental evidence from coherent three-dimensional electronic spectroscopy that the E to Z photoisomerization of phytochrome Cph1 is strongly nonadiabatic, and the simulations reproduce the measured features only when the photoisomerization proceeds nonadiabatically near, but not through, a conical intersection. The results broaden the general understanding of photoisomerization mechanisms and motivate future studies of nonadiabatic processes with multiple outcomes arising from branching on excited-state potential energy surfaces.

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.

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

Phytochromes are widespread photoreceptors responsive to red and far-red light that exist in two photochromic forms Pr (inactive) and Pfr (active). The Pr --> Pfr conversion proceeds through a series of events initiated by Z-->E photoisomerization of the tetrapyrrole chromophore, believed to occur at C15 of the methine bridge between rings C and D. Recent crystal structures show that ring D in Pr is less tightly packed by the protein than rings A, B and C, which should favor the C15 reaction over reactions at C4 (AB methine bridge) and C10 (BC). In the present work, quantum chemical methods are used to establish the intrinsic reactivity of the chromophore towards all three possible Z-->E isomerization events in the absence of steric effects and specific interactions with the protein. Using a level of theory that reproduces spectroscopic data with an accuracy of approximately 0.2 eV, it is demonstrated that isolated conditions allow the C10 photoreaction to substantially dominate. This finding suggests that the different degrees of ring-packing observed in the protein are crucial not only to facilitate a reaction at C15, but also to prevent an intrinsically more favorable reaction at C10 from taking place.

Which factors determine the acidity of the phytochromobilin chromophore of plant phytochrome?

Quantum chemical calculations aimed at identifying the factors controlling the acidity of phytochromobilin, the tetrapyrrole chromophore of the plant photoreceptor phytochrome, are reported. Phytochrome is converted from an inactive (Pr) to an active form (Pfr) through a series of events initiated by a Z --> E photoisomerization of phytochromobilin, forming the Lumi-R intermediate, and much controversy exists as to whether the protonation state of the chromophore (cationic in Pr with all nitrogens protonated) changes during the photoactivation. Here, relative ground (S0) and excited-state (S1) pKa s of all four pyrrole moieties of phytochromobilin in all 64 possible configurations with respect to the three methine bridges are calculated in a protein-like environment, using a recently benchmarked level of theory. Accordingly, the relationships between acidity and chromophore geometry and charge distribution, hydrogen bonding, and light absorption are investigated in some detail, and discussed in terms of possible mechanisms making a proton transfer reaction more probable along the Pr --> Pfr reaction than in the parent cationic Pr state. It is found that charge distribution in the cationic species, intra-molecular hydrogen bonding in the neutral, and hydrogen bonding with two highly conserved aspartate and histidine residues have a significant effect on the acidity, while overall chromophore geometry and electronic state are less important factors. Furthermore, based on the calculations, two processes that may facilitate a proton transfer by substantially lowering the pKa s relative to their Pr values are identified: (i) a thermal Z,anti --> Z,syn isomerization at C5, occurring after formation of Lumi-R; (ii) a perturbation of the hydrogen bonding network which in Pr comprises the nitrogens of pyrroles A, B and C and the two aspartate and histidine residues.

Relative Ground and Excited-State pKa Values of Phytochromobilin in the Photoactivation of Phytochrome: A Computational Study

The conversion of the plant photoreceptor phytochrome from an inactive (Pr) to an active form (Pfr) is accomplished by a red-light induced Z --> E photoisomerization of its phytochromobilin chromophore. In recent years, the question whether the photoactivation involves a change in chromophore protonation state has been the subject of many experimental studies. Here, we have used quantum chemical methods to calculate relative ground and excited-state pKa values of the different pyrrole moieties of phytochromobilin in a protein-like environment. Assuming (based on experimental data) a Pr ZaZsZa chromophore and considering isomerizations at C15 and C5, it is found that moieties B and C are the strongest acids both in the ground state and in the bright first singlet excited state, which is rationalized in simple geometric and electronic terms. It is also shown that neither light absorption nor isomerization increases the acidity of phytochromobilin relative to the reference Pr state with all pyrrolenic nitrogens protonated. Hence, provided that the subset of chromophore geometries under investigation is biologically relevant, there appears to be no intrinsic driving force for a proton-transfer event. In a series of benchmark calculations, the performance of ab initio and time-dependent density functional theory methods for excited-state studies of phytochromobilin is evaluated in light of available experimental data.

Protein-bound chromophores astaxanthin and phytochromobilin: excited state quantum chemical studies

We present an overview of excited state quantum chemical calculations aimed at elucidating controversial issues regarding the photochemistry of the protein-bound chromophores astaxanthin and phytochromobilin. In particular, we show how the application of time-dependent density functional theory and other single-reference quantum chemical excited state methods have contributed to shed new light on the origin of the >0.5 eV bathochromic shift of the electronic absorption by the carotenoid astaxanthin in the protein macromolecular complex crustacyanin, and the mechanism for C15-Z,syn --> C15-E,anti isomerization of the tetrapyrrole phytochromobilin that underlies the photoactivation of the plant photoreceptor phytochrome. Within the approximation that exciton coupling is neglected, the calculations on astaxanthin provide support for the notion that the bathochromic shift, which is responsible for the slate-blue coloration of lobster shell, is due to polarization rather than a conformational change of the chromophore in the protein-bound state. Furthermore, the polarization is attributed to a hydrogen-bonded protonated histidine residue. The calculations on phytochromobilin, in turn, suggest that a stepwise C15-Z,syn --> C15-E,syn (photochemical), C15-E,syn --> C15-E,anti (thermal) mechanism is much more favorable than a concerted, fully photochemical mechanism, and that neutral forms of the chromophore are much less likely to photoisomerize than the parent, protonated form. Accordingly, the calculations indirectly support the view that the photoactivation of phytochrome does not involve a proton transfer from the chromophore to the surrounding protein.