Chemical dynamics in proteins: the photoisomerization of retinal in bacteriorhodopsin

Mechanism and dynamics of photoswitchable flavoprotein charge-transfer complexes

Because of their diverse uses in biological science and engineering, continued effort has been made to expand the pool of photoswitchable protein systems. A recent study demonstrated that in monomeric sarcosine oxidase (MSOX), photoexcitation of a charge-transfer (CT) complex formed by a flavin cofactor and a nonreactive ligand (e.g., methylthioacetate) induces the ligand to reversibly change conformation, with implications for the development of flavin-dependent fast photochromic proteins. However, the factors that control the underlying switching mechanism and dynamics remain largely unexplored. Herein, combining extensive protein mutagenesis, ultrafast laser spectroscopic measurements and classical and quantum computational approaches, we assess those factors in a range of protein variants, including those of MSOX and another flavoenzyme, N-methyltryptophan oxidase (MTOX), where we find that a similar photoswitching cycle can occur. We demonstrate that (1) the kinetic behaviors of the photoswitching cycle are protein- and ligand-dependent; (2) the photoswitching and backward thermal recovery rates can be tuned by mutation of a specific active-site residue (Met245 and Thr239 in MSOX and MTOX, respectively), with recovery rates spanning over an order of magnitude, and (3) modifications of the protein environment alter the conformational energy landscape of the ligand-flavin complex, consequently regulating the photocycle. Taken together, these findings highlight the versatility of such photoswitchable systems, providing a molecular basis for fine-tuning their photophysical properties.

Mapping photoisomerization dynamics on a three-state model potential energy surface in bacteriorhodopsin using femtosecond stimulated Raman spectroscopy

The process of proton translocation in Halobacterium salinarum, triggered by light, is powered by the photoisomerization of all-trans-retinal in bacteriorhodopsin (bR). The primary events in bR involving rapid structural changes upon light absorption occur within subpicoseconds to picoseconds. While the three-state model has received extensive support in describing the primary events between the H and K states, precise characterization of each excited state in the three-state model during photoisomerization remains elusive. In this study, we investigate the ultrafast structural dynamics of all-trans-retinal in bR using femtosecond stimulated Raman spectroscopy. We report Raman modes at 1820 cm-1 which arise from C[double bond, length as m-dash]C stretch vibronic coupling and provide direct experimental evidence for the involvement of the I and J states with 2A- g symmetric character in the three-state model. The detection of the C[double bond, length as m-dash]C vibronic coupling mode, C[double bond, length as m-dash]N stretching mode (1700 cm-1), and hydrogen out-of-plane (HOOP) mode (954 cm-1) further supports the three-state model that elucidates the initial charge translocation along the conjugated chain accompanied by trans-to-cis photoisomerization dynamics through H(1B+ u) → I(2A- g) → J(2A- g) → K(13-cis ground state) transitions in all-trans-retinal in bR.

Reversible Photochromic Reactions of Bacteriorhodopsin from Halobacterium salinarum at Femto- and Picosecond Times

The operation of bacteriorhodopsin (BR) from the archaeon Halobacterium salinarum is based on the photochromic reaction of isomerization of the chromophore group (the retinal protonated Schiff base, RPSB) from the all-trans to the 13-cis form. The ultrafast dynamics of the reverse 13-cis → all-trans photoreaction was studied using femtosecond transient absorption spectroscopy in comparison with the forward photoreaction. The forward photoreaction was initiated by photoexcitation of BR by pulse I (540 nm). The reverse photoreaction was initiated by photoexcitation of the product K590 at an early stage of its formation (5 ps) by pulse II (660 nm). The conversion of the excited K590 to the ground state proceeds at times of 0.19, 1.1, and 16 ps with the relative contributions of ~20/60/20, respectively. All these decay channels lead to the formation of the initial state of BR as a product with a quantum yield of ~1. This state is preceded by vibrationally excited intermediates, the relaxation of which occurs in the 16 ps time range. Likely, the heterogeneity of the excited state of K590 is determined by the heterogeneity of its chromophore center. The forward photoreaction includes two components-0.52 and 3.5 ps, with the relative contributions of 91/9, respectively. The reverse photoreaction initiated from K590 proceeds more efficiently in the conical intersection (CI) region but on the whole at a lower rate compared to the forward photoreaction, due to significant heterogeneity of the potential energy surface.

