Chemical dynamics in proteins: the photoisomerization of retinal in bacteriorhodopsin

Photoexcited breathers in conjugated polyenes: an excited-state molecular dynamics study

pi-conjugated polymers have become an important class of materials for electronic devices. Design of these devices requires understanding such processes as photochemical reactions, spatial dynamics of photoexcitations, and energy and charge transport, which in turn involve complex coupled electron-vibrational dynamics. Here we study nonlinear photoexcitation dynamics in the polyene oligomers by using a quantum-chemical method suitable for the simulation of excited-state molecular dynamics in extended molecular systems with sizes up to hundreds of atoms. The method is based on the adiabatic propagation of the ground-state and transition single-electron density matrices along the trajectory. The simulations reveal formation of a self-localized vibronic excitation ("breather" or multiquanta bound state) with a typical period of 34 fs and allows us to identify specific slow and fast nuclear motions strongly coupled to the electronic degrees of freedom. The effect of chain imperfections and chemical defects on the dynamics is also investigated. A complementary two-dimensional analysis of corresponding transition density matrices provides an efficient way to monitor time-dependent real-space localization of the photoexcitation by identifying the underlying changes in charge densities and bond orders. Possible correlated electronic and vibrational spectroscopic signatures of photoexcited breathers are predicted, and generalizations to energy localization in complex macromolecules are discussed.

Comparison of the dynamics of the primary events of bacteriorhodopsin in its trimeric and monomeric states

In this paper, femtosecond pump-probe spectroscopy in the visible region of the spectrum has been used to examine the ultrafast dynamics of the retinal excited state in both the native trimeric state and the monomeric state of bacteriorhodopsin (bR). It is found that the excited state lifetime (probed at 490 nm) increases only slightly upon the monomerization of bR. No significant kinetic difference is observed in the recovery process of the bR ground state probed at 570 nm nor in the fluorescent state observed at 850 nm. However, an increase in the relative amplitude of the slow component of bR excited state decay is observed in the monomer, which is due to the increase in the concentration of the 13-cis retinal isomer in the ground state of the light-adapted bR monomer. Our data indicate that when the protein packing around the retinal is changed upon bR monomerization, there is only a subtle change in the retinal potential surface, which is dependent on the charge distribution and the dipoles within the retinal-binding cavity. In addition, our results show that 40% of the excited state bR molecules return to the ground state on three different time scales: one-half-picosecond component during the relaxation of the excited state and the formation of the J intermediate, a 3-ps component as the J changes to the K intermediate where retinal photoisomerization occurs, and a subnanosecond component during the photocycle.

Conformational dynamics of photoexcited conjugated molecules

Time-dependent photoexcitation and optical spectroscopy of pi-conjugated molecules is described using a new method for the simulation of excited state molecular dynamics in extended molecular systems with sizes up to hundreds of atoms. Applications are made to poly(p-phenylene vinylene) oligomers. Our analysis shows self-trapping of excitations on about six repeat units in the course of photoexcitation relaxation, identifies specific slow (torsion) and fast (bond-stretch) nuclear motions strongly coupled to the electronic degrees of freedom, and predicts spectroscopic signatures of molecular conformations.

Femtosecond infrared spectroscopy of bacteriorhodopsin chromophore isomerization

The vibrational dynamics of the retinal chromophore all-trans-to-13-cis photoisomerization in bacteriorhodopsin has been studied with mid-infrared absorption spectroscopy at high time resolution (about 200 femtoseconds). After photoexcitation of light-adapted bacteriorhodopsin, the transient infrared absorption was probed in a broad spectral region, including vibrations with dominant C-C, C=C, and C=NH stretching mode amplitude. All photoproduct modes, especially those around 1190 reciprocal-centimeters that are indicative for a 13-cis configuration of the chromophore, rise with a time constant of approximately 0.5 picosecond. The results presented give direct vibrational-spectroscopic evidence for the isomerization taking place within 0.5 picosecond, as has been suggested by previous optical femtosecond time-resolved experiments but questioned recently by picosecond time-resolved vibrational spectroscopy experiments.

Femtosecond and picosecond fluorescence of native bacteriorhodopsin and a nonisomerizing analog

