Signal transduction in the photoactive yellow protein. I. Photon absorption and the isomerization of the chromophore

Molecular dynamics simulation techniques together with time-dependent density functional theory calculations have been used to investigate the effect of photon absorption by a 4-hydroxy-cinnamic acid chromophore on the structural properties of the photoactive yellow protein (PYP) from Ectothiorodospira halophila. The calculations suggest that the protein not only modifies the absorption spectrum of the chromophore but also regulates the subsequent isomerization of the chromophore by stabilizing the isomerization transition state. Although signaling from PYP is thought to involve partial unfolding of the protein, the mechanical effects accompanying isomerization do not appear to directly destabilize the protein.     

Ultrafast spectroscopy of biological photoreceptors

We review recent new insights on reaction dynamics of photoreceptors proteins gained from ultrafast spectroscopy. In Blue Light sensing Using FAD (BLUF) domains, a hydrogen-bond rearrangement around the flavin chromophore proceeds through a radical-pair mechanism, by which light-induced electron and proton transfer from the protein to flavin result in rotation of a conserved glutamine that switches the hydrogen bond network. Femtosecond infrared spectroscopy has shown that in photoactive yellow protein (PYP), breaking of a hydrogen bond that connects the p-coumaric acid chromophore to the backbone is crucial for trans-cis isomerization and successful entry into the photocycle. Furthermore, isomerization reactions of phycocyanobilin in phytochrome and retinal in the rhodopsins have been revealed in detail through application of femtosecond infrared and femtosecond-stimulated Raman spectroscopy.     

Function and structure of GFP-like proteins in the protein data bank

The RCSB protein databank contains 266 crystal structures of green fluorescent proteins (GFP) and GFP-like proteins. This is the first systematic analysis of all the GFP-like structures in the pdb. We have used the pdb to examine the function of fluorescent proteins (FP) in nature, aspects of excited state proton transfer (ESPT) in FPs, deformation from planarity of the chromophore and chromophore maturation. The conclusions reached in this review are that (1) The lid residues are highly conserved, particularly those on the "top" of the β-barrel. They are important to the function of GFP-like proteins, perhaps in protecting the chromophore or in β-barrel formation. (2) The primary/ancestral function of GFP-like proteins may well be to aid in light induced electron transfer. (3) The structural prerequisites for light activated proton pumps exist in many structures and it's possible that like bioluminescence, proton pumps are secondary functions of GFP-like proteins. (4) In most GFP-like proteins the protein matrix exerts a significant strain on planar chromophores forcing most GFP-like proteins to adopt non-planar chromophores. These chromophoric deviations from planarity play an important role in determining the fluorescence quantum yield. (5) The chemospatial characteristics of the chromophore cavity determine the isomerization state of the chromophore. The cavities of highlighter proteins that can undergo cis/trans isomerization have chemospatial properties that are common to both cis and trans GFP-like proteins.     

Molecular modeling of green fluorescent protein: structural effects of chromophore deprotonation

Molecular dynamics (MD) simulations were carried out to study the conformational rearrangement induced by deprotonation of the fluorescent chromophore in GFP, as well as the associated changes in the hydrogen-bonding network. For both the structures with either a neutral or an anionic chromophore, it was found that the beta-barrel was stable and rigid, and the conformation of the chromophore was consistent with the available x-ray structure. The conformational change in Thr203 due to deprotonation was also found to be consistent with the three-state isomerization model. Although GFP is highly fluorescent, denatured-GFP is nonfluorescent, indicating that the environment of the protein plays an important role in its fluorescence behavior. Our MD simulations, which explore the effect of the protein shell on the conformation of the chromophore, find the flexibility of the central chromophore to be significantly restricted due to the rigid nature of the protein shell. The hydrogen-bonding between the chromophore and neighboring residues was also shown to contribute to the chromophore rigidity. In addition to the MD studies, quantum mechanics/molecular mechanics (QM/MM) ONIOM calculations were carried out to investigate the effect of the beta-barrel on the internal rotation in the chromophore. Along with providing quantitative values for torsional rotation barriers about the bridging bond in the chromophore, the ONIOM calculations also validate our MD force field parameters.     