Optical control of ultrafast structural dynamics in a fluorescent protein

The photoisomerization reaction of a fluorescent protein chromophore occurs on the ultrafast timescale. The structural dynamics that result from femtosecond optical excitation have contributions from vibrational and electronic processes and from reaction dynamics that involve the crossing through a conical intersection. The creation and progression of the ultrafast structural dynamics strongly depends on optical and molecular parameters. When using X-ray crystallography as a probe of ultrafast dynamics, the origin of the observed nuclear motions is not known. Now, high-resolution pump-probe X-ray crystallography reveals complex sub-ångström, ultrafast motions and hydrogen-bonding rearrangements in the active site of a fluorescent protein. However, we demonstrate that the measured motions are not part of the photoisomerization reaction but instead arise from impulsively driven coherent vibrational processes in the electronic ground state. A coherent-control experiment using a two-colour and two-pulse optical excitation strongly amplifies the X-ray crystallographic difference density, while it fully depletes the photoisomerization process. A coherent control mechanism was tested and confirmed the wave packets assignment.

Similarities and Differences in Photochemistry of Type I and Type II Rhodopsins

The diversity of the retinal-containing proteins (rhodopsins) in nature is extremely large. Fundamental similarity of the structure and photochemical properties unites them into one family. However, there is still a debate about the origin of retinal-containing proteins: divergent or convergent evolution? In this review, based on the results of our own and literature data, a comparative analysis of the similarities and differences in the photoconversion of the rhodopsin of types I and II is carried out. The results of experimental studies of the forward and reverse photoreactions of the bacteriorhodopsin (type I) and visual rhodopsin (type II) rhodopsins in the femto- and picosecond time scale, photo-reversible reaction of the octopus rhodopsin (type II), photovoltaic reactions, as well as quantum chemical calculations of the forward photoreactions of bacteriorhodopsin and visual rhodopsin are presented. The issue of probable convergent evolution of type I and type II rhodopsins is discussed.

Extended Single-Electron Transfer Model and Dynamically Associated Energy Transfer Event in a Dual-Functional Catalyst System

Organic photocatalysis has been developed flourishingly to rely on bimolecular energy transfer (EnT) or oxidative/reductive electron transfer (ET), promoting a variety of synthetic transformations. However, there are rare examples to merge EnT and ET processes rationally within one chemical system, of which the mechanistic investigation still remains in its infancy. Herein, the first mechanistic illustration and kinetic assessments of the dynamically associated EnT and ET paths were conducted for realizing the C-H functionalization in a cascade photochemical transformation of isomerization and cyclization by using the dual-functional organic photocatalyst of riboflavin. An extended single-electron transfer model of transition-state-coupled dual-nonadiabatic crossings was explored to analyze the dynamic behaviors in the proton transfer-coupled cyclization. This can also be used to clarify the dynamic correlation with the EnT-driven E → Z photoisomerization that has been kinetically evaluated by using Fermi's golden rule with the Dexter model. The present computational results of electron structures and kinetic data contribute to a fundamental basis for understanding the photocatalytic mechanism of the combined operation of EnT and ET strategies, which will guide the design and manipulation for the implementation of multiple activation modes based on a single photosensitizer.

From cis-Lobeline to trans-Lobeline: Study on the Pharmacodynamics and Isomerization Factors

Lobeline is an alkaloid derived from the leaves of an Indian tobacco plant (Lobelia inflata), which has been prepared by chemical synthesis. It is classified as a partial nicotinic agonist and has a long history of therapeutic usage ranging from emetic and respiratory stimulant to tobacco smoking cessation agent. The presence of both cis and trans isomers in lobeline is well known, and many studies on the relationship between the structure and pharmacological activity of lobeline and its analogs have been reported. However, it is a remarkable fact that no studies have reported the differences in pharmacological activities between the two isomers. In this article, we found that different degrees of isomerization of lobeline injection have significant differences in respiratory excitatory effects in pentobarbital sodium anesthetized rats. Compared with cis-lobeline injections, the respiratory excitatory effect was significantly reduced by 50.2% after administration of injections which contained 36.9% trans-lobeline. The study on the influencing factors of isomerization between two isomers shown that this isomerization was a one-way isomerism and only converted from cis to trans, where temperature was the catalytic factor and pH was the key factor. This study reports a new discovery. Despite the widespread use of ventilators, first-aid medicines such as nikethamide and lobeline has retired to second line, but as a nonselective antagonist with high affinity for a4b2 and a3b2 nicotinic acetylcholine receptors (nAChRs). In recent years, lobeline has shown great promise as a therapeutic drug for mental addiction and nervous system disorders, such as depression, Alzheimer disease and Parkinson disease. Therefore, we suggest that the differences between two isomers should be concerned in subsequent research papers and applications.

Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency

