Discriminating Clockwise and Counterclockwise Photoisomerization Paths in Achiral Photoswitches by Excited-State Electronic Circular Dichroism

Despite the numerous investigations of photoisomerization reactions from both the computational and experimental points of view, even in complex environments, to date there is no direct demonstration of the direction of rotation of the retinal chromophore, initiating the vision process in several organisms, occurring upon light irradiation. In the literature, many proposals have been formulated to shed light on the details of this process, most of which are extracted from semiclassical simulations. Although high hopes are held in the development of time-resolved X-ray spectroscopy, I argue in this work that simpler but less known techniques can be used to unravel the details of this fascinating photochemical process. In fact, chiroptical spectroscopy would unambiguously prove the direction of the rotatory motion of the chromophore during the photoisomerization process by probing excited state chirality, a piece of information that, so far, has been exclusively extracted from atomistic simulations. I demonstrate this statement by computing the expected chiroptical response along photoisomerization pathways for several models of the retinal chromophores that are found in nature bound to rhodopsins, including nuclear ensemble spectra from semiclassical dynamics simulations, that can be compared with time-resolved experiments.

Bidirectional Photochemistry of Antarctic Microbial Rhodopsin: Emerging Trend of Ballistic Photoisomerization from the 13-cis Resting State

The decades-long ultrafast examination of nearly a dozen microbial retinal proteins, ion pumps, and sensory photoreceptors has not identified structure-function indicators which predict photoisomerization dynamics, whether it will be sub-picosecond and ballistic or drawn out with complex curve-crossing kinetics. Herein, we report the emergence of such an indicator. Using pH control over retinal isomer ratios, photoinduced transient absorption is recorded in an inward proton pumping Antarctic microbial rhodopsin (AntR) for 13-cis and all-trans retinal resting states. The all-trans fluorescent state decays with 1 ps exponential kinetics. In contrast, in 13-cis it decays within ∼300 fs accompanied by continuous spectral evolution, indicating ballistic internal conversion. The coherent wave packet nature of 13-cis isomerization in AntR matches published results for bacteriorhodopsin (BR) and Anabaena sensory rhodopsin (ASR), which also accommodate both all-trans and 13-cis retinal resting states, marking the emergence of a first structure-photodynamics indicator which holds for all three tested pigments.

Rhodopsins: An Excitingly Versatile Protein Species for Research, Development and Creative Engineering

The first member and eponym of the rhodopsin family was identified in the 1930s as the visual pigment of the rod photoreceptor cell in the animal retina. It was found to be a membrane protein, owing its photosensitivity to the presence of a covalently bound chromophoric group. This group, derived from vitamin A, was appropriately dubbed retinal. In the 1970s a microbial counterpart of this species was discovered in an archaeon, being a membrane protein also harbouring retinal as a chromophore, and named bacteriorhodopsin. Since their discovery a photogenic panorama unfolded, where up to date new members and subspecies with a variety of light-driven functionality have been added to this family. The animal branch, meanwhile categorized as type-2 rhodopsins, turned out to form a large subclass in the superfamily of G protein-coupled receptors and are essential to multiple elements of light-dependent animal sensory physiology. The microbial branch, the type-1 rhodopsins, largely function as light-driven ion pumps or channels, but also contain sensory-active and enzyme-sustaining subspecies. In this review we will follow the development of this exciting membrane protein panorama in a representative number of highlights and will present a prospect of their extraordinary future potential.

Torsionally broken symmetry assists infrared excitation of biomimetic charge-coupled nuclear motions in the electronic ground state

The concerted interplay between reactive nuclear and electronic motions in molecules actuates chemistry. Here, we demonstrate that out-of-plane torsional deformation and vibrational excitation of stretching motions in the electronic ground...

Velocity Adjustment in Surface Hopping: Ethylene as a Case Study of the Maximum Error Caused by Direction Choice

The most common surface hopping dynamics algorithms require velocity adjustment after hopping to ensure total-energy conservation. Based on the semiclassical analysis, this adjustment must be made parallel to the nonadiabatic coupling vector's direction. Nevertheless, this direction is not always known, and the common practice has been to adjust the velocity in either the linear momentum or velocity directions. This paper benchmarks surface hopping dynamics of photoexcited ethylene with velocity adjustment in several directions, including those of the nonadiabatic coupling vector, the momentum, and the energy gradient difference. It is shown that differences in time constants and structural evolution fall within the statistical uncertainty of the method considering up to 500 trajectories in each dynamics set, rendering the three approaches statistically equivalent. For larger ensembles beyond 1000 trajectories, significant differences between the results arise, limiting the validity of adjustment in alternative directions. Other possible adjustment directions (velocity, single-state gradients, angular momentum) are evaluated as well. Given the small size of ethylene, the results reported in this paper should be considered an upper limit for the error caused by the choice of the velocity-adjustment direction on surface hopping dynamics.

Probing Structure and Reaction Dynamics of Proteins Using Time-Resolved Resonance Raman Spectroscopy

The mechanistic understanding of protein functions requires insight into the structural and reaction dynamics. To elucidate these processes, a variety of experimental approaches are employed. Among them, time-resolved (TR) resonance Raman (RR) is a particularly versatile tool to probe processes of proteins harboring cofactors with electronic transitions in the visible range, such as retinal or heme proteins. TR RR spectroscopy offers the advantage of simultaneously providing molecular structure and kinetic information. The various TR RR spectroscopic methods can cover a wide dynamic range down to the femtosecond time regime and have been employed in monitoring photoinduced reaction cascades, ligand binding and dissociation, electron transfer, enzymatic reactions, and protein un- and refolding. In this account, we review the achievements of TR RR spectroscopy of nearly 50 years of research in this field, which also illustrates how the role of TR RR spectroscopy in molecular life science has changed from the beginning until now. We outline the various methodological approaches and developments and point out current limitations and potential perspectives.

