Femtosecond infrared spectroscopy of bacteriorhodopsin chromophore isomerization

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.

SpaiNN: equivariant message passing for excited-state nonadiabatic molecular dynamics

Excited-state molecular dynamics simulations are crucial for understanding processes like photosynthesis, vision, and radiation damage. However, the computational complexity of quantum chemical calculations restricts their scope. Machine learning offers a solution by delivering high-accuracy properties at lower computational costs. We present SpaiNN, an open-source Python software for ML-driven surface hopping nonadiabatic molecular dynamics simulations. SpaiNN combines the invariant and equivariant neural network architectures of SchNetPack with SHARC for surface hopping dynamics. Its modular design allows users to implement and adapt modules easily. We compare rotationally-invariant and equivariant representations in fitting potential energy surfaces of multiple electronic states and properties arising from the interaction of two electronic states. Simulations of the methyleneimmonium cation and various alkenes demonstrate the superior performance of equivariant SpaiNN models, improving accuracy, generalization, and efficiency in both training and inference.

Multiple retinal isomerizations during the early phase of the bestrhodopsin photoreaction

Bestrhodopsins constitute a class of light-regulated pentameric ion channels that consist of one or two rhodopsins in tandem fused with bestrophin ion channel domains. Here, we report on the isomerization dynamics in the rhodopsin tandem domains of Phaeocystis antarctica bestrhodopsin, which binds all-trans retinal Schiff-base (RSB) absorbing at 661 nm and, upon illumination, converts to the meta-stable P540 state with an unusual 11-cis RSB. The primary photoproduct P682 corresponds to a mixture of highly distorted 11-cis and 13-cis RSB directly formed from the excited state in 1.4 ps. P673 evolves from P682 in 500 ps and contains highly distorted 13-cis RSB, indicating that the 11-cis fraction in P682 converts to 13-cis. Next, P673 establishes an equilibrium with P595 in 1.2 µs, during which RSB converts to 11-cis and then further proceeds to P560 in 48 µs and P540 in 1.0 ms while remaining 11-cis. Hence, extensive isomeric switching occurs on the early ground state potential energy surface (PES) on the hundreds of ps to µs timescale before finally settling on a metastable 11-cis photoproduct. We propose that P682 and P673 are trapped high up on the ground-state PES after passing through either of two closely located conical intersections that result in 11-cis and 13-cis RSB. Co-rotation of C11=C12 and C13=C14 bonds results in a constricted conformational landscape that allows thermal switching between 11-cis and 13-cis species of highly strained RSB chromophores. Protein relaxation may release RSB strain, allowing it to evolve to a stable 11-cis isomeric configuration in microseconds.

The Functionality of the DC Pair in a Rhodopsin Guanylyl Cyclase from Catenaria anguillulae

Rhodopsin guanylyl cyclases (RGCs) belong to the class of enzymerhodopsins catalyzing the transition from GTP into the second messenger cGMP, whereas light-regulation of enzyme activity is mediated by a membrane-bound microbial rhodopsin domain, that holds the catalytic center inactive in the dark. Structural determinants for activation of the rhodopsin moiety eventually leading to catalytic activity are largely unknown. Here, we investigate the mechanistic role of the D283-C259 (DC) pair that is hydrogen bonded via a water molecule as a crucial functional motif in the homodimeric C. anguillulae RGC. Based on a structural model of the DC pair in the retinal binding pocket obtained by MD simulation, we analyzed formation and kinetics of early and late photocycle intermediates of the rhodopsin domain wild type and specific DC pair mutants by combined UV-Vis and FTIR spectroscopy at ambient and cryo-temperatures. By assigning specific infrared bands to S-H vibrations of C259 we are able to show that the DC pair residues are tightly coupled. We show that deprotonation of D283 occurs already in the inactive L state as a prerequisite for M state formation, whereas structural changes of C259 occur in the active M state and early cryo-trapped intermediates. We propose a comprehensive molecular model for formation of the M state that activates the catalytic moiety. It involves light induced changes in bond strength and hydrogen bonding of the DC pair residues from the early J state to the active M state and explains the retarding effect of C259 mutants.

Voltage Imaging with Engineered Proton-Pumping Rhodopsins: Insights from the Proton Transfer Pathway