The activation of rhodopsin, the light-sensitive G-protein coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. Here, we use hundreds of quantum-classical trajectories to show that, 15 femtoseconds after light absorption, a degeneracy between the reactive excited state and a neighboring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show that such splitting is modulated by protein electrostatics, thus linking amino-acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies, would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.

A Unified View on Varied Ultrafast Dynamics of the Primary Process in Microbial Rhodopsins

All-trans to 13-cis photoisomerization of the protonated retinal Schiff base (PRSB) chromophore is the primary step that triggers various biological functions of microbial rhodopsins. While this ultrafast primary process has been extensively studied, it has been recognized that the relevant excited-state relaxation dynamics differ significantly from one rhodopsin to another. To elucidate the origin of the complicated ultrafast dynamics of the primary process in microbial rhodopsins, we studied the excited-state dynamics of proteorhodopsin, its D97N mutant, and bacteriorhodopsin by femtosecond time-resolved absorption (TA) spectroscopy in a wide pH range. The TA data showed that their excited-state relaxation dynamics drastically change when pH approaches the pKa of the counterion residue of the PRSB chromophore in the ground state. This result reveals that the varied excited-state relaxation dynamics in different rhodopsins mainly originate from the difference of the ground-state heterogeneity (i.e., protonation/deprotonation of the PRSB counterion).

Control of Protonated Schiff Base Excited State Decay within Visual Protein Mimics: A Unified Model for Retinal Chromophores

Artificial biomimetic chromophore-protein complexes inspired by natural visual pigments can feature color tunability across the full visible spectrum. However, control of excited state dynamics of the retinal chromophore, which is of paramount importance for technological applications, is lacking due to its complex and subtle photophysics/photochemistry. Here, ultrafast transient absorption spectroscopy and quantum mechanics/molecular mechanics simulations are combined for the study of highly tunable rhodopsin mimics, as compared to retinal chromophores in solution. Conical intersections and transient fluorescent intermediates are identified with atomistic resolution, providing unambiguous assignment of their ultrafast excited state absorption features. The results point out that the electrostatic environment of the chromophore, modified by protein point mutations, affects its excited state properties allowing control of its photophysics with same power of chemical modifications of the chromophore. The complex nature of such fine control is a fundamental knowledge for the design of bio-mimetic opto-electronic and photonic devices.

Comparative Femtosecond Spectroscopy of Primary Photoreactions of Exiguobacterium sibiricum Rhodopsin and Halobacterium salinarum Bacteriorhodopsin

The primary stages of the Exiguobacterium sibiricum rhodopsin (ESR) photocycle were investigated by femtosecond absorption laser spectroscopy in the spectral range of 400-900 nm with a time resolution of 25 fs. The dynamics of the ESR photoreaction were compared with the reactions of bacteriorhodopsin (bR) in purple membranes (bR_PM) and in recombinant form (bR_rec). The primary intermediates of the ESR photocycle were similar to intermediates I, J, and K in bacteriorhodopsin photoconversion. The CONTIN program was applied to analyze the characteristic times of the observed processes and to clarify the reaction scheme. A similar photoreaction pattern was observed for all studied retinal proteins, including two consecutive dynamic Stokes shift phases lasting ∼0.05 and ∼0.15 ps. The excited state decays through a femtosecond reactive pathway, leading to retinal isomerization and formation of product J, and a picosecond nonreactive pathway that leads only to the initial state. Retinal photoisomerization in ESR takes 0.69 ps, compared with 0.48 ps in bR_PM and 0.74 ps in bR_rec. The nonreactive excited state decay takes 5 ps in ESR and ∼3 ps in bR. We discuss the similarity of the primary reactions of ESR and other retinal proteins.