Mapping the ultrafast vibrational dynamics of all-trans and 13-cis retinal isomerization in Anabaena Sensory Rhodopsin

Discrepancies in the isomerization dynamics and quantum yields of the trans and cis retinal protonated Schiff base is a well-known issue in the context of retinal photochemistry. Anabaena Sensory Rhodopsin (ASR) is a microbial retinal protein that comprises a retinal chromophore in two ground state (GS) conformations: all-trans, 15-anti (AT) and 13-cis, 15-syn (13C). In this study, we applied impulsive vibrational spectroscopic techniques (DFWM, pump-DFWM and pump-IVS) to ASR to shed more light on how the structural changes take place in the excited state within the same protein environment. Our findings point to distinct features in the ground state structural conformations as well as to drastically different evolutions in the excited state manifold. The ground state vibrational spectra show stronger Raman activity of the C14-H out-of-plane wag (at about 805 cm-1) for the 13C isomer than that for the AT isomer, which hints at a pre-distortion of 13C in the ground state. Evolution of the Raman frequency after interaction with the actinic pulse shows a blue-shift for the C[double bond, length as m-dash]C stretching and CH3 rocking mode for both isomers. For AT, however, the blue-shift is not instantaneous as observed for the 13C isomer, rather it takes more than 200 fs to reach the maximum frequency shift. This frequency blue-shift is rationalized by a decrease in the effective conjugation length during the isomerization reaction, which further confirms a slower formation of the twisted state for the AT isomer and corroborates the presence of a barrier in the excited state trajectory previously predicted by quantum chemical calculations.

Two-dimensional electronic spectroscopy as a tool for tracking molecular conformations in DNA/RNA aggregates

A computational strategy to simulate two-dimensional electronic spectra (2DES) is introduced, which allows characterising ground state conformations of flexible nucleobase aggregates that play a crucial role in nucleic acid photochemistry.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

Coupled HOOP signature correlates with quantum yield of isorhodopsin and analog pigments

With a quantum yield of 0.66±0.03 the photoisomerization efficiency of the visual pigment rhodopsin (11-cis⇒all-trans chromophore) is exceptionally high. This is currently explained by coherent coupling of the excited state electronic wavepacket with local vibrational nuclear modes, facilitating efficient cross-over at a conical intersection onto the photoproduct energy surface. The 9-cis counterpart of rhodopsin, dubbed isorhodopsin, has a much lower quantum yield (0.26±0.03), which, however, can be markedly enhanced by modification of the retinal chromophore (7,8-dihydro and 9-cyclopropyl derivatives). The coherent coupling in the excited state is promoted by torsional skeletal and coupled HOOP vibrational modes, in combination with a twisted conformation around the isomerization region. Since such torsion will strongly enhance the infrared intensity of coupled HOOP modes, we investigated FTIR difference spectra of rhodopsin, isorhodopsin and several analog pigments in the spectral range of isolated and coupled HCCH wags. As a result we propose that the coupled HOOP signature in these retinal pigments correlates with the distribution of torsion over counteracting segments in the retinylidene polyene chain. As such the HOOP signature can act as an indicator for the photoisomerization efficiency, and can explain the higher quantum yield of the 7,8-dihydro and 9-cyclopropyl-isorhodopsin analogs.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Vibrational relaxation as the driving force for wavelength conversion in the peridinin-chlorophyll a-protein

We present a computationally derived energy transfer model for the peridinin-chlorophyll a-protein (PCP), which invokes vibrational relaxation in the two lowest singlet excited states rather than internal conversion between them. The model allows an understanding of the photoinduced processes without assuming further electronic states or a dependence of the 2Ag state character on the vibrational sub-state. We report molecular dynamics simulations (CHARMM22 force field) and quantum mechanics/molecular mechanics (QM/MM) calculations on PCP. In the latter, the QM region containing a single peridinin (Per) chromophore or a Per-Chl a (chlorophyll a) pair is treated by density functional theory (DFT, CAM-B3LYP) for geometries and by DFT-based multireference configuration interaction (DFT/MRCI) for excitation energies. The calculations show that Per has a bright, green light absorbing 2Ag state, in addition to the blue light absorbing 1Bu state found in other carotenoids. Both states undergo a strong energy lowering upon relaxation, leading to emission in the red, while absorbing in the blue or green. The orientation of their transition dipole moments indicates that both states are capable of excited-state energy transfer to Chl a, without preference for either 1Bu or 2Ag as donor state. We propose that the commonly postulated partial intramolecular charge transfer (ICT) character of a donating Per state can be assigned to the relaxed 1Bu state, which takes on ICT character. By assuming that both 1Bu and 2Ag are able to donate to the Chl a Q band, one can explain why different chlorophyll species in PCP exhibit different acceptor capabilities.

100 fs photo-isomerization with vibrational coherences but low quantum yield in Anabaena Sensory Rhodopsin

Counter-intuitive photochemistry: in Anabaena Sensory Rhodopsin, the retinal 13-cis isomer isomerizes much faster than all-trans ASR, but with a 3-times lower quantum yield.

Spectral properties and isomerisation path of retinal in C1C2 channelrhodopsin

Computed torsion profiles along the reactive coordinate in S₁reveal a two-path deactivation mechanism for retinal in C1C2 channelrhodopsin.