Voltage imaging using genetically encoded voltage indicators (GEVIs) has taken the field of neuroscience by storm in the past decade. Its ability to create subcellular and network level readouts of electrical dynamics depends critically on the kinetics of the response to voltage of the indicator used. Engineered microbial rhodopsins form a GEVI subclass known for their high voltage sensitivity and fast response kinetics. Here we review the essential aspects of microbial rhodopsin photocycles that are critical to understanding the mechanisms of voltage sensitivity in these proteins and link them to insights from efforts to create faster, brighter and more sensitive microbial rhodopsin-based GEVIs.

Versatile Femtosecond Laser Synchronization for Multiple-Timescale Transient Infrared Spectroscopy

Several ways to electronically synchronize different types of amplified femtosecond laser systems are presented based on a single freely programmable electronics hardware: arbitrary-detuning asynchronous optical sampling (ADASOPS), as well as actively locking two femtosecond laser oscillators, albeit not necessarily to the same round-trip frequency. They allow us to rapidly probe a very wide range of timescales, from picoseconds to potentially seconds, in a single transient absorption experiment without the need to move any delay stage. Experiments become possible that address a largely unexplored aspect of many photochemical reactions, in particular in the context of photo-catalysis as well as photoactive proteins, where an initial femtosecond trigger very often initiates a long-lasting cascade of follow-up processes. The approach is very versatile and allows us to synchronize very different lasers, such as a Ti:Sa amplifier and a 100 kHz Yb-laser system. The jitter of the synchronization, and therewith the time-resolution in the transient experiment, lies in the range from 1 to 3 ps, depending on the method. For illustration, transient IR measurements of the excited state solvation and decay of a metal carbonyl complex as well as the full reaction cycle of bacteriorhodopsin are shown. The pros and cons of the various methods are discussed, with regard to the scientific question one might want to address, and also with regard to the laser systems that might be already existent in a laser lab.

Swellable Array Strategy Based on Designed Flexible Double Hypercross-linked Polymers for Synergistic Adsorption of Toluene and Formaldehyde

High-capacity adsorption and removal of complex volatile organic compounds (VOCs) from real-world environments is a tough challenge for researchers. Herein, a swellable array adsorption strategy was proposed to realize the synergistic adsorption of toluene and formaldehyde on the flexible double hypercross-linked polymers (FD-HCPs). FD-HCPs exhibited multiple adsorption sites awarded by a hydrophobic benzene ring/pyrrole ring and a hydrophilic hydroxyl structural unit. The array benzene ring, hydroxyl, and pyrrole N sites in FD-HCPs effectively captured toluene and formaldehyde molecules through π-π conjugation and electrostatic interaction and weakened their mutual competitive adsorption. Interestingly, the strong binding force of toluene molecules to the skeleton deformed the pore structure of FD-HCPs and generated new adsorption microenvironments for the other adsorbate. This behavior significantly improved the adsorption capacity of FD-HCPs for toluene and formaldehyde by 20% under multiple VOCs. Moreover, the pyrrole group in FD-HCPs greatly hindered H₂O molecule diffusion in the pore, thus efficiently weakening the competitive adsorption of H₂O toward VOCs. These fascinating properties enabled FD-HCPs to achieve synergistic adsorption for multicomponent VOC vapor under a highly humid environment and overcame single-species VOC adsorption properties on state-of-the-art porous adsorbents. This work provides the practical feasibility of synergistic adsorption to remove complex VOCs in real-world environments.

Convergent evolution of animal and microbial rhodopsins

Rhodopsins, a family of photoreceptive membrane proteins, contain retinal as a chromophore and were firstly identified as reddish pigments from frog retina in 1876. Since then, rhodopsin-like proteins have been identified mainly from animal eyes. In 1971, a rhodopsin-like pigment was discovered from the archaeon Halobacterium salinarum and named bacteriorhodopsin. While it was believed that rhodopsin- and bacteriorhodopsin-like proteins were expressed only in animal eyes and archaea, respectively, before the 1990s, a variety of rhodopsin-like proteins (called animal rhodopsins or opsins) and bacteriorhodopsin-like proteins (called microbial rhodopsins) have been progressively identified from various tissues of animals and microorganisms, respectively. Here, we comprehensively introduce the research conducted on animal and microbial rhodopsins. Recent analysis has revealed that the two rhodopsin families have common molecular properties, such as the protein structure (i.e., 7-transmembrane structure), retinal structure (i.e., binding ability to cis- and trans-retinal), color sensitivity (i.e., UV- and visible-light sensitivities), and photoreaction (i.e., triggering structural changes by light and heat), more than what was expected at the early stages of rhodopsin research. Contrastingly, their molecular functions are distinctively different (e.g., G protein-coupled receptors and photoisomerases for animal rhodopsins and ion transporters and phototaxis sensors for microbial rhodopsins). Therefore, based on their similarities and dissimilarities, we propose that animal and microbial rhodopsins have convergently evolved from their distinctive origins as multi-colored retinal-binding membrane proteins whose activities are regulated by light and heat but independently evolved for different molecular and physiological functions in the cognate organism.