The Origin of Ultrafast Multiphasic Dynamics in Photoisomerization of Bacteriophytochrome

Red-light bacteriophytochromes regulate many physiological functions through photoisomerization of a linear tetrapyrrole chromophore. In this work, we mapped out femtosecond-resolved fluorescence spectra of the excited Pr state and observed unique active-site relaxations on the picosecond time scale with unusual spectral tuning of rises on the blue side and decays on the red side of the emission. We also observed initial wavepacket dynamics in femtoseconds with two low-frequency modes of 38 and 181 cm^-1 as well as the intermediate product formation after isomerization in hundreds of picoseconds. With critical mutations at the active site, we observed similar dynamic patterns with different times for both relaxation and isomerization, consistent with the structural and chemical changes induced by the mutations. The observed multiphasic dynamics clearly represents the active-site relaxation, not different intermediate reactions or excitation of heterogeneous ground states. The active-site relaxation must be considered in understanding overall isomerization reactions in phytochromes, and such a molecular mechanism should be general.

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.

Acid-base equilibrium of the chromophore counterion results in distinct photoisomerization reactivity in the primary event of proteorhodopsin

Proteorhodopsin (PR) is a proton-pumping rhodopsin, and it is known to exhibit a multi-phasic decay of the excited-state population in the primary process. So far, this complex excited-state decay has been attributed to the branching of the relaxation pathway on the excited-state potential energy surface. However, a recent ultrafast spectroscopic study on a sodium-pumping rhodopsin suggested that such a complex decay may originate from the heterogeneity in the ground state due to the acid-base equilibrium of the counterion of the protonated retinal Schiff base (PRSB). In this study, we studied the excited-state dynamics of PR at pH 11 and 4, in which the counterion of the PRSB, Asp97, is completely deprotonated and protonated, respectively. The obtained time-resolved absorption data revealed that the excited-state lifetime is decisively governed by the protonation state of Asp97, and the photoisomerization of the PRSB chromophore proceeds faster and more efficiently when Asp97 is deprotonated. This conclusion was further supported by high similarity of the excited-state dynamics between PR at pH 4 and the D97N mutant in which Asp97 is replaced with neutral Asn. The results of this study suggest that the protonation state of the PRSB counterion plays a decisive role in determining the excited-state dynamics and the photoisomerization reactivity of rhodopsins in general, by making a significant influence on the exited-state potential energy surface of the PRSB chromophore.

Dependence on the Initial Configuration of Strong Field-Driven Isomerization of C2H2 Cations and Anions

We have investigated the femtosecond laser-induced fragmentation of C₂H₂ ^q ion beam targets in various initial configurations, including acetylene (linear HCCH), vinylidene (H₂CC), and cis/ trans. The initial configuration is shown to have a tremendous impact on the branching ratio of acetylene-like (CH ^q₁ + CH ^q₂) and vinylidene-like (C ^q₁' + CH₂ ^q₂') dissociation of a specific C₂H₂ ^q molecular ion. In particular, whereas C₂H₂⁺ generated from C₂H₂, a linear HCCH target, exhibits comparable levels of acetylene-like and vinylidene-like fragmentation, vinylidene or cis/ trans configuration ion beams preferably undergo vinylidene-like fragmentation, with an acetylene branching ratio ranging from 13.9% to zero.

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.

Retinal isomerization in bacteriorhodopsin captured by a femtosecond x-ray laser

Ultrafast isomerization of retinal is the primary step in photoresponsive biological functions including vision in humans and ion transport across bacterial membranes. We used an x-ray laser to study the subpicosecond structural dynamics of retinal isomerization in the light-driven proton pump bacteriorhodopsin. A series of structural snapshots with near-atomic spatial resolution and temporal resolution in the femtosecond regime show how the excited all-trans retinal samples conformational states within the protein binding pocket before passing through a twisted geometry and emerging in the 13-cis conformation. Our findings suggest ultrafast collective motions of aspartic acid residues and functional water molecules in the proximity of the retinal Schiff base as a key facet of this stereoselective and efficient photochemical reaction.

How Methylation Modifies the Photophysics of the Native All- trans-Retinal Protonated Schiff Base: A CASPT2/MD Study in Gas Phase and in Methanol