A molecular movie at 1.8 A resolution displays the photocycle of photoactive yellow protein, a eubacterial blue-light receptor, from nanoseconds to seconds.

The photocycle of the bacterial blue-light photoreceptor, photoactive yellow protein, was stimulated by illumination of single crystals by a 7 ns laser pulse. The molecular events were recorded at high resolution by time-resolved X-ray Laue diffraction as they evolved in real time, from 1 ns to seconds after the laser pulse. The complex structural changes during the photocycle at ambient temperature are displayed in a movie of difference electron density maps relative to the dark state. The step critical to entry into the photocycle is identified as flipping of the carbonyl group of the 4-hydroxycinnamic acid chromophore into an adjacent, hydrophobic environment rather than the concomitant isomerization about the double bond of the chromophore tail. The structural perturbation generated at the chromophore propagates throughout the entire protein as a light-induced "protein quake" with its "epicenter" at the carbonyl moiety of the chromophore.     

A structural basis for reversible photoswitching of absorbance spectra in red fluorescent protein rsTagRFP

rsTagRFP is the first monomeric red fluorescent protein (FP) with reversibly photoswitchable absorbance spectra. The switching is realized by irradiation of rsTagRFP with blue (440 nm) and yellow (567 nm) light, turning the protein fluorescence ON and OFF, respectively. It is perhaps the most useful probe in this color class that has yet been reported. Because of the photoswitchable absorbance, rsTagRFP can be used as an acceptor in photochromic Förster resonance energy transfer. Yellow FPs, YPet and mVenus, are demonstrated to be excellent photochromic Förster resonance energy transfer donors for the rsTagRFP acceptor in its fusion constructs. Analysis of X-ray structures has shown that photoswitching of rsTagRFP is accompanied by cis–trans isomerization and protonation/deprotonation of the chromophore, with the deprotonated cis- and protonated trans-isomers corresponding to its ON and OFF states, respectively. Unlike in other photoswitchable FPs, both conformers of rsTagRFP chromophore are essentially coplanar. Two other peculiarities of the rsTagRFP chromophore are an essentially hydrophobic environment of its p-hydroxyphenyl site and the absence of direct hydrogen bonding between this moiety and the protein scaffold. The influence of the immediate environment on rsTagRFP chromophore was probed by site-directed mutagenesis. Residues Glu145 and His197 were found to participate in protonation/deprotonation of the chromophore accompanying the photoswitching of rsTagRFP fluorescence, whereas residues Met160 and Leu174 were shown to spatially restrict chromophore isomerization, favoring its radiative decay.     

Deactivation of rhodopsin in the transition from the signaling state meta II to meta III involves a thermal isomerization of the retinal chromophore C[double bond]D

Light-induced isomerization of rhodopsin's retinal chromophore to the activating all-trans geometry initializes the formation of the active receptor state, Meta II. In the absence of peripheral regulatory proteins, the activity of Meta II is switched off spontaneously by two independent pathways: either by hydrolysis of the retinal Schiff base and dissociation of the light receptor into apoprotein opsin plus free retinal or by formation of Meta III, an inactive species with intact retinal protonated Schiff base absorbing at 470 nm. By FTIR spectroscopy on rhodopsin reconstituted with isotopically labeled chromophores in combination with quantum mechanical DFT calculations, we show that the deactivating step during formation of Meta III involves a thermal isomerization of the chromophore C[double bond]N, such that the chromophore in Meta III is all-trans-15-syn. This isomerization step is catalyzed by the protein environment and proceeds via Meta I, as suggested by its dependence on pH and on properties of the lipid/detergent environment of the protein. In the long term, Meta III decays likewise to opsin and free retinal by slow hydrolysis of the Schiff base.     

Mechanism of rhodopsin activation as examined with ring-constrained retinal analogs and the crystal structure of the ground state protein