Direct structural observation of ultrafast photoisomerization dynamics in sinapate esters

Sinapate esters have been extensively studied for their potential application in 'nature-inspired' photoprotection. There is general consensus that the relaxation mechanism of sinapate esters following photoexcitation with ultraviolet radiation is mediated by geometric isomerization. This has been largely inferred through indirect studies involving transient electronic absorption spectroscopy in conjunction with steady-state spectroscopies. However, to-date, there is no direct experimental evidence tracking the formation of the photoisomer in real-time. Using transient vibrational absorption spectroscopy, we report on the direct structural changes that occur upon photoexcitation, resulting in the photoisomer formation. Our mechanistic analysis predicts that, from the photoprepared ππ* state, internal conversion takes place through a conical intersection (CI) near the geometry of the initial isomer. Our calculations suggest that different CI topographies at relevant points on the seam of intersection may influence the isomerization yield. Altogether, we provide compelling evidence suggesting that a sinapate ester's geometric isomerization can be a more complex dynamical process than originally thought.

Photoinduced isomerization sampling of retinal in bacteriorhodopsin

Photoisomerization of retinoids inside a confined protein pocket represents a critical chemical event in many important biological processes from animal vision, nonvisual light effects, to bacterial light sensing and harvesting. Light-driven proton pumping in bacteriorhodopsin entails exquisite electronic and conformational reconfigurations during its photocycle. However, it has been a major challenge to delineate transient molecular events preceding and following the photoisomerization of the retinal from noisy electron density maps when varying populations of intermediates coexist and evolve as a function of time. Here, I report several distinct early photoproducts deconvoluted from the recently observed mixtures in time-resolved serial crystallography. This deconvolution substantially improves the quality of the electron density maps, hence demonstrates that the all-trans retinal undergoes extensive isomerization sampling before it proceeds to the productive 13-cis configuration. Upon light absorption, the chromophore attempts to perform trans-to-cis isomerization at every double bond together with the stalled anti-to-syn rotations at multiple single bonds along its polyene chain. Such isomerization sampling pushes all seven transmembrane helices to bend outward, resulting in a transient expansion of the retinal binding pocket, and later, a contraction due to recoiling. These ultrafast responses observed at the atomic resolution support that the productive photoreaction in bacteriorhodopsin is initiated by light-induced charge separation in the prosthetic chromophore yet governed by stereoselectivity of its protein pocket. The method of a numerical resolution of concurrent events from mixed observations is also generally applicable.

Simple approach to broadband mid-infrared pulse generation with a mode-locked Yb-doped fiber laser

Broadband mid-infrared (MIR) molecular spectroscopy demands a bright and broadband light source in the molecular fingerprint region. To this end, intra-pulse difference frequency generation (IDFG) has shown excellent properties among various techniques. Although IDFG systems pumped with 1.5- or 2-µm ultrashort pulsed lasers have been extensively developed, few systems have been demonstrated with 1-µm lasers, which use bulky 100-W-class high-power Yb thin-disk lasers. In this work, we demonstrate a simple and robust approach of 1-µm-pumped broadband IDFG with a conventional mode-locked Yb-doped fiber laser. We first generate 3.3-W, 12.1-fs ultrashort pulses at 50 MHz by a simple combination of spectral broadening with a short single-mode fiber and pulse compression with chirped mirrors. Then, we use them for pumping a thin orientation-patterned gallium phosphide crystal, generating 1.2-mW broadband MIR pulses with the -20-dB bandwidth of 480 cm^-1 in the fingerprint region (760-1240 cm^-1, 8.1-13.1 µm). The 1-µm-based IDFG system allows for additional generations of ultrashort pulses in the ultraviolet and visible regions, enabling, for example, 50-MHz-level high-repetition-rate vibrational sum-frequency generation spectroscopy or pump-probe spectroscopy.

Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─Quo Vadis?

Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.