The spectrally and temporally resolved fluorescence properties of native bacteriorhodopsin (bR) and bR reconstituted with a nonisomerizing analog of the retinal Schiff base (bR5.12) are examined. The first attempt to experimentally monitor the excited state relaxation processes in both type of pigments using ultrafast fluorescence spectroscopy is reported. The fluorescence is emitted from retinal molecules in an all-trans configuration. Substantial energy relaxation involves very fast intramolecular and intermolecular vibrational modes and these are shown to occur on a time scale faster than isomerization. The possible contribution of dielectric interaction between the retinal Schiff base and the protein environment for the excited state energy relaxation is discussed.

A semiempirical study of the optimized ground and excited state potential energy surfaces of retinal and its protonated Schiff base

The electronic ground and first excited states of retinal and its Schiff base are optimized for the first time using the semiempirical AM1 Hamiltonian. The barrier for rotation about the C(11)-C(12) double bond is characterized by variation of both the twist angle delta(C(10)-C(11)-C(12)-C(13)) and the bond length d(C(11)-C(12)). The potential energy surface is obtained by varying these two parameters. The calculated ground state rotational barrier is equal to 15.6 kcal/mol for retinal and 20.5 kcal/mol for its Schiff base. The all-trans conformation is more stable by 3.7 kcal/mol than the 11-cis geometry. For the first excited state, S(1,) the 90 degrees twisted geometry represents a saddle point for retinal with the rotational barrier of 14.6 kcal/mol. In contrast, this conformation is an energy minimum for the Schiff base. It can be easily reached at room temperature from the planar minima since it is separated from them by a barrier of only 0.6 kcal/mol. The 90 degrees minimum conformation is more stable than the all-trans by 8.6 kcal/mol. We are thus able to present a reaction path on the S(1) surface of the retinal Schiff base with an almost barrier-less geometrical relaxation into a twisted minimum geometry, as observed experimentally. The character of the ground and first excited singlet states underscores the need for the inclusion of double excitations in the calculations.

Real-time spectroscopy of transition states in bacteriorhodopsin during retinal isomerization

Real-time investigations of the rearrangement of bonds during chemical transformations require femtosecond temporal resolution, so that the atomic vibrations within the reacting molecules can be observed. Following the development of lasers capable of emitting ultrashort laser flashes on this timescale, chemical reactions involving relatively simple molecules have been monitored in detail, revealing the transient existence of intermediate species as reactants are transformed into products. Here we report the direct observation of nuclear motion in a complex biological system, the retinal chromophore of bacteriorhodopsin (bR568), as it undergoes the trans-cis photoisomerization that is fundamental to the vision process. By using visible-light pulses of less than 5 femtosecond in duration, we are able to monitor changes in the vibrational spectra of the transition state and thus show that despite photoexcitation of the anti-bonding molecular orbital involved, isomerization does not occur instantly, but involves transient formation of a so-called 'tumbling state'. Our observations thus agree with growing experimental and ab initio evidence for a three-state photoisomerization model and firmly discount the initially suggested two-state model for this process.

The relaxation dynamics of the excited electronic states of retinal in bacteriorhodopsin by two-pump-probe femtosecond studies

