5122
J. Phys. Chem. B 1999, 103, 5122-5130
On the Nature of the Primary Light-Induced Events in Bacteriorhodopsin: Ultrafast Spectroscopy of Native and C13)C14 Locked Pigments T. Ye,† N. Friedman,‡ Y. Gat,‡ G. H. Atkinson,§ M. Sheves,*,‡ M. Ottolenghi,*,† and S. Ruhman*,† Department of Physical Chemistry, The Hebrew UniVersity, Jerusalem 91904, Israel, Department of Organic Chemistry, The Weizmann Institute of Science, RehoVot 76100, Israel, and Department of Chemistry and Optical Science Center, UniVersity of Arizona, Tucson, Arizona 85721 ReceiVed: December 4, 1998
The primary light-induced events in the photosynthetic retinal protein bacteriorhodopsin (bR) are investigated by ultrafast optical spectroscopy over the 440-1000 nm spectral range. The study compares the early dynamics of the native all-trans pigment bR570 with those of two synthetic analogues, bR5.12 and bR5.13, in which isomerization around the critical C13dC14 bond is blocked by a five-membered ring into all-trans and 13-cis configurations, respectively. Nearly identical spectral evolution is observed in both native and artificial systems over the first 100-200 fs of probe delay. During this period stimulated near-IR (∼900 nm) emission, and intense ∼460 nm absorption bands, due to analogous fluorescent I states (denoted as I460, I5.12 and I5.13, respectively), appear concurrently within 30 fs. In all systems continuous spectral shifting over tens of femtoseconds is observed in the 500-700 nm range. Native bR goes on to produce the J625 absorption band within ∼1 ps, which is accompanied by disappearance of the I460 emission and absorption features. In bR5.12 and bR5.13, aside from minor spectral modifications, the analogous dramatic changes associated with I5.12 and I5.13 are arrested beyond the first ∼100 fs, reverting uniformly to the initial ground state with exponential time constants of 19 ps and 11 ps, respectively. Analysis of the data calls for a major revision of models previously put forward for the primary events in bacteriorhodopsin. The close likeness of initial transient spectral evolution in both native and artificial pigments, despite the locking of the active isomerization coordinate in the synthetic chromophores, demonstrates that in bR570 the ultrafast changes in transmission leading to I460, previously believed to involve C13dC14 torsion, must be associated with other modes. The detailed comparison conducted here also identifies which of the later spectral changes in the native system requires torsional flexibility in C13dC14. Similarity of 660 nm probing data in both synthetic and native chromophores demonstrates that the sub-picosecond dynamic features uncovered at this probing wavelength commonly attributed to the evolution of J625, are not, as previously thought, reliable measures of all-trans S 13-cis isomerization dynamics.
I. Introduction Bacteriorhodopsin (bR) is the photosynthetic representative of a series of retinal proteins, also including visual pigments, the chloride ion pump halorhodopsin (hR), and sensory rhodopsins. (For a series of comprehensive reviews on retinal proteins, see ref 1.) bR acts as a proton pump, driven by the absorption of visible light by the all-trans retinal chromophore (Structure A) which is attached to the protein via a protonated Schiff base linkage with Lys-216. Light absorption by all-trans bR, denoted as bR570, is associated with the schematic photocycle shown in Scheme 1:
SCHEME 1 hν
2 µs
3 ps
bR570 98 [S1FC(H), I460, J625] 98 K610 98 50 µs
2 ms
8 ms
L550 {\} M412 {\} (N540, O640) 98 bR570 Subscripts refer to the corresponding absorption maxima, S1FC represents the primarily generated Franck-Condon excited state, also denoted as H, and time values are approximate transition times between the various photointermediates. * Authors to whom correspondence should be addressed. † The Hebrew University. ‡ The Weizmann Institute of Science. § University of Arizona.
Of special interest are the early light-induced ultrafast dynamic events, most commonly investigated by time-resolved pump-probe spectroscopies (for recent reviews, see refs 2 and 3). These events terminate with the formation of the K610 intermediate, which is the first phototransient which may be trapped at cryogenic temperatures. A major problem is the characterization of the early protein and chromophore dynamics, with special emphasis on determining the exact time scale of the C13dC14 (all-trans) f 13-cis isomerization, which is a basic feature of the bR570 photocycle. These ultrafast dynamic features are closely related to key questions such as the causes for the high quantum yield and selectivity of the photoreaction, as well as the molecular mechanism by which light energy is initially stored and subsequently used by the protein in the pump process. Early absorption experiments with picosecond time resolution4 indicated that K610 is preceded by a further red-shifted precursor denoted5 as J625. An analysis5 of the first subpicosecond experiments of Ippen et al.6 indicated that J625 is not the primary photointermediate, suggesting that it grows-in over a 100
19
2-4
19
Rh(Oct)
700 nm absorption range of J625 (features 2 and 3) and those of the decay of I460 at ∼460 and ∼900 nm. Moreover, both features 2 and 3 are assumed to reflect a substantial trans f 13-cis tortional motion. Both assumptions are invalidated by our present data (vide supra). (b) As outlined above, our data do not exclude a mechanism involving two initial molecular pathways, initiating from the S1FC state. This will not be consistent with the Class 1 model of Figure 9. (c) Most Class 1 and Class 2 (but not Class 3) models contradict our major conclusion, namely, that I460 cannot be substantially twisted around C13dC14. (d) The conclusion that the molecular changes mirrored in the ∼500 fs decay of I460 must require substantial torsional motion around C13dC14 contradicts the Class 3 model of Figure 9, in which I460 is on the nonreactive pathway. Our present data do not lead to an unequivocal mechanistic alternative to the above models, but they do provide feasible options. Thus, the suggestion of two distinct pathways (see above), would imply a Class 2 mechanism in which two I states (II and III) are generated from S1FC. However, in variance with the Class 2 models of Figure 9, neither II nor III can represent a C13dC14 isomerized chromophore. This conclusion, which is valid independently of any specific single or double pathway mechanism, implies that other coordinates associated with the chromphore and/or with the surrounding protein matrix should play a role in the early events in bR, hR, and possibly visual pigments. This is in keeping with the Fourier transform analysis of Akyama et al. of the optical absorption spectra of bR570 and bR5.12.28 However, neither the latter study nor our present experimental data can account on a molecular level for the nature of these early processes. A study of the primary reorganization of the excited state, also accounting for the excitation wavelength dependence of the bR spontaneous fluorescence, was recently carried out by Wexler et al.34 A multi-mode analysis identified vibrational degrees of freedom of the retinal chromophore, as well as protein (dielectric) relaxation coordinates, which might be associated with the primary relaxation out of the FranckCondon state. Specific information in this respect was provided by the general ab initio model for protonated Schiff bases of Garavelli et al.,29 which leads to the conclusion that, in retinal proteins, the initial excited-state motion out of the S1FC state is dominated by CdC stretching, rather than by twisting modes. This relaxation results in elongation of a central double bond, associated with change in bond order of the excited state. A recent time-resolved resonance Raman study in bR57030 is indeed indicative of changes in the CdC stretching frequency in the fluorescent state. Such a process may well be associated with the formation of the I460 states. Evidently, if branching into reactive and nonreactive pathways occurs from S1FC, it implies that additional relaxation of the polyene and/or the surrounding protein, should play a role in these early stages of the photoreaction. In other words, changes along one or more additional coordinates are required for differentiating between different fluorescent states such as II and III. It is important to emphasize that while our data are compatible with a Class 2 mechanism in which the (nonisomerized) I460 t II state is on
0.2
τ2
2-4
τ3
τ4
ref
the reactive pathway, we cannot exclude the alternative approach19,20 in which I460 is on the nonreactive pathway. In such a case I460 t III, and II will be identified with the 40-80 fs decay in the 520-560 nm range. Discrimination between these two alternatives will be the subject of future studies. We finally consider the “phase 3” features which in the native system has been interpreted in some models as representing the growing-in of J625 in several hundreds femtoseconds. In our previous report it was suggested22 that an analogous “J-like” species, denoted as T5.12, was formed in the photolysis of bR5.12. In fact, in view of the close similarity in the dynamic features at 660 nm between the native and locked molecules (Figures 6 and 8c), it is tempting to consider that even in the locked systems a species akin to J625 is formed. This would imply that J625 still retains the original all-trans configuration. However, in the context of the broad spectral coverage of the current results, the similarities of appearance in the data collected at 660 nm may be coincidental. In the locked species, excess absorption at 660 nm decays concurrently with I5.12, while in the native system, J625 persists after the decay of I460. Moreover, in the locked system, we do not observe an absorbance decay that would match the 300-500 fs evolution at 660 nm. Thus, an alternative plausible assignment of T5.12 is that it reflects a red band of the I5.12 state, as suggested by the analyses of refs 14-16. This raises the question of why the slower rise in absorption near 660 nm in bR5.12 has no observable counterpart in the near-IR and blue bands assigned to I5.12. A plausible explanation is that this very slight rise in OD is magnified in appearance due to the proposed cancellation of emission and absorption in this spectral range. A conclusive answer will have to await further investigation. V. Conclusions In the present study we have addressed the early light-induced events in bacteriorhodopsin. We have shown that the primary dynamic process following light absorption does not involve the C13dC14 torsional coordinate and that other chromophore and/or protein coordinates must be involved in the very early branching from the Franck-Condon excited state into reactive and nonreactive pathways. Establishing the nature of such pathways which is required from an identification of early species such as I460, and J625 is still an open question, and should be the subject of future investigations. These conclusions most probably also apply to the closely related chloride pump halorhodopsin and possibly to visual pigments. Acknowledgment. This research was supported by grants from the Israel National Science Foundation, founded by the Israel Academy of Sciences and Humanities, Centers of Excellence Program, and from the James-Franck and L. Farkas (HUMinerva) Research Centers. References and Notes (1) The Photophysics and Photochemistry of Retinal Proteins; Ottolenghi, M., Sheves, M., Eds.; Isr. J. Chem. 1995, 35 (3/4), pp 193-515.
5130 J. Phys. Chem. B, Vol. 103, No. 24, 1999 (2) Kochendoerfer, G. G.; Mathies, R. A. Isr. J. Chem. 1995, 35, 211226. (3) Stuart, J. A.; Birge, R. R. Biomembranes 1996, 2A, 33-139. (4) Applebury, M. L.; Peters, K. S.; Rentzepis, P. M. Biophys. J. 1978, 23, 375-382. (5) (a) Dinur, U.; Honig, B.; Ottolenghi, M. In DeVelopments in Biophysical Research; Borsellino, A., Ed.; Plenum: New YorksLondon, 1980; pp. 209-221. (b) Dinur, U.; Honig, B.; Ottolenghi, M. Photochem. Photobiol. 1981, 33, 523-527. (6) Ippen, E. P.; Shank, A.; Lewis, A.; Marcus, M. A. Science 1978, 200, 1279-1281. (7) (a) Polland, H. J.; Zinth, W.; Kaiser, W. In Ultrafast Phenomena; Auston, D. H., Esienthal, K. B., Eds.; Springer: Berlin, 1984; pp 456459. (b) Polland, H. J.; Franz, M. A.; Zinth, W.; Kaiser, W.; Kolling, E.; Oesterhelt, D. Biochim. Biophys. Acta 1984, 767, 635-639. (8) Sharkov, A. V.; Pakulev, A. V.; Chekalin, S. V.; Matveets, YV. A. Biochim. Biophys. Acta 1985, 808, 94-102. (9) Nuss, M. C.; Zinth, W.; Kaiser, W.; Koling, E.; Oesterhelt, D. Chem. Phys. Lett. 1986, 117, 1-7. (10) Petrich, J. W.; Breton, J.; Martin, J. L. In Primary Processes in Photobiology; Kobayashi, T., Ed.; Springer-Verlag: Berlin, 1986, pp 5260. (11) Mathies, R.; Brito Cruz, C. H., Pollard, T. W.; Shank, C. V. Science 1988, 240, 777-779. (12) Dobler, J.; Zinth, W.; Kaiser, K.; Oesterhelt, D. Chem. Phys. 1988, 144, 215-220. (13) Kobayashi, T.; Terauchi, M.; Kouyama, T.; Yoshizawa, M.; Taiji, M. SPIE 1990, 1403, 407-416. (14) Haran, G.; Wynne, K.; Xie, A.; He, Q.; Chance, M.; H, R. M. Chem. Phys. Lett. 1996, 26, 389-395. (15) Hasson, K. C.; Gai, F.; Anfinrud, P. A. Proc. Natl. Acad. Sci. U.S.A. 1996, 93, 15124-15129. (16) Gai, F.; McDonald, P. A.; Anfinrud, P. A. J. Am. Chem. Soc. 1997, 119, 6201-6202. (17) Rosenfeld, T.; Honig, B.; Ottolenghi, M.; Hurley, J.; Ebrey, T. G. Pure Appl. Chem. 1977, 49, 341-351. (18) Hurley, J.; Ebrey, T. G.; Honig, B.; Ottolenghi, M. Nature (London) 1977, 270, 540-542.
Ye et al. (19) Kobayashi, T.; Kim, M.; Taiji, M.; Iwasa, T.; Nakagawa, M.; Tsuda, M. J. Phys. Chem. B. 1998, 102, 272-280. (20) Gai, F.; Hasson, K. C.; McDonald Kooper, J.; Anfinrud, P. A. Science 1998, 279, 1886-1891. (21) Arlt, T.; Schmidt, S.; Zinth, W.; Haupts, U.; Oesterhelt, D. Chem. Phys. Lett. 1995, 241, 559-565. (22) Delaney, J. K.; Brack, T. L.; Atkinson, G. H.; Ottolenghi, M.; Steinberg, G.; Sheves, M. Proc. Natl. Acad. Sci. U.S.A., 1995, 92, 21012105. (23) Zhong, Q.; Ruhman, S.; Ottolenghi, M.; Sheves, M.; Friedman, N.; Atkinson, G. H.; Delaney, J. K. J. Am. Chem. Soc. 1996, 118, 1282812829. (24) Waldman, A.; Banin, U.; Rabani, E.; Ruhman, S. J. Phys. Chem. 1992, 96, 10824. (25) Polland, H. J.; Franz, M. A.; Zinth, W.; Kaiser, W.; Hegemann, P.; Oesterheld, D. Biophys. J. 1985, 47, 55-59. (26) Kandori, H.; Yoshihara, K.; Tomioka, H.; Sasabe, H. J. Phys. Chem. 1992, 96, 6066-6071. (27) Kandori, H.; Yoshihara, K.; Tomioka, H.; Sasabe, H.; Schichida, Y. Chem. Phys. Lett. 1993, 211, 385-391. (28) Akiyama, R.; Yoshimori, A.; Kakitani, T.; Imamoto, Y.; Schichida, Y.; Hatano, Y. J. Phys. Chem. A 1997, 101, 417. (29) Garavelli, M.; Celani, P.; Bernardi, F.; Robb, M. A.; Olivucci, M. J. Am. Chem. Soc. 1997, 119, 6891-6901. (30) Song, L.; El-Sayed, M. A. J. Am. Chem. Soc. 1998, 120, 88898890. (31) Xu, D.; Martin, C.; Schulten, K. Biophys. J. 1996, 70, 453-458. (32) Fang, J.; Carriker, J.; Balogh-Nair, V.; Nakanishi, K. J. Am. Chem. Soc. 1983, 105, 5162-5164. (33) Rousso, I.; Friedman, N.; Sheves, M.; Ottolenghi, M. Biochemistry 1995, 34, 12059-12065. (34) Wexler, G; Kochendoerfer, G.; Mathies, R. A. In Femtochemistry and Phemptobiology: Ultrafast Reaction Dynamics at Atomic Scale Resolution, Nobel Symposium 101; Sundrosto¨m, V., Ed.; Imperial College Press: London, 1997; pp. 724-752. (35) Du, M.; Fleming, G. Biophys. Chem. 1993, 48, 101-111. (36) Freedman, K.; Becker, R. S. J. Am. Chem. Soc. 1986, 108, 12451251.