Impact of Ultrafast Electric Field Changes on Photoreceptor Protein Dynamics

Studies on photoreceptors provide a wealth of information on cofactor and protein dynamics on the microsecond to seconds time-scale. Up to now, ultrafast dynamics addresses mainly the cofactor or chromophore, but ultrafast protein dynamics are poorly understood. Increasing evidence show that protein responses can occur even faster than the cofactor dynamics. The causal reason for the ultrafast protein response cannot be explained by the localized cofactor excitation or its excited-state decay, alone. We propose a Coulomb interaction mechanism started by a shock wave and stabilized by a dipole moment change at least partially responsible for coherent oscillations in proteins, protonation changes, water dislocations, and protein changes prior to and beyond chromophore's excited-state decay. Photoexcitation changes the electron density distribution of the chromophore within a few femtoseconds: The Coulomb shock wave affects polar groups, hydrogen bonds, and protein bound water molecules. The process occurs on a time-scale even faster than excited-state decay of the chromophore. We discuss studies on selected photoreceptors in light of this mechanism and its impact on a detailed understanding of protein dynamics.

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).

Efficient Degradation of Iopromide by Using Sulfite Activated with Mackinawite

Iopromide (IOP), an iodinated X-ray contrast medium (ICM), is identified as a precursor to iodide disinfection byproducts that have high genotoxicity and cytotoxicity to mammals. ICM remains persistent through typical wastewater treatment processes and even through some hydroxyl radical-based advanced oxidation processes. The development of new technologies to remove ICMs is needed. In this work, mackinawite (FeS)-activated sulfite autoxidation was employed for the degradation of IOP-containing water. The experiment was performed in a 500 mL self-made temperature-controlled reactor with online monitoring pH and dissolved oxygen in the laboratory. The effects of various parameters, such as initial pH values, sulfite dosages, FeS dosages, dissolved oxygen, and inorganic anions on the performance of the treatment process have been investigated. Eighty percent of IOP could be degraded in 15 min with 1 g L-1 FeS, 400 μmol L-1 sulfite at pH 8, and high efficiency on the removal of total organic carbon (TOC) was achieved, which is 71.8% via a reaction for 1 h. The generated hydroxyl and oxysulfur radicals, which contributed to the oxidation process, were identified through radical quenching experiments. The dissolved oxygen was essential for the degradation of IOP. The presence of Cl- could facilitate IOP degradation, while NO3- and CO32- could inhibit the degradation process. The reaction pathway involving H-abstraction and oxidative decarboxylation was proposed, based on product identification. The current system shows good applicability for the degradation of IOP and may help in developing a new approach for the treatment of ICM-containing water.

A computational methodology to diagnose sequence-variant dynamic perturbations by comparing atomic protein structures

MOTIVATION

The objective is to diagnose dynamics perturbations caused by amino-acid mutations as prerequisite to assess protein functional health or drug failure, simply using network models of protein X-ray structures.

RESULTS

We find that the differences in the allocation of the atomic interactions of each amino acid to 1D, 2D, 3D, 4D structural levels between variants structurally robust, recover experimental dynamic perturbations. The allocation measure validated on two B-pentamers variants of AB5 toxins having 17 mutations, also distinguishes dynamic perturbations of pathogenic and non-pathogenic Transthyretin single-mutants. Finally, the main proteases of the coronaviruses SARS-CoV and SARS-CoV-2 exhibit changes in the allocation measure, raising the possibility of drug failure despite the main proteases structural similarity.

AVAILABILITY AND IMPLEMENTATION

The Python code used for the production of the results is available at github.com/lorpac/protein_partitioning_atomic_contacts. The authors will run the analysis on any PDB structures of protein variants upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Multiphoton Control of 6π Photocyclization via State-Dependent Reactant-Product Correlations

