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

A dual experimental-theoretical perspective on ESPT photoacids and their challenges ahead

Photoacids undergo an increase in acidity upon electronic excitation, enabling excited-state proton transfer (ESPT) reactions. A multitude of compounds that allow ESPT has been identified and integrated in numerous applications, as is outlined by reviewing the rich history of photoacid research reaching back more than 90 years. In particular, achievements together with ambitions and challenges are highlighted from a combined experimental and theoretical perspective. Besides explicating the spectral signatures, transient ion-pair species, and electronic states involved in an ESPT, special emphasis is put on the diversity of methods used for studying photoacids as well as on the effects of the environment on the ESPT, illustrated in detail for 8-hydroxypyrene-1,3,6-trisulfonate (HPTS) and the naphthols as examples of prototypical photoacids. The development of exceptionally acidic super-photoacids and magic photoacids is subsequently discussed, which opens the way to applications even in aprotic solvents and provides additional insight into the mechanisms underlying ESPT. In the overview of highlights from theory, a comprehensive picture of the scope of studies on HPTS is presented, along with the general conceptualization of the electronic structure of photoacids and approaches for the quantification of excited-state acidity. We conclude with a juxtaposition of established applications of photoacids together with potential open questions and prospective research directions.

Long-Distance Protonation-Conformation Coupling in Phytochrome Species

Phytochromes are biological red/far-red light sensors found in many organisms. The connection between photoconversion and the cellular output signal involves light-mediated global structural changes in the interaction between the photosensory module (PAS-GAF-PHY, PGP) and the C-terminal transmitter (output) module. We recently showed a direct correlation of chromophore deprotonation with pH-dependent conformational changes in the various domains of the prototypical phytochrome Cph1 PGP. These results suggested that the transient phycocyanobilin (PCB) chromophore deprotonation is closely associated with a higher protein mobility both in proximal and distal protein sites, implying a causal relationship that might be important for the global large-scale protein rearrangements. Here, we investigate the prototypical biliverdin (BV)-binding phytochrome Agp1. The structural changes at various positions in Agp1 PGP were investigated as a function of pH using picosecond time-resolved fluorescence anisotropy and site-directed fluorescence labeling of cysteine variants of Agp1 PGP. We show that the direct correlation of chromophore deprotonation with pH-dependent conformational changes does not occur in Agp1. Together with the absence of long-range effects between the PHY domain and chromophore pK_a, in contrast to the findings in Cph1, our results imply phytochrome species-specific correlations between transient chromophore deprotonation and intramolecular signal transduction.

Theoretical Study on the Photoacidity of Hydroxypyrene Derivatives in DMSO Using ADC(2) and CC2

This work applies the thermodynamic Förster cycle to theoretically investigate the pK_a^*, i.e., excited-state pK_a values of pyranine-derived superphotoacids developed by Jung and co-workers. The latter photoacids are strong enough to transfer a proton to the aprotic solvent dimethyl sulfoxide (DMSO). The Förster cycle provides access to pK_a^* via the ground-state pK_a and the electronic excitation energies. We use the conductor-like screening model for real solvents (COSMO-RS) to compute the ground-state pK_a and the correlated wavefunction-based methods ADC(2) and CC2 with the continuum solvation model COSMO to calculate the pK_a change upon excitation. A comparison of the calculated UV/Vis absorption and fluorescence emission energies to the experimental results leads us to infer that this approach allows for a proper description of the electronic excitations. In particular, implicit solvation by means of the COSMO model appears to be sufficient for the treatment of these photoacids in DMSO. The calculations confirm the presumption that a charge redistribution from the hydroxy group to the aromatic ring and the electron-withdrawing substituents is the origin of photoacidity for these photoacids. Moreover, the calculations with the continuum solvation model predict that the pK_a jump upon excitation decreases with increasing solvent polarity, as rationalized based on the Förster cycle.

Photobase-Driven Excited-State Intramolecular Proton Transfer (ESIPT) in a Strapped π-Electron System