The guanine nucleotide-binding protein (G-protein)-coupled receptor superfamily (GPCR) is comprised of a large group of membrane proteins involved in a wide range of physiological signaling processes. The functional switch from a quiescent to an active conformation is at the heart of GPCR action. The GPCR rhodopsin has been studied extensively because of its key role in scotopic vision. The ground state chromophore, 11-cis-retinal, holds the transmembrane region of the protein in the inactive conformation. Light induces cis-trans isomerization and rhodopsin activation. Here we show that rhodopsin regenerated with a ring-constrained 11-cis-retinal analog undergoes photoisomerization; however, it remains marginally active because isomerization occurs without the chromophore-induced conformational change of the opsin moiety. Modeling the locked chromophore analogs in the active site of rhodopsin suggests that the beta-ionone ring rotates but is largely confined within the binding site of the natural 11-cis-retinal chromophore. This constraint is a result of the geometry of the stable 11-cis-locked configuration of the chromophore analogs. These results suggest that the native chromophore cis-trans isomerization is merely a mechanism for repositioning of the beta-ionone ring which ultimately leads to helix movements and determines receptor activation.     

Role of protein in the primary step of the photoreaction of yellow protein

We show the unexpectedly important role of the protein environment in the primary step of the photoreaction of the yellow protein after light illumination. The driving force of the trans-to-cis isomerization reaction was analyzed by a computational method. The force was separated into two different components: the term due to the protein-chromophore interaction and the intrinsic term of the chromophore itself. As a result, we found that the contribution from the interaction term was much greater than that coming from the intrinsic term. This accounts for the efficiency of the isomerization reaction in the protein environment in contrast to that in solution environments. We then analyzed the relaxation process of the chromophore on the excited-state energy surface and compared the process in the protein environment and that in a vacuum. Based on this analysis, we found that the bond-selectivity of the isomerization reaction also comes from the interaction between the chromophore and the protein environment.     

Red fluorescent protein DsRed: Parametrization of its chromophore as an amino acid residue for computer modeling in the OPLS-AA force field

Topology of the neutral form of the DsRed fluorescent protein chromophore as a residue of [(4-cis)-2-[(1-cis)-4-amino-4-oxobutanimidoyl]-4-(4-hydroxybenzylidene)-5-oxo-4,5-dihydro-1H-imidazol-1-yl]acetic acid was calculated with OPLS-AA force field. Use of this topology and molecular dynamics simulation allows calculating the parameters of proteins that contain such residue in their polypeptide chains. The chromophore parameters were obtained by ab initio (RHF/6-31G**) quantum chemical calculations applying density functional theory (B3LYP). Using this chromophore, we have calculated the molecular dynamics trajectory of tetrameric fluorescent protein DsRed in solution at 300 K (4 nsec). Correctness of the chromophore parametrization was revealed by comparison of quantitative characteristics of the chromophore structure obtained from the molecular dynamic simulations of DsRed protein with the quantitative characteristics of the chromophore based on the crystallographic X-ray data of fluorescent protein DsRed (PDB ID: 1ZGO, 1G7K, and 1GGX), and also with the quantitative characteristics of the chromophore obtained by quantum chemical calculations. Inclusion of the neutral form of DsRed protein chromophore topology into the OPLS-AA force field yielded the extended force field OPLS-AA/DsRed. This force field can be used for molecular dynamics calculations of proteins containing the DsRed chromophore. The parameter set presented in this study can be applied for similar extension in any other force fields.     

Chromophore protonation state controls photoswitching of the fluoroprotein asFP595