A comparison between the free-energy surfaces of the all- trans-retinal protonated Schiff base (RPSB) and its 10-methylated derivative in gas phase and methanol solution is performed at CASSCF//CASSCF and CASPT2//CASSCF levels. Solvent effects were included using the average solvent electrostatic potential from molecular dynamics method. This is a QM/MM (quantum mechanics/molecular mechanics) method that makes use of the mean field approximation. It is found that the methyl group bonded to C10 produces noticeable changes in the solution free-energy profile of the S₁ excited state, mainly in the relative stability of the minimum energy conical intersections (MECIs) with respect to the Franck-Condon (FC) point. The conical intersections yielding the 9- cis and 11- cis isomers are stabilized while that yielding the 13- cis isomer is destabilized; in fact, it becomes inaccessible by excitation to S₁. Furthermore, the planar S₁ minimum is not present in the methylated compound. The solvent notably stabilizes the S₂ excited state at the FC geometry. Therefore, if the S₂ state has an effect on the photoisomerization dynamics, it must be because it permits the RPSB population to branch around the FC point. All these changes combine to speed up the photoisomerization in the 10-methylated compound with respect to the native compound.

Blue-light-promoted carbon-carbon double bond isomerization and its application in the syntheses of quinolines

A blue-light-promoted carbon-carbon double bond isomerization in the absence of any photoredox catalyst is reported. It provides rapid access to a series of quinolines in good to excellent yields under simple aerobic conditions. The protocol is direct, catalyst-free and operationally convenient.

The origin of absorptive features in the two-dimensional electronic spectra of rhodopsin

A three-state three-mode model Hamiltonian reveals the origin of the absorptive features in the two-dimensional electronic spectra of rhodopsin.

Impact of Electronic State Mixing on the Photoisomerization Time Scale of the Retinal Chromophore

Spectral data show that the photoisomerization of retinal protonated Schiff base (rPSB) chromophores occurs on a 100 fs time scale or less in vertebrate rhodopsins, it is several times slower in microbial rhodopsins and it is between one and 2 orders of magnitude slower in solution. These time scale variations have been attributed to specific modifications of the topography of the first excited state potential energy surface of the chromophore. However, it is presently not clear which specific environment effects (e.g., electrostatic, electronic, or steric) are responsible for changing the surface topography. Here, we use QM/MM models and excited state trajectory computations to provide evidence for an increase in electronic mixing between the first and the second excited state of the chromophore when going from vertebrate rhodopsin to the solution environments. Ultimately, we argue that a correlation between the lifetime of the first excited state and electronic mixing between such state and its higher neighbor, may have been exploited to evolve rhodopsins toward faster isomerization and, possibly, light-sensitivity.

Fine Tuning of Retinal Photoinduced Decay in Solution

Single methylation at position C₁₀ of the all-trans retinal protonated Schiff base switches its excited-state decay in methanol from a slower picosecond into an ultrafast, protein-like subpicosecond process. QM/MM modeling in conjunction with on-the-fly excited-state dynamics provides fundamental understanding of the fine-tuning mechanics that "catalyzes" the photoinduced decay of solvated retinals. Methylation alters the interplay between the ionic S₁ and covalent S₂ states, reducing the excited-state lifetime by favoring the formation of a S₁ transient fluorescent state with fully inverted bond lengths that accounts for the recorded transient spectroscopy and from which a space-saving conical intersection seam is quickly (<1 ps) reached. Minimal and apparently innocent chemical modifications thus affect the characteristic intramolecular charge-transfer of the S₁ state as well as the interaction with the covalent S₂ excited state, eventually providing the high tunability of retinal photophysics and photochemistry and delivering a new concept for the rational design of retinal-based photoactive molecular devices.

Schiff Base Proton Acceptor Assists Photoisomerization of Retinal Chromophores in Bacteriorhodopsin