We report a new design strategy for an excited-state intramolecular proton transfer (ESIPT) fluorophore that can be used in acidic media. A photobasic pyridine-centered donor-acceptor-donor-type fluorophore is combined with a basic trialkylamine "strap". In the presence of an acid, protonation occurs predominantly at the amine moiety in the ground state. A single-crystal X-ray diffraction analysis confirmed the formation of a pre-organized intramolecular hydrogen-bonded structure between the resulting ammonium moiety and the pyridine ring. Upon excitation, the intramolecular charge-transfer transition increases the basicity of the pyridine moiety in the excited state, resulting in proton transfer from the amine to the pyridine moiety. Consequently, the fluorophore takes on a polymethine-dye character in the ESIPT state, which gives rise to significantly red-shifted emission with an increased fluorescence quantum yield.

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.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

Unusual spectral properties of bacteriophytochrome Agp2 result from a deprotonation of the chromophore in the red-absorbing form Pr

Phytochromes are widely distributed photoreceptors with a bilin chromophore that undergo a typical reversible photoconversion between the two spectrally different forms, Pr and Pfr. The phytochrome Agp2 from Agrobacterium tumefaciens belongs to the group of bathy phytochromes that have a Pfr ground state as a result of the Pr to Pfr dark conversion. Agp2 has untypical spectral properties in the Pr form reminiscent of a deprotonated chromophore as confirmed by resonance Raman spectroscopy. UV/visible absorption spectroscopy showed that the pKa is >11 in the Pfr form and ∼7.6 in the Pr form. Unlike other phytochromes, photoconversion thus results in a pKa shift of more than 3 units. The Pr/Pfr ratio after saturating irradiation with monochromatic light is strongly pH-dependent. This is partially due to a back-reaction of the deprotonated Pr chromophore at pH 9 after photoexcitation as found by flash photolysis. The chromophore protonation and dark conversion were affected by domain swapping and site-directed mutagenesis. A replacement of the PAS or GAF domain by the respective domain of the prototypical phytochrome Agp1 resulted in a protonated Pr chromophore; the GAF domain replacement afforded an inversion of the dark conversion. A reversion was also obtained with the triple mutant N12S/Q190L/H248Q, whereas each single point mutant is characterized by decelerated Pr to Pfr dark conversion.

Fluorescence quantum yield and photochemistry of bacteriophytochrome constructs

Bacteriophytochromes (Bphs) are red-light photoreceptor proteins with a photosensory core that consists of three distinct domains, PAS, GAF and PHY, and covalently binds biliverdin (BV) to a conserved cysteine in the PAS domain. In a recent development, PAS-GAF variants were engineered for use as a near-infrared fluorescent marker in mammalian tissues (Tsien and co-workers, Science, 2009, 324, 804-807). Here, we report the fluorescence quantum yield and photochemistry of two highly-related Bphs from Rps. palustris, RpBphP2 (P2) and RpBphP3 (P3) with distinct photoconversion and fluorescence properties. We applied ultrafast spectroscopy to wild type P3 and P2 PAS-GAF proteins and their P3 D216A, Y272F and P2 D202A PAS-GAF-PHY mutant proteins. In these mutants hydrogen-bond interactions between a conserved aspartate (Asp) which connects the BV chromophore with the PHY domains are disrupted. The excited-state lifetime of the truncated P3 and P2 PAS-GAF proteins was significantly longer than in their PAS-GAF-PHY counterparts that constitute the full photosensory core. Mutation of the conserved Asp to Ala in the PAS-GAF-PHY protein had a similar but larger effect. The fluorescence quantum yields of the P3 D216A and Y272F mutants were 0.066, higher than that of wild type P3 (0.043) and similar to the engineered Bph of Tsien and co-workers. We conclude that elimination of a key hydrogen-bond interaction between Asp and a conserved Arg in the PHY domain is responsible for the excited-state lifetime increase in all Bph variants studied here. H/D exchange resulted in a 1.4-1.7 fold increase of excited-state lifetime. The results support a reaction model in which deactivation of the BV chromophore proceeds via excited-state proton transfer from the BV pyrrole nitrogens to the backbone of the conserved Asp or to a bound water. This work may aid in rational structure- and mechanism-based conversion of constructs based on P3 and other BPhs into efficient near-IR, deep tissue, fluorescent markers.

Proton-transfer and hydrogen-bond interactions determine fluorescence quantum yield and photochemical efficiency of bacteriophytochrome