We present the results of two-pump and probe femtosecond experiments designed to follow the relaxation dynamics of the lowest excited state (S(1)) populated by different modes. In the first mode, a direct (S(0) --> S(1)) radiative excitation of the ground state is used. In the second mode, an indirect excitation is used where the S(1) state is populated by the use of two femtosecond laser pulses with different colors and delay times between them. The first pulse excites the S(0) --> S(1) transition whereas the second pulse excites the S(1) --> S(n) transition. The nonradiative relaxation from the S(n) state populates the lowest excited state. Our results suggest that the S(1) state relaxes faster when populated nonradiatively from the S(n) state than when pumped directly by the S(0) --> S(1) excitation. Additionally, the S(n) --> S(1) nonradiative relaxation time is found to change by varying the delay time between the two pump pulses. The observed dependence of the lowest excited state population as well as its dependence on the delay between the two pump pulses are found to fit a kinetic model in which the S(n) state populates a different surface (called S'(1)) than the one being directly excited (S(1)). The possible involvement of the A(g) type states, the J intermediate, and the conical intersection leading to the S(0) or to the isomerization product (K intermediate) are discussed in the framework of the proposed model.

An unusual pathway of excitation energy deactivation in carotenoids: singlet-to-triplet conversion on an ultrafast timescale in a photosynthetic antenna

Carotenoids are important biomolecules that are ubiquitous in nature and find widespread application in medicine. In photosynthesis, they have a large role in light harvesting (LH) and photoprotection. They exert their LH function by donating their excited singlet state to nearby (bacterio)chlorophyll molecules. In photosynthetic bacteria, the efficiency of this energy transfer process can be as low as 30%. Here, we present evidence that an unusual pathway of excited state relaxation in carotenoids underlies this poor LH function, by which carotenoid triplet states are generated directly from carotenoid singlet states. This pathway, operative on a femtosecond and picosecond timescale, involves an intermediate state, which we identify as a new, hitherto uncharacterized carotenoid singlet excited state. In LH complex-bound carotenoids, this state is the precursor on the reaction pathway to the triplet state, whereas in extracted carotenoids in solution, this state returns to the singlet ground state without forming any triplets. We discuss the possible identity of this excited state and argue that fission of the singlet state into a pair of triplet states on individual carotenoid molecules constitutes the mechanism by which the triplets are generated. This is, to our knowledge, the first ever direct observation of a singlet-to-triplet conversion process on an ultrafast timescale in a photosynthetic antenna.

Computational evidence in favor of a two-state, two-mode model of the retinal chromophore photoisomerization

In this paper we use ab initio multiconfigurational second-order perturbation theory to establish the intrinsic photoisomerization path model of retinal chromophores. This is accomplished by computing the ground state (S(0)) and the first two singlet excited-state (S(1), S(2)) energies along the rigorously determined photoisomerization coordinate of the rhodopsin chromophore model 4-cis-gamma-methylnona-2,4,6,8-tetraeniminium cation and the bacteriorhodopsin chromophore model all-trans-hepta-2,4, 6-trieniminium cation in isolated conditions. The computed S(2) and S(1) energy profiles do not show any avoided crossing feature along the S(1) reaction path and maintain an energy gap >20 kcal small middle dotmol(-1). In addition, the analysis of the charge distribution shows that there is no qualitative change in the S(2) and S(1) electronic structure along the path. Thus, the S(1) state maintains a prevalent ionic (hole-pair) character whereas the S(2) state maintains a covalent (dot-dot) character. These results, together with the analysis of the S(1) reaction coordinate, support a two-state, two-mode model of the photoisomerization that constitutes a substantial revision of the previously proposed models.

Proton transfer reactions across bacteriorhodopsin and along the membrane

Bacteriorhodopsin is probably the best understood proton pump so far and is considered to be a model system for proton translocating membrane proteins. The basis of a molecular description of proton translocation is set by having the luxury of six highly resolved structural models at hand. Details of the mechanism and reaction dynamics were elucidated by a whole variety of biophysical techniques. The current molecular picture of catalysis by BR will be presented with examples from time-resolved spectroscopy. FT-IR spectroscopy monitors single proton transfer events within bacteriorhodopsin and judiciously positioned pH indicators detect proton migration at the membrane surface. Emerging properties are briefly outlined that underlie the efficient proton transfer across and along biological membranes.

Protein-assisted pericyclic reactions: an alternate hypothesis for the action of quantal receptors

The rules for allowable pericyclic reactions indicate that the photoisomerizations of retinals in rhodopsins can be formally analogous to thermally promoted Diels-Alder condensations of monoenes with retinols. With little change in the seven-transmembrane helical environment these latter reactions could mimic the retinal isomerization while providing highly sensitive chemical reception. In this way archaic progenitors of G-protein-coupled chemical quantal receptors such as those for pheromones might have been evolutionarily plagiarized from the photon quantal receptor, rhodopsin, or vice versa. We investigated whether the known structure of bacteriorhodopsin exhibited any similarity in its active site with those of the two known antibody catalysts of Diels-Alder reactions and that of the photoactive yellow protein. A remarkable three-dimensional motif of aromatic side chains emerged in all four proteins despite the drastic differences in backbone structure. Molecular orbital calculations supported the possibility of transient pericyclic reactions as part of the isomerization-signal transduction mechanisms in both bacteriorhodopsin and the photoactive yellow protein. It appears that reactions in all four of the proteins investigated may be biological analogs of the organic chemists' chiral auxiliary-aided Diels-Alder reactions. Thus the light receptor and the chemical receptor subfamilies of the heptahelical receptor family may have been unified at one time by underlying pericyclic chemistry.

Femtochemistry

The topic of femtochemistry is surveyed from both theoretical and experimental points of view. A time-dependent wave packet description of the photodissociation of the O---C---S molecule reveals vibrational motion in the transition-state region and suggests targets for direct experimental observation. Theoretical approaches for treating femtosecond chemical phenomena in condensed phases are featured along with prospects for laser-controlled chemical reactions by using tailored ultrashort chirped pulses. An experimental study of the photoisomerization of retinal in the protein bacteriorhodopsin is discussed with an aim to gain insight into the potential energy surfaces on which this remarkably efficient and selective reactions proceeds. Finally, a prospective view of new frontiers in femtochemistry is given.

Closing in on bacteriorhodopsin: progress in understanding the molecule

Bacteriorhodopsin is the best understood ion transport protein and has become a paradigm for membrane proteins in general and transporters in particular. Models up to 2.5 A resolution of bacteriorhodopsin's structure have been published during the last three years and are basic for understanding its function. Thus one focus of this review is to summarize and to compare these models in detail. Another focus is to follow the protein through its catalytic cycle in summarizing more recent developments. We focus on literature published since 1995; a comprehensive series of reviews was published in 1995 (112).

Shedding light on the dark and weakly fluorescent states of green fluorescent proteins

Recent experiments on various similar green fluorescent protein (GFP) mutants at the single-molecule level and in solution provide evidence of previously unknown short- and long-lived "dark" states and of related excited-state decay channels. Here, we present quantum chemical calculations on cis-trans photoisomerization paths of neutral, anionic, and zwitterionic GFP chromophores in their ground and first singlet excited states that explain the observed behaviors from a common perspective. The results suggest that favorable radiationless decay channels can exist for the different protonation states along these isomerizations, which apparently proceed via conical intersections. These channels are suggested to rationalize the observed dramatic reduction of fluorescence in solution. The observed single-molecule fast blinking is attributed to conversions between the fluorescent anionic and the dark zwitterionic forms whereas slow switching is attributed to conversions between the anionic and the neutral forms. The predicted nonadiabatic crossings are seen to rationalize the origins of a variety of experimental observations on a common basis and may have broad implications for photobiophysical mechanisms in GFP.

Three electronic state model of the primary phototransformation of bacteriorhodopsin

The primary all-trans --> 13-cis photoisomerization of retinal in bacteriorhodopsin has been investigated by means of quantum chemical and combined classical/quantum mechanical simulations employing the density matrix evolution method. Ab initio calculations on an analog of a protonated Schiff base of retinal in vacuo reveal two excited states S1 and S2, the potential surfaces of which intersect along the reaction coordinate through an avoided crossing, and then exhibit a second, weakly avoided, crossing or a conical intersection with the ground state surface. The dynamics governed by the three potential surfaces, scaled to match the in situ level spacings and represented through analytical functions, are described by a combined classical/quantum mechanical simulation. For a choice of nonadiabatic coupling constants close to the quantum chemistry calculation results, the simulations reproduce the observed photoisomerization quantum yield and predict the time needed to pass the avoided crossing region between S1 and S2 states at tau1 = 330 fs and the S1 --> ground state crossing at tau2 = 460 fs after light absorption. The first crossing follows after a 30 degrees torsion on a flat S1 surface, and the second crossing follows after a rapid torsion by a further 60 degrees. tau1 matches the observed fluorescence lifetime of S1. Adjusting the three energy levels to the spectral shift of D85N and D212N mutants of bacteriorhodospin changes the crossing region of S1 and S2 and leads to an increase in tau1 by factors 17 and 10, respectively, in qualitative agreement with the observed increase in fluorescent lifetimes.

The structure and mechanism of the family of retinal proteins from halophilic archaea

Retinal proteins from halophilic archaea provide a unique opportunity to analyze vectorial ion translocation. Studies on its structure, conformational changes, proton conduction and electrogenic steps have helped to elucidate the catalytic cycle of bacteriorhodopsin in increasing detail. Experimental modulation of the vectoriality and ion specificity by altering the substrate availability, point mutations and light conditions for the different retinal proteins allows the proposal of a general model of ion transport for this protein family.