Fluorescent proteins have been widely used as genetically encodable fusion tags for biological imaging. Recently, a new class of fluorescent proteins was discovered that can be reversibly light-switched between a fluorescent and a non-fluorescent state. Such proteins can not only provide nanoscale resolution in far-field fluorescence optical microscopy much below the diffraction limit, but also hold promise for other nanotechnological applications, such as optical data storage. To systematically exploit the potential of such photoswitchable proteins and to enable rational improvements to their properties requires a detailed understanding of the molecular switching mechanism, which is currently unknown. Here, we have studied the photoswitching mechanism of the reversibly switchable fluoroprotein asFP595 at the atomic level by multiconfigurational ab initio (CASSCF) calculations and QM/MM excited state molecular dynamics simulations with explicit surface hopping. Our simulations explain measured quantum yields and excited state lifetimes, and also predict the structures of the hitherto unknown intermediates and of the irreversibly fluorescent state. Further, we find that the proton distribution in the active site of the asFP595 controls the photochemical conversion pathways of the chromophore in the protein matrix. Accordingly, changes in the protonation state of the chromophore and some proximal amino acids lead to different photochemical states, which all turn out to be essential for the photoswitching mechanism. These photochemical states are (i) a neutral chromophore, which can trans-cis photoisomerize, (ii) an anionic chromophore, which rapidly undergoes radiationless decay after excitation, and (iii) a putative fluorescent zwitterionic chromophore. The overall stability of the different protonation states is controlled by the isomeric state of the chromophore. We finally propose that radiation-induced decarboxylation of the glutamic acid Glu215 blocks the proton transfer pathways that enable the deactivation of the zwitterionic chromophore and thus leads to irreversible fluorescence. We have identified the tight coupling of trans-cis isomerization and proton transfers in photoswitchable proteins to be essential for their function and propose a detailed underlying mechanism, which provides a comprehensive picture that explains the available experimental data. The structural similarity between asFP595 and other fluoroproteins of interest for imaging suggests that this coupling is a quite general mechanism for photoswitchable proteins. These insights can guide the rational design and optimization of photoswitchable proteins.     

X-ray structure of Cerulean GFP: a tryptophan-based chromophore useful for fluorescence lifetime imaging

The crystal structure of the cyan-fluorescent Cerulean green fluorescent protein (GFP), a variant of enhanced cyan fluorescent protein (ECFP), has been determined to 2.0 A. Cerulean bears an internal fluorophore composed of an indole moiety derived from Y66W, conjugated to the GFP-like imidazolinone ring via a methylene bridge. Cerulean undergoes highly efficient fluorescence resonance energy transfer (FRET) to yellow acceptor molecules and exhibits significantly reduced excited-state heterogeneity. This feature was rationally engineered in ECFP by substituting His148 with an aspartic acid [Rizzo et al. (2004) Nat. Biotechnol. 22, 445], rendering Cerulean useful for fluorescence lifetime imaging microscopy (FLIM). The X-ray structure is consistent with a single conformation of the chromophore and surrounding residues and may therefore provide a structural rationale for the previously described monoexponential fluorescence decay. Unexpectedly, the carboxyl group of H148D is found in a buried position, directly contacting the indole nitrogen of the chromophore via a bifurcated hydrogen bond. Compared to the similarly constructed ECFP chromophore, the indole group of Cerulean is rotated around the methylene bridge to adopt a cis-coplanar conformation with respect to the imidazolinone ring, resulting in a close edge-to-edge contact of the two ring systems. The double-humped absorbance spectrum persists in single-crystal absorbance measurements, casting doubt on the idea that ground state conformational heterogeneity forms the basis of the two overlapping transitions. At low pH, a blue shift in absorbance of 10-15 nm suggests a pH-induced structural transition that proceeds with a time constant of 47 (+/-2) min and is reversible. Possible interpretations in terms of chromophore isomerization are presented.     

Initial steps of signal generation in photoactive yellow protein revealed with femtosecond mid-infrared spectroscopy

Photoactive yellow protein (PYP) is a bacterial blue light sensor that induces Halorhodospira halophila to swim away from intense blue light. Light absorption by PYP's intrinsic chromophore, p-coumaric acid, leads to the initiation of a photocycle that comprises several distinct intermediates. Here we describe the initial structural changes of the chromophore and its nearby amino acids, using visible pump/mid-infrared probe spectroscopy. Upon photoexcitation, the trans bands of the chromophore are bleached, and shifts of the phenol ring bands occur. The latter are ascribed to charge translocation, which probably plays an essential role in driving the trans to cis isomerization process. We conclude that breaking of the hydrogen bond of the chromophore's C=O group with amino acid Cys69 and formation of a stable cis ground state occur in approximately 2 ps. Dynamic changes also include rearrangements of the hydrogen-bonding network of the amino acids around the chromophore. Relaxation of the coumaryl tail of the chromophore occurs in 0.9-1 ns, which event we identify with the I(0) to I(1) transition observed in visible spectroscopy.     