In this study, we investigated the ultrafast dynamics of bacteriorhodopsins (BRs) from Haloquadratum walsbyi (HwBR) and Haloarcula marismortui (HmBRI and HmBRII). First, the ultrafast dynamics were studied for three HwBR samples: wild-type, D93N mutation, and D104N mutation. The residues of the D93 and D104 mutants correspond to the control by the Schiff base proton acceptor and donor of the proton translocation subchannels. Measurements indicated that the negative charge from the Schiff base proton acceptor residue D93 interacts with the ultrafast and substantial change of the electrostatic potential associated with chromophore isomerization. By contrast, the Schiff base proton donor assists the restructuring of the chromophore cavity hydrogen-bond network during the thermalization of the vibrational hot state. Second, the ultrafast dynamics of the wild-types of HwBR, HmBRI, and HmBRII were compared. Measurements demonstrated that the hydrogen-bond network in the extracellular region in HwBR and HmBRII slows the photoisomerization of retinal chromophores, and the negatively charged helices on the cytoplasmic side of HwBR and HmBRII accelerate the thermalization of the vibrational hot state of retinal chromophores. The similarity of the correlation spectra of the wild-type HmBRI and D104N mutant of HwBR indicates that inactivation of the Schiff base proton donor induces a positive charge on the helices of the cytoplasmic side.

An On-the-Fly Surface-Hopping Program JADE for Nonadiabatic Molecular Dynamics of Polyatomic Systems: Implementation and Applications

Nonadiabatic dynamics simulations have rapidly become an indispensable tool for understanding ultrafast photochemical processes in complex systems. Here, we present our recently developed on-the-fly nonadiabatic dynamics package, JADE, which allows researchers to perform nonadiabatic excited-state dynamics simulations of polyatomic systems at an all-atomic level. The nonadiabatic dynamics is based on Tully's surface-hopping approach. Currently, several electronic structure methods (CIS, TDHF, TDDFT(RPA/TDA), and ADC(2)) are supported, especially TDDFT, aiming at performing nonadiabatic dynamics on medium- to large-sized molecules. The JADE package has been interfaced with several quantum chemistry codes, including Turbomole, Gaussian, and Gamess (US). To consider environmental effects, the Langevin dynamics was introduced as an easy-to-use scheme into the standard surface-hopping dynamics. The JADE package is mainly written in Fortran for greater numerical performance and Python for flexible interface construction, with the intent of providing open-source, easy-to-use, well-modularized, and intuitive software in the field of simulations of photochemical and photophysical processes. To illustrate the possible applications of the JADE package, we present a few applications of excited-state dynamics for various polyatomic systems, such as the methaniminium cation, fullerene (C20), p-dimethylaminobenzonitrile (DMABN) and its primary amino derivative aminobenzonitrile (ABN), and 10-hydroxybenzo[h]quinoline (10-HBQ).

Wavelength-dependent photocycle activity of xanthorhodopsin in the visible region

Xanthorhodopsin (xR) is a dual-chromophore proton-pump photosynthetic protein comprising one retinal Schiff base and one light-harvesting antenna salinixanthin (SX). The excitation wavelength-dependent transient population of the intermediate M demonstrates that the excitation of the retinal at 570 nm leads to the highest photocycle activity and the excitations of SX at 460 and 430 nm reduce the activity to ca. 37% relatively, suggesting an energy transfer pathway from the S₂ state of the SX to the S₁ state of the retinal and a quick internal vibrational relaxation in the S₂ state of SX prior to the energy transfer from SX to retinal.

•

Energy transfer efficiency from the salinixanthin (SX) to the retinal is ca. 37%.

•

Energy transfer efficiency is not dependent on wavelength at 486–430 nm.

•

Energy transfer from the S2 state of SX to the S2 state of retinal is less accessible.

Light-Activated Reversible Imine Isomerization: Towards a Photochromic Protein Switch

Mutants of cellular retinoic acid-binding protein II (CRABPII), engineered to bind all-trans-retinal as an iminium species, demonstrate photochromism upon irradiation with light at different wavelengths. UV light irradiation populates the cis-imine geometry, which has a high pKa , leading to protonation of the imine and subsequent "turn-on" of color. Yellow light irradiation yields the trans-imine isomer, which has a depressed pKa , leading to loss of color because the imine is not protonated. The protein-bound retinylidene chromophore undergoes photoinduced reversible interconversion between the colored and uncolored species, with excellent fatigue resistance.

Molecular bases for the selection of the chromophore of animal rhodopsins