Phytochromes are red-light photoreceptor proteins that regulate a variety of responses and cellular processes in plants, bacteria, and fungi. The phytochrome light activation mechanism involves isomerization around the C15 horizontal lineC16 double bond of an open-chain tetrapyrrole chromophore, resulting in a flip of its D-ring. In an important new development, bacteriophytochrome (Bph) has been engineered for use as a fluorescent marker in mammalian tissues. Here we report that an unusual Bph, RpBphP3 from Rhodopseudomonas palustris, denoted P3, is fluorescent. This Bph modulates synthesis of light-harvesting complex in combination with a second Bph exhibiting classical photochemistry, RpBphP2, denoted P2. We identify the factors that determine the fluorescence and isomerization quantum yields through the application of ultrafast spectroscopy to wild-type and mutants of P2 and P3. The excited-state lifetime of the biliverdin chromophore in P3 was significantly longer at 330-500 ps than in P2 and other classical phytochromes and accompanied by a significantly reduced isomerization quantum yield. H/D exchange reduces the rate of decay from the excited state of biliverdin by a factor of 1.4 and increases the isomerization quantum yield. Comparison of the properties of the P2 and P3 variants shows that the quantum yields of fluorescence and isomerization are determined by excited-state deprotonation of biliverdin at the pyrrole rings, in competition with hydrogen-bond rupture between the D-ring and the apoprotein. This work provides a basis for structure-based conversion of Bph into an efficient near-IR fluorescent marker.

Ultrafast excited-state isomerization in phytochrome revealed by femtosecond stimulated Raman spectroscopy

Photochemical interconversion between the red-absorbing (P(r)) and the far-red-absorbing (P(fr)) forms of the photosensory protein phytochrome initiates signal transduction in bacteria and higher plants. The P(r)-to-P(fr) transition commences with a rapid Z-to-E photoisomerization at the C(15)=C(16) methine bridge of the bilin prosthetic group. Here, we use femtosecond stimulated Raman spectroscopy to probe the structural changes of the phycocyanobilin chromophore within phytochrome Cph1 on the ultrafast time scale. The enhanced intensity of the C(15)-H hydrogen out-of-plane (HOOP) mode, together with the appearance of red-shifted C=C stretch and N-H in-plane rocking modes within 500 fs, reveal that initial distortion of the C(15)=C(16) bond occurs in the electronically excited I* intermediate. From I*, 85% of the excited population relaxes back to P(r) in 3 ps, whereas the rest goes on to the Lumi-R photoproduct consistent with the 15% photochemical quantum yield. The C(15)-H HOOP and skeletal modes evolve to a Lumi-R-like pattern after 3 ps, thereby indicating that the C(15)=C(16) Z-to-E isomerization occurs on the excited-state surface.

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.

A second conserved GAF domain cysteine is required for the blue/green photoreversibility of cyanobacteriochrome Tlr0924 from Thermosynechococcus elongatus

Phytochromes are widely occurring red/far-red photoreceptors that utilize a linear tetrapyrrole (bilin) chromophore covalently bound within a knotted PAS-GAF domain pair. Cyanobacteria also contain more distant relatives of phytochromes that lack this knot, such as the phytochrome-related cyanobacteriochromes implicated to function as blue/green switchable photoreceptors. In this study, we characterize the cyanobacteriochrome Tlr0924 from the thermophilic cyanobacterium Thermosynechococcus elongatus. Full-length Tlr0924 exhibits blue/green photoconversion across a broad range of temperatures, including physiologically relevant temperatures for this organism. Spectroscopic characterization of Tlr0924 demonstrates that its green-absorbing state is in equilibrium with a labile, spectrally distinct blue-absorbing species. The photochemically generated blue-absorbing state is in equilibrium with another species absorbing at longer wavelengths, giving a total of 4 states. Cys499 is essential for this behavior, because mutagenesis of this residue results in red-absorbing mutant biliproteins. Characterization of the C 499D mutant protein by absorbance and CD spectroscopy supports the conclusion that its bilin chromophore adopts a similar conformation to the red-light-absorbing P r form of phytochrome. We propose a model photocycle in which Z/ E photoisomerization of the 15/16 bond modulates formation of a reversible thioether linkage between Cys499 and C10 of the chromophore, providing the basis for the blue/green switching of cyanobacteriochromes.

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.