Multiphoton excitation promises opportunities for opening new photochemical reaction pathways and controlling photoproduct distributions. We demonstrate photonic control of the 6π photocyclization of ortho-terphenyl to make 4a,4b-dihydrotriphenylene (DHT). Using pump-repump-probe spectroscopy we show that 1 + 1' excitation to a high-lying reactant electronic state generates a metastable species characterized by a red absorption feature that accompanies a repump-induced depletion in the one-photon trans-dihydro product (trans-DHT); signatures of the new photoproduct are clearer for a structural analogue of the reactant that is sterically inhibited against one-photon cyclization. Quantum-chemical computations support assignment of this species to cis-DHT, which is accessible photochemically along a disrotatory coordinate from high-lying electronic states reached by 1 + 1' excitation. We use time-resolved spectroscopy to track photochemical dynamics producing cis-DHT. In total, we demonstrate that selective multiphoton excitation opens a new photoreaction channel in these photocyclizing reactants by taking advantage of state-dependent correlations between reactant and product electronic states.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Ultrafast proton release reaction and primary photochemistry of phycocyanobilin in solution observed with fs-time-resolved mid-IR and UV/Vis spectroscopy

Deactivation processes of photoexcited (λ_ex = 580 nm) phycocyanobilin (PCB) in methanol were investigated by means of UV/Vis and mid-IR femtosecond (fs) transient absorption (TA) as well as static fluorescence spectroscopy, supported by density-functional-theory calculations of three relevant ground state conformers, PCB_A, PCB_B and PCB_C, their relative electronic state energies and normal mode vibrational analysis. UV/Vis fs-TA reveals time constants of 2.0, 18 and 67 ps, describing decay of PCB_B*, of PCB_A* and thermal re-equilibration of PCB_A, PCB_B and PCB_C, respectively, in line with the model by Dietzek et al. (Chem Phys Lett 515:163, 2011) and predecessors. Significant substantiation and extension of this model is achieved first via mid-IR fs-TA, i.e. identification of molecular structures and their dynamics, with time constants of 2.6, 21 and 40 ps, respectively. Second, transient IR continuum absorption (CA) is observed in the region above 1755 cm^-1 (CA1) and between 1550 and 1450 cm^-1 (CA2), indicative for the IR absorption of highly polarizable protons in hydrogen bonding networks (X-H^…Y). This allows to characterize chromophore protonation/deprotonation processes, associated with the electronic and structural dynamics, on a molecular level. The PCB photocycle is suggested to be closed via a long living (> 1 ns), PCB_C-like (i.e. deprotonated), fluorescent species.

Broadband second-harmonic generation from artificial optical nonlinearity

In this Letter, we present a mechanism for effectively broadening the bandwidth of second-harmonic generation (SHG) with the metamaterial-based artificial optical nonlinearity. As the nonlinear response of the artificial nonlinearity arising from the magnetoelectric coupling constructed by the meta-molecule (MM) structure, the broadband second-order nonlinearity can be built by simply combining the MMs with different geometrical sizes together. The physical model and the numerical simulation fully support the artificial generation and modulation of the broadband second harmonics. Our work suggests a new route for realizing the on-chip custom-designed nonlinear optical devices with broadband operation.

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.

Analysis of Nanoconfined Protein Dielectric Signals Using Charged Amino Acid Network Models

Protein slow motions involving collective molecular fluctuations on the timescale of microseconds to seconds are difficult to measure and not well understood despite being essential to sustain protein folding and protein function. Broadband dielectric spectroscopy (BDS) is one of the most powerful experimental techniques to monitor, over a broad frequency and temperature range, the molecular dynamics of soft matter through the orientational polarisation of permanent dipole moments that are generated by the chemical structure and morphological organisation of matter. Its typical frequency range goes from 107 Hz down to 10−3 Hz, being thus suitable for investigations on slow motions in proteins. Moreover, BDS has the advantage of providing direct experimental access to molecular fluctuations taking place on different length-scales, from local to cooperative dipolar motions. The unfolding of the cholera toxin B pentamer (CtxB5) after thermal treatment for 3h at 80°C is investigated by BDS under nanoconfined and dehydrated conditions. From the X-ray structure of the toxin pentamer, network-based models are used to infer the toxin dipoles present in the native state and to compute their stability and dielectric properties. Network analyses highlight three domains with distinct dielectric and stability properties that support a model where the toxin unfolds into three conformations after the treatment at 80°C. This novel integrative approach offers some perspective into the investigation of the relation between local perturbations (e.g. mutation, thermal treatment) and larger scale protein conformational changes. It might help ranking protein sequence variants according to their respective scale of dynamics perturbations.

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.

Experimental Protein Molecular Dynamics: Broadband Dielectric Spectroscopy coupled with nanoconfinement