Crystallographic structures of Discosoma red fluorescent protein with immature and mature chromophores: linking peptide bond trans-cis isomerization and acylimine formation in chromophore maturation

The mature self-synthesizing p-hydroxybenzylideneimidazolinone-like fluorophores of Discosoma red fluorescent protein (DsRed) and Aequorea victoria green fluorescent protein (GFP) are extensively studied as powerful biological markers. Yet, the spontaneous formation of these fluorophores by cyclization, oxidation, and dehydration reactions of tripeptides within their protein environment remains incompletely understood. The mature DsRed fluorophore (Gln 66, Tyr 67, and Gly 68) differs from the GFP fluorophore by an acylimine that results in Gln 66 Calpha planar geometry and by a Phe 65-Gln 66 cis peptide bond. DsRed green-to-red maturation includes a green-fluorescing immature chromophore and requires a chromophore peptide bond trans-cis isomerization that is slow and incomplete. To clarify the unique structural chemistry for the individual immature "green" and mature "red" chromophores of DsRed, we report here the determination and analysis of crystal structures for the wild-type protein (1.4 A resolution), the entirely green DsRed K70M mutant protein (1.9 A resolution), and the DsRed designed mutant Q66M (1.9 A resolution), which shows increased red chromophore relative to the wild-type DsRed. Whereas the mature, red-fluorescing chromophore has the expected cis peptide bond and a sp(2)-hybridized Gln 66 Calpha with planar geometry, the crystal structure of the immature green-fluorescing chromophore of DsRed, presented here for the first time, reveals a trans peptide bond and a sp(3)-hybridized Gln 66 Calpha with tetrahedral geometry. These results characterize a GFP-like immature green DsRed chromophore structure, reveal distinct mature and immature chromophore environments, and furthermore provide evidence for the coupling of acylimine formation with trans-cis isomerization.     

Kinetics of and intermediates in a photocycle branching reaction of the photoactive yellow protein from Ectothiorhodospira halophila.

We have studied the kinetics of the blue light-induced branching reaction in the photocycle of photoactive yellow protein (PYP) from Ectothiorhodospira halophila, by nanosecond time-resolved UV/Vis spectroscopy. As compared to the parallel dark recovery reaction of the presumed blue-shifted signaling state pB, the light-induced branching reaction showed a 1000-fold higher rate. In addition, a new intermediate was detected in this branching pathway, which, compared to pB, showed a larger extinction coefficient and a blue-shifted absorption maximum. This substantiates the conclusion that isomerization of the chromophore is the rate-controlling step in the thermal photocycle reactions of PYP and implies that absorption of a blue photon leads to cis-->trans isomerization of the 4-hydroxy-cinnamyl chromophore of PYP in its pB state.     

Correlation between absorption maxima and thermal isomerization rates in bacteriorhodopsin

The reported rates of thermal 13-cis to all-trans isomerization of the protonated Schiff base of retinal (PSBR) in solution and in bacteriorhodopsin (BR) are shown to be correlated with the red shift in the absorption maximum of the chromophore, though the linear fit is different for BR and for a model PSBR in solution. Because the red shift in the absorption has been previously shown to be correlated with π-electron delocalization in the chromophore, this suggests that the thermal isomerization rate is largely regulated by the amount of double bond character in the chromophore. Because the linear fit of isomerization rates with absorption maxima is different for BR and the model PSBR, specific interactions of the protein with the chromophore must also be a factor in determining thermal isomerization rates in BR. A model of the later steps in the photocycle of BR is presented in which the 13-cis to all-trans thermal isomerization occurs during the O intermediate.     

Resonance Raman analysis of the Pr and Pfr forms of phytochrome