The functions of microbial and animal rhodopsins are triggered by the isomerization of their all-trans and 11-cis retinal chromophores, respectively. To lay the molecular basis driving the evolutionary transition from the all-trans to the 11-cis chromophore, multiconfigurational quantum chemistry is used to compare the isomerization mechanisms of the sensory rhodopsin from the cyanobacterium Anabaena PCC 7120 (ASR) and of the bovine rhodopsin (Rh). It is found that, despite their evolutionary distance, these eubacterial and vertebrate rhodopsins start to isomerize via distinct implementations of the same bicycle-pedal mechanism originally proposed by Warshel [Warshel A (1976) Nature 260:678-683]. However, by following the electronic structure changes of ASR (featuring the all-trans chromophore) during the isomerization, we find that ASR enters a region of degeneracy between the first and second excited states not found in Rh (featuring the 11-cis chromophore). We show that such degeneracy is modulated by the preorganized structure of the chromophore and by the position of the reactive double bond. It is argued that the optimization of the electronic properties of the chromophore, which affects the photoisomerization efficiency and the thermal isomerization barrier, provided a key factor for the emergence of the striking amino acid sequence divergence observed between the microbial and animal rhodopsins.

A Bio-Inspired, Catalytic E → Z Isomerization of Activated Olefins

Herein, Nature's flavin-mediated activation of complex (poly)enes has been translated to a small molecule paradigm culminating in a highly (Z)-selective, catalytic isomerization of activated olefins using (-)-riboflavin (up to 99:1 Z/E). In contrast to the prominent Z → E isomerization of the natural system, it was possible to invert the directionality of the isomerization (E → Z) by simultaneously truncating the retinal scaffold, and introducing a third olefin substituent to augment A1,3-strain upon isomerization. Consequently, conjugation is reduced in the product chromophore leading to a substrate/product combination with discrete photophysical signatures. The operationally simple isomerization protocol has been applied to a variety of enone-derived substrates and showcased in the preparation of the medically relevant 4-substituted coumarin scaffold. A correlation of sensitizer triplet energy (ET) and reaction efficiency, together with the study of additive effects and mechanistic probes, is consistent with a triplet energy transfer mechanism.

An artificial molecular pump

Carrier proteins consume fuel in order to pump ions or molecules across cell membranes, creating concentration gradients. Their control over diffusion pathways, effected entirely through noncovalent bonding interactions, has inspired chemists to devise artificial systems that mimic their function. Here, we report a wholly artificial compound that acts on small molecules to create a gradient in their local concentration. It does so by using redox energy and precisely organized noncovalent bonding interactions to pump positively charged rings from solution and ensnare them around an oligomethylene chain, as part of a kinetically trapped entanglement. A redox-active viologen unit at the heart of a dumbbell-shaped molecular pump plays a dual role, first attracting and then repelling the rings during redox cycling, thereby enacting a flashing energy ratchet mechanism with a minimalistic design. Our artificial molecular pump performs work repetitively for two cycles of operation and drives rings away from equilibrium toward a higher local concentration.

Kinetics of exchange between zero-, one-, and two-hydrogen-bonded states of methyl and ethyl acetate in methanol

It has recently been shown that the ester carbonyl stretching vibration can be used as a sensitive probe of local electrostatic field in molecular systems. To further characterize this vibrational probe and extend its potential applications, we studied the kinetics of chemical exchange between differently hydrogen-bonded (H-bonded) ester carbonyl groups of methyl acetate (MA) and ethyl acetate (EA) in methanol. We found that, while both MA and EA can form zero, one, or two H-bonds with the solvent, the population of the 2hb state in MA is significantly smaller than that in EA. Using a combination of linear and nonlinear infrared measurements and numerical simulations, we further determined the rate constants for the exchange between these differently H-bonded states. We found that for MA the chemical exchange reaction between the two dominant states (i.e., 0hb and 1hb states) has a relaxation rate constant of 0.14 ps(-1), whereas for EA the three-state chemical exchange reaction occurs in a predominantly sequential manner with the following relaxation rate constants: 0.11 ps(-1) for exchange between 0hb and 1hb states and 0.12 ps(-1) for exchange between 1hb and 2hb states.

Investigation of the nonlinear absorption spectrum of all-trans retinoic acid by using the steady and transient two-photon absorption spectroscopy

This paper reports for the first time a complete study on the steady and transient excited state dynamics induced by 2PA for ATRA.