Protein dynamics covers multiple spatiotemporal scale processes, among which slow motions, not much understood even though they are underlying protein folding and protein functions. Protein slow motions are associated with structural heterogeneity, short-lived and poorly populated conformations, hard to detect individually. In addition, they involve collective motions of many atoms, not easily tracked by simulation and experimental devices. Here we propose a biophysical approach, coupling geometrical nanoconfinement and broadband dielectric spectroscopy (BDS), which distinguishes protein conformations by their respective molecular dynamics. In particular, protein-unfolding intermediates, usually poorly populated in macroscopic solutions are detected. The protein dynamics is observed under unusual conditions (sample nanoconfinement and dehydration) highlighting the robustness of protein structure and protein dynamics to a variety of conditions consistent with protein sustainability. The protein dielectric signals evolve with the temperature of thermal treatments indicating sensitivity to atomic and molecular interaction changes triggered by the protein thermal unfolding. As dipole fluctuations depend on both collective large-scale motions and local motions, the approach offers a prospect to track in-depth unfolding events.

Machine learning enables long time scale molecular photodynamics simulations

Photo-induced processes are fundamental in nature but accurate simulations of their dynamics are seriously limited by the cost of the underlying quantum chemical calculations, hampering their application for long time scales. Here we introduce a method based on machine learning to overcome this bottleneck and enable accurate photodynamics on nanosecond time scales, which are otherwise out of reach with contemporary approaches. Instead of expensive quantum chemistry during molecular dynamics simulations, we use deep neural networks to learn the relationship between a molecular geometry and its high-dimensional electronic properties. As an example, the time evolution of the methylenimmonium cation for one nanosecond is used to demonstrate that machine learning algorithms can outperform standard excited-state molecular dynamics approaches in their computational efficiency while delivering the same accuracy.

Protein Dynamics Preceding Photoisomerization of the Retinal Chromophore in Bacteriorhodopsin Revealed by Deep-UV Femtosecond Stimulated Raman Spectroscopy

Bacteriorhodopsin is a prototypical photoreceptor protein that functions as a light-driven proton pump. The retinal chromophore of bacteriorhodopsin undergoes C₁₃═C₁₄ trans-to-cis isomerization upon photoexcitation, and it has been believed to be the first event that triggers the cascaded structural changes in bacteriorhodopsin. We investigated the protein dynamics of bacteriorhodopsin using deep-ultraviolet resonance femtosecond stimulated Raman spectroscopy. It was found that the stimulated Raman signals of tryptophan and tyrosine residues exhibit significant changes within 0.2 ps after photoexcitation while they do not noticeably change during the isomerization process. This result implies that the protein environment changes first, and its change is small during isomerization. The obtained femtosecond stimulated Raman data indicate that ultrafast change is induced in the protein part by the sudden creation of the large dipole of the excited-state chromophore, providing an environment that realizes efficient and selective isomerization.

Three-dimensional view of ultrafast dynamics in photoexcited bacteriorhodopsin

Bacteriorhodopsin (bR) is a light-driven proton pump. The primary photochemical event upon light absorption is isomerization of the retinal chromophore. Here we used time-resolved crystallography at an X-ray free-electron laser to follow the structural changes in multiphoton-excited bR from 250 femtoseconds to 10 picoseconds. Quantum chemistry and ultrafast spectroscopy were used to identify a sequential two-photon absorption process, leading to excitation of a tryptophan residue flanking the retinal chromophore, as a first manifestation of multiphoton effects. We resolve distinct stages in the structural dynamics of the all-trans retinal in photoexcited bR to a highly twisted 13-cis conformation. Other active site sub-picosecond rearrangements include correlated vibrational motions of the electronically excited retinal chromophore, the surrounding amino acids and water molecules as well as their hydrogen bonding network. These results show that this extended photo-active network forms an electronically and vibrationally coupled system in bR, and most likely in all retinal proteins.

Bacteriorhodopsin (bR) is a light-driven proton pump. Here the authors combine time-resolved crystallography at a free-electron laser, ultrafast spectroscopy and quantum chemistry to study the structural changes following multiphoton photoexcitation of bR and find that they occur within 300 fs not only in the light-absorbing chromophore but also in the surrounding protein.

Carrier-to-envelope phase-stable, mid-infrared, ultrashort pulses from a hybrid parametric generator: Cr:ZnSe laser amplifier system

Our Cr:ZnSe laser amplifier, seeded by parametric difference mixing, produces 72fs long pulses at the central wavelength of ~2.37µm. The stability of the carrier-to-envelope phase of the amplified seed pulses, attained at the stage of their parametric generation, is preserved through 6 orders of magnitude of laser amplification.