Resonance Raman vibrational spectra of the Pr and Pfr forms of oat phytochrome have been obtained at room temperature. When Pr is converted to Pfr, new bands appear in the C = C and C = N stretching region at 1622, 1599, and 1552 cm-1, indicating that a major structural change of the chromophore has occurred. The Pr to Pfr conversion results in an 11 cm-1 lowering of the N-H rocking band from 1323 to 1312 cm-1. Normal mode calculations correlate this frequency drop with a Z----E isomerization about the C15 = C16 bond. A line at 803 cm-1 in Pr is replaced by an unusually intense mode at 814 cm-1 in Pfr. Calculations on model tetrapyrrole chromophores suggest that these low-wavenumber modes are hydrogen out-of-plane (HOOP) wagging vibrations of the bridging C15 methine hydrogen and that both the intensity and frequency of the C15 HOOP mode are sensitive to the geometry around the C14-C15 and C15 = C16 bonds. The large intensity of the 814-cm-1 mode in Pfr indicates that the chromophore is highly distorted from planarity around the C15 methine bridge. If the Pr----Pfr conversion does involve a C15 = C16 Z----E isomerization, then the intensity of the C15 HOOP mode in Pfr argues that the chromophore has an E,anti conformation. On the basis of a comparison with the vibrational calculations, the low frequency (803 cm-1) and the reduced intensity of the C15 HOOP mode in Pr suggest that the chromophore in Pr adopts the C15-Z,syn conformation.     

Molecular dynamics simulations of retinal in rhodopsin: from the dark-adapted state towards lumirhodopsin

The formation of photointermediates and conformational changes observed in the retinal chromophore of bilayer-embedded rhodopsin during the early steps of the protein activation have been studied by molecular dynamics (MD) simulation. In particular, the lysine-bound retinal has been examined, focusing on its conformation in the dark-adapted state (10 ns) and on the early steps after the isomerization of the 11-cis bond to trans (up to 10 ns). The parametrization for the chromophore is based on a recent quantum study [Sugihara, M., Buss, V., Entel, P., Elstner, M., and Frauenheim, T. (2002) Biochemistry 41, 15259-15266] and shows good conformational agreement with recent experimental results. The isomerization, induced by switching the function governing the dihedral angle for the C11=C12 bond, was repeated with several different starting conformations. From the repeated simulations, it is shown that the retinal model exhibits a conserved activation pattern. The conformational changes are sequential and propagate outward from the C11=C12 bond, starting with isomerization of the C11=C12 bond, then a rotation of methyl group C20, and followed by increased fluctuations at the beta-ionone ring. The dynamics of these changes suggest that they are linked with photointermediates observed by spectroscopy. The exact moment when these events occur after the isomerization is modulated by the starting conformation, suggesting that retinal isomerizes through multiple pathways that are slightly different. The amplitudes of the structural fluctuations observed for the protein in the dark-adapted state and after isomerization of the retinal are similar, suggesting a subtle mechanism for the transmission of information from the chromophore to the protein.     

Conformation and stability of alpha-helical membrane proteins. 2. Influence of pH and salts on stability and unfolding of rhodopsin

We studied the stability and pH-induced denaturation of rhodopsin and its photoproducts as a model for alpha-helical membrane proteins. The increased stability of the dark state of rhodopsin as compared to its photoproduct states allows the initiation of unfolding of the protein by light-dependent isomerization of the chromophore. We could therefore characterize the transition from the native to either acid or alkaline denatured states by light-induced Fourier transform infrared difference spectroscopy, UV-visible spectroscopy, and intrinsic tryptophan fluorescence spectroscopy. The results indicate a loss of important tertiary interactions within the protein and between the protein and the retinal chromophore in the denatured state, despite that the secondary structure of the protein is almost fully retained during the transition. We therefore propose that in this denatured state the protein adopts the conformation of a loose bundle of preserved, but only weakly interacting, transmembrane helices with a largely des-oriented and partly solvent-exposed chromophore. We further characterized the influence of salts on the stability of the rhodopsin helix bundle, which was found to follow the Hofmeister series. We found that the effect of sodium chloride may be stabilizing or destabilizing, depending on the intrinsic stability of the examined protein conformation and on salt concentration. In particular, sodium chloride is shown to counteract the formation of the denatured loose bundle state presumably by increasing the lateral pressure on the helix bundle, thereby stabilizing native-like tertiary contacts within the protein.