Bacteriorhodopsin: Structural Insights Revealed Using X-Ray Lasers and Synchrotron Radiation

Directional transport of protons across an energy transducing membrane—proton pumping—is ubiquitous in biology. Bacteriorhodopsin (bR) is a light-driven proton pump that is activated by a buried all- trans retinal chromophore being photoisomerized to a 13- cis conformation. The mechanism by which photoisomerization initiates directional proton transport against a proton concentration gradient has been studied by a myriad of biochemical, biophysical, and structural techniques. X-ray free electron lasers (XFELs) have created new opportunities to probe the structural dynamics of bR at room temperature on timescales from femtoseconds to milliseconds using time-resolved serial femtosecond crystallography (TR-SFX). Wereview these recent developments and highlight where XFEL studies reveal new details concerning the structural mechanism of retinal photoisomerization and proton pumping. We also discuss the extent to which these insights were anticipated by earlier intermediate trapping studies using synchrotron radiation. TR-SFX will open up the field for dynamical studies of other proteins that are not naturally light-sensitive.

Effect of point mutations on the ultrafast photo-isomerization of Anabaena sensory rhodopsin

Anabaena sensory rhodopsin (ASR) is a particular microbial retinal protein for which light-adaptation leads to the ability to bind both the all-trans, 15-anti (AT) and the 13-cis, 15-syn (13C) isomers of the protonated Schiff base of retinal (PSBR). In the context of obtaining insight into the mechanisms by which retinal proteins catalyse the PSBR photo-isomerization reaction, ASR is a model system allowing to study, within the same protein, the protein-PSBR interactions for two different PSBR conformers at the same time. A detailed analysis of the vibrational spectra of AT and 13C, and their photo-products in wild-type ASR obtained through femtosecond (pump-) four-wave-mixing is reported for the first time, and compared to bacterio- and channelrhodopsin. As part of an extensive study of ASR mutants with blue-shifted absorption spectra, we present here a detailed computational analysis of the origin of the mutation-induced blue-shift of the absorption spectra, and identify electrostatic interactions as dominating steric effects that would entail a red-shift. The excited state lifetimes and isomerization reaction times (IRT) for the three mutants V112N, W76F, and L83Q are studied experimentally by femtosecond broadband transient absorption spectroscopy. Interestingly, in all three mutants, isomerization is accelerated for AT with respect to wild-type ASR, and this the more, the shorter the wavelength of maximum absorption. On the contrary, the 13C photo-reaction is slightly slowed down, leading to an inversion of the ESLs of AT and 13C, with respect to wt-ASR, in the blue-most absorbing mutant L83Q. Possible mechanisms for these mutation effects, and their steric and electrostatic origins are discussed.

Ultrafast photoisomerisation of an isolated retinoid

The photoinduced excited state dynamics of gas-phase trans-retinoate (deprotonated trans-retinoic acid, trans-RA^-) are studied using tandem ion mobility spectrometry coupled with laser spectroscopy, and frequency-, angle- and time-resolved photoelectron imaging. Photoexcitation of the bright S₃(ππ*) ← S₀ transition leads to internal conversion to the S₁(ππ*) state on a ≈80 fs timescale followed by recovery of S₀ and concomitant isomerisation to give the 13-cis (major) and 9-cis (minor) photoisomers on a ≈180 fs timescale. The sub-200 fs stereoselective photoisomerisation parallels that for the retinal protonated Schiff base chromophore in bacteriorhodopsin. Measurements on trans-RA^- in methanol using the solution photoisomerisation action spectroscopy technique show that 13-cis-RA^- is also the principal photoisomer, although the 13-cis and 9-cis photoisomers are formed with an inverted branching ratio with photon energy in methanol when compared with the gas phase, presumably due to solvent-induced modification of potential energy surfaces and inhibition of electron detachment processes. Comparison of the gas-phase time-resolved data with transient absorption spectroscopy measurements on retinoic acid in methanol suggest that photoisomerisation is roughly six times slower in solution. This work provides clear evidence that solvation significantly affects the photoisomerisation dynamics of retinoid molecules.

Excited State Vibrational Spectra of All- trans Retinal Derivatives in Solution Revealed By Pump-DFWM Experiments

The ultrafast structural changes during the photoinduced isomerization of the retinal-protonated Schiff base (RPSB) is still a poorly understood aspect in the retinal's photochemistry. In this work, we apply pump-degenerate four-wave mixing (pump-DFWM) to all- trans retinal (ATR) and retinal Schiff bases (RSB) to resolve coherent high- and low-frequency vibrational signatures from excited electronic states. We show that the vibrational spectra of excited singlet states in these samples exhibit pronounced differences compared to the relaxed ground state. Pump-DFWM results indicate three major features for ATR and RSB. (i) Excited state vibrational spectra of ATR and RSB consist predominately of low-frequency modes in the energetic range 100-500 cm^-1. (ii) Excited state vibrational spectra show distinct differences for excitation in specific regions of electronic transitions of excited state absorption and emission. (iii) Low-frequency modes in ATR and RSB are inducible during the entire lifetime of the excited electronic states. This latter effect points to a transient molecular structure that, following initial relaxation between different excited electronic states, does not change anymore over the lifetime of the finally populated excited electronic state.

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.

Spiers Memorial Lecture. Introductory lecture: the impact of structure on photoinduced processes in nucleic acids and proteins

Light is an important environmental variable and most organisms have evolved means to sense, exploit or avoid it and to repair detrimental effects on their genome. In general, light absorption is the task of specific chromophores, however other biomolecules such as oligonucleotides also do so which can result in undesired outcomes such as mutations and cancer. Given the biological importance of light-induced processes and applications for imaging, optogenetics, photodynamic therapy or photovoltaics, there is a great interest in understanding the detailed molecular mechanisms of photoinduced processes in proteins and nucleic acids. The processes are typically characterized by time-resolved spectroscopic approaches or computation, inferring structural information on transient species from stable ground state structures. Recently, however, structure determination of excited states or other short-lived species has become possible with the advent of X-ray free-electron lasers. This review gives an overview of the impact of structure on the understanding of photoinduced processes in macromolecules, focusing on systems presented at this Faraday Discussion meeting.

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.

Protein dynamics observed by tunable mid-IR quantum cascade lasers across the time range from 10ns to 1s

We have developed a spectrometer based on tunable quantum cascade lasers (QCLs) for recording time-resolved absorption spectra of proteins in the mid-infrared range. We illustrate its performance by recording time-resolved difference spectra of bacteriorhodopsin in the carboxylic range (1800-1700cm^-1) and on the CO rebinding reaction of myoglobin (1960-1840cm^-1), at a spectral resolution of 1cm^-1. The spectrometric setup covers the time range from 4ns to nearly a second with a response time of 10-15ns. Absorption changes as low as 1×10^-4 are detected in single-shot experiments at t>1μs, and of 5×10^-6 in kinetics obtained after averaging 100 shots. While previous time-resolved IR experiments have mostly been conducted on hydrated films of proteins, we demonstrate here that the brilliance of tunable quantum cascade lasers is superior to perform ns time-resolved experiments even in aqueous solution (H₂O).

Different hydrogen bonding environments of the retinal protonated Schiff base control the photoisomerization in channelrhodopsin-2

The first event of the channelrhodopsin-2 (ChR2) photocycle, i.e. trans-to-cis photoisomerization, is studied by means of quantum mechanics/molecular mechanics, taking into account the flexible retinal environment in the ground state.

Shedding Light on the Photoisomerization Pathway of Donor-Acceptor Stenhouse Adducts

Donor-acceptor Stenhouse adducts (DASAs) are negative photochromes that hold great promise for a variety of applications. Key to optimizing their switching properties is a detailed understanding of the photoswitching mechanism, which, as yet, is absent. Here we characterize the actinic step of DASA-photoswitching and its key intermediate, which was studied using a combination of ultrafast visible and IR pump-probe spectroscopies and TD-DFT calculations. Comparison of the time-resolved IR spectra with DFT computations allowed to unambiguously identify the structure of the intermediate, confirming that light absorption induces a sequential reaction path in which a Z-E photoisomerization of C2-C3 is followed by a rotation around C3-C4 and a subsequent thermal cyclization step. First and second-generation DASAs share a common photoisomerization mechanism in chlorinated solvents with notable differences in kinetics and lifetimes of the excited states. The photogenerated intermediate of the second-generation DASA was photo-accumulated at low temperature and probed with time-resolved spectroscopy, demonstrating the photoreversibility of the isomerization process. Taken together, these results provide a detailed picture of the DASA isomerization pathway on a molecular level.